CN104948163B - A kind of shale gas well deliverability assay method - Google Patents
A kind of shale gas well deliverability assay method Download PDFInfo
- Publication number
- CN104948163B CN104948163B CN201410112840.7A CN201410112840A CN104948163B CN 104948163 B CN104948163 B CN 104948163B CN 201410112840 A CN201410112840 A CN 201410112840A CN 104948163 B CN104948163 B CN 104948163B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mfrac
- area
- fracturing reform
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to a kind of shale gas well deliverability computational methods.The characteristics of there is fracturing reform area and non-fracturing reform area according to shale gas reservoir Area of a well in the present invention, consider adsorbed gas desorption diffusion and the influence of fracturing reform area and non-fracturing reform area to shale gas well deliverability, set up and consider that the two-region of adsorbed gas is combined shale material balance equation for gas reservoir, and combine the channelling equation between Gas Well Productivity and Liang Ge areas, it is established that a kind of fast and accurately shale gas well deliverability computational methods.This method does not need complicated shale gas reservoir geological model when calculating, it is not required that many geology and pressure break dynamic parameter.Result of calculation is suitable for a variety of applications such as shale gas well rational proration, developing value and development plan optimization.
Description
Technical field
The present invention relates to geophysical exploration technology, especially with regard to a kind of shale gas well deliverability assay method.
Background technology
Shale gas refers to that main body is located in the furvous mud shale rich in organic matter, using absorption or free state as main tax
The oil-gas accumulation of mode is deposited, shale is both reservoir, is source bed and capping layer again.
The dual media gas reservoir that shale gas reservoir is made up of the intrinsic fracture of special hypotonic matrix pores and Natural closure, base
Matter is the main preservation space of free gas and adsorbed gas.Shale gas well natural production capacity is very low, mainly passes through horizontal well multistage pressure
Split volume transform to form " artificial gas reservoir " after just have industrial production capacity, now fracture network is the main thoroughfare of seepage flow, and is being opened
There are the desorption diffusion phenomena of adsorbed gas during adopting.
Gas well deliverability, which is determined, to be referred to assess and predict the production rate-maintenance capability of gas well, the lapse rate of depletion stage and final tires out
Count gas production.Assess and predict that exactly gas well deliverability is the basis of gas well rational proration, be also optimized well pattern spacing, establishment gas
Hide the important evidence of development plan.
Normal gas pools are usually completed on the basis of geological knowledge by Geologic modeling, numerical simulation.Shale gas reservoir
Geologic feature and development scheme determine that its production capacity measuring method is more much more complex than normal gas pools.Main pass through external at present passes through
The methods such as relational expression, Production Decline Analysis and single well numerical simulation estimation shale gas well deliverability is tested, the country mainly follows conventional lines
The open-flow capacity that the method estimation shale gas well of gas reservoir flow-after-flow test test is gone into operation initial stage.Mode is bled off pressure different from general use to give birth to
The external shale gas well of production, domestic shale gas well requires stable yields in operation and supplied with keeping stable initial stage.Due to domestic and international shale
Gas well liquid loading mode is different, domestic to count the rule-of-thumb relation to be formed to calculate page based on a large amount of creation datas using foreign countries
Rock gas well deliverability.Production Decline Analysis is required could when bottom pressure changes less and gas well enters depletion stage more than half a year
Use, and domestic gas well is mainly produced initial stage in operation in pressure control stable yields mode, therefore Production Decline Analysis method can not yet
Measuring and calculating for the gas well deliverability at the initial stage of going into operation.Given many strict geology and the pressure break ginsengs of single well numerical simulation method requirement
Number, and these parameter values are difficult often accurate setting, the shale gas well for causing single well numerical simulation method to be calculated in early stage
Production capacity resultant error is big.Flow-after-flow test test can only reflect the maximum production potential of gas well at operation initial stage, it is impossible to directly and accurately
Obtain stable production period, production decline and the final cumulative gas of gas well.
In summary, although existing shale gas reservoir production capacity measuring method has his own strong points, but is due to the visitors such as the mode of production
The limitation of sight factor, difficulty is big when using, and is difficult to obtain effectively accurate shale gas Well Productivity Analysis result, especially exists
The early stage of shale gas reservoir exploitation.
The content of the invention
In view of the above-mentioned problems, the present invention proposes a kind of new fast and accurately shale gas well deliverability assay method.This hair
It is bright to have considered adsorbed gas desorption diffusion phenomena, and fracturing reform area and non-fracturing reform area are to the shadow of shale gas well deliverability
Ring, propose a kind of based on shale gas reservoir matter balance equation and to combine Gas Well Productivity and calculate the side of shale gas well deliverability
Method.
This method, comprises the following steps：
S10, test, collection and setting gas reservoir engineering parameter, and set up Gas Well Productivity；
S20, consideration adsorbed gas desorption diffusion, set up the shale gas reservoir in gas well fracturing transformation area and non-fracturing reform area respectively
Matter balance equation；
S30, according to Gas Well Productivity, calculate the initil output of initial time gas well；
S40, setting time step step-length, calculate next time step corresponding time, and update current time step；
S50, the shale gas reservoir matter balance equation according to fracturing reform area, iterate to calculate current time step fracturing transformation area
Mean reservoir pressure；
S60, according to current time step fracturing transform area mean reservoir pressure, and non-fracturing reform area shale gas reservoir
Matter balance equation, iterates to calculate the mean reservoir pressure in current time step Xia Wei fracturing reforms area；
S70, the mean reservoir pressure according to current time step fracturing transformation area and non-fracturing reform area, calculate current time step
Under shale gas well deliverability；
S80, judge whether current time is more than and given maximum evaluate number of days：
If it is not, regarding the mean reservoir pressure value of current time step as the initial value of next time step iteration, return to step S40；
If so, output shale gas well deliverability result of calculation.
Embodiments in accordance with the present invention, above-mentioned Gas Well Productivity is：
In above formula, A is Gas Well Productivity Monomial coefficient；B is the secondary term coefficient of Gas Well Productivity；piFor primitively
Stressor layer；Q is that gas well produces tolerance, 10 daily4m3/d；pwfStablize flowing bottomhole pressure (FBHP), MPa for gas well is corresponding in stable yields yield q.
In embodiments in accordance with the present invention, above-mentioned steps S20,
The shale gas reservoir matter balance equation in fracturing reform area is：
The shale gas reservoir matter balance equation in non-fracturing reform area is：
In above formula, Gp1For the cumulative gas of gas well, 104m3；Gp2For fracturing reform area adding up to non-fracturing reform area
Channelling amount, 104m3；G1For fracturing reform area reserves；G2For non-fracturing reform area reserves；p1For the fifty-fifty lamination in fracturing reform area
Power, MPa；p2For the mean reservoir pressure in non-fracturing reform area, MPa；piFor original formation pressure；Z1 a、Z2 aAnd Zi aDefine respectively
It is as follows：
In above formula, z is equation of gas state deviation factors, is the function of gas reservoir pressure；pscFor standard atmospheric pressure, MPa；
zscFor the equation of gas state deviation factors under standard state, ρBFor density, t/m3, VLFor Lan Shi volumes, m3/ t, PLFor Lan Shi
Pressure, MPa, Ф is effecive porosity, t0For formation temperature, sw0For water saturation.
Embodiments in accordance with the present invention, above-mentioned steps S50 further comprises the following steps：
S51, the fracturing reform area mean reservoir pressure made under current time step initial value p1Pressure break for a upper time step changes
Make the mean reservoir pressure value p in area10；
S52, the shale gas reservoir matter balance equation according to non-fracturing reform area, iterate to calculate under current time step, work as pressure break
The mean reservoir pressure value for transforming area is p1When corresponding non-fracturing reform area mean reservoir pressure p2；
S53, according to the current mean reservoir pressure value p in fracturing reform area and non-fracturing reform area1And p2Calculate current time step
Under gas well gas production and Liang Ge areas between channelling amount；
S54, according to time step step-length and yield, calculate the always accumulative channelling between total cumulative gas of gas well and Liang Ge areas
Amount；
S55, calculating fracturing reform area matter balance equation residual error；
S56, when residual error absolute value be less than assigned error when, exit iteration；Otherwise continue the new pressure break of iterative calculation to change
Make the mean reservoir pressure value in area, and return to step S52.
In embodiments in accordance with the present invention, above-mentioned steps S53, the average stratum in area is transformed according to current time step fracturing
Pressure value p1, and Gas Well Productivity, according to stable yields yield qgsCalculate flowing bottomhole pressure (FBHP) square value p2 wf：If the square value is small
In 0, then gas well is by pwLLevel pressure is produced；If pwfFlowing bottomhole pressure (FBHP) p minimum higher than the stable yields end of termwL, then gas well gas production is qgw=
qgs, otherwise according to the following formula with minimum flowing bottomhole pressure (FBHP) pwLProduce to calculate the gas production q of gas wellgw：
Embodiments in accordance with the present invention, calculate fracturing reform area matter balance equation residual error rsd according to the following formula：
In above formula, p1For the fracturing reform area mean reservoir pressure under current time step；piFor original formation pressure；Gp1For gas
The cumulative gas of well, 104m3；Gp2For accumulative channelling amount from fracturing reform area to non-fracturing reform area, 104m3；G1Change for pressure break
Make area's reserves；Z1 aAnd Zi aIt is defined respectively as：
In above formula, z is equation of gas state deviation factors, is the function of gas reservoir pressure；pscFor standard atmospheric pressure, MPa；
zscFor the equation of gas state deviation factors under standard state, ρBFor density, t/m3, VLFor Lan Shi volumes, m3/ t, PLFor Lan Shi
Pressure, MPa, Ф is effecive porosity, t0For formation temperature, sw0For water saturation.
Embodiments in accordance with the present invention, above-mentioned steps S60 further comprises the following steps：
S61, the non-fracturing reform area mean reservoir pressure made under current time step initial value p2For not pressing for a upper time step
Split the mean reservoir pressure value p in transformation area20；
S62, the mean reservoir pressure value p for obtaining the corresponding fracturing reform area of current time step iteration1,
S63, according to the current mean reservoir pressure value p in fracturing reform area and non-fracturing reform area1And p2Calculate current time step
Channelling amount between Xia Liangge areas；
S64, according to time step step-length and yield, calculate the always accumulative channelling amount between Liang Ge areas；
S65, the non-fracturing reform area matter balance equation residual error of calculating；
S66, when residual error absolute value be less than assigned error when, exit iteration；Otherwise continue to iterate to calculate non-fracturing reform
The mean reservoir pressure value in area, and return to step S62.
Embodiments in accordance with the present invention, calculate non-fracturing reform area matter balance equation residual error rsd according to the following formula：
In above formula, p2For the non-fracturing reform area mean reservoir pressure under current time step；piFor original formation pressure；Gp2For
Accumulative channelling amount from fracturing reform area to non-fracturing reform area, 104m3；G2For non-fracturing reform area reserves；Z2 aAnd Zi aIt is fixed respectively
Justice is as follows：
In above formula, z is equation of gas state deviation factors, is the function of gas reservoir pressure；pscFor standard atmospheric pressure, MPa；
zscFor the equation of gas state deviation factors under standard state, ρBFor density, t/m3, VLFor Lan Shi volumes, m3/ t, PLFor Lan Shi
Pressure, MPa, Ф is effecive porosity, t0For formation temperature, sw0For water saturation.
Embodiments in accordance with the present invention, above-mentioned shale gas well deliverability results of measuring can include the gas production of gas well, always tire out
Count the channelling amount between gas production, and fracturing reform area and non-fracturing reform area, always add up channelling amount.
Embodiments in accordance with the present invention, can be iterated calculating using Newton Raphson methods.
Compared with prior art, shale gas well deliverability assay method proposed by the present invention has advantages below：
1) the shale gas reservoir geological model of complexity need not be set up, it is not required that obtain many geology and pressure break dynamic ginseng
Number, can be greatly shortened operation time during implementation, simpler quickly to obtain results of measuring；
2) influence to shale gas well deliverability such as adsorbed gas desorption, fracturing reform area and non-fracturing reform area is considered,
Results of measuring is more accurate；
3) Gas Well Productivity needed for calculating typically can easily be obtained in Productivity Testing of Gas Well test, be produced
The restriction of the objective factors such as condition is smaller.
This method, which is particularly suitable for use in analyze, calculates the shale gas well production that initial stage of development does not set up shale gas reservoir geological model
Energy.Other features and advantages of the present invention will be illustrated in the following description, and partly become from specification it is aobvious and
It is clear to, or is understood by implementing the present invention.The purpose of the present invention and other advantages can be by specification, claims
And specifically noted structure is realized and obtained in accompanying drawing.
Brief description of the drawings
Accompanying drawing is used for providing a further understanding of the present invention, and constitutes a part for specification, the reality with the present invention
Apply example to be provided commonly for explaining the present invention, be not construed as limiting the invention.In the accompanying drawings：
Fig. 1 is certain shale gas reservoir multistage pressure break horizontal well Area of a well schematic diagram in the embodiment of the present invention；
Fig. 2 is the flow chart of shale gas well deliverability assay method proposed by the present invention；
Fig. 3 is to test the Gas Well Productivity set up in the embodiment of the present invention according to burnt page 1-HF wells stable productivity；
Fig. 4 is the schematic diagram of measuring and calculating shale gas well deliverability in the embodiment of the present invention；
Fig. 5 is the iterative process schematic diagram of solution fracturing reform area mean reservoir pressure in the embodiment of the present invention；
Fig. 6 is the iterative process schematic diagram that non-fracturing reform area mean reservoir pressure is solved in the embodiment of the present invention；
Fig. 7 is Jiao page 1-HF of embodiment of the present invention wells respectively by the day calculated when matching somebody with somebody production in 80,000 sides/day and 6.5 ten thousand sides/day
Gas production and cumulative gas；
Fig. 8 be the burnt page 1-HF wells of the embodiment of the present invention by 80,000 sides/day with flowing bottomhole pressure (FBHP) and fracturing reform area when producing production and
The variation diagram of non-fracturing reform area strata pressure.
Wherein, reference：
The non-fracturing reform areas of 100- fracturing reforms area 200-
Choked flow border between the 400- wells of 300- fracturing reforms area border
Embodiment
In order that disclosed technology contents are more detailed and complete, it is detailed with specific embodiment with reference to the accompanying drawings
Carefully illustrate the purpose of the present invention, technical scheme and technique effect.Although it should be strongly noted that being said for accompanying drawing
It is bright, it will be appreciated by those of skill in the art that drawings and Examples are not used for limiting the scope that the present invention is covered.
As shown in figure 1, the present invention is shale gas reservoir multistage pressure break horizontal well Area of a well schematic diagram.Shale gas reservoir list
Well control area can be divided into fracturing reform area 100 and non-two parts of fracturing reform area 200 according to well pattern spacing.Fracturing reform
Area 100 is the region for forming complex fracture network in reservoir after volume fracturing is transformed；Non- fracturing reform area 200 is adjacent
Not by the region of fracturing reform in the middle of two mouthfuls of wells.With the growth of recovery time, strata pressure in fracturing reform area 100 by
Step is reduced, and the pressure difference between non-fracturing reform area 200 gradually amplifies and forms channelling, the ground between channelling Liang Yuliangge areas
Stressor layer difference is relevant.After strata pressure is reduced to desorption pressures, desorption diffusion occurs in the adsorbed gas on matrix granule surface,
And enter Fracture System.
The characteristics of present invention has fracturing reform area and non-fracturing reform area according to shale gas reservoir Area of a well, sets up
Consider that the two-region of adsorbed gas is combined shale material balance equation for gas reservoir, and combine binomial potential curve and equation and the Liang Ge areas of gas well
Between channelling equation, it is established that a kind of fast and accurately shale gas well deliverability assay method.
Need given using method proposed by the present invention measuring and calculating shale gas well deliverability or test following parameter：
The manufacturing parameter of given gas well, such as stable production period including gas well is with production qgs(104m3/ d), the stable yields end of term minimum well
P is pressed in underflowwL(MPa) it is, maximum to evaluate number of days tmax(d)；
Given shale gas reservoir geologic parameter, such as including original formation pressure pi(MPa), formation temperature t0(DEG C), pressure break change
Make the interporosity flow coefficient λ (10 between area and non-fracturing reform area4m3/(d.MPa2)), fracturing reform area reserves G1(104m3), do not press
Split transformation area reserves G2(104m3)；
Given shale gas PVT physical parameters, such as including gas relative density γ；
Shale core parameter is measured, such as including effecive porosity Ф, density pB(t/m3), water saturation sw0, stratum
Under the conditions of the corresponding Lan Shi volumes V of shale adsorption isothermal curveL(m3/ t) and Lan Shi pressure PL(MPa)；
Gas well parameter is measured, for example, the binomial potential curve and equation coefficient of gas well is obtained including being tested by stable productivity well testing
A and B.
As shown in Fig. 2 shale gas well deliverability assay method proposed by the present invention is mainly included the following steps that：
S10, test, collection and setting gas reservoir engineering parameter, and set up Gas Well Productivity.
Collect and the manufacturing parameter of gas well, geologic parameter are set；Test and shale gas PVT parameters are set, shale core is joined
Number etc..In the present embodiment, the gas reservoir engineering parameter of burnt page 1-HF wells may refer to table 1：
Table 1
According to People's Republic of China (PRC) oil and gas industry standard SY/T 5440-2000, carry out shale gas well and go into operation just
Phase stable productivity well testing is tested.4~5 working systems can be selected during test, measuring point yield is ascending, be progressively incremented by.When
After gas well gas production and flowing bottomhole pressure (FBHP) stabilization under each working system, corresponding gas well gas production and flowing bottomhole pressure (FBHP) number are recorded
According to.After the completion of test, Gas Well Productivity (binomial equation) is set up by following formula recurrence：
In above formula, A is Gas Well Productivity Monomial coefficient；B is the secondary term coefficient of Gas Well Productivity；Q is gas well day
Gas production, 104m3/d；piFor original formation pressure, pwfStablize flowing bottomhole pressure (FBHP), MPa for gas well is corresponding in stable yields yield q.Such as
Shown in Fig. 3, in the present embodiment, by taking burnt page 1-HF wells as an example, gas well binomial potential curve and equation coefficient is respectively A=21.107, B
=2.9494.
S20, point fracturing reform area and non-fracturing reform area, set up the two-region composed page for considering that adsorbed gas desorbs diffusion phenomena
Rock material balance equation for gas reservoir：
1. the shale gas reservoir matter balance equation in fracturing reform area：
The shale gas reservoir matter balance equation in ②Wei fracturing reforms area：
In above formula, Gp1For the cumulative gas of gas well, 104m3；Gp2For fracturing reform area adding up to non-fracturing reform area
Channelling amount, 104m3；p1For the mean reservoir pressure in fracturing reform area, MPa；p2For the mean reservoir pressure in non-fracturing reform area,
MPa；Z1 a、Z2 aAnd Zi aThe influence of free gas and adsorbed gas to shale gas reservoir matter balance equation has been considered respectively, and it is defined
It is as follows：
In above formula, z is equation of gas state deviation factors, is the function of gas reservoir pressure；pscFor standard atmospheric pressure, MPa；
zscFor the equation of gas state deviation factors under standard state.
S30, according to Gas Well Productivity, calculate the initil output of t=0 moment (also referred to as initial time) gas well.
According to the gas well binomial potential curve and equation (shown in such as formula (1)) under original formation pressure, according to previously given
Gas well stable production period is with production qgsCalculate flowing bottomhole pressure (FBHP) square value p2 wf.If the square value is less than 0, illustrate that gas well can not be with yield
qgsStable yields is, it is necessary to the stable yields end of term minimum bottom pressure pwLCarry out level pressure production；If flowing bottomhole pressure (FBHP) is higher than the stable yields end of term minimum well
P is pressed in underflowwL, then gas well is q according to initil outputgsPressure control stable yields mode produce, otherwise gas well is according to the stable yields end of term minimum well
P is pressed in underflowwLProduction.
When gas well can be with stable yields, the flowing bottomhole pressure (FBHP) of gas well can be calculated according to Gas Well Productivity, otherwise need to be according to most
Low flowing bottomhole pressure (FBHP) pwLTo calculate the gas production of gas well.Circular is as follows：
1. gas well is according to stable yields yield qgsFlowing bottomhole pressure (FBHP) during production：
2. gas well is according to minimum flowing bottomhole pressure (FBHP) pwLYield during production：
Wherein, there is no producing pressure differential between initial time fracturing reform area and non-fracturing reform area, channelling amount is 0.
Further, it is preferable to the t=0 moment is set, namely the cumulative gas of gas well is G under original statep1=0, non-pressure break
Transform accumulative channelling amount from area to fracturing reform area be Gp2=0.
S40, setting time step step delta t, calculate next time step corresponding time：T=t+ Δ t, and as current
Time step is to perform following steps S50 and S60 (namely update current time step).
S50, the shale gas reservoir matter balance equation according to fracturing reform area, change according to Newton method (Newton Raphson)
In generation, calculates the mean reservoir pressure p in the fracturing reform area under current time step1。
Understand that in the specific implementation, step S50 can further comprise the following steps with reference to Fig. 4 and Fig. 5：
S51, the mean reservoir pressure value p by the fracturing reform area of a upper time step10Change as the pressure break under current time step
The initial value of area's mean reservoir pressure is made, even p1=p10。
S52, the shale gas reservoir matter balance equation according to non-fracturing reform area, current time step is calculated according to Newton iteration method
Under, when fracturing reform area mean reservoir pressure value be p1When corresponding non-fracturing reform area mean reservoir pressure p2。
S53, the gas well gas production calculated under current time step.
The mean reservoir pressure value p in area is transformed according to current time step fracturing1, and Gas Well Productivity, according to stable yields
Yield qgsCalculate flowing bottomhole pressure (FBHP) square value p2 wf.If the square value is less than 0, gas well presses pwLLevel pressure is produced；If pwfIt is higher than
The stable yields end of term minimum flowing bottomhole pressure (FBHP) pwL, then gas well gas production is qgw=qgs, otherwise according to the following formula with minimum flowing bottomhole pressure (FBHP) pwLProduction
To calculate the gas production q of gas wellgw：
And according to the current mean reservoir pressure value p in fracturing reform area and non-fracturing reform area1And p2, altered according to following
Flow equation calculates the channelling amount q between current time step Xia Liangge areas2：
S54, according to time step step delta t and yield, calculate total cumulative gas G of gas wellp1, it is total accumulative between Liang Ge areas
Channelling amount Gp2：
The accumulative aerogenesis of gas well in current step：
Total cumulative gas of gas well：Gp1=Gp10+ΔGp1 (10)
Accumulative channelling amount between current step Nei Liangge areas：
Always accumulative channelling amount between Liang Ge areas：Gp2=Gp20+ΔGp2 (12)
In above formula, qgw0For the gas well gas production of last time step, 104m3/d；q20For the channelling between last time step Liang Ge areas
Amount, 104m3/d；Gp10For the total cumulative gas of gas well of last time step, 104m3；Gp20To be total between the Liang Ge areas of last time step
Accumulative channelling amount, 104m3。
S55, according to the fracturing reform area mean reservoir pressure p under current time step1, calculate equation of gas state deviation factors
Z and matter balance equation coefficient Z1 a, according to original formation pressure piCalculate Zi a, bring following formula into and calculate fracturing reform area material balance
Equation residual error rsd：
S56, when residual error rsd absolute value be less than assigned error ε 0 when, exit iteration；Otherwise, according to Newton iteration method meter
New p1Value, and return to step 52), continue iteration.
S60, the p according to the current time step of iterative calculation1Value, according to the shale gas reservoir material balance side in non-fracturing reform area
Journey, iterates to calculate the mean reservoir pressure p in current time step Xia Wei fracturing reforms area2。
It should be noted that the process of step iterative calculation is as above-mentioned steps S52 iterative process.Wherein,
S52 is to calculate p1The pilot process called is needed during value, is to work as p herein1After value is finally determined through iteration, based on the p1Value meter
Calculate corresponding p2Value.
Understand that in the specific implementation, step S60 can further comprise the following steps with reference to Fig. 4 and Fig. 6：
S61, the mean reservoir pressure value p by the non-fracturing reform area of a upper time step20It is used as not pressing under current time step
The initial value of transformation area's mean reservoir pressure is split, even p2=p20。
S62, the mean reservoir pressure value p for obtaining the corresponding fracturing reform area of current time step iteration1。
S63, according to the current mean reservoir pressure value p in fracturing reform area and non-fracturing reform area1And p2, altered according to following
Flow equation (same to formula (8)) calculates the channelling amount q between current time step Xia Liangge areas2：
S64, according to time step step delta t and yield, calculate the always accumulative channelling amount G between Liang Ge areasp2：
Accumulative channelling amount between current step Nei Liangge areas：
Always accumulative channelling amount between Liang Ge areas：Gp2=Gp20+ΔGp2 (16)
In above formula, q2For the channelling amount between current time step Liang Ge areas, 104m3/d；q20For between last time step Liang Ge areas
Channelling amount, 104m3/d；Gp20For the always accumulative channelling amount between the Liang Ge areas of last time step, 104m3。
S65, according to the non-fracturing reform area mean reservoir pressure p under current time step2, calculate equation of gas state deviation because
Sub- z and matter balance equation coefficient Z2 a, according to original formation pressure piCalculate Zi a, bring following formula into and calculate non-fracturing reform area material
Equilibrium equation residual error rsd：
S66, when residual error rsd absolute value be less than assigned error ε 0 when, exit iteration；Otherwise calculated according to Newton iteration method
New p2Value, and return to step 62), continue iteration.
S70, the current time step fracturing obtained according to step S50 and S60 transform area and non-fracturing reform area fifty-fifty
Stressor layer, calculates the shale gas well deliverability under current time step, includes gas production, the total cumulative gas of gas well, and pressure break changes
Make the channelling amount between area and non-fracturing reform area, always add up channelling amount.
According to the fracturing reform area mean reservoir pressure value p of current time step1, gas well gas production is calculated according to step S53；
According to the mean reservoir pressure in fracturing reform area and non-fracturing reform area, according between formula (8) calculating Liang Ge areas
Channelling amount；
According to total cumulative gas of current time step and the gas production of last time step, and last time step, according to formula
(10) total cumulative gas of current time step is calculated, and formula (12) calculates the always accumulative channelling amount of current time step.
S80, judge whether current time t is more than and given maximum evaluate number of days tmax：
If not, return to step S40, continues to calculate.
If it is, output shale gas well deliverability result of calculation, such as gas well gas production, total cumulative gas, and pressure break
Transform the channelling amount between area and non-fracturing reform area, always add up the result of calculations such as channelling amount.
As shown in fig. 7, in the present embodiment, burnt page 1-HF wells are matched somebody with somebody according to the stable production period in 80,000 sides/day and 6.5 ten thousand sides/day respectively
Production, is estimated gas well deliverability using the above method.Wherein, gas well is in stable production period, the cumulative gas in the stable yields end of term and pre-
The estimation result of total cumulative gas at the end of the survey phase may refer to table 2.Because the gas well requires that stable production period reaches 1.5 years,
Therefore 6.5 ten thousand sides/day is not to be exceeded with production in the AOF calculation result obtained according to the present invention, the well.
| Index | With production scheme one | With production scheme two |
| Stable production period yield (104m3/d) | 8 | 6.5 |
| Stable production period (year) | 1.08 | 1.66 |
| Stable yields end of term cumulative gas (104m3) | 2935.6 | 3815.2 |
| Finally tire out gas production (104m3) | 9084.5 | 9079 |
Table 2
Accordingly, Fig. 8 shows the flowing bottomhole pressure (FBHP) of the gas well, fracturing reform area and non-fracturing reform area strata pressure at any time
Between situation of change.
Although presently disclosed embodiment is as above, described content is only to facilitate understanding the present invention and adopting
Embodiment, is not limited to the present invention.Any those skilled in the art to which this invention pertains, are not departing from this
On the premise of the disclosed spirit and scope of invention, any modification made in the formal and details of implementation and change, all
Should be in the scope of patent protection of the present invention.
Claims (10)
1. a kind of shale gas well deliverability assay method, comprises the following steps：
S10, test, collection and setting gas reservoir engineering parameter, and set up Gas Well Productivity；
S20, consideration adsorbed gas desorption diffusion, set up the shale gas reservoir material in gas well fracturing transformation area and non-fracturing reform area respectively
Equilibrium equation；
S30, according to Gas Well Productivity, calculate the initil output of initial time gas well；
S40, setting time step step-length, calculate next time step corresponding time, and update current time step；
S50, the shale gas reservoir matter balance equation according to fracturing reform area, iterate to calculate the flat of current time step fracturing transformation area
Equal strata pressure；
S60, according to current time step fracturing transform area mean reservoir pressure, and non-fracturing reform area shale gas reservoir material
Equilibrium equation, iterates to calculate the mean reservoir pressure in current time step Xia Wei fracturing reforms area；
S70, the mean reservoir pressure according to current time step fracturing transformation area and non-fracturing reform area, are calculated under current time step
Shale gas well deliverability；
S80, judge whether current time is more than and given maximum evaluate number of days：
If it is not, regarding the mean reservoir pressure value of current time step as the initial value of next time step iteration, return to step S40；
If so, output shale gas well deliverability result of calculation.
2. shale gas well deliverability assay method as claimed in claim 1, it is characterised in that the Gas Well Productivity is：
<mrow>
<mfrac>
<mrow>
<msubsup>
<mi>p</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>wf</mi>
<mn>2</mn>
</msubsup>
</mrow>
<mi>q</mi>
</mfrac>
<mo>=</mo>
<mi>A</mi>
<mo>+</mo>
<mi>Bq</mi>
</mrow>
In above formula, A is Gas Well Productivity Monomial coefficient；B is the secondary term coefficient of Gas Well Productivity；piTo be primitively laminated
Power；Q is that gas well produces tolerance, 10 daily4m3/d；pwfStablize flowing bottomhole pressure (FBHP), MPa for gas well is corresponding in stable yields yield q.
3. shale gas well deliverability assay method as claimed in claim 1, it is characterised in that in the step S20,
The shale gas reservoir matter balance equation in fracturing reform area is：
<mrow>
<mfrac>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>Z</mi>
<mn>1</mn>
<mi>a</mi>
</msubsup>
</mfrac>
<mo>=</mo>
<mfrac>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>Z</mi>
<mi>i</mi>
<mi>a</mi>
</msubsup>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>p</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>p</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>G</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
The shale gas reservoir matter balance equation in non-fracturing reform area is：
<mrow>
<mfrac>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>Z</mi>
<mn>2</mn>
<mi>a</mi>
</msubsup>
</mfrac>
<mo>=</mo>
<mfrac>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>Z</mi>
<mi>i</mi>
<mi>a</mi>
</msubsup>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>G</mi>
<mrow>
<mi>p</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>G</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
In above formula, Gp1For the cumulative gas of gas well, 104m3；Gp2For accumulative channelling from fracturing reform area to non-fracturing reform area
Amount, 104m3；G1For fracturing reform area reserves；G2For non-fracturing reform area reserves；p1For the mean reservoir pressure in fracturing reform area,
MPa；p2For the mean reservoir pressure in non-fracturing reform area, MPa；piFor original formation pressure；Z1 a、Z2 aAnd Zi aDefine respectively such as
Under：
<mrow>
<msubsup>
<mi>Z</mi>
<mi>j</mi>
<mi>a</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mi>z</mi>
<mrow>
<mo>[</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>w</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>]</mo>
<mo>+</mo>
<mfrac>
<mrow>
<mi>z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mn>273.15</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>p</mi>
<mi>SC</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mi>B</mi>
</msub>
<msub>
<mi>V</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>293.15</mn>
<msub>
<mi>Z</mi>
<mi>SC</mi>
</msub>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>L</mi>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mfrac>
<mo>,</mo>
<mo>=</mo>
<mn>1,2</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
In above formula, z is equation of gas state deviation factors, is the function of gas reservoir pressure；pscFor standard atmospheric pressure, MPa；zscFor mark
Equation of gas state deviation factors under quasi- state, ρBFor density, t/m3, VLFor Lan Shi volumes, m3/t, PLFor Lan Shi pressure,
MPa, Ф is effecive porosity, t0For formation temperature, sw0For water saturation.
4. shale gas well deliverability assay method as claimed in claim 1, it is characterised in that the step S50 is further wrapped
Include following steps：
S51, the fracturing reform area mean reservoir pressure made under current time step initial value p1For the fracturing reform area of a upper time step
Mean reservoir pressure value p10；
S52, the shale gas reservoir matter balance equation according to non-fracturing reform area, iterate to calculate under current time step, work as fracturing reform
The mean reservoir pressure value in area is p1When corresponding non-fracturing reform area mean reservoir pressure p2；
S53, according to the current mean reservoir pressure value p in fracturing reform area and non-fracturing reform area1And p2Calculate under current time step
Channelling amount between gas well gas production and Liang Ge areas；
S54, according to time step step-length and yield, calculate the always accumulative channelling amount between total cumulative gas of gas well and Liang Ge areas；
S55, calculating fracturing reform area matter balance equation residual error；
S56, when residual error absolute value be less than assigned error when, exit iteration；Otherwise continue to iterate to calculate new fracturing reform area
Mean reservoir pressure value, and return to step S52.
5. shale gas well deliverability assay method as claimed in claim 4, it is characterised in that in the step S53, according to current
Time step fracturing transforms the mean reservoir pressure value p in area1, and Gas Well Productivity, according to stable yields yield qgsCalculate shaft bottom stream
Press square value p2 wf：If the square value is less than 0, gas well presses pwLLevel pressure is produced；If pwfHigher than the stable yields end of term minimum shaft bottom
Stream pressure pwL, then gas well gas production is qgw=qgs, otherwise according to the following formula with minimum flowing bottomhole pressure (FBHP) pwLProduce to calculate the aerogenesis of gas well
Measure qgw：
<mrow>
<msub>
<mi>q</mi>
<mi>gw</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<mi>A</mi>
<mo>+</mo>
<msqrt>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>4</mn>
<mi>B</mi>
<mo>[</mo>
<msubsup>
<mi>p</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>p</mi>
<mi>wL</mi>
<mn>2</mn>
</msubsup>
</msqrt>
</mrow>
<mrow>
<mn>2</mn>
<mi>B</mi>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
6. shale gas well deliverability assay method as claimed in claim 4, it is characterised in that calculate fracturing reform area according to the following formula
Matter balance equation residual error rsd：
<mrow>
<mi>rsd</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>Z</mi>
<mn>1</mn>
<mi>a</mi>
</msubsup>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>Z</mi>
<mi>i</mi>
<mi>a</mi>
</msubsup>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>p</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>p</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<msub>
<mi>G</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
In above formula, p1For the fracturing reform area mean reservoir pressure under current time step；piFor original formation pressure；Gp1For gas well
Cumulative gas, 104m3；Gp2For accumulative channelling amount from fracturing reform area to non-fracturing reform area, 104m3；G1For fracturing reform area
Reserves；Z1 aAnd Zi aIt is defined respectively as：
<mrow>
<msubsup>
<mi>Z</mi>
<mi>j</mi>
<mi>a</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mi>z</mi>
<mrow>
<mo>[</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>w</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>]</mo>
<mo>+</mo>
<mfrac>
<mrow>
<mi>z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mn>273.15</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>p</mi>
<mi>SC</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mi>B</mi>
</msub>
<msub>
<mi>V</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>293.15</mn>
<msub>
<mi>Z</mi>
<mi>SC</mi>
</msub>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>L</mi>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mfrac>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
In above formula, z is equation of gas state deviation factors, is the function of gas reservoir pressure；pscFor standard atmospheric pressure, MPa；zscFor mark
Equation of gas state deviation factors under quasi- state, ρBFor density, t/m3, VLFor Lan Shi volumes, m3/ t, PLFor Lan Shi pressure,
MPa, Ф are effecive porosity, t0For formation temperature, sw0For water saturation.
7. the shale gas well deliverability assay method as described in claim 1~6 any one, it is characterised in that the step S60
Further comprise the following steps：
S61, the non-fracturing reform area mean reservoir pressure made under current time step initial value p2For the non-fracturing reform of a upper time step
The mean reservoir pressure value p in area20；
S62, the mean reservoir pressure value p for obtaining the corresponding fracturing reform area of current time step iteration1,
S63, according to the current mean reservoir pressure value p in fracturing reform area and non-fracturing reform area1And p2Calculate two under current time step
Channelling amount between individual area；
S64, according to time step step-length and yield, calculate the always accumulative channelling amount between Liang Ge areas；
S65, the non-fracturing reform area matter balance equation residual error of calculating；
S66, when residual error absolute value be less than assigned error when, exit iteration；Otherwise continue to iterate to calculate non-fracturing reform area
Mean reservoir pressure value, and return to step S62.
8. shale gas well deliverability assay method as claimed in claim 7, it is characterised in that calculate non-fracturing reform according to the following formula
Area matter balance equation residual error rsd：
<mrow>
<mi>rsd</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>Z</mi>
<mn>2</mn>
<mi>a</mi>
</msubsup>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>Z</mi>
<mi>i</mi>
<mi>a</mi>
</msubsup>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>G</mi>
<mrow>
<mi>p</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>G</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
In above formula, p2For the non-fracturing reform area mean reservoir pressure under current time step；piFor original formation pressure；Gp2For pressure break
Transform accumulative channelling amount of the area to non-fracturing reform area, 104m3；G2For non-fracturing reform area reserves；Z2 aAnd Zi aDefine respectively such as
Under：
<mrow>
<msubsup>
<mi>Z</mi>
<mi>j</mi>
<mi>a</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mi>z</mi>
<mrow>
<mo>[</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>s</mi>
<mrow>
<mi>w</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>]</mo>
<mo>+</mo>
<mfrac>
<mrow>
<mi>z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mn>273.15</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>p</mi>
<mi>SC</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mi>B</mi>
</msub>
<msub>
<mi>V</mi>
<mi>L</mi>
</msub>
</mrow>
<mrow>
<mn>293.15</mn>
<msub>
<mi>Z</mi>
<mi>SC</mi>
</msub>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>L</mi>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mfrac>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>2</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
In above formula, z is equation of gas state deviation factors, is the function of gas reservoir pressure；pscFor standard atmospheric pressure, MPa；zscFor mark
Equation of gas state deviation factors under quasi- state, ρBFor density, t/m3, VLFor Lan Shi volumes, m3/ t, PLFor Lan Shi pressure,
MPa, Ф are effecive porosity, t0For formation temperature, sw0For water saturation.
9. shale gas well deliverability assay method as claimed in claim 1, it is characterised in that the shale gas well deliverability calculates knot
Fruit includes gas production, the total cumulative gas of gas well, and channelling amount between fracturing reform area and non-fracturing reform area, always tires out
Count channelling amount.
10. shale gas well deliverability assay method as claimed in claim 1, it is characterised in that using Newton Raphson sides
Method is iterated calculating.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410112840.7A CN104948163B (en) | 2014-03-24 | 2014-03-24 | A kind of shale gas well deliverability assay method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410112840.7A CN104948163B (en) | 2014-03-24 | 2014-03-24 | A kind of shale gas well deliverability assay method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104948163A CN104948163A (en) | 2015-09-30 |
| CN104948163B true CN104948163B (en) | 2017-09-15 |
Family
ID=54163171
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410112840.7A Active CN104948163B (en) | 2014-03-24 | 2014-03-24 | A kind of shale gas well deliverability assay method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104948163B (en) |
Families Citing this family (34)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107239648B (en) * | 2016-03-25 | 2020-07-10 | 中国石油化工股份有限公司 | Shale gas well yield composition determination method and device |
| CN107435528B (en) * | 2016-05-20 | 2020-08-07 | 中国石油天然气股份有限公司 | Volcanic gas reservoir gas well production allocation method |
| CN107563899B (en) * | 2016-06-30 | 2021-04-30 | 中国石油天然气股份有限公司 | Oil-gas well productivity prediction method and device |
| CN106227995B (en) * | 2016-07-16 | 2019-01-11 | 中石化重庆涪陵页岩气勘探开发有限公司 | Shale gas reservoir single well productivity calculation method |
| CN106326549B (en) * | 2016-08-23 | 2019-11-08 | 中国石油天然气股份有限公司 | Method and system for determining formation parameters of vertical well shale gas reservoir |
| CN106285568B (en) * | 2016-08-31 | 2018-08-21 | 延安大学 | A kind of branch horizontal well recovery method based on Green's function |
| CN106285569B (en) * | 2016-08-31 | 2018-08-21 | 延安大学 | A kind of branch horizontal well recovery method based on slippage coefficient |
| CN106437629B (en) * | 2016-10-31 | 2018-08-21 | 延安大学 | A kind of branch horizontal well recovery method based on starting pressure gradient |
| CN106437630B (en) * | 2016-10-31 | 2018-08-21 | 延安大学 | A kind of branch horizontal well recovery method based on diffusion coefficient |
| CN106351612B (en) * | 2016-10-31 | 2018-08-21 | 延安大学 | A kind of branch horizontal well recovery method based on fractal dimension |
| CN106640042A (en) * | 2016-12-23 | 2017-05-10 | 中国石油天然气股份有限公司 | Device and method for evaluating single-well energy production of gas well |
| CN107066674B (en) * | 2017-01-18 | 2021-08-10 | 西南石油大学 | Method for calculating unsteady state yield of shale gas reservoir volume fracturing horizontal well |
| CN106761677B (en) * | 2017-03-09 | 2020-10-09 | 长江大学 | Logging prediction method for single-well productivity of shale gas horizontal well |
| CN108625823A (en) * | 2017-03-20 | 2018-10-09 | 中国石油天然气股份有限公司 | Method for determining the capacity of a gas reservoir |
| CN107120106B (en) * | 2017-05-25 | 2020-06-23 | 中石化石油工程技术服务有限公司 | Shale quality evaluation method based on organic porosity and total organic carbon content |
| CN107622165B (en) * | 2017-09-25 | 2018-06-05 | 西南石油大学 | A kind of shale gas horizontal well refracturing Productivity |
| CN109554303B (en) * | 2017-09-27 | 2022-04-12 | 中国石油化工股份有限公司 | Microbacterium aurantiacus and preparation and application thereof |
| CN109594968B (en) * | 2017-09-28 | 2022-04-12 | 中国石油化工股份有限公司 | Shale gas multistage fracturing horizontal well post-fracturing fracture parameter evaluation method and system |
| CN108180008B (en) * | 2017-12-27 | 2021-08-31 | 中国石油化工股份有限公司 | Method and device for predicting collectable reserve, electronic device and storage medium |
| CN110080743B (en) * | 2018-01-24 | 2022-03-29 | 中国石油天然气股份有限公司 | Oil well potential detection method |
| CN108547610B (en) * | 2018-02-07 | 2021-06-01 | 中国石油天然气股份有限公司 | Method and device for determining horizontal well productivity under volume fracturing |
| CN108442911A (en) * | 2018-02-28 | 2018-08-24 | 西南石油大学 | A kind of shale gas horizontal well refracturing hydraulic fracture method for optimally designing parameters |
| CN108547612B (en) | 2018-03-02 | 2022-03-29 | 中国石油天然气股份有限公司 | Method and device for determining oil yield and gas yield of shale oil in-situ conversion |
| CN109236258B (en) * | 2018-10-27 | 2019-07-23 | 中国地质大学（北京） | A kind of compact oil reservoir pressure break horizontal well optimization method based on Adaptive proxy model |
| CN111199010B (en) * | 2018-11-20 | 2023-03-31 | 中国石油化工股份有限公司 | Correction isochronous well testing-based tight gas reservoir gas well productivity calculation method and device |
| CN109657299B (en) * | 2018-11-28 | 2023-01-17 | 郭建林 | Shale gas reservoir mining method |
| CN109441422B (en) | 2018-12-03 | 2020-06-26 | 北京科技大学 | Shale gas well spacing optimization mining method |
| CN110135122B (en) * | 2019-06-17 | 2021-04-23 | 中国科学技术大学 | Low-permeability natural gas reservoir yield prediction method and system |
| CN111506865B (en) * | 2020-07-01 | 2020-09-15 | 西南石油大学 | Shale gas well production fitting and predicting method |
| CN111911115B (en) * | 2020-08-12 | 2021-09-17 | 西南石油大学 | Shale gas well dynamic production allocation method |
| CN112392473B (en) * | 2020-11-17 | 2023-11-28 | 中国石油天然气股份有限公司 | Method for evaluating injection and production capacity of low-permeability gas reservoir gas well |
| CN112964597B (en) * | 2021-03-09 | 2022-09-16 | 中国石油大学(华东) | Multi-scale rock core adsorption-analysis-exploitation simulation experiment device and method |
| US11578595B2 (en) | 2021-04-06 | 2023-02-14 | Saudi Arabian Oil Company | Systems and methods for selecting and performing gas deliverability tests |
| CN115126473B (en) * | 2022-06-29 | 2023-08-15 | 重庆科技学院 | Calculation method for standardized test yield of shale gas well |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7369979B1 (en) * | 2005-09-12 | 2008-05-06 | John Paul Spivey | Method for characterizing and forecasting performance of wells in multilayer reservoirs having commingled production |
| CN102509179B (en) * | 2011-12-26 | 2014-09-24 | 中国石油天然气股份有限公司 | Gas well dynamic capacity prediction method |
| US20130262069A1 (en) * | 2012-03-29 | 2013-10-03 | Platte River Associates, Inc. | Targeted site selection within shale gas basins |
| CN102707333A (en) * | 2012-06-06 | 2012-10-03 | 陕西延长石油(集团)有限责任公司研究院 | Shale gas resource/reserve measurement method |
| CN103266881B (en) * | 2013-05-22 | 2015-05-20 | 中国石化集团华北石油局 | Method for predicting yield of compact hypotonic gas field multistage fracturing horizontal well |
-
2014
- 2014-03-24 CN CN201410112840.7A patent/CN104948163B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104948163A (en) | 2015-09-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104948163B (en) | A kind of shale gas well deliverability assay method | |
| CN106484933B (en) | It is a kind of for determining the method and system of shale gas well well control dynamic holdup | |
| CN104018829B (en) | A kind of method that utilization coal bed gas well creation data measures air water phase percolation curve | |
| CN107044277B (en) | Low permeable and heterogeneity reservoir horizontal well refracturing yield potential evaluation method | |
| Kikani et al. | Perturbation analysis of stress-sensitive reservoirs | |
| CN105507893A (en) | Method for predicting production effect of shale reservoir after size modification | |
| CN107301306A (en) | Dynamic open-flow capacity Forecasting Methodology for DAMAGE OF TIGHT SAND GAS RESERVOIRS pressure break horizontal well | |
| CN104318052B (en) | A kind of evaluation method of low-permeability oil deposit water drive sweep efficiency | |
| CN110147561B (en) | Method for predicting volume fracture network of tight oil and gas reservoir containing natural fracture | |
| CN104750896B (en) | A kind of fractured-cavernous carbonate reservoir method for numerical simulation | |
| CN107038313B (en) | Layered crustal stress fine description method based on numerical value core | |
| CN105484741A (en) | Prediction method for yield of low-permeability, heterogeneous and stress-sensitive reservoir fractured horizontal well | |
| CN106481332A (en) | Method for determining area's dynamic holdup inside and outside shale gas multistage pressure break horizontal well | |
| Kalantari-Dahaghi et al. | Coupling numerical simulation and machine learning to model shale gas production at different time resolutions | |
| CN105672973A (en) | Coalbed methane overall volume fracturing optimization design method | |
| Zheng et al. | Numerical prediction of the decline of the shale gas production rate with considering the geomechanical effects based on the two-part Hooke’s model | |
| CN105089566A (en) | Gas well system production allocation method | |
| CN105404735A (en) | Quantitative evaluation method for contribution rate of fractures and matrix to single well yield of ultra-low permeability reservoir | |
| Young | Computer modeling and simulation of coalbed methane resources | |
| CN105649616B (en) | Method for evaluating dynamic reserve of gas well under underground throttling condition of low-permeability gas reservoir | |
| CN108729912B (en) | Yield splitting method suitable for numerical reservoir simulation | |
| Chen et al. | A workflow based on a semianalytical model to estimate the properties of stimulated reservoir volume of tight-oil wells | |
| Sun et al. | The optimum pressure drawdown for production from a shale gas reservoir: A numerical study with a coupled geomechanics and reservoir model | |
| CN107704646B (en) | Modeling method after compact reservoir volume modification | |
| CN105134187B (en) | Method and device for assisting oil deposit sedimentary facies division and connectivity evaluation |
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 |