CN104594872B - A kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy - Google Patents

A kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy Download PDF

Info

Publication number
CN104594872B
CN104594872B CN201510001929.0A CN201510001929A CN104594872B CN 104594872 B CN104594872 B CN 104594872B CN 201510001929 A CN201510001929 A CN 201510001929A CN 104594872 B CN104594872 B CN 104594872B
Authority
CN
China
Prior art keywords
crack
horizontal well
reservoir
congruence
pressure break
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510001929.0A
Other languages
Chinese (zh)
Other versions
CN104594872A (en
Inventor
曾凡辉
郭建春
龙川
蒋豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN201510001929.0A priority Critical patent/CN104594872B/en
Publication of CN104594872A publication Critical patent/CN104594872A/en
Application granted granted Critical
Publication of CN104594872B publication Critical patent/CN104594872B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/25Methods for stimulating production
    • E21B43/26Methods for stimulating production by forming crevices or fractures
    • E21B41/0092

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The present invention relates to the fracturing reform in oil-gas field development field, the method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy is more particularly to, is mainly comprised the following steps:(1) reservoir, fluid properties, horizontal well pit shaft basic parameter are collected;(2) pressure break horizontal well crack basic parameter is collected;(3) crack of pressure break horizontal well is averagely divided into the congruence of equal length along seam length direction;(4) the Reservoir Seepage model of tight gas reservoir pressure break horizontal well Fracture System is set up;(5) Pressure Drop Model that gas flows in crack is set up;(6) flow model of the coupled gas in Reservoir Seepage and crack is set up, the Production rate model of tight gas reservoir pressure break horizontal well is formed;(7) fracture condudtiviy of tight gas reservoir pressure break horizontal well is optimized.The method provided using the present invention can overcome the deficiencies in the prior art, effectively solve optimization problem of the pressure break horizontal well along the seam non-constant fracture condudtiviy of length direction, so as to provide rational basis for the optimization design of reservoir reconstruction, improve reservoir reconstruction effect.

Description

A kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy
Technical field
The present invention relates to the fracturing reform in oil-gas field development field, it is more particularly to for optimizing tight gas reservoir frac water The method of horizontal well fracture condudtiviy.
Background technology
Staged fracturing of horizontal well is the effective measures for developing tight gas reservoir.Hydraulic fracturing is a kind of by injecting high to stratum Pressure fluid makes formation fracture in the case of more than stratum imbibition ability;The mulling liquid containing proppant is subsequently injected into, promotes crack Further extend into reservoir deep;After construction terminates, due to the supporting role of proppant fracture wall so that fracture surface will not Close completely and remain in that certain flow conductivity.These have the crack of flow conductivity ensure reservoir hydrocarbons can Channel Group, So that Oil/gas Well is increased production and improved recovery ratio;By the flow conductivity energy for optimizing tight gas reservoir pressure break horizontal well crack Enough significantly improve the yield of tight gas reservoir pressure break horizontal well.
The method of conventional optimization tight gas reservoir pressure break horizontal well fracture condudtiviy is at present:Initially set up pressure break level Well production computation model;The pressure break horizontal well Production rate model set up is then based on, is joined in given reservoir basic parameter, gas Under conditions of number, horizontal well parameter and fracture parameters, the pressure break horizontal well production under different fracture condudtiviies is simulated;Finally Carry out preferred fracture condudtiviy according to cumulative production under the certain time of simulation.Accurate preferably tight gas reservoir pressure break horizontal well crack The key of flow conductivity is to set up correct pressure break horizontal well Production rate model.At present, relevant pressure break horizontal well Production rate Model be broadly divided into two classes:One class is the production to fractured horizontal well formation many cracks using commercial numerical simulation software Amount prediction carries out calculating analysis, and (finite element method [J] hydrodynamicses of the new Horizontal Well Reservoir Numerical Simulations of Zhou Zhumei, Lang Zhao are ground Study carefully and A volumes of .1996,11 (3) of progress:261-270;Zhang Xuewen, Fang Hongchang, Qiu Yinan wait low-permeability oil deposit pressure breaks horizontal well to produce Can influence factor [J] petroleum journals, 1999,20 (4):51-54).Another kind of is to use Analytic Calculation Method, is converged according to source Thought, is divided into two stages of flowing in Reservoir Seepage and crack, respectively by the flow process of tight gas reservoir pressure break horizontal well gas Set up corresponding mathematical modeling and describe its flow process;Then the equation that coupled gas flows in Reservoir Seepage and crack, Solve the yield of tight gas reservoir pressure break horizontal well and the flow conductivity in preferred crack.Its basic process solved includes five substantially Step:
The first step:Pressure break horizontal well crack single-blade is divided into the congruence of equal length, the production of each congruence can With approximate processing into the production of a bite straight well, principle of stacking is utilized, it is established that in crack tip production multiple congruences are produced simultaneously when Raw drop of pressure (Aissa Zerzar.Interpretation of Multiple Hydraulically Fractured Horizontal Wells in Closed Systems[R]SPE 84888,2003;Zeng Fanhui, Guo Jianchun, Zhao Jinzhou, wait Influence factor analysis [J] Petroleum finances of fractured horizontal well's productivity, 2007,34 (4):474-477;Sun Hai, Yao Jun, Lian Peiqing, waits to consider pressure break horizontal well unsteady Model [J] the petroleum journals of basement rock to pit shaft feed flow, 2012,33 (1): 117-122)。
Second step:Solve the pressure drop equation that gas flows in crack.Flowing for gas in crack at present, or It is to be processed into along the invariable miniature planar radial flow of seam length direction flow conductivity that (Zeng Fanhui, Guo Jianchun, Zhao Jinzhou wait shadows Ring factor analysis [J] Petroleum finances .2007,34 (4) of fractured horizontal well's productivity:474-477;Sun Hai, Yao Jun, it is honest and clean Training celebrating, waits to consider pressure break horizontal well unsteady Model [J] the petroleum journals of basement rock to pit shaft feed flow, 2012,33 (1):117- 122), otherwise be processed into along the invariable sheet flow model of seam length direction flow conductivity (Gringarten, A.C., Ramey,H.I.Jr.,Raghavan,R.Unsteady State Pressure Distributions Created by a Well with a Single Infinite-Conductivity Vertical Fracture[J]:SPEJ,1974:347- 60)。
3rd step:Meet that pressure is equal, flow is continuously substantially former when being flowed using gas in Reservoir Seepage and crack Then, it is established that containing using each congruence yield as variable and the production of equal with the several numbers of total congruence tight gas reservoir pressure break horizontal well Equation is measured, because the number of equation and the variable number of congruence yield are equal, the yield for solving and obtaining each congruence can be closed.
4th step:The yield of each congruence is added up, the yield of whole pressure break horizontal well is obtained.
5th step:It is distributed by the fracture condudtiviy for changing pressure break horizontal well, contrasts different fracture condudtiviy distributions Under pressure break horizontal well crack output distribotion and cumulative production, the fracture condudtiviy of Optimum Fracturing horizontal well.
The method of current these methods optimization tight gas reservoir pressure break horizontal well fracture condudtiviy has following significant drawback:
Assuming that fracture condudtiviy is invariable along seam length, the result for causing optimization can only be average crack flow conductivity (Zeng Fanhui, Guo Jianchun, Zhao Jinzhou, wait to influence factor analysis [J] Petroleum finance .2007 of fractured horizontal well's productivity, 34(4):474-477;Sun Hai, Yao Jun, Lian Peiqing, wait to consider pressure break horizontal well unsteady Model of the basement rock to pit shaft feed flow [J] petroleum journals, 2012,33 (1):117-122;Wang Xiaodong, Zhang Yitang, the kind group vertical fractures well capacity of Liu and flow conductivity Optimizing research [J] Petroleum finance, 2004,31 (6):78-81), this with tight gas reservoir horizontal well in hydraulic fracturing process In, because the non-constant flow conductivity of length direction is stitched on the edge that the factors such as laying concentration difference, the fracturing fluid residue of proppant are caused There is marked difference in distribution.
It is of the invention by tight gas in order to realize that the fracture condudtiviy of tight gas reservoir pressure break horizontal well under real conditions optimizes The gas flowing for hiding pressure break horizontal well is divided into Reservoir Seepage and crack in the processes that flowing two is coupled, it is considered to which gas is along fracture surface On non-homogeneous inflow and crack in become the actual conditions of fracture condudtiviy, using room and time discrete technology, utilize Instantaneous point source function and barrier layer such as add at the general principle, have initially set up the unstable state Production rate mould of tight gas reservoir pressure break horizontal well Type;And further by contrasting the output condition under the distribution of different fracture condudtiviies, preferred tight gas reservoir pressure break horizontal well Fracture condudtiviy.
The content of the invention
It is an object of the invention to provide a kind of method of quantitative optimization tight gas reservoir pressure break horizontal well fracture condudtiviy, It can overcome the deficiencies in the prior art using this method, effectively solve pressure break horizontal well along the seam non-constant fracture guide energy of length direction The optimization problem of power, so as to provide rational basis for the optimization design of reservoir reconstruction, improves reservoir reconstruction effect.
A kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy, is mainly included the following steps that:
1) reservoir, fluid properties, horizontal well pit shaft basic parameter are collected;
2) pressure break horizontal well crack basic parameter is collected;
3) crack of pressure break horizontal well is averagely divided into the congruence of equal length along seam length direction;
4) the Reservoir Seepage model of tight gas reservoir pressure break horizontal well Fracture System is set up;
5) Pressure Drop Model that gas flows in crack is set up;
6) flow model of the coupled gas in Reservoir Seepage and crack is set up, the production of tight gas reservoir pressure break horizontal well is formed Measure computation model;
7) fracture condudtiviy of tight gas reservoir pressure break horizontal well is optimized.
In the present invention, the step 1) in collect reservoir, fluid properties, horizontal well pit shaft basic parameter, specifically include Original strata stress direction, reservoir thickness, porosity, permeability, gas viscosity, gas critical pressure, deviation factor for gas, reservoir temperature Degree, gas critical-temperature, horizontal well orientation, Horizontal Well tube length degree.
In the present invention, the step 2) in collect pressure break horizontal well crack basic parameter and include:Hydraulic fracture orientation, split Stitch bar number, fracture interval, crack location, fracture length.
In the present invention, the step 3) in the crack of pressure break horizontal well along seam length direction is averagely divided into equal length The congruence.
In the present invention, the step 4) set up the Reservoir Seepage mould for calculating tight gas reservoir pressure break horizontal well Fracture System Type, including following basic step:
(1) a bite pressure break level in physical model such as Fig. 1 of pressure break horizontal well, the infinitely great reservoir of upper top bottom closing is set up Well, the fracturing fracture that horizontal well is completely penetrated reservoir is divided into some sections, and reservoir gas flows first to crack and flowed again by crack To horizontal wellbore.Horizontal well radius is rw, length is L, and well is (x in reservoir center position coordinates0,y1,z0)-(x0,y2,z0), And it is parallel with y-axis.Reservoir homogeneous isotropism, thickness is that h, porosity φ, permeability K are constant;Horizontal wellbore is infinitely to lead Stream, reservoir initial pressure is constant pi.Because fracturing fracture runs through oil reservoir, therefore infinitely great formation breakdown horizontal well whole system Flowing can be reduced to the Radial Flow in plane oil reservoir, crack can be reduced to a congruence, and (Li Jun poem pressure breaks horizontal well is dynamically point Analysis research [D] Beijing:China University of Geosciences, 2005).
(2) for the ease of solving, kth crack single-blade is divided into ns sections, is Δ x per segment lengthfk, Δ xfk=xfk/ns; Each congruence can be processed into the production of a bite straight well and consider (Fig. 2).Yield is continually changing straight well in process of production, but if Time interval is obtained into very little, it is definite value that can be approximately considered the yield within this time, by the whole production rank of pressure break horizontal well Section t is divided into equal time interval Δ t, the t=m Δ t of m spacing.
Being located on kth crack has a congruence i (yield is qfk,i) it is located at M (xfk,i,yfk) place, have one on the k+1 of crack Point of observation O (xfk+1,j,yfk+1), in the case where considering fluid volume coefficient, congruence i is produced after production time Δ t in O points Drop of pressure be (Fig. 2):
In formula:piFor original formation pressure, MPa;pfk+1,jFor the pressure at the crack jth congruence sides of fracture of kth+1, MPa;qfk,iFor the yield of i-th congruence on kth crack, m3/s;μ is viscosity of crude, mPa.s;B is volume factor, zero dimension; (xfk,i, yfk) for the coordinate (m, m) of i-th congruence on kth crack;(xfk+1,j, yfk+1) it is the jth congruence on the crack of kth+1 Coordinate (m, m);K is reservoir permeability, 10‐3μm2;H is reservoir thickness, m;η is piezometric conductivity, μm2/(mPa·s×MPa-1), η =K/ μ c φ;C is stratum system compressibility, MPa‐1;φ is reservoir porosity, zero dimension;Δ t is production time, s;K is to split Loop bonding number;I, j number for the crack congruence.
(3) when containing multiple continuous congruences on the k of crack, each congruence of k cracks can obtain according to principle of stacking simultaneously raw The drop of pressure produced during production at point of observation O;After the same method, pressure break horizontal well formation N cracks can be also obtained, always In the drop of pressure of O points the common N × 2ns crack congruence is produced simultaneously when:
In formula:N is pressure break horizontal well Number of Fractures, bar;Ns is crack single-blade congruence number.
(4) when reservoir fluid is gas, according to the definition of pressure function and actual gas state equation, and by underground yield (Ning Zhengfu, Han Shugang, Cheng Linsong wait Low permeable oil and gas reservoirs pressure break horizontal wells to the gas yield for being converted into the case of ground standard Productivity [J] petroleum journals, 2002,23 (2):68-71), formula (2) can be written as:
In formula:μgFor gas viscosity, mPas;pscFor gas critical pressure, MPa;Z is deviation factor for gas, zero dimension; T is reservoir temperature, DEG C;TscFor gas critical-temperature, DEG C.
Formula (3) is exactly all slits congruence while the drop of pressure equation produced after production time Δ t in O points, the equation Consider and interfere effect between each congruence in crack and crack.
In the present invention, the step 5) calculate the Pressure Drop Model that gas flows in crack, including following basic step:
(1) consider that hydraulic fracture width becomes narrow gradually along seam length direction, by crack section processing in echelon, kth+1 is split Seam root width is wfk+1,max, end portion width is wfk+1,min, as shown in Figure 3.For the ease of calculating, each congruence is handled Rectangular to be calculated, line taking converges central point width for congruence width.Using the intersection point of crack and horizontal wellbore as origin, apart from water Horizontal well cylinder xjThe width of j-th of congruence at place is:
wfk+1,j=wfk+1,min+(wfk+1,max-wfk+1,min)×(xfk+1-xfk+1,j)/xfk+1 (4)
The average crack width w of j-th of congruencefk+1,averjFor:
wfk+1,averj=(wfk+1,j+wfk+1,j-1)/2 (5)
In formula:wfk+1,jFor the fracture width of the crack jth congruence of kth+1, m;wfk+1,maxIt is wide for the crack of kth+1 root Degree, m;wfk+1,minFor the crack end portion width of kth+1, m;xfk+1Long, the m for the crack single slot of kth+1;xfk+1,jFor kth+1 The crack jth congruence is in the x-direction apart from the distance of pit shaft, m;wfk+1,averjFor the crack jth congruence mean breadth of kth+1, m.
(2) consider that gas meets linear flow in crack, congruence j (points on the crack of kth+1 are calculated according to Darcy's law Ofk+1,j) pressure loss that flow to horizontal wellbore is:
In formula:pfk+1,0For the crack of kth+1 and the pressure of horizontal wellbore intersection, MPa;kfk+1For oozing for the crack of kth+1 Saturating rate, 10‐3μm2
In the present invention, the step 6) calculate the pressure break horizontal well production of tight gas reservoir, including following basic step:
(1) coupling model flowed in Reservoir Seepage-crack is set up
The process that gas flows to horizontal wellbore from reservoir can be divided into two processes of flowing, formula in gas reservoir seepage flow and crack (3), formula (6) has had been built up corresponding descriptive equation;Flowed in view of gas in storage intrastratal flow and fluid in crack When pressure at the wall of crack it is equal, pressure continuity equation can be set up;As it is assumed that pit shaft infinite fluid diversion, each crack exists With horizontal wellbore intersection O0Pressure it is equal, stable bottom hole pressure (pwf) production when boundary condition be:
pfk+1,0=pwf (7)
In formula:pwfFor horizontal wellbore flowing bottomhole pressure (FBHP), MPa;
Simultaneous equations (3), (6), (7) arrange the flow equation for obtaining the crack of kth+1 jth congruence gas reservoir-crack coupling:
Variable is the yield of each congruence in formula (8), thus establishes pressure break horizontal well crack and is coupled with gas reservoir Transient Flow mathematical modeling.Due to bottom pressure pwf, it is known that the line that the N × 2ns linear equation set up according to formula (8) is constituted Property equation group, can try to achieve the yield of each congruence.
(2) the transient state yield of tight gas reservoir pressure break horizontal well is calculated
Due to considering that gas is only flowed into horizontal wellbore by crack, so the total output of pressure break horizontal well is a frac water The sum of all congruence yield on horizontal well crack:
The yield of tight gas reservoir pressure break horizontal well can be tried to achieve according to formula (9).
(3) the unstable state yield of tight gas reservoir pressure break horizontal well is calculated
Equation (8), (9) establish the computational methods of each congruence and pressure break horizontal well production under the Δ t times;For each line For remittance, because yield can change with the production time, at this moment need to try to achieve the value of any time using time superposition.If Time step is Δ t, as t=m Δ t (m-th of period ending), then can obtain (honest and clean to the crack jth point sink of kth+1 Training celebrating, Cheng Linsong, Cao Renyi, etc. unsteady Model [J] calculating that low-permeability oil deposit pressure break horizontal well pit shafts are coupled with oil reservoir Physics, 2010,27 (2):203‐210):
In formula:
As producing the Production rate of Δ t time spans, constrained according to bottom pressure constraint and total output, composition closing Solving equations, since the 1st time Δ t step-length, circulation is solved, until calculating the Production rate under m-th of Δ t step-length As a result.
In the present invention, the step 7) optimization tight gas reservoir pressure break horizontal well fracture condudtiviy, including following base This step:
(1) different fracture condudtiviy distribution schemes are set along seam length direction;
(2) the crack output distribotion situation under different fracture condudtiviy distributions is calculated;
(3) cumulative production that the pressure break horizontal well under different fracture condudtiviy distributions produces 360 days is calculated, with reference to accumulative Yield result simultaneously considers engineering technology condition, preferably optimal fracture condudtiviy distribution scheme;
Compared with prior art, beneficial effects of the present invention:
Technical scheme can realize the fracture condudtiviy quantitative optimization to tight gas reservoir pressure break horizontal well, profit It can be realized fine and close according to underlying parameters such as the physical parameter of reservoir, fluid properties, horizontal wellbore and fracturing fractures with this method Quantitative optimization of the gas reservoir pressure break horizontal well fracture condudtiviy along seam length direction.It can only be realized in the prior art pair so as to overcome Average crack flow conductivity optimization of the fracture condudtiviy along whole seam length direction, improves the transformation of tight gas reservoir fractured horizontal well Validity and effect.
Brief description of the drawings
Fig. 1 is the infinitely great gas reservoir frac water horizontal well physical model schematic diagram of upper top bottom closing;
Fig. 2 is that two dimensional surface crack is segmented and discretization schematic diagram;
Fig. 3 is gas flow schematic diagram in gas reservoir pressure break horizontal well single-blade crack;
Fig. 4 is yield spatial distribution map when the 1st crack is produced 1 day under different flow conductivities are distributed;
Fig. 5 is the cumulative production comparison diagram of production 360 days under different fracture condudtiviy distribution schemes;
Fig. 6 is the fracture condudtiviy optimum results of example pressure break horizontal well.
Embodiment
Illustrate further below using how the present invention specifically realizes that a certain tight gas reservoir pressure break horizontal well crack is led The Optimization Design of stream ability, it is specific as follows:
It is a certain to treat that fracturing stratum reservoir, fluid and crack basic parameter are as follows:
Certain tight gas reservoir horizontal well buried depth 2800m, horizontal section length is 500m, and effective thickness is 25m, and porosity is 12.0%, permeability is 0.75 × 10-3μm2, gas-bearing formation temperature is 68 DEG C, and gas viscosity is 0.035mPas, deviation factor for gas For 0.89, gas critical pressure 4.64MPa, 88 DEG C of gas critical-temperature;Stratum system compressibility is 3.5 × 10-4MPa-1, Strata pressure is that 30MPa, flowing bottomhole pressure (FBHP) are 25MPa.Uniformly press off to form 4 cracks, fracture length along horizontal well pit shaft 75.0m, fracture permeabgility is 200 μm2
1) fracture condudtiviy distribution scheme is designed
Under the conditions of the fracture condudtiviy of pressure break horizontal well is fracture width after pressure break and reservoir buried depth closure stress The product of proppant permeability (Li Yingchuan is edited, petroleum production engineering [J] petroleum industry publishing house, 2009).Under reservoir buried depth Proppant permeability determines that therefore the key that fracture condudtiviy optimizes is exactly that hydraulic fracture is wide by selected proppant type The optimization of degree.
For the ease of optimizing the fracture width of tight gas reservoir pressure break horizontal well, it is assumed that fracture width is linear gradually along seam length direction Become:The fracture width of (crack root) is maximum i.e. at pressure break horizontal well pit shaft, tapered into along seam length direction fracture width, it is assumed that 4 kinds of fracture width change programmes, scheme I:wmax=wmin=2.75mm;Scheme II:wmax=3mm, wmin=2.5mm;Scheme III:wmax=4mm, wmin=1.5mm;Scheme IV:wmax=5mm, wmin=0.5mm.
2) the 1st crack output distribotion under different fracture condudtiviy distribution schemes is calculated
Crack output distribotion curve when Fig. 4 is the crack of pressure break horizontal well the 1st production 1d.As can be seen that in crack root There is local peaking in portion's yield, is caused because pressure is lower in root, crack;From contrast fracture condudtiviy Distribution scheme I, II, III, IV crack volume analysis can be seen that with the increase of crack root width, the yield peak value of root Reduction;" double U " types distribution characteristics show as typical non-homogeneous aerogenesis feature, gas for output distribotion presentation on whole crack Mainly pass through crack root and end output;With the increase of crack root width, interlude yield increase in crack is intended to along whole The feature that individual crack section is uniformly produced, is conducive to well production of improving the standard.
3) cumulative production that the pressure break horizontal well under different fracture condudtiviy distributions produces 360 days is calculated
Fig. 5 is the cumulative production comparison diagram that pressure break horizontal well produces 360d along under seam 4 kinds of width distribution schemes of length direction.Can To find out, under the conditions of the average slit width (being 2.75mm) of identical, due to the distributional difference of fracture width, 360d is caused to add up There is larger difference (scheme I in yield:1534.50×104m3, scheme IV:1564.06×104m3), this shows can be by excellent Change slit width distribution to improve pressure break horizontal well production.From the point of view of influence of the changes of slit length to yield, with root and end portion width Difference increases, and pressure break horizontal well increases;But after difference increase to a certain extent, yield increasing degree reduces, and exists preferred Fracture width is distributed.In this example, it is contemplated that during long-term production, be crushed due to proppant, the factor such as fracturing fluid residue Crack effective width is caused to narrow, preferably crack root width is that 4mm, fracture tip width are 1.5mm, that is to say crack root Flow conductivity is 200 μm2× 4mm=80 μm2.cm, fracture tip width is 200 μm2× 1.5mm=30 μm2.cm。
Fig. 6 is the fracture condudtiviy optimum results of example well.
4) 10 mouthfuls of water horizontal wells have carried out field conduct to the optimisation technique at the scene.Present invention optimization tight gas reservoir pressure break level The method of well fracture condudtiviy makes crack inner support agent distribution more be intended to the laying of optimization, obtains preferably stable life Production capacity power:3.52 × l0 of average temperature yield is obtained after 10 mouthfuls of well transformations4m3/ d, higher than 3.30 × l0 before optimization4m3/ d, increases Produce effect obvious.Illustrate the applicability and reliability of the present invention.

