CN106481332A - Method for determining area's dynamic holdup inside and outside shale gas multistage pressure break horizontal well - Google Patents
Method for determining area's dynamic holdup inside and outside shale gas multistage pressure break horizontal well Download PDFInfo
- Publication number
- CN106481332A CN106481332A CN201510549414.4A CN201510549414A CN106481332A CN 106481332 A CN106481332 A CN 106481332A CN 201510549414 A CN201510549414 A CN 201510549414A CN 106481332 A CN106481332 A CN 106481332A
- Authority
- CN
- China
- Prior art keywords
- area
- fracturing reform
- pressure
- reform area
- gas
- 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.)
- Granted
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of method for determining area's dynamic holdup inside and outside shale gas multistage pressure break horizontal well, including：Gas well binomial potential curve and equation is set up by gas testing data, calculates actual flowing bottomhole pressure (FBHP) by creation data and gas well structural parameters；Consider adsorbed gas desorbing and abnormal pressure impact, set up respectively and be directed to the fracturing reform area of multistage pressure break horizontal well and the matter balance equation in non-fracturing reform area；Give the dynamic holdup in fracturing reform area and non-fracturing reform area respectively, and the channelling equation between the matter balance equation with reference to fracturing reform area and non-fracturing reform area and twoth area, binomial potential curve and equation adjustment dynamic holdup carrys out the actual flowing bottomhole pressure (FBHP) of matching, so that it is determined that area's dynamic holdup inside and outside pressure break horizontal well.The present invention does not need closing well to test mean reservoir pressure it is possible to consider adsorbed gas desorbing and the impact of abnormal pressure and pressure break recombination region in use, can determine fracturing reform area and periphery non-fracturing reform area dynamic holdup respectively.
Description
Technical field
The invention belongs to shale gas exploration and development technical field, specifically, it is related to one kind and is used for determining shale gas
The method of area's dynamic holdup inside and outside multistage pressure break horizontal well.
Background technology
The evaluation methodology of shale Estimation of Gas Well Dynamic Reserves is widely different with conventional gas well, is mainly reflected in：(1) gas
Hide geologic feature and development technique is different.Shale gas reservoir is certainly to be conigenous gas storage to hide, and permeability spy is low, needs to pass through
Volume fracturing transformation comes " man Tibetan ", and gas well yield and recoverable reserves are limited by pressure break, individual well dynamic holdup
Relevant with volume fracturing transformation.(2) seepage flow characteristics are different.Due to shale gas reservoir, matrix permeability spy is low, gas
Tibetan is difficult to enter the boundary control stream stage, and dynamic holdup evaluation requires to reach the boundary control stream stage, otherwise comments
The dynamic holdup of valency be only pressure involve in the range of mining-employed reserves.(3) adsorbed gas desorbing diffusion.Free gas
All affect shale gas well deliverability with adsorbed gas, dynamic holdup evaluation needs to consider adsorbed gas desorbing impact.
North America shale gas well produces essentially according to bleeding off pressure production decline mode at present, mainly adopts production rate decline curve
To analyze the ultimate recoverable reserves (EUR) of gas well, mainly have improved Arps method, power law index method,
Diffusion index method and Duong method etc..Produced come matching by adjusting the coefficient in these decline curve models
Amount data, then prediction gas well yield and recoverable reserves.
Decline curve analysis method does not require nothing more than gas well in the change of production phase flowing bottomhole pressure (FBHP) less, and also requires that
Keep constant in forecast period flowing bottomhole pressure (FBHP).Additionally, the method also requires gas well to reach the boundary control stream stage,
I.e. pressure wave reaches reservoir physical boundary or choked flow border, and the yield otherwise predicted and ultimate recoverable reserves are higher.
Domestic gas well generally requires the stable production period of 2-3 to ensure supply of stabilizing the market at present, afterwards just according to
Level pressure production decline mode produces, for example Fuling reef dam shale gas field development program design gas well stable production period be
2 years.It can thus be seen that due to the difference of the mode of production, domestic shale gas well is difficult to adopt in commitment
Decline curve analysis method is calculating Estimation of Gas Well Dynamic Reserves.
Material balance method is for determining the conventional method of Estimation of Gas Well Dynamic Reserves in normal gas pools, and the method needs flat
All strata pressuresWith tired gas production GPData.Calculated by tired gas production and mean reservoir pressure valueValue, then
Rectangular plots are drawn out a series ofWith GPData point, fitting a straight line is extrapolated to x-axis and is gas well waterout
Reserves.Wherein, mean reservoir pressure is mainly obtained by gas well shut-in pressure recovery well testing test interpretation.
The method is applied to mainly there is problems with when calculating shale Estimation of Gas Well Dynamic Reserves：One is to move in calculating gas well
The impact of shale substrate adsorbed gas desorbing diffusion cannot be considered during state reserves；Two is that shale matrix permeability spy is low,
It is difficult to recover transient well test test to explain mean reservoir pressure by shut-in pressure；Three is to cannot be distinguished by pressure break
Transformation area's dynamic holdup and periphery non-fracturing reform area dynamic holdup.
From the point of view of comprehensive shale gas well waterout Reserve Estimation Method both at home and abroad, lack effectively accurate shale gas well at present
Dynamic holdup computational methods are it is impossible to distinguish fracturing reform area dynamic holdup and non-fracturing reform area dynamic holdup.
Content of the invention
For solving problem above, the invention provides one kind is used for determining area inside and outside shale gas multistage pressure break horizontal well
The method of dynamic holdup.
According to one embodiment of present invention, there is provided one kind is used for determining inside and outside shale gas multistage pressure break horizontal well
The method of area's dynamic holdup, including：
Obtain shale gas reservoir gas testing data, creation data and gas well structural parameters, and built by described gas testing data
Vertical gas well binomial potential curve and equation, calculates actual flowing bottomhole pressure (FBHP) by described creation data and gas well structural parameters；
Consider adsorbed gas desorbing and abnormal pressure impact, set up the fracturing reform for multistage pressure break horizontal well respectively
Area and the matter balance equation in non-fracturing reform area；
Give the dynamic holdup in described fracturing reform area and described non-fracturing reform area respectively, and combine described pressure break
Channelling equation between the matter balance equation in transformation area and described non-fracturing reform area and twoth area, described binomial
Deliverability equation adjusts described dynamic holdup and carrys out actual flowing bottomhole pressure (FBHP) described in matching, so that it is determined that inside and outside pressure break horizontal well
Area's dynamic holdup.
According to one embodiment of present invention, set up fracturing reform area and the non-pressure break for multistage pressure break horizontal well
The step of the matter balance equation in transformation area further includes：
Consider adsorbed gas desorbing and abnormal pressure impact, set up shale gas reservoir matter balance equation；
Multistage pressure break horizontal well Area of a well is divided into by fracturing reform area and non-pressure break according to well pattern spacing
Transformation area；
Described fracturing reform area is set up respectively based on described shale gas reservoir matter balance equation and described non-pressure break changes
Make the matter balance equation in area.
According to one embodiment of present invention, the step setting up described shale gas reservoir matter balance equation is wrapped further
Include：
Based on pore space compressibility of rock, calculate during strata pressure change due to rock matrix compression and fluid expansion
The underground pore volume decrement causing；
Drop in formation pressure is calculated based on described underground pore volume decrement and Langmuir isothermal adsorpting equation
Shale gas reservoir residual free gas reserves afterwards and remaining adsorbed gas reserves；
According to the law of indestructibility of matter：Original free gas reserves+original adsorbed gas reserves=residual free gas reserves+surplus
Remaining adsorbed gas yield+tired gas production, sets up described shale gas reservoir matter balance equation.
According to one embodiment of present invention, described shale gas reservoir matter balance equation is：
Wherein, Za=Za(p),
Work as p=piWhen,P is strata pressure, piFor original formation pressure, GpFor cumulative production,
G is dynamic holdup, cfFor pore space compressibility of rock, cwFor stratum water coefficient of compressibility, SwiContain for gas reservoir is original
Water saturation, SgiFor the original gas saturation of gas reservoir, φ is effecive porosity, ρBFor shale density, VLFor orchid
Ge Miaoer volume, PLFor Lan Shi pressure, pscFor gas pressure under standard state, T is temperature, TscFor standard shape
Temperature under state, z is compressibility factor, zscFor the compressibility factor under standard state.
According to one embodiment of present invention, the matter balance equation in described fracturing reform area is：
Wherein, p1For the mean reservoir pressure in fracturing reform area,G1Dynamic storage for fracturing reform area
Amount, Gp1For gas well cumulative production, Gp2For accumulative channelling amount from non-fracturing reform area to fracturing reform area.
According to one embodiment of present invention, the matter balance equation in described non-fracturing reform area is：
Wherein, p2For the mean reservoir pressure in non-fracturing reform area,Gp2For non-fracturing reform area to
The accumulative channelling amount in fracturing reform area, G2Dynamic holdup for non-fracturing reform area.
According to one embodiment of present invention, determine that the step of area's dynamic holdup inside and outside pressure break horizontal well is wrapped further
Include：
The fracturing reform area dynamic holdup G of given shale gas well1, non-fracturing reform area dynamic holdup G2；
Based on fracturing reform area dynamic holdup G1, non-fracturing reform area dynamic holdup G2, the material in fracturing reform area
Channelling equation between equilibrium equation, the matter balance equation in non-fracturing reform area, Gas Well Productivity and two-region,
The average stratum that tolerance is iterated to calculate fracturing reform area by fracturing reform area matter balance equation is produced daily according to gas well
Pressure；
Mean reservoir pressure p based on described fracturing reform area1With described binomial potential curve and equation, by daily output tolerance
Prediction flowing bottomhole pressure (FBHP)；
The flowing bottomhole pressure (FBHP) of matching prediction and described actual flowing bottomhole pressure (FBHP) are to determine fracturing reform area described in horizontal well
Dynamic holdup and the dynamic holdup in described non-fracturing reform area.
According to one embodiment of present invention, the step of the mean reservoir pressure in iterative calculation fracturing reform area enters one
Step includes：
By upper time step ti-1Fracturing reform area mean reservoir pressure value p10As current tiTime step fracturing reform area
The iterative initial value p of mean reservoir pressure1If, ti=0, then iterative initial value is original formation pressure value；
According to fracturing reform area mean reservoir pressure p1, by the matter balance equation iterative calculation in non-fracturing reform area
Non- fracturing reform area mean reservoir pressure p2；
According to fracturing reform area mean reservoir pressure p1, non-fracturing reform area mean reservoir pressure p2And interval channelling
Coefficient calculates current tiTime step non-fracturing reform area is to the channelling amount in fracturing reform area and accumulative channelling amount；
Based on fracturing reform area mean reservoir pressure value p1, original formation pressure pi, gas well cumulative gas Gp1、
With described accumulative channelling amount Gp2Calculate the matter balance equation residual error in current time step fracturing reform area：
When the absolute value of residual error is less than assigned error, iteration convergence, iteration exits；Otherwise with current average
Value formation pressure p1As initial value, continue the fifty-fifty lamination in the new fracturing reform area of iterative calculation according to Newton method
Force value is till meeting error requirements.
According to one embodiment of present invention, the step iterating to calculate the mean reservoir pressure in non-fracturing reform area is entered
One step includes：
Take time step ti-1Non- fracturing reform area mean reservoir pressure p20As this time step tiNon- fracturing reform area
Mean reservoir pressure initial value p2If, ti=0, then it is taken as original formation pressure；
According to current tiCurrent mean reservoir pressure value p in time step fracturing reform area and non-fracturing reform area1And p2,
Calculate channelling amount q between Liang Ge area2And accumulative channelling amount Gp2,
Channelling amount q between current time step Xia Liangge area2Calculated by following formula：
Always accumulative channelling amount G between Liang Ge areap2Calculated by following formula：
Gp2=Gp20+ΔGp2,
Wherein, λ is interporosity flow coefficient, Δ Gp2For the accumulative channelling amount between current step Nei Liangge area,Δ t is time step step-length, q20For the channelling amount between last time step Liang Ge area, Gp20For
Always accumulative channelling amount between the Liang Ge area of last time step；
Based on current tiTime step non-fracturing reform area mean reservoir pressure value p2, original formation pressure piTire out with interval
Meter channelling yield Gp2Calculate the matter balance equation residual error in non-fracturing reform area,
When the absolute value of residual error is less than assigned error, exit iteration, otherwise with current mean reservoir pressure value
p2As initial value, the mean reservoir pressure value continuing to iterate to calculate non-fracturing reform area until meeting error requirements is
Only.
According to one embodiment of present invention, the flowing bottomhole pressure (FBHP) of matching prediction and described actual flowing bottomhole pressure (FBHP) are to determine
The step of the dynamic holdup in the dynamic holdup in fracturing reform area described in horizontal well and described non-fracturing reform area enters one
Step includes：
Mean reservoir pressure based on described fracturing reform area and described binomial potential curve and equation calculate flowing bottomhole pressure (FBHP)；
Adjust the dynamic holdup G in described fracturing reform area by Optimizing Algorithm for Fitting1With described non-fracturing reform area
Dynamic holdup G2, so that calculated flowing bottomhole pressure (FBHP) is being set with the square-error value preset of actual flowing bottomhole pressure (FBHP)
Determine in scope, so that it is determined that the dynamic holdup G in described fracturing reform area1Dynamic with described non-fracturing reform area
Reserves G2.
Beneficial effects of the present invention：
The present invention does not need closing well test mean reservoir pressure in use it is possible to consider adsorbed gas desorbing and different
Often the impact of high pressure and physical property recombination region, can determine that fracturing reform area and periphery non-fracturing reform area are dynamic respectively
Reserves.
Other features and advantages of the present invention will illustrate in the following description, and, partly from description
In become apparent, or by implement the present invention and understand.The purpose of the present invention and other advantages can be passed through
In description, claims and accompanying drawing, specifically noted structure is realizing and to obtain.
Brief description
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment
Or required accompanying drawing does simple introduction in description of the prior art：
Fig. 1 is shale gas multistage pressure break horizontal well Area of a well schematic diagram；
Fig. 2 is method flow diagram according to an embodiment of the invention；
Fig. 3 is algorithm flowchart according to an embodiment of the invention；
Fig. 4 is the flow chart of iterative inner region stratum according to an embodiment of the invention average pressure；
Fig. 5 is the flow chart of iterative outskirt stratum according to an embodiment of the invention average pressure；
Fig. 6 is daily output tolerance and the test flowing bottomhole pressure (FBHP) schematic diagram of well A according to an embodiment of the invention；
Fig. 7 is the binomial deliverability curve schematic diagram of well A according to an embodiment of the invention；
Fig. 8 be well A according to an embodiment of the invention history matching before calculate flowing bottomhole pressure (FBHP) and measured value
Correlation curve schematic diagram；And
Fig. 9 be well A according to an embodiment of the invention history matching after calculate flowing bottomhole pressure (FBHP) and measured value
Correlation curve schematic diagram.
Specific embodiment
To describe embodiments of the present invention below with reference to drawings and Examples in detail, whereby to the present invention how
Application technology means are solving technical problem, and reach realizing process and fully understanding and real according to this of technique effect
Apply.As long as it should be noted that not constituting conflict, in each embodiment in the present invention and each embodiment
Each feature can be combined with each other, and the technical scheme being formed is all within protection scope of the present invention.
Shale gas multistage pressure break horizontal well Area of a well can be divided into fracturing reform area according to well pattern spacing
(inner region) and non-fracturing reform area (outskirt) two parts.It is illustrated in figure 1 a horizontal well Area of a well
Schematic diagram, its control area includes fracturing reform area 100 and non-fracturing reform area 200, and the boundary in twoth area is pressure
Split transformation area border, the outer boundary in non-fracturing reform area 200 is choked flow border between well.Inner region is through multistage
The complex fracture network area being made up of main fracturing fracture and dry being formed after volume fracturing transformation, outskirt
For between adjacent well not by the region of fracturing reform.
Shale gas pressure break horizontal well production phase early stage aerogenesis, essentially from inner region, drops with inner region strata pressure
Low, the pressure reduction between interior outskirt gradually amplifies and forms interval channelling.Channelling yield is poor with interior outskirt strata pressure
And interval interporosity flow coefficient is relevant.After inner region drop in formation pressure is to critical desorption pressures, adsorbed gas starts desorbing
Diffusion simultaneously produces impact to gas well deliverability.
Therefore, the invention provides a kind of for determining area's dynamic holdup inside and outside shale gas multistage pressure break horizontal well
Method.The present invention according to shale gas pressure break horizontal well control area have the characteristics that in outskirt, by Liang Ge area respectively
Set up and consider adsorbed gas desorbing, the shale gas reservoir matter balance equation of abnormal pressure impact, and combine the two of gas well
Channelling equation between item formula deliverability equation and Liang Ge area, sets up a kind of shale gas well gas production, flowing bottomhole pressure (FBHP)
Computation model, by adjusting outskirt Reserve Fitting flowing bottomhole pressure (FBHP) in gas well, inside and outside final determination shale gas well, area is moved
State reserves.
It is illustrated in figure 2 method flow diagram according to an embodiment of the invention, be according to this as shown in Figure 3
The algorithm flowchart of one embodiment of invention, the present invention is carried out in detail below with reference to Fig. 2 and Fig. 3
Explanation.
First, in step s 110, shale gas reservoir gas testing data, creation data and gas well structural parameters are obtained,
And gas well binomial potential curve and equation is set up by gas testing data, calculated real by creation data and gas well structural parameters
Border flowing bottomhole pressure (FBHP).
Specifically, the shale gas reservoir gas testing data of acquisition, creation data and gas well structural parameters are as shown in table 1.
According to People's Republic of China (PRC) oil and gas industry standard SY/T 5440-2000, built by matching gas testing data
Vertical gas well binomial potential curve and equation is expressed as follows：
Wherein, A is gas well binomial potential curve and equation Monomial coefficient；B is that gas well binomial potential curve and equation is secondary
Term coefficient；Q is gas well daily output tolerance, 104m3/d；piFor original formation pressure, pwfFor gas well in stable yields yield
Stablize flowing bottomhole pressure (FBHP), MPa for corresponding during q.From formula (1-1), if given mean reservoir pressure and
It is possible to be calculated the prediction flowing bottomhole pressure (FBHP) of gas well according to this formula after gas well daily output tolerance.Fifty-fifty it is laminated
Power can be calculated in below step of the present invention.
Next, in the step s 120 it is considered to adsorbed gas desorbing and abnormal pressure impact, setting up respectively and be directed to
The fracturing reform area of multistage pressure break horizontal well and the matter balance equation in non-fracturing reform area.This step specifically includes
Following several steps.
It is considered to adsorbed gas desorbing and abnormal pressure impact in step S1201, set up shale gas reservoir material balance
Equation.In this step, when setting up shale gas reservoir matter balance equation, need to consider free gas, absorption
The impact of gas desorbing diffusion.Additionally, for abnormal pressure shale gas reservoir in addition it is also necessary to consider that rock elasticity can affect.
Because Abnormal High Pressure Gas Reservoirs initial stage of development pore space compressibility of rock approximate with gas compressive coefficient same
The individual order of magnitude, its flexible drive can be can not ignore, and the dynamic holdup otherwise evaluated is higher.Different below by way of analyzing
The often high pressure and adsorbed gas desorbing impact to shale gas reservoir matter balance equation, according to the principle of mass conservation, sets up
Consider abnormal pressure and the shale gas reservoir matter balance equation of adsorbed gas desorbing impact.
Hypothesis gas reservoir original free gas reserves are Gf, when strata pressure is by PiWhen being reduced to P, due to rock matrix
Compression and the underground pore volume decrement that causes of fluid expansion are expressed as：
Wherein, GfFor gas reservoir original free gas reserves, 104m3；BgiFor the natural-gas under original formation pressure
Volume factor；SgiFor the original gas saturation of gas reservoir；SwiFor the original water saturation of gas reservoir；cfFor active porosity
Coefficient of compressibility, MPa-1；cwFor stratum water coefficient of compressibility, MPa-1；ΔVeDeclined by strata pressure and cause
Underground pore volume change, 104m3.
Due to when x → 0, ex≈ 1+x, formula (2-1) can be reduced to：
Then formula (2-2) represents that Abnormal High Pressure Gas Reservoirs strata pressure declines the underground pore volume change causing.
Because the natural gas in shale gas reservoir is with free gas and adsorbed gas mode preservation, wherein adsorbed gas accounting reaches
20～85%.Shale adsorbed gas air content can be described with following Langmuir isothermal adsorpting equation：
Wherein, V (p) represents that pressure is shale saturation adsorbed gas air content during p, m3/t；VLFor Langmuir body
Long-pending, represent the maximum saturation adsorbed gas air content of the shale when strata pressure tends to infinity, m3/t；PLFor orchid
Family name's pressure, represents 50%V in Lan Shi adsorption isothermal curveLCorresponding pressure, MPa.Lan Shi pressure is lower, then
Adsorbed gas is less susceptible to desorbing in recovery process.
Based on formula (2-3), under any pressure, the total absorption tolerance of shale gas reservoir can be expressed as：
Wherein, ρBFor shale density, t/m3；VBFor shale cumulative volume, m3.
Based on formula (2-4), original formation pressure is piWhen free gas and adsorbed gas gross reserves be represented by：
When strata pressure is by piBeing reduced to tired gas production during p is Gp, shale gas reservoir pore volume decrement can lead to
Cross formula (2-2) to be calculated：
Residual free gas reserves can be calculated by following formula：
Remaining adsorbed gas reserves：
According to the law of indestructibility of matter：Original free gas reserves+original adsorbed gas reserves=residual free gas reserves+residue is inhaled
Attached gas yield+tired gas production, can obtain：
By the B in formula (2-9)gWith the conversion of compressibility factor z, then formula (2-9) can be arranged and be：
As order：
Then formula (2-10) can arrange and be：
Formula (2-12) is the shale gas reservoir matter balance equation considering adsorbed gas desorbing and abnormal pressure impact.
In step S1202, the Area of a well of multistage pressure break horizontal well is divided into by pressure according to well pattern spacing
Split transformation area and non-fracturing reform area, its subregion is as shown in Figure 1.
In step S1203, set up respectively for fracturing reform area based on shale gas reservoir matter balance equation and not
The matter balance equation in fracturing reform area.Specifically, changed based on the pressure break that matter balance equation (2-12) is set up
The matter balance equation making area is：
Wherein, p1For the mean reservoir pressure in fracturing reform area, G1For the dynamic holdup in fracturing reform area, Gp1For
Gas well cumulative production, Gp2For accumulative channelling amount from non-fracturing reform area to fracturing reform area.
The matter balance equation in the non-fracturing reform area being set up based on matter balance equation (2-12) is：
Wherein, p2For the mean reservoir pressure in non-fracturing reform area, G2Dynamic holdup for non-fracturing reform area.
Finally, in step s 130, give the dynamic holdup in fracturing reform area and non-fracturing reform area respectively,
And the channelling equation between the matter balance equation with reference to fracturing reform area and non-fracturing reform area and twoth area, binomial
Formula deliverability equation adjustment dynamic holdup carrys out the actual flowing bottomhole pressure (FBHP) of matching, so that it is determined that area is dynamic inside and outside pressure break horizontal well
Reserves.
In this step, determine that in horizontal well, outskirt dynamic holdup specifically includes following several steps.First,
In step S1301, given shale gas well fracturing reform area dynamic holdup is G1, non-fracturing reform area dynamically stores up
Measure as G2, this two dynamic holdups are estimated value, are not end value.
In step S1302, based on fracturing reform area dynamic holdup G1, non-fracturing reform area dynamic holdup G2、
The matter balance equation in fracturing reform area, the matter balance equation in non-fracturing reform area, Gas Well Productivity and double
Channelling equation between area, produces, according to gas well, the mean reservoir pressure that tolerance iterates to calculate fracturing reform area daily.
Specifically, the mean reservoir pressure in iterative calculation fracturing reform area can be by the several steps shown in Fig. 4
Realize.First, by upper time step ti-1Mean reservoir pressure value p in interior fracturing reform area10As current time step
tiInterior fracturing reform area mean reservoir pressure p1Iterative initial value.If ti=0, iterative initial value is then primitively to be laminated
Force value pi.
Then, according to fracturing reform area mean reservoir pressure p1, changed by the matter balance equation in non-fracturing reform area
In generation, calculates non-fracturing reform area mean reservoir pressure p2.
Iterate to calculate non-fracturing reform area mean reservoir pressure p2When include the several steps shown in Fig. 5.First,
Take time step ti-1Nei Wei fracturing reform area mean reservoir pressure p20As this time step tiNei Wei fracturing reform area
Mean reservoir pressure initial value p2If, ti=0, iterative initial value is then taken as original formation pressure pi.
Then, based on current tiThe current mean reservoir pressure in fracturing reform area in time step and non-fracturing reform area
Value p1And p2, channelling amount q between current time step Xia Liangge area is calculated according to following formula2And accumulative channelling amount Gp2.
Wherein, calculate channelling amount q between current time step Xia Liangge area2Expression formula：
Calculate always accumulative channelling amount G between Liang Ge areap2Expression formula be：
Gp2=Gp20+ΔGp2, (2-16)
Wherein,λ is interporosity flow coefficient, and Δ t is time step step-length, q20For last time step two
Channelling amount between individual area, Gp20For the always accumulative channelling amount between the Liang Ge area of last time step, Δ Gp2For current
Accumulative channelling amount between step-length Nei Liangge area.
Then, based on current tiTime step Nei Wei fracturing reform area mean reservoir pressure value p2, original formation pressure pi
With interval accumulative channelling yield Gp2Calculate the matter balance equation residual error in non-fracturing reform area.This residual error is：
When the absolute value of residual error is less than assigned error, exit iteration；Otherwise, with current mean reservoir pressure
Value p2As iterative initial value, continue to iterate to calculate the mean reservoir pressure value in non-fracturing reform area, until this is average
Till value formation pressure meets error requirements.
Mean reservoir pressure value p in non-fracturing reform area can be obtained based on above calculating2.Next, according to pressure
Split transformation area mean reservoir pressure p1, non-fracturing reform area mean reservoir pressure p2And interval interporosity flow coefficient calculates not
Fracturing reform area is to the current time step channelling amount in fracturing reform area and accumulative channelling amount.Current time step channelling amount herein
And accumulative channelling amount can using ask for above non-fracturing reform area mean reservoir pressure value when calculated
Value.
Next, being based on fracturing reform area mean reservoir pressure value p1, original formation pressure pi, gas well add up produce
Tolerance Gp1With accumulative channelling amount Gp2Calculate the matter balance equation residual error in fracturing reform area in current time step.This is residual
Difference is：
When the absolute value of residual error is less than assigned error, iteration convergence, iteration exits；Otherwise, put down with current
All value formation pressure p1As iterative initial value, continue the average of the new fracturing reform area of iterative calculation according to Newton method
Value formation pressure, till this mean reservoir pressure value meets error requirements.
Mean reservoir pressure value p in step S1303, based on the above fracturing reform area obtaining1And binomial
Deliverability equation, by daily output tolerance prediction flowing bottomhole pressure (FBHP).Specifically, under other specification known case, as long as giving
Determine fracturing reform area dynamic holdup G1Non- fracturing reform area dynamic holdup G2It is possible to be obtained by creation data
Take gas well daily output tolerance, produce tolerance daily further according to gas well and deliverability equation calculates corresponding flowing bottomhole pressure (FBHP).
In step S1304, the flowing bottomhole pressure (FBHP) of matching prediction and actual flowing bottomhole pressure (FBHP) are changed with determining fractured horizontal well
Make the dynamic holdup in area and the dynamic holdup in non-fracturing reform area.Specifically, in this step, changed based on pressure break
Make the mean reservoir pressure in area and binomial potential curve and equation formula (1-1) calculates flowing bottomhole pressure (FBHP), calculated by matching optimization
Method adjusts the dynamic holdup G in fracturing reform area1The dynamic holdup G in non-fracturing reform area2, so that calculating
The square-error value preset of the flowing bottomhole pressure (FBHP) arriving and actual flowing bottomhole pressure (FBHP) is in set point, so that it is determined that fracturing reform
The dynamic holdup G in area1The dynamic holdup G in non-fracturing reform area2.
Table 1
Below by way of a specific embodiment, the present invention will be described.This sentences Sichuan Basin Lower Silurian Series
As a example Longma small stream group shale gas well A, this well major parameter is shown in Table 1.
This well adopts 5 cun of half-covering tube pilot productions 9 months after terminating from December, 2012 gas testing, then enters down 2 cun
Half oil pipe produces, and period adds up to carry out 9 bottom pressures and thermograde test altogether.Fig. 6 is that this well gas testing tries
Daily output tolerance and test flowing bottomhole pressure (FBHP) during adopting.
Before dynamic holdup calculates, first according to this well binomial potential curve and equation of gas testing data evaluation.According to middle Chinese
People republic oil and gas industry standard SY/T 5440-2000, sets up gas well binomial by matching gas testing data
Formula deliverability equation, result is as shown in Figure 7.Wherein, well A binomial potential curve and equation coefficient is respectively A=21.107,
B=2.9494.
According to method of the present invention, set up algorithm routine, input above-mentioned table 1 parameter and Fig. 6 daily output tolerance
Deng creation data, and give inside and outside dynamic holdup G respectively1、G2Initial value, G1=2.5 hundred million sides, G2=3 hundred million sides.
The flowing bottomhole pressure (FBHP) calculating is as shown in Figure 8.Dynamic holdup G is adjusted by Optimizing Algorithm for Fitting1And G2, matching meter
The flowing bottomhole pressure (FBHP) calculated and actual flowing bottomhole pressure (FBHP), final fitting result is as shown in Figure 9.The well A's of matching is interior at present
The dynamic holdup in area is 1.83 hundred million sides, and the mining-employed reserves of outskirt is 0.4 hundred million sides, then the development degree of outskirt is low.
The matter balance equation set up in the present invention considers adsorbed gas desorbing and pore space compressibility of rock impact,
And comprehensive matter balance equation and quasi-stable state Gas Well Productivity establish based on Production development number on this basis
According to dynamic holdup evaluation methodology.The method can consider shale gas adsorbed gas desorbing and blowhole when calculating
The impact to shale gas reservoir matter balance equation for the coefficient of compressibility, and do not need shut-in pressure to recover to comment in the calculation
The mean reservoir pressure of valency.Shale gas well rational proration, developing value and exploitation that result of calculation is suitable for
Multiple application such as scheme optimization.
While it is disclosed that embodiment as above, but described content is only to facilitate understand the present invention
And the embodiment adopting, it is not limited to the present invention.Technology people in any the technical field of the invention
Member, on the premise of without departing from spirit and scope disclosed in this invention, can be in the formal and details implemented
On make any modification and change, but the scope of patent protection of the present invention, still must be with appending claims institute
The scope defining is defined.
Claims (10)
1. a kind of method for determining area's dynamic holdup inside and outside shale gas multistage pressure break horizontal well, including：
Obtain shale gas reservoir gas testing data, creation data and gas well structural parameters, and built by described gas testing data
Vertical gas well binomial potential curve and equation, calculates actual flowing bottomhole pressure (FBHP) by described creation data and gas well structural parameters；
Consider adsorbed gas desorbing and abnormal pressure impact, set up the fracturing reform for multistage pressure break horizontal well respectively
Area and the matter balance equation in non-fracturing reform area；
Give the dynamic holdup in described fracturing reform area and described non-fracturing reform area respectively, and combine described pressure break
Channelling equation between the matter balance equation in transformation area and described non-fracturing reform area and twoth area, described binomial
Deliverability equation adjusts described dynamic holdup and carrys out actual flowing bottomhole pressure (FBHP) described in matching, so that it is determined that inside and outside pressure break horizontal well
Area's dynamic holdup.
2. method according to claim 1 is it is characterised in that set up for multistage pressure break horizontal well
The step of the matter balance equation in fracturing reform area and non-fracturing reform area further includes：
Consider adsorbed gas desorbing and abnormal pressure impact, set up shale gas reservoir matter balance equation；
Multistage pressure break horizontal well Area of a well is divided into by fracturing reform area and non-pressure break according to well pattern spacing
Transformation area；
Described fracturing reform area is set up respectively based on described shale gas reservoir matter balance equation and described non-pressure break changes
Make the matter balance equation in area.
3. method according to claim 2 is it is characterised in that set up described shale gas reservoir material balance
The step of equation further includes：
Based on pore space compressibility of rock, calculate during strata pressure change due to rock matrix compression and fluid expansion
The underground pore volume decrement causing；
Drop in formation pressure is calculated based on described underground pore volume decrement and Langmuir isothermal adsorpting equation
Shale gas reservoir residual free gas reserves afterwards and remaining adsorbed gas reserves；
According to the law of indestructibility of matter：Original free gas reserves+original adsorbed gas reserves=residual free gas reserves+surplus
Remaining adsorbed gas yield+tired gas production, sets up described shale gas reservoir matter balance equation.
4. according to the method in claim 2 or 3 it is characterised in that described shale gas reservoir material balance
Equation is：
Wherein, Za=Za(p),
Work as p=piWhen,P is strata pressure, piFor original formation pressure, GpFor cumulative production,
G is dynamic holdup, cfFor pore space compressibility of rock, cwFor stratum water coefficient of compressibility, SwiContain for gas reservoir is original
Water saturation, SgiFor the original gas saturation of gas reservoir, φ is effecive porosity, ρBFor shale density, VLFor orchid
Ge Miaoer volume, PLFor Lan Shi pressure, pscFor gas pressure under standard state, T is temperature, TscFor standard shape
Temperature under state, z is compressibility factor, zscFor the compressibility factor under standard state.
5. method according to claim 4 is it is characterised in that the material balance in described fracturing reform area
Equation is：
Wherein, p1For the mean reservoir pressure in fracturing reform area,G1Dynamic storage for fracturing reform area
Amount, Gp1For gas well cumulative production, Gp2For accumulative channelling amount from non-fracturing reform area to fracturing reform area.
6. method according to claim 4 is it is characterised in that the material in described non-fracturing reform area is put down
Weighing apparatus equation be：
Wherein, p2For the mean reservoir pressure in non-fracturing reform area,Gp2For non-fracturing reform area to
The accumulative channelling amount in fracturing reform area, G2Dynamic holdup for non-fracturing reform area.
7. the method according to claim 5 or 6 is it is characterised in that determine area inside and outside pressure break horizontal well
The step of dynamic holdup further includes：
The fracturing reform area dynamic holdup G of given shale gas well1, non-fracturing reform area dynamic holdup G2；
Based on fracturing reform area dynamic holdup G1, non-fracturing reform area dynamic holdup G2, the thing in fracturing reform area
Channelling side between matter equilibrium equation, the matter balance equation in non-fracturing reform area, Gas Well Productivity and two-region
Journey, produces tolerance daily according to gas well and iterates to calculate fracturing reform area fifty-fifty by fracturing reform area matter balance equation
Stressor layer；
Mean reservoir pressure p based on described fracturing reform area1With described binomial potential curve and equation, by daily output tolerance
Prediction flowing bottomhole pressure (FBHP)；
The flowing bottomhole pressure (FBHP) of matching prediction and described actual flowing bottomhole pressure (FBHP) are to determine fracturing reform area described in horizontal well
Dynamic holdup and the dynamic holdup in described non-fracturing reform area.
8. method according to claim 7 is it is characterised in that iterate to calculate the average of fracturing reform area
The step of strata pressure further includes：
By upper time step ti-1Fracturing reform area mean reservoir pressure value p10As current tiTime step fracturing reform area
The iterative initial value p of mean reservoir pressure1If, ti=0, then iterative initial value is original formation pressure value；
According to fracturing reform area mean reservoir pressure p1, by the matter balance equation iterative calculation in non-fracturing reform area
Non- fracturing reform area mean reservoir pressure p2；
According to fracturing reform area mean reservoir pressure p1, non-fracturing reform area mean reservoir pressure p2And interval channelling
Coefficient calculates current tiTime step non-fracturing reform area is to the channelling amount in fracturing reform area and accumulative channelling amount；
Based on fracturing reform area mean reservoir pressure value p1, original formation pressure pi, gas well cumulative gas Gp1、
With described accumulative channelling amount Gp2Calculate the matter balance equation residual error in current time step fracturing reform area：
When the absolute value of residual error is less than assigned error, iteration convergence, iteration exits, otherwise with current average
Value formation pressure p1As initial value, continue the fifty-fifty lamination in the new fracturing reform area of iterative calculation according to Newton method
Force value is till meeting error requirements.
9. method according to claim 8 is it is characterised in that iterate to calculate the flat of non-fracturing reform area
The step of all strata pressures further includes：
Take time step ti-1Non- fracturing reform area mean reservoir pressure p20As this time step tiNon- fracturing reform area
Mean reservoir pressure initial value p2If, ti=0, then it is taken as original formation pressure；
According to current tiCurrent mean reservoir pressure value p in time step fracturing reform area and non-fracturing reform area1And p2,
Calculate channelling amount q between Liang Ge area2And accumulative channelling amount Gp2,
Channelling amount q between current time step Xia Liangge area2Calculated by following formula：
Always accumulative channelling amount G between Liang Ge areap2Calculated by following formula：
Gp2=Gp20+ΔGp2,
Wherein, λ is interporosity flow coefficient, Δ Gp2For the accumulative channelling amount between current step Nei Liangge area,Δ t is time step step-length, q20For the channelling amount between last time step Liang Ge area, Gp20For
Always accumulative channelling amount between the Liang Ge area of last time step；
Based on current tiTime step non-fracturing reform area mean reservoir pressure value p2, original formation pressure piTire out with interval
Meter channelling yield Gp2Calculate the matter balance equation residual error in non-fracturing reform area：
When the absolute value of residual error is less than assigned error, exit iteration, otherwise with current mean reservoir pressure value
p2As initial value, the mean reservoir pressure value continuing to iterate to calculate non-fracturing reform area until meeting error requirements is
Only.
10. method according to claim 9 is it is characterised in that the flowing bottomhole pressure (FBHP) and described of matching prediction
Actual flowing bottomhole pressure (FBHP) is to determine the dynamic holdup in fracturing reform area described in horizontal well and described non-fracturing reform area
The step of dynamic holdup further includes：
Mean reservoir pressure based on described fracturing reform area and described binomial potential curve and equation calculate flowing bottomhole pressure (FBHP)；
Adjust the dynamic holdup G in described fracturing reform area by Optimizing Algorithm for Fitting1With described non-fracturing reform area
Dynamic holdup G2, so that calculated flowing bottomhole pressure (FBHP) is being set with the square-error value preset of actual flowing bottomhole pressure (FBHP)
Determine in scope, so that it is determined that the dynamic holdup G in described fracturing reform area1Dynamic with described non-fracturing reform area
Reserves G2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510549414.4A CN106481332B (en) | 2015-08-31 | 2015-08-31 | Method for determining area's dynamic holdup inside and outside shale gas multistage pressure break horizontal well |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510549414.4A CN106481332B (en) | 2015-08-31 | 2015-08-31 | Method for determining area's dynamic holdup inside and outside shale gas multistage pressure break horizontal well |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106481332A true CN106481332A (en) | 2017-03-08 |
| CN106481332B CN106481332B (en) | 2019-06-18 |
Family
ID=58235493
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510549414.4A Active CN106481332B (en) | 2015-08-31 | 2015-08-31 | Method for determining area's dynamic holdup inside and outside shale gas multistage pressure break horizontal well |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106481332B (en) |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107806916A (en) * | 2017-09-11 | 2018-03-16 | 中国石油天然气股份有限公司 | The determination method and apparatus of reservoir reconstruction volume |
| CN108166963A (en) * | 2017-12-13 | 2018-06-15 | 中国海洋石油集团有限公司 | A kind of offshore oil gas well evaluation of Fracturing Effect on Compact Sandstone method |
| CN108612525A (en) * | 2018-04-19 | 2018-10-02 | 重庆科技学院 | A kind of gas reservoir protection Reserve Estimation Method |
| CN108694254A (en) * | 2017-04-06 | 2018-10-23 | 中国石油化工股份有限公司 | One kind selling of one's property transformation production gas well yield experience decline curve analysis method |
| CN108804819A (en) * | 2018-06-10 | 2018-11-13 | 西南石油大学 | A kind of low permeability gas reservoirs dynamic holdup evaluation method |
| CN110656915A (en) * | 2018-06-13 | 2020-01-07 | 中国石油化工股份有限公司 | Shale gas multi-section fracturing horizontal well multi-working-system productivity prediction method |
| CN110778306A (en) * | 2019-10-22 | 2020-02-11 | 中国海洋石油集团有限公司 | Method for correcting abnormity of gas well productivity equation |
| CN110863818A (en) * | 2018-08-08 | 2020-03-06 | 中国石油化工股份有限公司 | Description method and device for residual oil/gas distribution |
| CN111927421A (en) * | 2020-08-18 | 2020-11-13 | 中国石油化工股份有限公司 | Repeated fracturing well selection method based on shale gas unstable linear flow theory |
| CN112065339A (en) * | 2020-09-02 | 2020-12-11 | 中国石油大学(北京) | Multi-reservoir gas reservoir commingled production capacity prediction method and device |
| CN112434426A (en) * | 2020-11-26 | 2021-03-02 | 北京科技大学 | Shale gas multistage fracturing horizontal well step gradient pressure drop development method and device |
| CN112528455A (en) * | 2019-09-03 | 2021-03-19 | 中国石油天然气股份有限公司 | Method and system for calculating reservoir pressure and dynamic reserve of adsorption unsaturated coal-bed gas well |
| CN112878987A (en) * | 2021-01-25 | 2021-06-01 | 长江大学 | Method for calculating shale gas well control reserves by using production data |
| CN112966422A (en) * | 2021-03-18 | 2021-06-15 | 中国石油大学(华东) | Flow simulation method suitable for shale gas reservoir thermal recovery |
| CN116050629A (en) * | 2023-01-18 | 2023-05-02 | 重庆科技学院 | Dynamic prediction method for storage capacity of gas storage in consideration of formation water evaporation salting-out |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2468184B (en) * | 2009-02-25 | 2011-08-24 | Logined Bv | Modeling a reservoir using a compartment model and a geomechanical model |
| CN103413030A (en) * | 2013-07-24 | 2013-11-27 | 中国石油天然气股份有限公司 | Dynamic analysis method and system of fracture-vug type carbonatite gas reservoir |
| CN104481520A (en) * | 2014-10-13 | 2015-04-01 | 成都创源油气技术开发有限公司 | Early evaluation method of recoverable reserves of shale gas wells |
| CN104632187A (en) * | 2013-11-14 | 2015-05-20 | 中国石油化工股份有限公司 | Method for determining dynamic reserve volume of water production coal seam gas well |
-
2015
- 2015-08-31 CN CN201510549414.4A patent/CN106481332B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2468184B (en) * | 2009-02-25 | 2011-08-24 | Logined Bv | Modeling a reservoir using a compartment model and a geomechanical model |
| CN103413030A (en) * | 2013-07-24 | 2013-11-27 | 中国石油天然气股份有限公司 | Dynamic analysis method and system of fracture-vug type carbonatite gas reservoir |
| CN104632187A (en) * | 2013-11-14 | 2015-05-20 | 中国石油化工股份有限公司 | Method for determining dynamic reserve volume of water production coal seam gas well |
| CN104481520A (en) * | 2014-10-13 | 2015-04-01 | 成都创源油气技术开发有限公司 | Early evaluation method of recoverable reserves of shale gas wells |
Non-Patent Citations (6)
| Title |
|---|
| 刘欢: "物质平衡方法计算页岩气储量", 《石油工程技术》 * |
| 张烈辉: "改进的页岩气藏物质平衡方程及储量计算方法", 《天然气工程》 * |
| 李艳丽: "页岩气储量计算方法探讨", 《天然气地球科学》 * |
| 王卫红 等: "超高压气藏渗流机理及气井生产动态特征", 《天然气地球科学》 * |
| 王德龙 等: "新吸附气藏物质平衡方程推导及储量计算", 《岩性油气藏》 * |
| 钟海全 等: "流动物质平衡法计算低渗透气藏单井动态储量", 《岩性油气藏》 * |
Cited By (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108694254A (en) * | 2017-04-06 | 2018-10-23 | 中国石油化工股份有限公司 | One kind selling of one's property transformation production gas well yield experience decline curve analysis method |
| CN108694254B (en) * | 2017-04-06 | 2021-10-08 | 中国石油化工股份有限公司 | Empirical decreasing curve analysis method for production of variable-production variable-pressure production gas well |
| CN107806916B (en) * | 2017-09-11 | 2019-12-10 | 中国石油天然气股份有限公司 | Method and device for determining reservoir reconstruction volume |
| CN107806916A (en) * | 2017-09-11 | 2018-03-16 | 中国石油天然气股份有限公司 | The determination method and apparatus of reservoir reconstruction volume |
| CN108166963B (en) * | 2017-12-13 | 2020-02-14 | 中国海洋石油集团有限公司 | Method for evaluating fracturing effect of offshore oil and gas well |
| CN108166963A (en) * | 2017-12-13 | 2018-06-15 | 中国海洋石油集团有限公司 | A kind of offshore oil gas well evaluation of Fracturing Effect on Compact Sandstone method |
| CN108612525A (en) * | 2018-04-19 | 2018-10-02 | 重庆科技学院 | A kind of gas reservoir protection Reserve Estimation Method |
| CN108804819A (en) * | 2018-06-10 | 2018-11-13 | 西南石油大学 | A kind of low permeability gas reservoirs dynamic holdup evaluation method |
| CN110656915A (en) * | 2018-06-13 | 2020-01-07 | 中国石油化工股份有限公司 | Shale gas multi-section fracturing horizontal well multi-working-system productivity prediction method |
| CN110656915B (en) * | 2018-06-13 | 2021-10-08 | 中国石油化工股份有限公司 | Shale gas multi-section fracturing horizontal well multi-working-system productivity prediction method |
| CN110863818A (en) * | 2018-08-08 | 2020-03-06 | 中国石油化工股份有限公司 | Description method and device for residual oil/gas distribution |
| CN110863818B (en) * | 2018-08-08 | 2023-08-29 | 中国石油化工股份有限公司 | Description method and device for residual oil/gas distribution |
| CN112528455B (en) * | 2019-09-03 | 2022-10-04 | 中国石油天然气股份有限公司 | Method and system for calculating reservoir pressure and dynamic reserve of adsorption unsaturated coal-bed gas well |
| CN112528455A (en) * | 2019-09-03 | 2021-03-19 | 中国石油天然气股份有限公司 | Method and system for calculating reservoir pressure and dynamic reserve of adsorption unsaturated coal-bed gas well |
| CN110778306A (en) * | 2019-10-22 | 2020-02-11 | 中国海洋石油集团有限公司 | Method for correcting abnormity of gas well productivity equation |
| CN111927421A (en) * | 2020-08-18 | 2020-11-13 | 中国石油化工股份有限公司 | Repeated fracturing well selection method based on shale gas unstable linear flow theory |
| CN112065339B (en) * | 2020-09-02 | 2021-10-26 | 中国石油大学(北京) | Multi-reservoir gas reservoir commingled production capacity prediction method and device |
| CN112065339A (en) * | 2020-09-02 | 2020-12-11 | 中国石油大学(北京) | Multi-reservoir gas reservoir commingled production capacity prediction method and device |
| CN112434426B (en) * | 2020-11-26 | 2022-09-27 | 北京科技大学 | Shale gas multistage fracturing horizontal well step gradient pressure drop development method and device |
| CN112434426A (en) * | 2020-11-26 | 2021-03-02 | 北京科技大学 | Shale gas multistage fracturing horizontal well step gradient pressure drop development method and device |
| CN112878987A (en) * | 2021-01-25 | 2021-06-01 | 长江大学 | Method for calculating shale gas well control reserves by using production data |
| CN112878987B (en) * | 2021-01-25 | 2024-01-19 | 长江大学 | Method for calculating shale gas well control reserves by using production data |
| CN112966422A (en) * | 2021-03-18 | 2021-06-15 | 中国石油大学(华东) | Flow simulation method suitable for shale gas reservoir thermal recovery |
| CN116050629A (en) * | 2023-01-18 | 2023-05-02 | 重庆科技学院 | Dynamic prediction method for storage capacity of gas storage in consideration of formation water evaporation salting-out |
| CN116050629B (en) * | 2023-01-18 | 2023-09-12 | 重庆科技学院 | Dynamic prediction method for storage capacity of gas storage in consideration of formation water evaporation salting-out |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106481332B (en) | 2019-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106481332A (en) | Method for determining area's dynamic holdup inside and outside shale gas multistage pressure break horizontal well | |
| CN106484933B (en) | It is a kind of for determining the method and system of shale gas well well control dynamic holdup | |
| CN104948163B (en) | A kind of shale gas well deliverability assay method | |
| CN104632187B (en) | A kind of method for determining production water coal bed gas well dynamic holdup | |
| CN106351651B (en) | The prediction technique and device of shale gas well deliverability | |
| Meng et al. | A mathematical model for gas and water production from overlapping fractured coalbed methane and tight gas reservoirs | |
| Xu et al. | A dynamic prediction model for gas–water effective permeability based on coalbed methane production data | |
| CN104612635B (en) | Coal seam group gas joint extraction pre- evaluation method up to standard | |
| Mazumder et al. | Permeability increase in Bowen Basin coal as a result of matrix shrinkage during primary depletion | |
| CN107301306A (en) | Dynamic open-flow capacity Forecasting Methodology for DAMAGE OF TIGHT SAND GAS RESERVOIRS pressure break horizontal well | |
| Liu et al. | Experimental evaluation of interlayer interference during commingled production in a tight sandstone gas reservoir with multi-pressure systems | |
| CN104750896B (en) | A kind of fractured-cavernous carbonate reservoir method for numerical simulation | |
| Xu et al. | Study on relative permeability characteristics affected by displacement pressure gradient: Experimental study and numerical simulation | |
| Liu et al. | Interactions between coal seam gas drainage boreholes and the impact of such on borehole patterns | |
| CN108133080B (en) | Heterogeneous fractured shale gas reservoir numerical simulation method considering unsteady adsorption | |
| CN105089649B (en) | A kind of inter well connectivity method for establishing model | |
| CN104695950A (en) | Prediction method for volcanic rock oil reservoir productivity | |
| CN105089566A (en) | Gas well system production allocation method | |
| Yang et al. | A semianalytical method for modeling two-phase flow in coalbed-methane reservoirs with complex fracture networks | |
| CN116306385B (en) | Oil reservoir fracturing imbibition energy increasing numerical simulation method, system, equipment and medium | |
| CN106202748B (en) | A kind of predictive methodology of coalbed methane recovery efficiency based on permeability and well spacing | |
| Li et al. | A new method for production data analysis in shale gas reservoirs | |
| Li et al. | Effect of adsorption and permeability correction on transient pressures in organic rich gas reservoirs: Vertical and hydraulically fractured horizontal wells | |
| Shi et al. | Pressure transient behavior of layered commingled reservoir with vertical inhomogeneous closed boundary | |
| Zhong et al. | Study on Co-production compatibility evaluation method of multilayer tight gas reservoir |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |