CN105422070B - Optimize the method for the heterogeneous gas reservoir pressure break horizontal well crack location of extra-low permeability - Google Patents
Optimize the method for the heterogeneous gas reservoir pressure break horizontal well crack location of extra-low permeability Download PDFInfo
- Publication number
- CN105422070B CN105422070B CN201510889494.8A CN201510889494A CN105422070B CN 105422070 B CN105422070 B CN 105422070B CN 201510889494 A CN201510889494 A CN 201510889494A CN 105422070 B CN105422070 B CN 105422070B
- Authority
- CN
- China
- Prior art keywords
- unit
- flow
- pressure
- crack
- permeable
- 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
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
Abstract
This application discloses a kind of methods for optimizing the heterogeneous gas reservoir pressure break horizontal well crack location of extra-low permeability, mainly include the following steps that：Collect the basic parameter of reservoir, fluid, horizontal wellbore；Heterogeneous gas reservoir is divided at least two vadose zones along the length direction of horizontal wellbore, vadose zone contains man-made fracture；Establish the permeable unit simulation model of each vadose zone；Establish the seepage simulation model of at least two permeable units；The seepage simulation model for coupling each permeable unit calculates the yield of heterogeneous gas reservoir pressure break horizontal well；The yield is compared, the heterogeneous highest crack location of the gas reservoir frac water amount of showing no increases in output of extra-low permeability is selected.The technical solution of the application has fully considered the influence of reservoir heterogeneity and starting pressure gradient, so that overcoming can only realize in the prior art to homogeneous gas reservoir, not consider to start the pressure break horizontal well crack location selection of pressure, to improve the validity and effect of the heterogeneous Horizontal Wells For Gas Reservoirs fracturing reform of extra-low permeability.
Description
Technical field
The application belongs to oil-gas field development field, specifically, being related to a kind of heterogeneous gas reservoir pressure break of optimization extra-low permeability
The method of horizontal well crack location.
Background technique
Ultralow permeability gas reservoir outstanding behaviours is reservoir permeability plane heterogeneity and there are the characteristics that starting pressure gradient,
By the way that yield can be significantly improved to heterogeneous Horizontal Wells For Gas Reservoirs staged fracturing.To ultralow permeability gas reservoir pressure break horizontal well crack
When position is preferred, it is necessary to comprehensively consider the influence of the factors such as permeability heterogeneity and starting pressure gradient.It is directed to pressure break at present
Horizontal well production prediction mainly has the methods of electrical analogue, conformal mapping, principle of stacking, the theory of complex functions.
When the above method does not consider that ultralow permeability gas reservoir heterogeneity and ultralow permeability gas reservoir fluid flow simultaneously
There are the combined influences of starting pressure gradient.
Summary of the invention
In view of this, technical problems to be solved in this application be not and meanwhile consider ultralow permeability gas reservoir heterogeneity with
And there are the combined influences of starting pressure gradient when the flowing of ultralow permeability gas reservoir fluid.
In order to solve the above-mentioned technical problem, this application discloses a kind of heterogeneous gas reservoir pressure break horizontal wells of optimization extra-low permeability
The method of crack location, mainly includes the following steps that：
1) basic parameter of reservoir, fluid, horizontal wellbore is collected；
2) heterogeneous gas reservoir is divided at least two vadose zones, the seepage flow along the length direction for being the horizontal wellbore
Band contains man-made fracture；
3) the permeable unit simulation model of each vadose zone is established；
4) the seepage simulation model of at least two permeable units is established；
5) simulation model for coupling the seepage flow of each permeable unit, calculates the production of heterogeneous gas reservoir pressure break horizontal well
Amount；
6) according to yield, optimize the crack location of the heterogeneous gas reservoir pressure break horizontal well of extra-low permeability.
Further, in the step 1), the basic parameter of the reservoir includes：Original strata stress direction, reservoir thickness, storage
Layer temperature, porosity, permeability；The basic parameter of the fluid includes：Gas viscosity, gas critical pressure, gas deviation because
Son, gas critical-temperature；The basic parameter of the horizontal wellbore includes：Horizontal well orientation, Horizontal Well tube length degree.
Further, the step 2) is specially：It, will according to the difference of reservoir permeability on horizontal wellbore length direction
The identical reservoir division of permeability is at the same vadose zone containing man-made fracture permeable belt.
Further, the basic flowing process of the permeable unit includes：Matrix linear flow, equivalent radius to stream and crack
Interior linear seepage flow；The resistance of each flow process in permeable unit is calculated by taking the 1st permeable unit in 1st article of carrying substrates as an example
Power.
Further, the step 3) specifically includes：
A, according to relational expression：The resistance of I area's linear flow is calculated, in formula：RuijFor
J-th of I territorial matrix linear flow resistance of permeable unit, unit MPa in i-th permeable belt2/((m3/d)·(mPa·s))；T is gas
Hide temperature, unit K；TscFor standard atmospheric temperature, unit K；pscFor standard atmospheric pressure, unit MPa；lijFor the i-th infiltration
J-th strip crack and upper crack the flowing distance, unit m in band；xfFor fracture half-length, unit m；kiFor the i-th infiltration tape base
Matter permeability, unit mD；H is reservoir thickness, unit m；
B, according to relational expression：Calculate II area's radial flow resistance；Formula
In：RnijIt is j-th of permeable unit equivalent radius in the i-th permeable belt to flow resistance power, unit MPa2/((m3/d)·(mPa·
s))；kfijFor fracture permeabgility, unit Dc；wfijFor fracture width, unit cm；
C, according to the resistance of I area's linear flow and II area's radial flow resistance, the permeable unit of the vadose zone is established
Simulation model.
Further, in the step a, the relational expressionIn, Ruij=ψuij-
ψdij,
According to the matrix linear flow flow and pseudopressure relationship under gas reservoir starting pressure gradient, formula is obtained：
It is defined according to pseudopressure, it is G that gas, which flows through the starting pressure consumed when distance d,ijD, by GijD is defined as opening
Dynamic pressure pB；Therefore following conversion can be done：
In formula：Ψuij、ΨdijThe outer pseudopressure of respectively I territorial matrix linear flow and interior pseudopressure, ΨBuij-ΨBdijFor fluid
Flowing through length is dijPressure difference, unit MPa are intended in the starting of I territorial matrix2/(mPa·s)；GijFor the starting pressure ladder of I territorial matrix
Degree, unit MPa/m；dijFor I territorial matrix linear flow length, unit m；qijIt is single for j-th strip crack flow in the i-th permeable belt
Position is m3/d。
Further, in the step b：
According in the same permeable unit, the pressure at the area I and the area II boundary position is equal, and the internal pressure of I area's linear flow is
For the external pressure of II area's radial flow, then radial flow flow is with pressure difference relational expression：
In formula：ΨwfijFor pseudopressure in II area's radial flow, unit MPa2/(mPa·s)；rwijIt is single for equivalent well radius
Position is m.
The radius r of the equivalent straight wellwIt is sought by equivalent radius model, specially：
By establishing crack yield and the equal relational expression of straight well yield to obtain between equivalent radius and fracture half-length
Relational expression establishes the yield for seeking equivalent diameter and crack separately below：
Straight well yield：Seek a bite straight well A in level pressure supply boundary band gas reservoir1Steady state production, according to mirror image
Principle, is translated into row's producing well and row's gas injection well yield stable state solves；
According to equivalent flow resistance method, a bite well production is obtained：
Limited fluid diversion crack yield：
In formula：Ψ is pseudopressure at x point in fractuer direction, unit MPa2/(mPa·s)。
Solving average pressure in crack is approximately pressure in crack, and fluid flows in crack from matrix meets linear flow rule
Rule, the outer pseudopressure of linear flow is Ψuij, matrix linear flow flow is found out, then crack yield is obtained by the conservation of mass：
The yield formula (8) for comparing the straight well yield formula (6) of equivalent and limited fluid diversion crack obtains equivalent well radius formula：
In formula, rwFor equivalent well radius, unit m；
Formula (8), which are substituted into equivalent radial flow resistance power near formula (5) available II area's equivalent well of permeable unit, is：
Further, the step 4) is specially：
It is basic according to meeting that pressure is continuous, flow is equal between different the flow domains after acquiring each area's flow resistance
Principle acquires fractured horizontal well's productivity formula using equivalent flow resistance method.
Further, the step 5) is specially：Fluid is first analyzed both sides permeable belt crack and is produced from both sides to intermediate flow
Energy equation analyzes intermediate cracking deliverability equation again, connects each resistance with equivalent flow resistance method between adjacent slits and obtainsThe system of linear equations of a equation composition；
Crack yield in the 1st permeable belt is calculated, formula is：
Crack yield in the 3rd permeable belt is calculated, formula is：
Since outermost two permeable units of the 2nd permeable belt are connected with 1,3 permeable belts respectively, then its linear flow resistance of I territorial matrix
Power is：
Assuming that the well row that intermediate distributary crack is equivalent for m crack in the 2nd permeable belt, well row is two flows of left and right
Meet, so N can be listed2+ 1 equation：
In formula：q2mThe flow of shunting well row, unit m are flowed to for the left side3/d；q2(m+1)Flow to shunting well row's for the right
Flow, unit m3/d；The actual flow for shunting crack is sum of the two, unit m3/d；
It is defined according to pseudopressure, the expression of the pressure difference of two squares can be used in the quasi- quasi- pressure difference of starting：
The pseudopressure function expression of any point pressure p can similarly be obtained：
According to the symmetry of gas reservoir, the total yield of pressure break horizontal well is：
In formula, Q is pressure break horizontal well total output, unit m3/d。
Further, the step 6) specifically includes：
A, the different crack location scheme of setting heterogeneous reservoir；
B, the yield of comparing calculation difference crack location scheme；
C, crack location is selected most preferably with post-fracturing yield.
Compared with prior art, the application can be obtained including following technical effect：
Fully consider the influence of reservoir heterogeneity and starting pressure gradient, it in the prior art can only be real to overcome
Now consider that the pressure break horizontal well crack location for starting pressure is selected to homogeneous gas reservoir, not, to improve the heterogeneous gas of extra-low permeability
Hide the validity and effect of fractured horizontal well transformation.
Certainly, any product for implementing the application must be not necessarily required to reach all the above technical effect simultaneously.
Detailed description of the invention
The drawings described herein are used to provide a further understanding of the present application, constitutes part of this application, this Shen
Illustrative embodiments and their description please are not constituted an undue limitation on the present application for explaining the application.In the accompanying drawings：
Fig. 1 is the heterogeneous reservoir pressure break horizontal well schematic diagram of the embodiment of the present application；
Fig. 2A is the practical equivalent schematic flowed of the permeable unit of the embodiment of the present application；
Fig. 2 B is the equivalent schematic of the equivalent flowing of permeable unit of the embodiment of the present application；
Fig. 3 A is the level pressure closed boundary mirror image schematic diagram of the embodiment of the present application；
Fig. 3 B is the circuit diagram after the level pressure closed boundary mirror image of the embodiment of the present application；
Fig. 4 is the gas reservoir pressure break horizontal well equivalent circuit diagram of the embodiment of the present application；
Fig. 5 is that the different permeabilities of the embodiment of the present application combine lower horizontal well production pattern；
Fig. 6 is the different crack location arrangement schematic diagrams of the embodiment of the present application；
1 indicates that constant-pressure boundary, 2 indicate that closed boundary, 3 indicate that horizontal wellbore, 4 indicate crack in Fig. 1；3 indicate water in Fig. 6
Horizontal well cylinder, 4 indicate crack.
Specific embodiment
Presently filed embodiment is described in detail below in conjunction with accompanying drawings and embodiments, thereby how the application is applied
Technological means solves technical problem and reaches the realization process of technical effect to fully understand and implement.
It is main to wrap this application discloses a kind of method for optimizing the heterogeneous gas reservoir pressure break horizontal well crack location of extra-low permeability
Include following steps：
1) basic parameter of reservoir, fluid, horizontal wellbore is collected；
2) heterogeneous gas reservoir is divided at least two vadose zones, the seepage flow along the length direction for being the horizontal wellbore
Band contains man-made fracture；
3) the permeable unit simulation model of each vadose zone is established；
4) the seepage simulation model of at least two permeable units is established；
5) simulation model for coupling the seepage flow of each permeable unit, calculates the production of heterogeneous gas reservoir pressure break horizontal well
Amount；
6) according to yield, optimize the crack location of the heterogeneous gas reservoir pressure break horizontal well of extra-low permeability.
Further, in the step 1), the basic parameter of the reservoir includes：Original strata stress direction, reservoir thickness, storage
Layer temperature, porosity, permeability；The basic parameter of the fluid includes：Gas viscosity, gas critical pressure, gas deviation because
Son, gas critical-temperature；The basic parameter of the horizontal wellbore includes：Horizontal well orientation, Horizontal Well tube length degree.
As shown in Figure 1, the step 2) is specially：According to the difference of reservoir permeability on horizontal wellbore length direction,
By the identical reservoir division of permeability at the same vadose zone containing man-made fracture permeable belt.
As shown in Figure 2 A and 2 B, the seepage flow of pressure break horizontal well can be regarded as by the unit group of many similar flow processes
At, referred to as permeable unit, each permeable unit contain the linear seepage flow of matrix linear flow, equivalent radius into stream and crack
Three basic flowing processes, the present embodiment calculate each in permeable unit by taking the 1st permeable unit in 1st article of carrying substrates as an example
The resistance of flow process；
The step 3) specifically includes：
A, the resistance of I area's linear flow is calculated：
Consider that the matrix linear flow flow under gas reservoir starting pressure gradient is with pseudopressure relational expression：
Integral term in above-mentioned (1) formulaIt is difficult to determine, pseudopressure definition can be used.Gas flows through distance d
When the starting pressure that consumes be GijD, by GijD is defined as starting pressure PB.Therefore following conversion can be done：
In formula：Ψuij、ΨdijThe outer pseudopressure of respectively I territorial matrix linear flow and interior pseudopressure, ΨBuij-ΨBuijFor fluid
Flowing through length is dijPressure difference, unit MPa are intended in the starting of I territorial matrix2/(mPa·s)；GijFor the starting pressure ladder of I territorial matrix
Degree, unit MPa/m；dijFor I territorial matrix linear flow length, unit m；qijFor flow, unit m3/d；xfFor crack half
It is long, unit m；kiFor the i-th permeable belt matrix permeability, unit mD；H is reservoir thickness, unit m；T is reservoir temperature,
Unit is K；TscFor standard atmospheric temperature, unit K；pscFor standard atmospheric pressure, unit MPa.
Since equivalent well is nearby there are equivalent radial punctured area, then I area linear flow length dijLess than fracture interval lij,
The radial flow perimeter that loses heart can be approximately fracture half-length, the line source that area approximation is crack of losing heart to interstitial rectangular area,
Then I area's linear flow length is：
So the resistance of I area's linear flow is：
In formula：RuijFor j-th of I territorial matrix linear flow resistance of permeable unit, unit MPa in the i-th permeable belt2/((m3/
d)·(mPa·s))；
B, II area's radial flow resistance is calculated：
According in the same permeable unit, the pressure at the area I and the area II boundary position is equal, and the internal pressure of I area's linear flow is
For the external pressure of II area's radial flow, then radial flow flow is with pressure difference relational expression：
In formula：ΨwfijFor pseudopressure in II area's radial flow, unit MPa2/(mPa·s)；rwijIt is single for equivalent well radius
Position is m；
The radius r of equivalent straight well is sought using equivalent radius modelw, basic thought be exactly by establish crack yield and
The equal relational expression of straight well yield is to obtain relational expression between equivalent radius and fracture half-length.It establishes separately below and seeks equivalent
The yield of diameter and crack.
As shown in Figure 3A and Figure 3B, straight well yield：In order to seek a bite straight well A in level pressure supply boundary band gas reservoir1's
Steady state production can be translated into row's producing well and row's gas injection well yield stable state is asked according to mirror image principle
Solution；
A bite well production is acquired using equivalent flow resistance method：
Limited fluid diversion crack yield：The Production rate in limited fluid diversion crack is described using following ODEs：
In formula：In formula：Ψ is pseudopressure at x point in fractuer direction, unit MPa2/(mPa·s)；kfijFor crack infiltration
Rate, unit Dc；wfijFor fracture width, unit cm.
Solving average pressure in crack is approximately pressure in crack, considers that fluid flows in crack from matrix and meets linearly
Flowing law, the outer pseudopressure of linear flow is Ψuij, matrix linear flow flow is found out, then crack yield is obtained by the conservation of mass：
The yield formula (8) for comparing the straight well yield formula (6) of equivalent and limited fluid diversion crack obtains equivalent well radius formula：
Formula (8), which are substituted into equivalent radial flow resistance power near formula (5) available II area's equivalent well of permeable unit, is：
In formula：RnijIt is j-th of permeable unit equivalent radius in the i-th permeable belt to flow resistance power, unit MPa2/((m3/
d)·(mPa·s))。
Further, the step 4) specifically includes：
Computational physics model：
As shown in figure 4, after acquiring each area's flow resistance, according to meeting between different the flow domains, pressure is continuous, flow
Equal basic principle acquires fractured horizontal well's productivity formula using equivalent flow resistance method；Here there are three not with gas reservoir
With permeability zones distribution for, N that the i-th permeable belt presses offiCrack ignores the influence of horizontal wellbore pressure drop, by difference
The combination distribution characteristics of permeable unit obtains gas reservoir equivalent circuit diagram；
Fluid first analyzes both sides permeable belt crack deliverability equation and analyzes intermediate cracking production capacity side again from both sides to intermediate flow
Journey connects each resistance with equivalent flow resistance method between adjacent slits and obtainsThe system of linear equations of a equation composition.
Further, the step 4) specifically includes：
Calculate yield：
Crack yield in the 1st permeable belt is calculated, formula is：
Crack yield in the 3rd permeable belt is calculated, formula is：
Since outermost two permeable units of the 2nd permeable belt are connected with 1,3 permeable belts respectively, then its linear flow resistance of I territorial matrix
Power is：
Assuming that the well row that intermediate distributary crack is equivalent for m crack in the 2nd permeable belt, well row is two flows of left and right
Meet, so N can be listed2+ 1 equation：
In formula：q2mThe flow of shunting well row, unit m are flowed to for the left side3/d；q2(m+1)Flow to shunting well row's for the right
Flow, unit m3/ d, the actual flow for shunting crack is sum of the two, unit m3/d；；
It is defined according to pseudopressure, the expression of the pressure difference of two squares can be used in the quasi- quasi- pressure difference of starting：
The pseudopressure function expression of any point pressure p can similarly be obtained：
According to the symmetry of gas reservoir, the total yield of pressure break horizontal well is：
In formula, Q is pressure break horizontal well total output, unit m3/d。
Further, the step 6) specifically includes：
A, the different crack location scheme of setting heterogeneous reservoir；
B, the yield of comparing calculation difference crack location scheme；
C, crack location is selected most preferably with post-fracturing yield.
Embodiment
The method of the heterogeneous gas reservoir pressure break horizontal well crack location of the present embodiment optimizing application extra-low permeability, it is specific as follows：
Certain heterogeneous Horizontal Wells For Gas Reservoirs underlying parameter is as follows：The heterogeneous a length of 800m of gas reservoir, width 120m, with a thickness of 15m,
Both sides terminal pressure be 25MPa, horizontal well cylinder pressure 22MPa, both sides permeable belt (the 1st and 3 permeable belt) length be 300m,
Matrix permeability is 2 0.6mD, equal pressure break equally distributed vertical fractures, and intermediate permeable belt (the 2nd permeable belt) is a length of
200m, matrix permeability 1.2mD, 1 vertical fracture, a length of 60m of all slits half, fracture condudtiviy 30D are pressed off
Cm, Natural Gas Viscosity 0.0215mPas, deviation factor 0.9218.
Table 1 is three permeable belt different substrates permeability combinations, in order to compare homogeneous gas reservoir and heterogeneous gas reservoir frac water
The difference of horizontal well Production rate result, while considering the combination combined II, combine the other two kinds of reservoir permeability parameters of III,
Its parameter is same as above.
Table 1
According to the relationship of matrix starting pressure gradient and permeability：
In formula：A is gas reservoir constant；N value is 1.
As shown in figure 5, existing, by A, successively value is 0,0.05,0.1,0.15,0.2, obtains permeating from different infiltration carrying substrates
Rate combines the combination of corresponding starting pressure gradient, in the daily output of different gas reservoir constant A fracturing horizontal wells.
As shown in Figure 5：As starting pressure gradient increases (increase of A value), same combination fracturing horizontal well production is gradually
It reduces；It is different that various combination reduces trend, and the downward trend for the combination (I) that both sides permeability is low, intermediate permeability is high is greater than
Low combination (III) the pressure break horizontal well production of high, the intermediate permeability of both sides permeability；Under same starting pressure gradient, both sides
Gas reservoir (I) the pressure break horizontal well low output that permeability is low, intermediate permeability is high low gas of high, intermediate permeability in both sides permeability
(III) pressure break horizontal well production is hidden, therefore for heterogeneous gas reservoir types I, according to the average weighted method of matrix permeability
Processing is homogeneous gas reservoir II, then predicts that production capacity is obviously higher；For heterogeneous gas reservoir types III, according to matrix permeability plus
The method processing of weight average is homogeneous gas reservoir II, then predicts that production capacity is obviously relatively low.
As shown in fig. 6, further compared the gas reservoir pressure break horizontal well 3 forms 5 cracks 4, in different 4 positions of crack
3 yield of pressure break horizontal well under arrangement.
Table 2 is the pressure break horizontal well production under the distribution of different crack locations, note:D1-2 indicate the 1st article and the 2nd crack it
Between distance；D2-3 indicates the distance between the 2nd article and the 3rd crack；D3-4 indicate the 3rd article between the 4th crack away from
From；D4-5 indicates the distance between the 4th article and the 5th crack.As shown in Table 2, the pressure break horizontal well production stitched using scheme I cloth
Highest, therefore preferred embodiment I carries out cloth seam.
Table 2
The application has fully considered the influence of reservoir heterogeneity and starting pressure gradient, to overcome in the prior art
It can only realize to homogeneous gas reservoir, not consider to start the pressure break horizontal well crack location selection of pressure, so that it is non-to improve extra-low permeability
The validity and effect of homogeneous gas reservoir fractured horizontal well transformation.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
Include, so that commodity or system including a series of elements not only include those elements, but also including not clear
The other element listed, or further include for this commodity or the intrinsic element of system.In the feelings not limited more
Under condition, the element that is limited by sentence "including a ...", it is not excluded that in the commodity or system for including the element also
There are other identical elements.
Above description shows and describes several preferred embodiments of the present application, but as previously described, it should be understood that the application
Be not limited to forms disclosed herein, should not be regarded as an exclusion of other examples, and can be used for various other combinations,
Modification and environment, and the above teachings or related fields of technology or knowledge can be passed through in application contemplated scope described herein
It is modified.And changes and modifications made by those skilled in the art do not depart from spirit and scope, then it all should be in this Shen
It please be in the protection scope of appended claims.
Claims (7)
1. a kind of method for optimizing the heterogeneous gas reservoir pressure break horizontal well crack location of extra-low permeability, which is characterized in that mainly include
Following steps：
1) basic parameter of reservoir, fluid, horizontal wellbore is collected；
2) heterogeneous gas reservoir is divided at least two vadose zones along the length direction of the horizontal wellbore, the vadose zone contains
Man-made fracture；
3) the seepage simulation model of the permeable unit of each vadose zone is established；
4) the seepage simulation model of at least two permeable units is established；
5) the seepage simulation model for coupling each permeable unit, calculates the yield of heterogeneous gas reservoir pressure break horizontal well；
6) according to yield, optimize the crack location of the heterogeneous gas reservoir pressure break horizontal well of extra-low permeability；
The step 3) specifically includes：
A, according to relational expression：The resistance of I area's linear flow is calculated, in formula：RuijIt is seeped for i-th
J-th of I territorial matrix linear flow resistance of permeable unit in band thoroughly, unit MPa2/((m3/d)·(mPa·s))；T is gas reservoir temperature
Degree, unit K；TscFor standard atmospheric temperature, unit K；pscFor standard atmospheric pressure, unit MPa；lijFor in the i-th permeable belt
J-th strip crack and upper crack the flowing distance, unit m；xfFor fracture half-length, unit m；kiIt is seeped for the i-th infiltration carrying substrates
Saturating rate, unit mD；H is reservoir thickness, unit m；
B, according to relational expression：Calculate II area's radial flow resistance；
In formula：RnijIt is j-th of permeable unit equivalent radius in the i-th permeable belt to flow resistance power, unit MPa2/((m3/d)·
(mPa·s))；kfijFor fracture permeabgility, unit Dc；wfijFor fracture width, unit cm；
C, according to the resistance of I area's linear flow and II area's radial flow resistance, the mould of the permeable unit of the vadose zone is established
Analog model.
2. the method as described in claim 1, which is characterized in that in the step 1), the basic parameter of the reservoir includes：It is former
Stress direction, reservoir thickness, reservoir temperature, porosity, permeability；The basic parameter of the fluid includes：Gas viscosity, gas
Body critical pressure, deviation factor for gas, gas critical-temperature；The basic parameter of the horizontal wellbore includes：Horizontal well orientation, water
Horizontal well pit shaft length.
3. method according to claim 2, which is characterized in that the step 2) is specially：According to reservoir in horizontal well tube length
The difference for spending permeability on direction, by the identical reservoir division of permeability at the same seepage flow containing man-made fracture permeable belt
Band.
4. method as claimed in claim 3, which is characterized in that the basic flowing process of the permeable unit includes：Matrix line
The linear seepage flow of property stream, equivalent radius into stream and crack；It is calculated by taking the 1st permeable unit in 1st article of carrying substrates as an example
The resistance of each flow process in permeable unit.
5. method as claimed in claim 4, which is characterized in that in the step a, the relational expressionIn, Ruij=ψuij-ψdij,
According to the matrix linear flow flow and pseudopressure relationship under gas reservoir starting pressure gradient, formula is obtained：
It is defined according to pseudopressure, it is G that gas, which flows through the starting pressure consumed when distance d,ijD, by GijD is defined as starting pressure
PB；Therefore following conversion can be done：
In formula：Ψuij、ΨdijThe outer pseudopressure of respectively I territorial matrix linear flow and interior pseudopressure, ΨBuij - Ψ Bdij For fluid
Flowing through length is dijI territorial matrix starting intend pressure difference, unit MPa2/(mPa·s)；GijFor the starting pressure ladder of I territorial matrix
Degree, unit MPa/m；dijFor I territorial matrix linear flow length, unit m；qijIt is single for j-th strip crack flow in the i-th permeable belt
Position is m3/d；xfFor fracture half-length, unit m；kiFor the i-th permeable belt matrix permeability, unit mD；H is reservoir thickness, single
Position is m；T is reservoir temperature, unit K；TscFor standard atmospheric temperature, unit K；pscFor standard atmospheric pressure, unit MPa.
6. method as claimed in claim 5, which is characterized in that in the step b：
According in the same permeable unit, the pressure at the area I and the area II boundary position is equal, and the internal pressure of I area's linear flow is II
The external pressure of area's radial flow, then radial flow flow be with pressure difference relational expression：
In formula：ΨwfijFor pseudopressure in II area's radial flow, unit MPa2/(mPa·s)；rwFor equivalent straight well radius, unit is
m；
The radius r of the equivalent straight wellwIt is sought by equivalent radius model, specially：
By establishing crack yield and the equal relational expression of straight well yield to obtain relationship between equivalent radius and fracture half-length
Formula establishes the yield for seeking equivalent diameter and crack separately below：
Straight well yield：Seek a bite straight well A in level pressure supply boundary band gas reservoir1Steady state production, according to mirror image principle,
It is translated into row's producing well and row's gas injection well yield stable state solves；
According to equivalent flow resistance method, a bite well production is obtained：
Limited fluid diversion crack yield：
In formula：Ψ is pseudopressure at x point in fractuer direction, unit MPa2/(mPa·s)；kfijFor fracture permeabgility, unit is
Dc；wfijFor fracture width, unit cm；
Solving average pressure in crack is approximately pressure in crack, and fluid flows in crack from matrix meets linear flow rule,
The outer pseudopressure of linear flow is Ψuij, matrix linear flow flow is found out, then crack yield is obtained by the conservation of mass：
The yield formula (8) for comparing the straight well yield formula (6) of equivalent and limited fluid diversion crack obtains equivalent well radius formula：
In formula, rwFor equivalent straight well radius, unit m；
Formula (8), which are substituted into equivalent radial flow resistance power near formula (5) available II area's equivalent well of permeable unit, is：
7. method as claimed in claim 6, which is characterized in that the step 4) is specially：
After acquiring each area's flow resistance, according to meeting the basic principle that pressure is continuous, flow is equal between different the flow domains,
Fractured horizontal well's productivity formula is acquired using equivalent flow resistance method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510889494.8A CN105422070B (en) | 2015-12-07 | 2015-12-07 | Optimize the method for the heterogeneous gas reservoir pressure break horizontal well crack location of extra-low permeability |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510889494.8A CN105422070B (en) | 2015-12-07 | 2015-12-07 | Optimize the method for the heterogeneous gas reservoir pressure break horizontal well crack location of extra-low permeability |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105422070A CN105422070A (en) | 2016-03-23 |
| CN105422070B true CN105422070B (en) | 2018-11-23 |
Family
ID=55500565
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510889494.8A Active CN105422070B (en) | 2015-12-07 | 2015-12-07 | Optimize the method for the heterogeneous gas reservoir pressure break horizontal well crack location of extra-low permeability |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105422070B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108256208B (en) * | 2018-01-15 | 2021-07-27 | 中国石油化工股份有限公司 | Method for calculating staged fracturing capacity of high-inclination quasi-horizontal well |
| CN109033674B (en) * | 2018-08-09 | 2022-06-10 | 西南石油大学 | Method for optimizing fracture parameters of target fracturing acidizing well |
| CN109236287A (en) * | 2018-11-08 | 2019-01-18 | 中国石油大学(北京) | One layer of shallow low permeability oil reservoir multilevel seam fracturing fracture parameter determination method and device |
| CN112765534A (en) * | 2021-01-05 | 2021-05-07 | 中国石油天然气股份有限公司 | Application method, device, terminal and storage medium of binomial capacity equation |
| CN116029232B (en) * | 2023-03-28 | 2023-07-28 | 西南石油大学 | Complex well type water invasion front prediction method for bottom water and gas reservoir |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103256035A (en) * | 2013-05-22 | 2013-08-21 | 中国石化集团华北石油局 | Fracturing fracture geological design method of horizontal well of dense gas field |
| CN104500017A (en) * | 2014-12-12 | 2015-04-08 | 中国石油天然气集团公司 | Method for optimizing staged fracturing position of horizontal well |
| CN104594872A (en) * | 2015-01-04 | 2015-05-06 | 西南石油大学 | Method for optimizing fracture conductivity of tight gas-reservoir fractured horizontal well |
| CN105019876A (en) * | 2014-04-24 | 2015-11-04 | 中国石油化工股份有限公司 | Staged fracturing horizontal well water-flooding fracture interval and well spacing determining method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120325462A1 (en) * | 2011-06-24 | 2012-12-27 | Roussel Nicolas P | Method for Determining Spacing of Hydraulic Fractures in a Rock Formation |
-
2015
- 2015-12-07 CN CN201510889494.8A patent/CN105422070B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103256035A (en) * | 2013-05-22 | 2013-08-21 | 中国石化集团华北石油局 | Fracturing fracture geological design method of horizontal well of dense gas field |
| CN105019876A (en) * | 2014-04-24 | 2015-11-04 | 中国石油化工股份有限公司 | Staged fracturing horizontal well water-flooding fracture interval and well spacing determining method |
| CN104500017A (en) * | 2014-12-12 | 2015-04-08 | 中国石油天然气集团公司 | Method for optimizing staged fracturing position of horizontal well |
| CN104594872A (en) * | 2015-01-04 | 2015-05-06 | 西南石油大学 | Method for optimizing fracture conductivity of tight gas-reservoir fractured horizontal well |
Non-Patent Citations (1)
| Title |
|---|
| "新马滚动区水平井分段压裂优化设计研究";于东;《中国优秀硕士学位论文全文数据库工程科技I辑》;20150915;第1页、第34页至第61页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105422070A (en) | 2016-03-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105422070B (en) | Optimize the method for the heterogeneous gas reservoir pressure break horizontal well crack location of extra-low permeability | |
| CN104594872B (en) | A kind of method for optimizing tight gas reservoir pressure break horizontal well fracture condudtiviy | |
| Economides et al. | How to decide between horizontal transverse, horizontal longitudinal and vertical fractured completion | |
| CN105484741A (en) | Prediction method for yield of low-permeability, heterogeneous and stress-sensitive reservoir fractured horizontal well | |
| CN109408859B (en) | Method for establishing two-dimensional triple medium numerical model of shale gas reservoir fractured horizontal well | |
| CN112434426B (en) | Shale gas multistage fracturing horizontal well step gradient pressure drop development method and device | |
| CN105404735B (en) | Quantitative evaluation method of the crack with matrix to extra-low-permeability reservoir well yield contribution rate | |
| RU2455471C1 (en) | System of solid low-productive zonary-heterogeneous oil formation development | |
| CN105350960B (en) | The method for determining low-permeable heterogeneous gas reservoir fracture parameters of fractured horizontal wells | |
| CN106321051B (en) | A method of for optimizing multistage pressure break horizontal well network fracture parameter | |
| CN109826609B (en) | It is a kind of for mostly every the accurate method for placing of proppant of interlayer oil-gas reservoir hydraulic fracturing | |
| CN105350961A (en) | Yield prediction method for volume fracturing horizontal well of low-permeability heterogeneous stress-sensitive reservoir stratum | |
| Economides et al. | Fracturing horizontal transverse, horizontal longitudinal and vertical wells: criteria for decision | |
| CN109101773A (en) | Fine and close grease horizontal well solid seam net cluster network pressure splits optimization method | |
| CN106469333B (en) | A kind of hypotonic horizontal wells in heavy oil reservoir thermal recovery pressure distribution forecasting method | |
| CA2744749A1 (en) | Basal planer gravity drainage | |
| CN105545275B (en) | Gas field gas well dispositions method and device | |
| CN106547930A (en) | Consider the gas drainage radius computational methods of tight gas reservoir seepage flow mechanism | |
| CN105205318A (en) | Method and device for determining total yield of multi-layer multi-section horizontal fractured producing well | |
| CN109146185A (en) | A kind of flooding pattern differentiation water drive optimizes and revises method | |
| WO2012145286A3 (en) | Well system with lateral main bore and strategically disposed lateral bores and method of forming | |
| Shen et al. | Evaluation of the vertical producing degree of commingled production via waterflooding for multilayer offshore heavy oil reservoirs | |
| Hough et al. | Impact of geological variation and completion type in the US Bakken oil shale play using decline curve analysis and transient flow character | |
| CN104156617A (en) | Six-stage modeling method for quality classification representation of multiple sandstone gas reservoir gas layers | |
| CN108520143A (en) | A kind of gas injection development oil reservoir gas-oil ratio climbing characterizing method |
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 |