Claims (4)

1. a kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy, is mainly included the following steps that:
1) reservoir, fluid properties, horizontal well pit shaft basic parameter are collected;
2) pressure break horizontal well crack basic parameter is collected;
3) crack of pressure break horizontal well is averagely divided into the congruence of equal length along seam length direction;
4) the Reservoir Seepage model of tight gas reservoir pressure break horizontal well Fracture System is set up;
5) Pressure Drop Model that gas flows in crack is set up;
6) flow model of the coupled gas in Reservoir Seepage and crack is set up, the yield meter of tight gas reservoir pressure break horizontal well is formed Calculate model;
7) fracture condudtiviy of tight gas reservoir pressure break horizontal well is optimized;
The step 1) in collect reservoir, fluid properties, horizontal well pit shaft basic parameter, specifically include Original strata stress direction, storage Thickness degree, porosity, permeability, gas viscosity, gas critical pressure, deviation factor for gas, reservoir temperature, gas stagnation temperature Degree, horizontal well orientation, Horizontal Well tube length degree;
The step 2) in collect pressure break horizontal well crack basic parameter and include:Between hydraulic fracture orientation, Number of Fractures, crack Away from, crack location, fracture length;
The step 3) in the crack of pressure break horizontal well is averagely divided into the congruence of equal length along seam length direction;
The step 4) set up the Reservoir Seepage model for calculating tight gas reservoir pressure break horizontal well Fracture System, including following basic step Suddenly:
(1) a bite pressure break horizontal well, horizontal well in the physical model of pressure break horizontal well, the infinitely great reservoir of upper top bottom closing are set up The fracturing fracture for being completely penetrated reservoir is divided into some sections, and reservoir gas flows first to crack and flows to horizontal well by crack again Cylinder;Horizontal well radius is rw, length is L, and well is (x in reservoir center position coordinates0,y1,z0)-(x0,y2,z0), and it is flat with y-axis OK;Reservoir homogeneous isotropism, thickness is that h, porosity φ, permeability K are constant;At the beginning of horizontal wellbore is infinite fluid diversion, reservoir Beginning pressure is constant pi;Because fracturing fracture runs through oil reservoir, therefore the flowing of infinitely great formation breakdown horizontal well whole system can The Radial Flow in plane oil reservoir is reduced to, crack can be reduced to a congruence;
(2) for the ease of solving, kth crack single-blade is divided into ns sections, is Δ x per segment lengthfk, Δ xfk=xfk/ns;It is each The individual congruence can be processed into the production consideration of a bite straight well;Yield is continually changing straight well in process of production, but if by between the time Every obtaining very little, it is definite value that can be approximately considered the yield within this time, and the whole production phase t of pressure break horizontal well is divided into Equal time interval Δ t, the t=m Δ t of m spacing;
Being located on kth crack has a congruence i (yield is qfk,i) it is located at M (xfk,i,yfk) place, there is a point of observation on the k+1 of crack O(xfk+1,j,yfk+1), in the case where considering fluid volume coefficient, the pressure that congruence i is produced after production time Δ t in O points Land and be:
In formula:piFor original formation pressure, MPa;pfk+1,jFor the pressure at the crack jth congruence sides of fracture of kth+1, MPa;qfk,i For the yield of i-th congruence on kth crack, m3/s;μ is viscosity of crude, mPa.s;B is volume factor, zero dimension;(xfk,i, yfk) for the coordinate (m, m) of i-th congruence on kth crack;(xfk+1,j, yfk+1) for the coordinate of the jth congruence on the crack of kth+1 (m,m);K is reservoir permeability, 10-3μm2;H is reservoir thickness, m;η is piezometric conductivity, μm2/(mPa·s×MPa-1), η=K/ μcφ;C is stratum system compressibility, MPa-1;φ is reservoir porosity, zero dimension;Δ t is production time, s;K compiles for crack Number;I, j number for the crack congruence;
(3) when containing multiple continuous congruences on the k of crack, each congruence of k cracks can obtain while when producing according to principle of stacking The drop of pressure produced at point of observation O;After the same method, pressure break horizontal well formation N cracks can be also obtained, altogether N In the drop of pressure of O points the × 2ns crack congruence is produced simultaneously when:
In formula:N is pressure break horizontal well Number of Fractures, bar;Ns is crack single-blade congruence number;
(4) when reservoir fluid is gas, changed according to the definition of pressure function and actual gas state equation, and by underground yield Into the gas yield in the case of ground standard, formula (2) can be written as:
In formula:μgFor gas viscosity, mPas;pscFor gas critical pressure, MPa;Z is deviation factor for gas, zero dimension;T is Reservoir temperature, DEG C;TscFor gas critical-temperature, DEG C;
Formula (3) is exactly all slits congruence while the drop of pressure equation produced after production time Δ t in O points, the equation considers Effect is interfered between crack and crack each congruence.
2. the method as described in claim 1, it is characterised in that the step 5) calculate the pressure drop mould that gas flows in crack Type, including following basic step:
(1) consider that hydraulic fracture width becomes narrow gradually along seam length direction, by crack section processing in echelon, the crack root of kth+1 Portion's width is wfk+1,max, end portion width is wfk+1,min;For the ease of calculate, by each congruence processing it is rectangular calculated, Line taking converges central point width for congruence width;Using the intersection point of crack and horizontal wellbore as origin, apart from horizontal wellbore xjThe jth at place The width of the individual congruence is:
wfk+1,j=wfk+1,min+(wfk+1,max-wfk+1,min)×(xfk+1-xfk+1,j)/xfk+1 (4)
The average crack width w of j-th of congruencefk+1,averjFor:
wfk+1,aver j=(wfk+1,j+wfk+1,j-1)/2 (5)
In formula:wfk+1,jFor the fracture width of the crack jth congruence of kth+1, m;wfk+1,maxFor the crack of kth+1 root width, m; wfk+1,minFor the crack end portion width of kth+1, m;xfk+1Long, the m for the crack single slot of kth+1;xfk+1,jFor the crack of kth+1 The jth congruence is in the x-direction apart from the distance of pit shaft, m;wfk+1,averjFor the crack jth congruence mean breadth of kth+1, m;
(2) consider that gas meets linear flow in crack, congruence j (point O on the crack of kth+1 are calculated according to Darcy's lawfk+1,j) The pressure loss for flowing to horizontal wellbore is:
In formula:pfk+1,0For the crack of kth+1 and the pressure of horizontal wellbore intersection, MPa;kfk+1For the permeability in the crack of kth+1, 10-3μm2
3. the method as described in claim 1, it is characterised in that the step 6) calculate the pressure break horizontal well production of tight gas reservoir Amount, including following basic step:
(1) coupling model flowed in Reservoir Seepage-crack is set up
The process that gas flows to horizontal wellbore from reservoir can be divided into two processes of flowing, formula (3), formula in gas reservoir seepage flow and crack (6) corresponding descriptive equation has been had been built up;Split in view of gas when storing up intrastratal flow and fluid flows in crack Pressure at slotted wall face is equal, can set up pressure continuity equation;As it is assumed that pit shaft infinite fluid diversion, each crack with level Pit shaft intersection O0Pressure it is equal, stable bottom hole pressure (pwf) production when boundary condition be:
pfk+1,0=pwf (7)
In formula:pwfFor horizontal wellbore flowing bottomhole pressure (FBHP), MPa;
Simultaneous equations (3), (6), (7) arrange the flow equation for obtaining the crack of kth+1 jth congruence gas reservoir-crack coupling:
Variable is the yield of each congruence in formula (8), thus establishes the wink that pressure break horizontal well crack is coupled with gas reservoir State seepage experiment;Due to bottom pressure pwf, it is known that the linear side that the N × 2ns linear equation set up according to formula (8) is constituted Journey group, can try to achieve the yield of each congruence;
(2) the transient state yield of tight gas reservoir pressure break horizontal well is calculated
Due to considering that gas is only flowed into horizontal wellbore by crack, so the total output of pressure break horizontal well is a pressure break horizontal well The sum of all congruence yield on crack:
The yield of tight gas reservoir pressure break horizontal well can be tried to achieve according to formula (9);
(3) the unstable state yield of tight gas reservoir pressure break horizontal well is calculated
Equation (8), (9) establish the computational methods of each congruence and pressure break horizontal well production under the Δ t times;For each congruence Speech, because yield can change with the production time, at this moment needs to try to achieve the value of any time using time superposition;If the time Step-length is Δ t, as t=m Δ t (m-th of period ending), then the crack jth point sink of kth+1 can be obtained:
In formula:
As producing the Production rate of Δ t time spans, constrained according to bottom pressure constraint and total output, constitute Closure equation Group is solved, since the 1st time Δ t step-length, and circulation is solved, until calculating the Production rate knot under m-th of Δ t step-length Really.
4. the method as described in claim 1, it is characterised in that the step 7) crack of optimization tight gas reservoir pressure break horizontal well Flow conductivity, including following basic step:
(1) different fracture condudtiviy distribution schemes are set along seam length direction;
(2) the crack output distribotion situation under different fracture condudtiviy distributions is calculated;
(3) cumulative production that the pressure break horizontal well under different fracture condudtiviy distributions produces 360 days is calculated, with reference to cumulative production As a result and engineering technology condition is considered, preferably optimal fracture condudtiviy distribution scheme.
CN201510001929.0A 2015-01-04 2015-01-04 A kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy Active CN104594872B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510001929.0A CN104594872B (en) 2015-01-04 2015-01-04 A kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510001929.0A CN104594872B (en) 2015-01-04 2015-01-04 A kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy

Publications (2)

Publication Number Publication Date
CN104594872A CN104594872A (en) 2015-05-06
CN104594872B true CN104594872B (en) 2017-08-15

Family

ID=53120881

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510001929.0A Active CN104594872B (en) 2015-01-04 2015-01-04 A kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy

Country Status (1)

Country Link
CN (1) CN104594872B (en)

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105089595B (en) * 2015-05-27 2019-04-12 中国石油天然气股份有限公司 Oil reservoir numerical simulation method and device under horizontal fracturing fracture diversion action
CN104895550B (en) * 2015-06-04 2018-03-13 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 A kind of tight gas pressure break horizontal well numerical well testing model establishes method for solving
WO2016192077A1 (en) * 2015-06-04 2016-12-08 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 Method for establishing and solving numerical well-testing model of horizontal well for tight gas hydraulic fracturing
CN106321051B (en) * 2015-07-01 2019-02-15 中国石油化工股份有限公司 A method of for optimizing multistage pressure break horizontal well network fracture parameter
CN106484930A (en) * 2015-08-26 2017-03-08 中国石油化工股份有限公司 For determining the method and system of heterogeneous reservoir labyrinth well production
CN105422071B (en) * 2015-12-07 2018-12-11 西南石油大学 Evaluate the rational method of low-permeable heterogeneous gas reservoir fracture parameters of fractured horizontal wells
CN105422070B (en) * 2015-12-07 2018-11-23 西南石油大学 Optimize the method for the heterogeneous gas reservoir pressure break horizontal well crack location of extra-low permeability
CN105350960B (en) * 2015-12-07 2018-11-23 西南石油大学 The method for determining low-permeable heterogeneous gas reservoir fracture parameters of fractured horizontal wells
CN105484741A (en) * 2015-12-07 2016-04-13 西南石油大学 Prediction method for yield of low-permeability, heterogeneous and stress-sensitive reservoir fractured horizontal well
CN105840187B (en) * 2016-06-03 2018-12-25 陕西延长石油(集团)有限责任公司研究院 Compactness reservoir horizontal well staged fracturing Productivity
CN106337677B (en) * 2016-10-17 2019-03-26 长江大学 Shale gas pressure-break net air water two-phase flow conductivity test macro and test method
CN106650100B (en) * 2016-12-23 2020-01-10 西南石油大学 Alternate volume fracturing method for horizontal well of experimental shale reservoir
CN108572401B (en) * 2017-03-08 2020-04-03 中国石油化工股份有限公司 Construction method of fracture-cavity combined model and method for detecting deformation of reservoir fracture-cavity
CN107832515B (en) * 2017-11-01 2020-08-28 中国石油大学(北京) Method and device for simulating coupling of oil reservoir and shaft
CN110318742B (en) * 2018-03-30 2022-07-15 中国石油化工股份有限公司 Method and system for determining fracture closure length based on fractured well production data
CN109033677B (en) * 2018-08-09 2022-05-03 西南石油大学 Fracture conductivity optimization method for fractured and acidized well
CN109025942B (en) * 2018-08-09 2021-08-17 西南石油大学 Yield calculation method for irregular and multi-crack fracturing of tight gas reservoir inclined shaft
CN108979612B (en) * 2018-08-09 2020-07-07 西南石油大学 Method for optimizing flow capacity of complex fracture of fracturing and acidizing of compact oil and gas reservoir
CN109033674B (en) * 2018-08-09 2022-06-10 西南石油大学 Method for optimizing fracture parameters of target fracturing acidizing well
CN111400853B (en) * 2019-01-03 2022-10-04 中国石油天然气股份有限公司 Method and device for predicting unsteady state capacity of closed boundary fractured horizontal well
CN110805436B (en) * 2019-10-09 2021-12-03 中国石油大学(北京) Single-segment fracture liquid production contribution rate evaluation method and equipment based on pressure drop data
CN111079341B (en) * 2020-01-19 2021-10-01 西安石油大学 Intelligent well completion and oil reservoir unsteady state coupling method based on iterative algorithm
GB2595062B (en) * 2020-04-24 2022-05-04 Univ China Mining Experimental method and system for simulating evolution of reservoir fracture stress field
CN112081583B (en) * 2020-09-25 2021-08-17 西南石油大学 Unconventional storage layer seam network fracturing multi-scale support yield calculation method and device
CN115828651B (en) * 2023-02-24 2023-05-05 中国石油大学(华东) Method, system, equipment and medium for determining reasonable diversion capability of hydraulic fracturing fracture
CN115994500B (en) * 2023-03-22 2023-06-02 北京科技大学 Method and system for evaluating dynamic change of fracture conductivity of shale gas well
CN116877067B (en) * 2023-07-18 2024-03-12 重庆地质矿产研究院 Method for predicting hydraulic fracturing generated cracks and swept area fluid pressure
CN117287150B (en) * 2023-08-31 2024-04-19 中国地质大学(北京) Method, device, terminal and storage medium for acquiring economic recoverable resource amount of coalbed methane
CN117307152B (en) * 2023-11-28 2024-02-09 四川省能源地质调查研究所 Shale gas well fracturing fracture conductivity dynamic evaluation method and device

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5515923A (en) * 1994-08-26 1996-05-14 Loree; Dwight N. Oil and gas well productivity
CN103256035A (en) * 2013-05-22 2013-08-21 中国石化集团华北石油局 Fracturing fracture geological design method of horizontal well of dense gas field
CN103266881A (en) * 2013-05-22 2013-08-28 中国石化集团华北石油局 Method for predicting yield of compact hypotonic gas field multistage fracturing horizontal well
CN103527163A (en) * 2013-09-24 2014-01-22 西南石油大学 Tight reservoir horizontal well volume fracturing process
CN103577886A (en) * 2012-08-06 2014-02-12 中国石油化工股份有限公司 Staged fracturing yield prediction method of low-permeability gas reservoir horizontal well

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9378310B2 (en) * 2011-10-13 2016-06-28 Los Alamos National Security, Llc Material point method modeling in oil and gas reservoirs

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5515923A (en) * 1994-08-26 1996-05-14 Loree; Dwight N. Oil and gas well productivity
CN103577886A (en) * 2012-08-06 2014-02-12 中国石油化工股份有限公司 Staged fracturing yield prediction method of low-permeability gas reservoir horizontal well
CN103256035A (en) * 2013-05-22 2013-08-21 中国石化集团华北石油局 Fracturing fracture geological design method of horizontal well of dense gas field
CN103266881A (en) * 2013-05-22 2013-08-28 中国石化集团华北石油局 Method for predicting yield of compact hypotonic gas field multistage fracturing horizontal well
CN103527163A (en) * 2013-09-24 2014-01-22 西南石油大学 Tight reservoir horizontal well volume fracturing process

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
低渗气藏压裂水平井渗流与井筒管流耦合模型;孙福街等;《西南石油学院学报》;20050228;第27卷(第1期);第32页-第36页 *

Also Published As

Publication number Publication date
CN104594872A (en) 2015-05-06

Similar Documents

Publication Publication Date Title
CN104594872B (en) A kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy
Jing et al. Mechanism simulation of oil displacement by imbibition in fractured reservoirs
CN106437674B (en) Imitative water injection of horizontal well well pattern adaptation method
CN103590824B (en) Capacity calculation method for compact gas reservoir horizontal well after multi-section fracturing modification
CN104989341B (en) Method for determining effective displacement injection-production well spacing of low-permeability oil reservoir
CN104632157A (en) Low permeability reservoir equilibrium displacement method
WO2021031543A1 (en) Method and apparatus for determining integrated development approach for shale and adjacent oil layers
CN105089612A (en) Determining method for distance of well-drain and length of pressure break of low penetration oil reservoir artificial fracture
CN102777157A (en) CO2 drive oil-gas-water separate well injecting oil reservoir mixing drive development method
CN107355206A (en) A kind of shale gas horizontal well refracturing temporarily blocks up critical pressure method of testing
CN109101773A (en) Fine and close grease horizontal well solid seam net cluster network pressure splits optimization method
CN106547930A (en) Consider the gas drainage radius computational methods of tight gas reservoir seepage flow mechanism
CN112541287A (en) Loose sandstone fracturing filling sand control production increase and profile control integrated design method
Zhao et al. Physical simulation of waterflooding development in large-scale fractured-vuggy reservoir considering filling characteristics
Ignatyev et al. Multistage hydraulic fracturing in horizontal wells as a method for the effective development of gas-condensate fields in the arctic region
Xue et al. Status and Prospect of Improved Oil Recovery Technology of High Water Cut Reservoirs
Pan et al. Combined technology of weak gel flooding assisting thermal huff and puff enhances oil recovery for offshore heavy oil field
Chen et al. Application of the unstructured grids in the numerical simulation of fractured horizontal wells in ultra-low permeability gas reservoirs
Hu et al. Quick assessment to ascertain technical rational well spacing density in artificial water flooding oilfield
CN106401549A (en) Method of building influence plate on horizontal well production effect by crack and well trajectory angle
CN111287741B (en) Rapid calculation method for permeability of compact oil reservoir volume fracturing transformation area
CN106777847A (en) Low-permeability oil deposit staged fracturing horizontal well unstable state productivity model
Guo et al. Visual experiments of bottom water and multi-well water and gas injection flooding for fault-controlled fractured-vuggy reservoirs
Cui et al. Improved oil recovery for the complex fault-block and multilayered reservoir with edge water
Fu et al. Method of Predicting the Location of Water Cresting for Horizontal Wells in a Water-Drive Reservoir for Early Prevention

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant