CN107463761A - gas storage gas injection well control dynamic evaluation method - Google Patents
gas storage gas injection well control dynamic evaluation method Download PDFInfo
- Publication number
- CN107463761A CN107463761A CN201610390275.XA CN201610390275A CN107463761A CN 107463761 A CN107463761 A CN 107463761A CN 201610390275 A CN201610390275 A CN 201610390275A CN 107463761 A CN107463761 A CN 107463761A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- represent
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention provides a kind of gas storage gas injection well control dynamic evaluation method, this method includes:Derive step:Derive zero dimension gas injection rate decreasing model;Modeling procedure:Gas injection rate theoretical model is established according to zero dimension gas injection rate decreasing model;Draw letter and build figure step:Regularization pseudopressure function and material balance pseudotime function are introduced, zero dimension gas injection rate and material balance pseudotime theory plate are established according to gas injection rate theoretical model;Curve fitting step:Obtain zero dimension gas injection rate and material balance pseudotime relation curve and match fitting nondimensional number gas injection relation curve and zero dimension gas injection rate and material balance pseudotime theory plate;Calculation procedure:Calculate the storage of gas well gas injection process and ooze parameter.The present invention can carry out dynamic evaluation to gas storage gas injection well well control.
Description
Technical field
The present invention relates to gas storage gas injection well control assessment technique field, in particular to a kind of gas storage gas injection well control
Dynamic evaluation method.
Background technology
The gas well modern times unstable analysis method of yield is widely used in gas reservoir development dynamic analysis, is also succeeded in recent years
Introduce and be applied in gas storage gas well gas production dynamic tracking evaluation, achieve preferable application effect, but it is only applicable to adopt
Gas process, the gas injection operation of gas storage are the processes of a high speed injection, and gas injection dynamic tracking evaluation technical research is less, lacks
Necessary theory and technology.
The unstable analysis method of gas storage gas well gas injection can provide and normal pressures transient well test analysis identical ginseng
Number information and required precision, a kind of new method is provided for the dynamic tracking evaluation during gas well injection gas, with cost it is low,
Data is wide, explains the advantages of simple, precision is high.
At present, the dynamic evaluation expansion in the determination of gas storage gas well rational productivity, especially gas production process is corresponding
Research, such as the patent " underground natural gas storage tank rational productivity forecast value revision method " that number of patent application is 201510157312.8 are open
One kind is related to underground natural gas storage tank gas well rational productivity forecast value revision method, by establishing gas well reasonable production pressure, so that it is right
Gas well is flowed into and is modified with elution curve joint, gas well ability of reasonable production after being corrected.Number of patent application is
201410809184.6 patent " Forecasting Methodology and device of gas storage air water interactive areas well capacity " discloses a kind of gas storage
The Forecasting Methodology and device of air water interactive areas well capacity, this method Seepage Experiment modified result gas well binomial potential curve and equation,
Obtain deducting in gas storage running and lose the research object well capacity of percolation ability because air water interacts displacement.
Nowadays, more ripe for gas storage productivity evaluation of gas well method, the dynamic monitoring in gas production process also has certain method, but all
Fail to carry out gas injection process the well control dynamic evaluation during system evaluation, especially high speed injection.
The content of the invention
It is a primary object of the present invention to provide a kind of gas storage gas injection well control dynamic evaluation method, to solve prior art
In the problem of can not evaluating gas injection well control degree.
To achieve these goals, the invention provides a kind of gas storage gas injection well control dynamic evaluation method, this method bag
Include:Derive step:Derive zero dimension gas injection rate decreasing model;Modeling procedure:Established and noted according to zero dimension gas injection rate decreasing model
Tolerance theoretical model;Draw letter and build figure step:Regularization pseudopressure function and material balance pseudotime function are introduced, according to gas injection rate
Theoretical model establishes zero dimension gas injection rate and material balance pseudotime theory plate;Curve fitting step:Obtain zero dimension gas injection
Amount and material balance pseudotime relation curve simultaneously match fitting nondimensional number gas injection relation curve and zero dimension gas injection rate and material
Balance pseudotime theory plate;Calculation procedure:Calculate the storage of gas well gas injection process and ooze parameter.
Further, the derivation formula of zero dimension gas injection rate decreasing model is:
Wherein, r represents distance of the gas generation border away from wellbore centre, and P represents pseudopressure, and μ represents gas viscosity, K generations
Table in-place permeability, Ct represent stratigraphic compression coefficient, rwWellbore radius is represented, φ represents formation porosity, and it is effective that h represents gas-bearing formation
Thickness, PiRepresent note just strata pressure, PwfRepresent flowing bottomhole pressure (FBHP), reWell control radius is represented, B represents gas volume factor, and q is represented
Gas injection rate.
Further, the gas injection rate theoretical model formula in modeling procedure is:
Wherein,Dimensionless production is represented, α represents the first dimensionless production coefficient, and θ represents the second dimensionless production system
Number, K1Represent the first Bessel function, K0Represent the second Bessel function, I1Represent the 3rd Bessel function, I0Represent the 4th shellfish
Sai Er functions, behalf skin factor, reDRepresent zero dimension well control radius.
Further, material balance pseudotime function is:
Wherein, q represents gas injection rate, CtRepresent gas compressibility factor, CtiIt is gas viscosity for original gas compressed coefficient μ, t
The time is represented, G represents well head reserves,Represent original formation pressure Regularization pseudopressure, ppRepresent strata pressure Regularization plan
Pressure.
Further, Regularization pseudopressure function is:
Wherein, ppStrata pressure Regularization pseudopressure is represented, μ is gas viscosity, and Z represents deviation factor for gas, and p, which is represented, to be intended
Pressure.
Further, curve fitting step includes stream backfin calculation step, and conversion stream pressure number is obtained by flowing backfin calculation step
According to, conversion stream pressure data are contrasted with the actual pressure data that flow, wherein, the calculation formula that step is calculated in stream backfin is:
Log δ=0.31-0.49 × Tr+0.18×Tr 2
S=0.03415 × γgH/(TavZav),
Wherein, α represents the first deviation factors design factor, and β represents the second deviation factors design factor, and it is inclined that γ represents the 3rd
Poor factor design factor, δ represent the 4th deviation factors design factor, and p represents pseudopressure, TrRepresent pseudoreduced temperature, PrRepresent and intend
Reduced pressure, Z represent deviation factor for gas, PwfRepresent bottom pressure, PwhWell head pressure, behalf skin factor are represented, λ is represented
Oil pipe resistance coefficient, ε represent conversion factor, and qi represents gas injection rate, TavRepresent post mean temperature of being taken offence in pit shaft, ZavRepresent well
Taken offence post coefficient of mean deviation in cylinder, d represents oil pipe interior diameter, γgRepresent natural gas relative density.
Further, curve fitting step also includes obtaining Regularization gas injection integration, Regularization gas injection integral derivative and thing
Matter equilibration time relation curve step, passes through formula
Obtain Regularization gas injection integration, Regularization gas injection integral derivative and material balance time curve;
Wherein, QiRepresent Regularization gas injection integration, tcThe gas well material balance time is represented, q represents gas injection rate, Δ ppRepresent
Regularization pseudopressure is poor,Bottom pressure Regularization pseudopressure is represented,Original formation pressure Regularization pseudopressure is represented,
QidRepresent Regularization gas injection integral derivative.
Further, parameter is oozed in storage includes gas injection well control scope, gas injection well control reserves, effective permeability, skin factor, pressure
Power fitting precision and yield fitting precision.
Further, gas injection well control scope formula is:
Wherein, reRepresent well control radius, BiRepresent original gas volume factor, SwRepresent water saturation, CtiRepresent original
Gas compressibility factor, tcRepresent gas well material balance time, tcDRepresent gas well zero dimension material balance time, QiRepresent Regularization
Gas injection integrates, and qDi represents zero dimension gas injection rate integration, and M represents match point footmark at M, and it is effective to represent gas-bearing formation by h at the M of the application
Thickness.
Further, gas injection well control reserves formula is:
Wherein, G represents gas injection well control reserves, CtiRepresent the original gas compressed coefficient, tcThe gas well material balance time is represented,
tcDRepresent gas well zero dimension material balance time, QiRepresent Regularization gas injection integration, qDiRepresent zero dimension gas injection rate integration.
Further, effective permeability formula is:
Wherein,Original formation pressure Regularization pseudopressure is represented,Represent bottom pressure Regularization pseudopressure, q generations
Table gas injection rate, CtiRepresent the original gas compressed coefficient, tcaThe average gas well material balance time is represented, G represents the storage of gas injection well control
Amount, μ represent gas viscosity, and B represents gas volume factor, and i represents initial data, and K represents in-place permeability, and h, which represents gas-bearing formation, to be had
Imitate thickness, CAForm factor is represented, γ represents gas relative density, rwWellbore radius is represented, A is represented, reDRepresent zero dimension well
Radius is controlled, M represents match point footmark at M.
Further, the calculation formula of skin factor is:
Wherein, rwaEffective hole diameter is represented, K represents in-place permeability, and φ represents formation porosity, and μ represents gas viscosity, Sw
Represent water saturation, CtiRepresent the original gas compressed coefficient, tcRepresent gas well material balance time, tcDRepresent gas well zero dimension
Material balance time, M represent match point footmark at M, reDRepresent zero dimension well control radius, S is pull-type variable, rwRepresent pit shaft half
Footpath.
Further, the fitting of pressure fitting precision and yield fitting precision is with reference to formula:
Wherein, qiGas injection rate is represented, K represents in-place permeability, and h represents gas-bearing net pay, TavRepresent in pit shaft and take offence
Post mean temperature, PwfRepresent bottom pressure, PrRepresent pseudoreduced pressure, μgRepresent, ZavRepresent post average deviation of being taken offence in pit shaft
Coefficient, PavRepresent post average pressure of being taken offence in pit shaft, TfRepresent reservoir temperature, S is pull-type variable, rwRepresent wellbore radius, reGeneration
Table well control radius.
Apply the technical scheme of the present invention, the present invention is noted by deriving zero dimension gas injection rate decreasing model according to zero dimension
Tolerance decreasing model establishes gas injection rate theoretical model, introduces Regularization pseudopressure function and material balance pseudotime function afterwards,
Zero dimension gas injection rate and material balance pseudotime theory plate are established according to gas injection rate theoretical model, and obtain zero dimension gas injection rate
With material balance pseudotime relation curve and match fitting nondimensional number gas injection relation curve put down with zero dimension gas injection rate with material
Weigh pseudotime theory plate, and finally oozing parameter to each storage calculates, and realizes to the well control dynamic in gas storage gas injection process
Evaluation.
Brief description of the drawings
The Figure of description for forming the part of the application is used for providing a further understanding of the present invention, and of the invention shows
Meaning property embodiment and its illustrate be used for explain the present invention, do not form inappropriate limitation of the present invention.In the accompanying drawings:
Fig. 1 shows the schematic diagram of the zero dimension gas injection plate of the gas storage gas injection well control dynamic evaluation method of the present invention;
Fig. 2 shows the contrast of the actual measurement stream pressure and fitted flow pressure of the gas storage gas injection well control dynamic evaluation method of the present invention
Figure;
Fig. 3 shows the actual gas injection fitted figure of zero dimension of the gas storage gas injection well control dynamic evaluation method of the present invention;
Fig. 4 shows the multicycle gas injection performance matching of the gas storage gas injection well control dynamic evaluation method of the present invention;And
Fig. 5 shows the flow chart of the gas storage gas injection well control dynamic evaluation method of the present invention.
Embodiment
It should be noted that in the case where not conflicting, the feature in embodiment and embodiment in the application can phase
Mutually combination.Describe the present invention in detail below with reference to the accompanying drawings and in conjunction with the embodiments.
It is noted that described further below is all exemplary, it is intended to provides further instruction to the application.It is unless another
Indicate, all technologies used herein and scientific terminology are with usual with the application person of an ordinary skill in the technical field
The identical meanings of understanding.
It should be noted that term used herein above is merely to describe embodiment, and be not intended to restricted root
According to the illustrative embodiments of the application.As used herein, unless the context clearly indicates otherwise, otherwise singulative
It is also intended to include plural form, additionally, it should be understood that, when in this manual using term "comprising" and/or " bag
Include " when, it indicates existing characteristics, step, operation, device, component and/or combinations thereof.
It should be noted that term " first " in the description and claims of this application and above-mentioned accompanying drawing, "
Two " etc. be for distinguishing similar object, without for describing specific order or precedence.It should be appreciated that so use
Data can exchange in the appropriate case, so that presently filed embodiment described herein for example can be with except herein
Order beyond those of diagram or description is implemented.
In addition, term " comprising " and " having " and their any deformation, it is intended that cover it is non-exclusive include, example
Such as, process, method, system, product or the equipment for containing series of steps or unit are not necessarily limited to those clearly listed
Step or unit, but may include not list clearly or for intrinsic other of these processes, method, product or equipment
Step or unit.
For the ease of description, space relative terms can be used herein, as " ... on ", " ... top ",
" ... upper surface ", " above " etc., for describing such as a device shown in the figure or feature and other devices or spy
The spatial relation of sign.It should be appreciated that space relative terms are intended to comprising the orientation except device described in figure
Outside different azimuth in use or operation.
For example, if the device in accompanying drawing is squeezed, it is described as " above other devices or construction " or " in other devices
On part or construction " device after will be positioned as or " in other devices or constructing it " below other devices or construction "
Under ".Thus, exemplary term " in ... top " can include " in ... top " and " in ... lower section " two kinds of orientation.Should
Device can also other different modes positioning (being rotated by 90 ° or in other orientation), it is and relative to space used herein above
Respective explanations are made in description.
Referring to shown in Fig. 1 to Fig. 5, the invention provides a kind of gas storage gas injection well control dynamic evaluation method, the gas storage
Gas injection well control dynamic evaluation method builds figure step, curve fitting step and calculating step including deriving step, modeling procedure, drawing letter
Suddenly.
Wherein, derive step to refer to derive zero dimension gas injection rate decreasing model, modeling procedure refers to according to zero dimension gas injection
Amount decreasing model establishes gas injection rate theoretical model, draws letter and builds figure step and refers to that introducing Regularization pseudopressure function and material balance intends
The function of time, zero dimension gas injection rate and material balance pseudotime theory plate, curve are established according to the gas injection rate theoretical model
Fit procedure refers to obtain zero dimension gas injection rate and material balance pseudotime relation curve and matches the fitting nondimensional number note
Gas relation curve and the zero dimension gas injection rate and material balance pseudotime theory plate, calculation procedure refer to calculate gas well gas injection
Parameter is oozed in process storage.The present invention can carry out dynamic evaluation to gas storage gas injection well well control.
As shown in figure 1, the present invention applies to the gas injection process dynamic evaluation method of gas reservoir type gas storage, this method is to build
The overall merit means stood on the basis of zero dimension gas injection plate.Abscissa is the material balance time in Fig. 1, ordinate be without because
Secondary gas injection integrates qDiWith zero dimension gas injection integral derivative qDid, its role is under the conditions of reacting different zero dimension well control radiuses
Set of curves, as direction of arrow dimensionless radius gradually increases.
Assuming that it is r in well control radiuse, outer sealed boundary circular boundary formation in, a bite gas well is carried out with constant gas injection rate q
Gas injection, the derivation formula of zero dimension gas injection rate decreasing model are:
Wherein, r represents distance of the gas generation border away from wellbore centre, unit m, and P represents pseudopressure, unit MPa,
μ represents gas viscosity, zero dimension, and K represents in-place permeability, and unit Md, Ct represent stratigraphic compression coefficient, zero dimension, rwGeneration
Table wellbore radius, unit m,Formation porosity, % are represented, h represents gas-bearing net pay, unit m, PiRepresent note just
Stressor layer, unit Mpa, PwfRepresent flowing bottomhole pressure (FBHP), unit Mpa, reWell control radius is represented, unit m, B represent gas body
Product coefficient, zero dimension, q represent gas injection rate, unit 104m3/d。
By above-mentioned Definite problem zero dimension, and Superposition Principle is utilized, obtain zero dimension gas injection rate under the conditions of level pressure
Laplace space expression formula, i.e., the present invention in modeling procedure in the gas injection rate theoretical model formula it is as follows:
Wherein,Dimensionless production, zero dimension are represented, α represents the first dimensionless production coefficient, zero dimension, and θ represents second
Dimensionless production coefficient, zero dimension, K1Represent the first Bessel function, K0Represent the second Bessel function, I1Represent the 3rd shellfish plug
That function, I0Represent the 4th Bessel function, behalf skin factor, zero dimension, reDRepresent zero dimension well control radius, zero dimension.
As shown in Fig. 2 the oil pressure in gas injection process needs to convert to shaft bottom, then calculate Regularization gas injection integration and Regularization
Gas injection integral derivative, realize such as the fitting effect in Fig. 3.
During flowing bottomhole pressure (FBHP) considers gas testing, stream extrudes the existing phenomenon easily fluctuated, the curve matching of the invention
Step include stream backfin calculate step, by it is described stream backfin calculate step obtain conversion stream pressure data, it is described conversion stream pressure data with
Actual stream pressure data are contrasted, wherein, the calculation formula that step is calculated in the stream backfin is:
Log δ=0.31-0.49 × Tr+0.18×Tr 2
S=0.03415 × γgH/(TavZav) (3)
Carry out flowing backfin calculation using formula (2) and formula (3), conversion is flowed and presses data and actual measurement to flow pressure and contrasted, is intended
92.5% can be reached by closing precision.
Wherein, α represents the first deviation factors design factor, zero dimension, and β represents the second deviation factors design factor, it is no because
Secondary, γ represents the 3rd deviation factors design factor, zero dimension, and δ represents the 4th deviation factors design factor, zero dimension, and p, which is represented, to be intended
Pressure, unit Mpa, TrRepresent pseudoreduced temperature, zero dimension, PrRepresent pseudoreduced pressure, zero dimension, Z represent gas deviation because
Son, zero dimension, PwfRepresent bottom pressure, unit Mpa, PwhRepresent well head pressure, unit Mpa, behalf skin factor, nothing
Dimension, λ represent oil pipe resistance coefficient, zero dimension, and ε represents conversion factor, zero dimension, qiRepresent gas injection rate, unit 104m3/ d,
TavRepresent post mean temperature of being taken offence in pit shaft, unit K, ZavRepresent post coefficient of mean deviation of being taken offence in pit shaft, zero dimension, d generations
Table oil pipe interior diameter, unit m, γgRepresent natural gas relative density, zero dimension.
As shown in figure 3, to each actual gas injection data point, substance for calculation equilibration time, Regularization gas injection integration and rule
Integralization gas injection integral derivative, for theoretical chart fitting, with determine zero dimension well control radius, effective permeability, effective hole diameter,
Skin factor, well control radius and well control reserves.
Specifically, the material balance pseudotime function is:
Wherein, q represents gas injection rate, unit, CtRepresent gas compressibility factor, unit MPa-1, CtiCompressed for original gas
Coefficient, unit Mpa-1, μ is gas viscosity, and unit Mpa/s, t represent the time, and unit d, G represent well control reserves, unit
For 108m3/ d,Represent original formation pressure Regularization pseudopressure, unit Mpa, ppRepresent strata pressure Regularization and intend pressure
Power, unit Mpa.
The Regularization pseudopressure function is:
(5)
Wherein, ppRepresent strata pressure Regularization pseudopressure, unit Mpa, μ are gas viscosity, unit Mpa/s, Z generation
Table deviation factor for gas, p represent pseudopressure, unit Mpa.
The curve fitting step of the present invention also include obtaining Regularization gas injection integration, Regularization gas injection integral derivative with
Material balance time curve step, passes through formula
Regularization gas injection integration, Regularization gas injection integral derivative and material balance time curve are obtained, wherein, QiGeneration
The gas injection of table Regularization integrates, tcThe gas well material balance time is represented, unit d, q represent gas injection rate, unit 104m3/ d, Δ pp
It is poor to represent Regularization pseudopressure, unit Mpa,Represent bottom pressure Regularization pseudopressure, unit Mpa,Represent original
Beginning strata pressure Regularization pseudopressure, unit Mpa, QidRepresent Regularization gas injection integral derivative.
The material balance pseudotime is iterated to calculate to error condition using setting well control reserves, takes a well control reserves G, according to
Formula (4), formula (5) calculate actual contents balance pseudotime, Regularization pseudopressure, drawStream material
The profile of equilibrium, G is determined according to regression straight line slope variance, iteration tentative calculation is until slope variance meets condition, then passes through formula (6)
It is double right with the material balance time that the Regularization gas injection integration having determined that under G, Regularization gas injection integral derivative are drawn with formula (7)
Number curve.
Parameter is oozed in the storage of the present invention includes gas injection well control scope, gas injection well control reserves, effective permeability, epidermis system
Number, pressure fitting precision and yield fitting precision.
Zero dimension well control radius is to be fitted using above-mentioned Regularization gas injection curve with theoretical plate, and fitting first determines
Zero dimension well control radius reD;
Wherein, gas reservoir effective permeability is that selection match point progress inverse obtains.Arbitrarily 1 match point of selection, record are real
Border match point (tc, Qi) and M and corresponding theory chart fitting point (tcD,qDi)M;According to yield match point, using monophasic fluid not
Stationary flow and pseudostable flow merge solution formula (8), if known gas viscosity, volume factor, gas reservoir thickness and zero dimension well control half
Ask for effective permeability and see formula (9) in footpath;
Wherein,Represent original formation pressure Regularization pseudopressure, unit Mpa,Represent bottom pressure Regularization
Pseudopressure, unit Mpa, q represent gas injection rate, unit 104m3/ d, CtiRepresent the original gas compressed coefficient, unit MPa-1,
tcaThe average gas well material balance time is represented, unit d, G represent gas injection well control reserves, unit 108m3/ d, μ represent gas
Viscosity, unit Mpa/s, B represent gas volume factor, and i represents initial data, and K represents in-place permeability, unit md, h generation
Table gas-bearing net pay, unit m, CAForm factor is represented, γ represents gas relative density, zero dimension, rwRepresent pit shaft half
Footpath, unit m, A represent well control area, reDZero dimension well control radius is represented, M represents match point footmark at M, wherein, position at M
Obtained to be random.
Effective hole diameter rwaIt is to react the wellbore radius that can be described under the actual pollution condition in shaft bottom.When being intended according to material balance
Between relation ask for effective hole diameter rwaSuch as formula (10);
Skin factor can reflect the pollution level in shaft bottom.In known wellbore radius rwUnder conditions of, calculate epidermis coefficient S
Such as formula (11);
Well control radius reThe control range that gas can feed through in the actual gas injection process of reaction, well control reserves G be
Well control radius reUnder circular bounded homogeneous formation straight well control reserve.Determine well control radius reScope and the storage of gas injection rate well control
Measure G such as formula (12) and formula (13);
Wherein, reRepresent well control radius, unit m, BiRepresent original gas volume factor, SwRepresentative contains water saturation
Degree, %, CtiRepresent the original gas compressed coefficient, unit MPa-1, tcRepresent gas well material balance time, unit d, tcDGeneration
Table gas well zero dimension material balance time, unit d, QiRepresent Regularization gas injection integration, qDiZero dimension gas injection rate integration is represented,
M represents match point footmark at M, and h represents gas-bearing net pay, unit m.
Wherein, G represents gas injection well control reserves, CtiRepresent the original gas compressed coefficient, unit MPa-1, tcRepresent gas well thing
Matter equilibration time, unit d, tcDRepresent gas well zero dimension material balance time, QiRepresent Regularization gas injection integration, qDiRepresent nothing
Dimension gas injection rate integrates.
As shown in figure 4, using the well control degree evaluation to each gas injection cycle of individual well, individual well each gas injection week is obtained
Phase gas reservoir effective permeability, skin factor and well control radius, and carry out being segmented yield and stream pressure fitting accordingly.
The fitting of the pressure fitting precision and yield fitting precision of the present invention is with reference to formula:
Wherein, qiRepresenting gas injection rate, K represents in-place permeability, unit md, and h represents gas-bearing net pay, unit m,
Tav represents post mean temperature of being taken offence in pit shaft, unit K, PwfBottom pressure is represented, unit Mpa, Pr, which is represented, intends reduced pressure
Power, unit Mpa, μ g represent gas viscosity, unit Mpa/s, ZavRepresent post coefficient of mean deviation of being taken offence in pit shaft, it is no because
It is secondary, PavRepresent post average pressure of being taken offence in pit shaft, unit Mpa, TfRepresent reservoir temperature, S is pull-type variable, rwRepresent pit shaft
Radius, unit m, reRepresent well control radius, unit m.
Piecewise fitting refers to that bringing Fig. 3 interpretive analysis results into note adopts in dynamic model, as shown in formula (14), it can be seen that
Be advantageous to improve individual well history matching precision, by taking individual well shown in Fig. 3 as an example, the well using this method evaluation gas injection well control degree
When 3 injection-production cycles are predicted using level pressure, yield fitting precision is up to 90.3%, is laid the foundation for later stage injection allocation prediction.
The preferred embodiments of the present invention are the foregoing is only, are not intended to limit the invention, for the skill of this area
For art personnel, the present invention can have various modifications and variations.Within the spirit and principles of the invention, that is made any repaiies
Change, equivalent substitution, improvement etc., should be included in the scope of the protection.
Claims (13)
- A kind of 1. gas storage gas injection well control dynamic evaluation method, it is characterised in that including:Derive step:Derive zero dimension gas injection rate decreasing model;Modeling procedure:Gas injection rate theoretical model is established according to the zero dimension gas injection rate decreasing model;Draw letter and build figure step:Regularization pseudopressure function and material balance pseudotime function are introduced, it is theoretical according to the gas injection rate Model establishes zero dimension gas injection rate and material balance pseudotime theory plate;Curve fitting step:Obtain zero dimension gas injection rate and material balance pseudotime relation curve and match and be fitted the zero dimension Gas injection rate relation curve and the zero dimension gas injection rate and material balance pseudotime theory plate;Calculation procedure:Calculate the storage of gas well gas injection process and ooze parameter.
- 2. gas storage gas injection well control dynamic evaluation method according to claim 1, it is characterised in that the zero dimension gas injection Amount decreasing model derivation formula be:<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <mi>r</mi> </mfrac> <mfrac> <mo>&part;</mo> <mrow> <mo>&part;</mo> <mi>r</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mi>r</mi> <mfrac> <mrow> <mo>&part;</mo> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <mi>r</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&phi;&mu;C</mi> <mi>t</mi> </msub> </mrow> <mi>K</mi> </mfrac> <mfrac> <mrow> <mo>&part;</mo> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <mi>t</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mrow> <mo>(</mo> <mi>r</mi> <mfrac> <mrow> <mo>&part;</mo> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <mi>r</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mi>r</mi> <mo>=</mo> <msub> <mi>r</mi> <mi>w</mi> </msub> </mrow> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mi>q</mi> <mi>&mu;</mi> <mi>B</mi> </mrow> <mrow> <mn>2</mn> <mi>&pi;</mi> <mi>K</mi> <mi>h</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mrow> <mo>&part;</mo> <mi>p</mi> </mrow> <mrow> <mo>&part;</mo> <mi>r</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>r</mi> <mo>=</mo> <msub> <mi>r</mi> <mi>e</mi> </msub> </mrow> </msub> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>p</mi> <mrow> <mo>(</mo> <mrow> <mi>r</mi> <mo>,</mo> <mn>0</mn> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>p</mi> <mrow> <mi>w</mi> <mi>f</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced>Wherein, r represents distance of the gas generation border away from wellbore centre, and P represents pseudopressure, and μ represents gas viscosity, and K represents ground Layer permeability, CtRepresent stratigraphic compression coefficient, rwWellbore radius is represented, φ represents formation porosity, and h represents gas-bearing net pay, PiRepresent note just strata pressure, PwfRepresent flowing bottomhole pressure (FBHP), reWell control radius is represented, B represents gas volume factor, and q represents gas injection Amount.
- 3. gas storage gas injection well control dynamic evaluation method according to claim 1, it is characterised in that in the modeling procedure The gas injection rate theoretical model formula be:<mrow> <msub> <mover> <mi>q</mi> <mo>&OverBar;</mo> </mover> <mi>D</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&alpha;</mi> <mi>&theta;</mi> </mrow> <msqrt> <mi>s</mi> </msqrt> </mfrac> <mrow> <mo>&lsqb;</mo> <mfrac> <mrow> <msub> <mi>K</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>&theta;</mi> <msqrt> <mi>s</mi> </msqrt> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>I</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>&theta;</mi> <msqrt> <mi>s</mi> </msqrt> <mo>)</mo> </mrow> <mfrac> <mrow> <msub> <mi>K</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>e</mi> <mi>D</mi> </mrow> </msub> <mi>&theta;</mi> <msqrt> <mi>s</mi> </msqrt> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>I</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>e</mi> <mi>D</mi> </mrow> </msub> <mi>&theta;</mi> <msqrt> <mi>s</mi> </msqrt> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow> <mrow> <msub> <mi>K</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>&theta;</mi> <msqrt> <mi>s</mi> </msqrt> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>I</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>&theta;</mi> <msqrt> <mi>s</mi> </msqrt> <mo>)</mo> </mrow> <mfrac> <mrow> <msub> <mi>K</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>e</mi> <mi>D</mi> </mrow> </msub> <mi>&theta;</mi> <msqrt> <mi>s</mi> </msqrt> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>I</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>e</mi> <mi>D</mi> </mrow> </msub> <mi>&theta;</mi> <msqrt> <mi>s</mi> </msqrt> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow> </mfrac> <mo>&rsqb;</mo> </mrow> <mo>,</mo> </mrow>Wherein,Dimensionless production is represented, α represents the first dimensionless production coefficient, and θ represents the second dimensionless production coefficient, K1Generation The Bessel function of table first, K0Represent the second Bessel function, I1Represent the 3rd Bessel function, I0Represent the 4th Bezier letter Number, behalf skin factor, reDRepresent zero dimension well control radius.
- 4. gas storage gas injection well control dynamic evaluation method according to claim 1, it is characterised in that the material balance is intended The function of time is:<mrow> <msub> <mi>t</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <msub> <mrow> <mo>(</mo> <msub> <mi>&mu;C</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> <mi>i</mi> </msub> <mi>q</mi> </mfrac> <msubsup> <mo>&Integral;</mo> <mn>0</mn> <mi>t</mi> </msubsup> <mfrac> <mi>q</mi> <mrow> <msub> <mi>&mu;C</mi> <mi>t</mi> </msub> </mrow> </mfrac> <mi>d</mi> <mi>t</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>GC</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> <mi>q</mi> </mfrac> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mi>p</mi> </msub> <mo>-</mo> <msub> <mi>p</mi> <mrow> <mi>p</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> </mrow>Wherein, q represents gas injection rate, CtRepresent gas compressibility factor, CtiIt is gas viscosity for original gas compressed coefficient μ, t is represented Time, G represent well head reserves,Represent original formation pressure Regularization pseudopressure, ppRepresent strata pressure Regularization pseudopressure.
- 5. gas storage gas injection well control dynamic evaluation method according to claim 1, it is characterised in that the Regularization intends pressure Force function is:<mrow> <msub> <mi>p</mi> <mi>p</mi> </msub> <mo>=</mo> <msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>&mu;</mi> <mi>Z</mi> </mrow> <mi>p</mi> </mfrac> <mo>)</mo> </mrow> <mi>i</mi> </msub> <msubsup> <mo>&Integral;</mo> <mn>0</mn> <mi>p</mi> </msubsup> <mfrac> <mi>p</mi> <mrow> <mi>&mu;</mi> <mi>Z</mi> </mrow> </mfrac> <mi>d</mi> <mi>p</mi> <mo>,</mo> </mrow> 1Wherein, ppStrata pressure Regularization pseudopressure is represented, μ is gas viscosity, and Z represents deviation factor for gas, and p represents pseudopressure.
- 6. gas storage gas injection well control dynamic evaluation method according to claim 1, it is characterised in that the curve matching step It is rapid to calculate step including stream backfin, step is calculated by the stream backfin and obtains converting stream pressure data, the conversion stream presses data and reality Border stream pressure data are contrasted, wherein, the calculation formula that step is calculated in the stream backfin is:<mrow> <mi>&alpha;</mi> <mo>=</mo> <mfrac> <mrow> <mn>0.32</mn> <mo>&times;</mo> <msubsup> <mi>P</mi> <mi>r</mi> <mn>6</mn> </msubsup> </mrow> <msup> <mn>10</mn> <mrow> <mn>9</mn> <mo>&times;</mo> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>r</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msup> </mfrac> <mo>+</mo> <mrow> <mo>(</mo> <mfrac> <mn>0.07</mn> <mrow> <msub> <mi>T</mi> <mi>r</mi> </msub> <mo>-</mo> <mn>0.86</mn> </mrow> </mfrac> <mo>-</mo> <mn>0.04</mn> <mo>)</mo> </mrow> <mo>&times;</mo> <msubsup> <mi>P</mi> <mi>r</mi> <mn>2</mn> </msubsup> <mo>+</mo> <mrow> <mo>(</mo> <mn>0.62</mn> <mo>-</mo> <mn>0.23</mn> <mo>&times;</mo> <msub> <mi>T</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> <mo>&times;</mo> <msub> <mi>P</mi> <mi>r</mi> </msub> </mrow><mrow> <mi>&beta;</mi> <mo>=</mo> <mn>0.89</mn> <mo>-</mo> <mn>1.39</mn> <mo>&times;</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>r</mi> </msub> <mo>-</mo> <mn>0.92</mn> <mo>)</mo> </mrow> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </msup> <mo>+</mo> <mn>0.36</mn> <mo>&times;</mo> <msub> <mi>T</mi> <mi>r</mi> </msub> </mrow><mrow> <mi>&gamma;</mi> <mo>=</mo> <mn>0.132</mn> <mo>-</mo> <mn>0.32</mn> <mo>&times;</mo> <msub> <mi>log</mi> <mn>10</mn> </msub> <msub> <mi>T</mi> <mi>r</mi> </msub> <mo>&times;</mo> <msubsup> <mi>P</mi> <mi>r</mi> <mi>&delta;</mi> </msubsup> </mrow>Log δ=0.31-0.49 × Tr+0.18×Tr 2<mrow> <mi>Z</mi> <mo>=</mo> <mn>1</mn> <mo>+</mo> <mfrac> <mrow> <mi>&beta;</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>e</mi> <mi>&alpha;</mi> </msup> <mo>)</mo> </mrow> </mrow> <msup> <mi>e</mi> <mi>&alpha;</mi> </msup> </mfrac> <mo>+</mo> <mi>&gamma;</mi> <mo>,</mo> </mrow><mrow> <msub> <mi>P</mi> <mrow> <mi>w</mi> <mi>f</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <msub> <mi>P</mi> <mrow> <mi>w</mi> <mi>h</mi> </mrow> </msub> <mn>2</mn> </msup> <msup> <mi>e</mi> <mrow> <mn>2</mn> <mi>s</mi> </mrow> </msup> <mo>-</mo> <mn>1.3243</mn> <mo>&times;</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mn>10</mn> </mrow> </msup> <msup> <msub> <mi>&lambda;&epsiv;q</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <msup> <msub> <mi>T</mi> <mrow> <mi>a</mi> <mi>v</mi> </mrow> </msub> <mn>2</mn> </msup> <msup> <msub> <mi>Z</mi> <mrow> <mi>a</mi> <mi>v</mi> </mrow> </msub> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msup> <mi>e</mi> <mrow> <mn>2</mn> <mi>s</mi> </mrow> </msup> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>/</mo> <msup> <mi>d</mi> <mn>5</mn> </msup> </mrow> </msqrt> </mrow>S=0.03415 × γgH/(TavZav),Wherein, α represents the first deviation factors design factor, and β represents the second deviation factors design factor, γ represent the 3rd deviation because Sub- design factor, δ represent the 4th deviation factors design factor, and p represents pseudopressure, TrPseudoreduced temperature is represented, Pr represents plan pair Than pressure, Z represents deviation factor for gas, PwfRepresent bottom pressure, PwhWell head pressure, behalf skin factor are represented, λ represents oil Pipe resistance coefficient, ε represent conversion factor, and qi represents gas injection rate, TavRepresent post mean temperature of being taken offence in pit shaft, ZavRepresent pit shaft Inside take offence post coefficient of mean deviation, d represents oil pipe interior diameter, γgRepresent natural gas relative density.
- 7. gas storage gas injection well control dynamic evaluation method according to claim 6, it is characterised in that the curve matching step Suddenly also include obtaining Regularization gas injection integration, Regularization gas injection integral derivative and material balance time curve step, pass through Formula<mrow> <msub> <mi>Q</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>t</mi> <mi>c</mi> </msub> </mfrac> <msubsup> <mo>&Integral;</mo> <mn>0</mn> <msub> <mi>t</mi> <mi>c</mi> </msub> </msubsup> <mfrac> <mi>q</mi> <mrow> <msub> <mi>&Delta;p</mi> <mi>p</mi> </msub> </mrow> </mfrac> <mi>d</mi> <mi>&tau;</mi> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>t</mi> <mi>c</mi> </msub> </mfrac> <msubsup> <mo>&Integral;</mo> <mn>0</mn> <msub> <mi>t</mi> <mi>c</mi> </msub> </msubsup> <mfrac> <mi>q</mi> <mrow> <msub> <mi>p</mi> <msub> <mi>p</mi> <mrow> <mi>w</mi> <mi>f</mi> </mrow> </msub> </msub> <mo>-</mo> <msub> <mi>p</mi> <msub> <mi>p</mi> <mi>i</mi> </msub> </msub> </mrow> </mfrac> <mi>d</mi> <mi>&tau;</mi> <mo>,</mo> </mrow><mrow> <msub> <mi>Q</mi> <mrow> <mi>i</mi> <mi>d</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>dQ</mi> <mi>i</mi> </msub> </mrow> <mrow> <msub> <mi>dlnt</mi> <mi>c</mi> </msub> </mrow> </mfrac> <mo>=</mo> <msub> <mi>t</mi> <mi>c</mi> </msub> <mfrac> <mrow> <msub> <mi>dQ</mi> <mi>i</mi> </msub> </mrow> <mrow> <msub> <mi>dt</mi> <mi>c</mi> </msub> </mrow> </mfrac> <mo>,</mo> </mrow>Obtain Regularization gas injection integration, Regularization gas injection integral derivative and material balance time curve;Wherein, QiRepresent Regularization gas injection integration, tcThe gas well material balance time is represented, q represents gas injection rate, Δ ppRepresent regular It is poor to change pseudopressure,Bottom pressure Regularization pseudopressure is represented,Represent original formation pressure Regularization pseudopressure, QidRepresent Regularization gas injection integral derivative.
- 8. gas storage gas injection well control dynamic evaluation method according to claim 6, it is characterised in that parameter bag is oozed in the storage Include gas injection well control scope, gas injection well control reserves, effective permeability, skin factor, pressure fitting precision and yield fitting precision.
- 9. gas storage gas injection well control dynamic evaluation method according to claim 8, it is characterised in that the gas injection well control model Enclosing formula is:<mrow> <msub> <mi>r</mi> <mi>e</mi> </msub> <mo>=</mo> <msqrt> <mfrac> <mrow> <mfrac> <msub> <mi>B</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>S</mi> <mi>w</mi> </msub> <mo>)</mo> <msub> <mi>C</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mfrac> <msub> <mrow> <mo>(</mo> <mfrac> <msub> <mi>t</mi> <mi>c</mi> </msub> <msub> <mi>t</mi> <mrow> <mi>c</mi> <mi>D</mi> </mrow> </msub> </mfrac> <mo>)</mo> </mrow> <mi>M</mi> </msub> <msub> <mrow> <mo>(</mo> <mfrac> <msub> <mi>Q</mi> <mi>i</mi> </msub> <msub> <mi>q</mi> <mrow> <mi>D</mi> <mi>i</mi> </mrow> </msub> </mfrac> <mo>)</mo> </mrow> <mi>M</mi> </msub> </mrow> <mrow> <mi>&pi;</mi> <mi>h</mi> <mi>&phi;</mi> </mrow> </mfrac> </msqrt> <mo>,</mo> </mrow>Wherein, reRepresent well control radius, BiRepresent original gas volume factor, SwRepresent water saturation, CtiRepresent original gas The compressed coefficient, tcRepresent gas well material balance time, tcDRepresent gas well zero dimension material balance time, QiRepresent Regularization gas injection Integration, qDi represent zero dimension gas injection rate integration, and M represents match point footmark at M, and h represents gas-bearing net pay.
- 10. gas storage gas injection well control dynamic evaluation method according to claim 8, it is characterised in that the gas injection well control Reserves formula is:<mrow> <mi>G</mi> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>C</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mfrac> <msub> <mrow> <mo>(</mo> <mfrac> <msub> <mi>t</mi> <mi>c</mi> </msub> <msub> <mi>t</mi> <mrow> <mi>c</mi> <mi>D</mi> </mrow> </msub> </mfrac> <mo>)</mo> </mrow> <mi>M</mi> </msub> <msub> <mrow> <mo>(</mo> <mfrac> <msub> <mi>Q</mi> <mi>i</mi> </msub> <msub> <mi>q</mi> <mrow> <mi>D</mi> <mi>i</mi> </mrow> </msub> </mfrac> <mo>)</mo> </mrow> <mi>M</mi> </msub> <mo>,</mo> </mrow>Wherein, G represents gas injection well control reserves, CtiRepresent the original gas compressed coefficient, tcRepresent gas well material balance time, tcDGeneration Table gas well zero dimension material balance time, QiRepresent Regularization gas injection integration, qDiRepresent zero dimension gas injection rate integration.
- 11. gas storage gas injection well control dynamic evaluation method according to claim 8, it is characterised in that the effectively infiltration Rate formula is:<mrow> <mfrac> <mrow> <msub> <mi>p</mi> <msub> <mi>p</mi> <mrow> <mi>w</mi> <mi>f</mi> </mrow> </msub> </msub> <mo>-</mo> <msub> <mi>p</mi> <msub> <mi>p</mi> <mi>i</mi> </msub> </msub> </mrow> <mi>q</mi> </mfrac> <mo>=</mo> <mfrac> <msub> <mi>t</mi> <mrow> <mi>c</mi> <mi>a</mi> </mrow> </msub> <mrow> <msub> <mi>GC</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> </mrow> </mfrac> <mo>+</mo> <mfrac> <msub> <mrow> <mo>(</mo> <mi>&mu;</mi> <mi>B</mi> <mo>)</mo> </mrow> <mi>i</mi> </msub> <mrow> <mn>2</mn> <mi>&pi;</mi> <mi>K</mi> <mi>h</mi> </mrow> </mfrac> <mo>&lsqb;</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>ln</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>4</mn> <mi>A</mi> </mrow> <mrow> <msub> <mi>C</mi> <mi>A</mi> </msub> <msup> <mi>e</mi> <mi>&gamma;</mi> </msup> <msubsup> <mi>r</mi> <mi>w</mi> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>&rsqb;</mo> <mo>,</mo> </mrow><mrow> <mi>K</mi> <mo>=</mo> <msub> <mrow> <mo>(</mo> <mfrac> <msub> <mi>Q</mi> <mi>i</mi> </msub> <msub> <mi>q</mi> <mrow> <mi>D</mi> <mi>i</mi> </mrow> </msub> </mfrac> <mo>)</mo> </mrow> <mi>M</mi> </msub> <mfrac> <mrow> <mi>&mu;</mi> <mi>B</mi> </mrow> <mrow> <mn>2</mn> <mi>&pi;</mi> <mi>h</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>lnr</mi> <mrow> <mi>e</mi> <mi>D</mi> </mrow> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> <mo>,</mo> </mrow>Wherein,Original formation pressure Regularization pseudopressure is represented,Bottom pressure Regularization pseudopressure is represented, q represents note Tolerance, CtiRepresent the original gas compressed coefficient, tcaThe average gas well material balance time is represented, G represents gas injection well control reserves, μ generations Table gas viscosity, B represent gas volume factor, and i represents initial data, and K represents in-place permeability, and h represents gas-bearing net pay, CAForm factor is represented, γ represents gas relative density, rwWellbore radius is represented, A represents well control area, reDRepresent zero dimension well Radius is controlled, M represents match point footmark at M.
- 12. gas storage gas injection well control dynamic evaluation method according to claim 8, it is characterised in that the skin factor Calculation formula be:<mrow> <msub> <mi>r</mi> <mrow> <mi>w</mi> <mi>a</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <mfrac> <mrow> <mn>2</mn> <mi>K</mi> <mo>/</mo> <mrow> <mo>(</mo> <mi>&phi;</mi> <mi>&mu;</mi> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>S</mi> <mi>w</mi> </msub> </mrow> <mo>)</mo> <msub> <mi>C</mi> <mrow> <mi>t</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>(</mo> <msubsup> <mi>r</mi> <mrow> <mi>e</mi> <mi>D</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <msub> <mi>lnr</mi> <mrow> <mi>e</mi> <mi>D</mi> </mrow> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> </mfrac> <msub> <mrow> <mo>(</mo> <mfrac> <msub> <mi>t</mi> <mi>c</mi> </msub> <msub> <mi>t</mi> <mrow> <mi>c</mi> <mi>D</mi> </mrow> </msub> </mfrac> <mo>)</mo> </mrow> <mi>M</mi> </msub> </mrow> </msqrt> <mo>,</mo> </mrow><mrow> <mi>s</mi> <mo>=</mo> <mi>ln</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>r</mi> <mi>w</mi> </msub> <msub> <mi>r</mi> <mrow> <mi>w</mi> <mi>a</mi> </mrow> </msub> </mfrac> <mo>)</mo> </mrow> <mo>,</mo> </mrow>Wherein, rwaEffective hole diameter is represented, K represents in-place permeability, and φ represents formation porosity, and μ represents gas viscosity, SwRepresent Water saturation, CtiRepresent the original gas compressed coefficient, tcRepresent gas well material balance time, tcDRepresent gas well zero dimension material Equilibration time, M represent match point footmark at M, reDRepresent zero dimension well control radius, S is pull-type variable, rwRepresent wellbore radius.
- 13. gas storage gas injection well control dynamic evaluation method according to claim 8, it is characterised in that the pressure fitting essence Spend and the fitting of yield fitting precision is with reference to formula:<mrow> <msub> <mi>q</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mn>2.714</mn> <mo>&times;</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mn>5</mn> </mrow> </msup> <msub> <mi>khT</mi> <mrow> <mi>a</mi> <mi>v</mi> </mrow> </msub> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mrow> <mi>w</mi> <mi>f</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>p</mi> <mi>r</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>&mu;</mi> <mi>g</mi> </msub> <msub> <mi>Z</mi> <mrow> <mi>a</mi> <mi>v</mi> </mrow> </msub> <msub> <mi>P</mi> <mrow> <mi>a</mi> <mi>v</mi> </mrow> </msub> <msub> <mi>T</mi> <mi>f</mi> </msub> <mrow> <mo>(</mo> <mi>l</mi> <mi>n</mi> <mfrac> <mrow> <mn>0.472</mn> <msub> <mi>r</mi> <mi>e</mi> </msub> </mrow> <msub> <mi>r</mi> <mi>w</mi> </msub> </mfrac> <mo>+</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>Wherein, qi represents gas injection rate, and K represents in-place permeability, and h represents gas-bearing net pay, TavPost of taking offence in pit shaft is represented to put down Equal temperature, PwfRepresent bottom pressure, PrRepresent pseudoreduced pressure, μgRepresent gas viscosity, ZavPost of taking offence in pit shaft is represented to be averaged Deviation factor, PavRepresent post average pressure of being taken offence in pit shaft, TfRepresent reservoir temperature, S is pull-type variable, rwRepresent pit shaft half Footpath, reRepresent well control radius.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610390275.XA CN107463761A (en) | 2016-06-02 | 2016-06-02 | gas storage gas injection well control dynamic evaluation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610390275.XA CN107463761A (en) | 2016-06-02 | 2016-06-02 | gas storage gas injection well control dynamic evaluation method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107463761A true CN107463761A (en) | 2017-12-12 |
Family
ID=60544794
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610390275.XA Pending CN107463761A (en) | 2016-06-02 | 2016-06-02 | gas storage gas injection well control dynamic evaluation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107463761A (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108804819A (en) * | 2018-06-10 | 2018-11-13 | 西南石油大学 | A kind of low permeability gas reservoirs dynamic holdup evaluation method |
CN111155989A (en) * | 2020-01-14 | 2020-05-15 | 成都理工大学 | Method for predicting drying range of reservoir around injection and production well of gas reservoir type gas storage |
CN111927444A (en) * | 2020-08-31 | 2020-11-13 | 中国石油集团渤海钻探工程有限公司 | Method for evaluating gas injection capacity of depleted oil-gas reservoir gas storage |
CN112395724A (en) * | 2019-07-30 | 2021-02-23 | 中国石油天然气股份有限公司 | Method and system for predicting hydraulic fracturing stratum effect |
CN113047815A (en) * | 2021-04-02 | 2021-06-29 | 西南石油大学 | Reasonable production allocation method for gas well |
CN113294149A (en) * | 2020-08-05 | 2021-08-24 | 中国石油天然气股份有限公司 | Method and device for determining correlation between bottom hole flowing pressure and injection and production amount of injection and production well |
CN113323648A (en) * | 2020-02-12 | 2021-08-31 | 中国石油天然气股份有限公司 | Method and device for determining unobstructed flow of gas well |
CN113567316A (en) * | 2020-04-29 | 2021-10-29 | 中国石油天然气股份有限公司 | Dynamic evaluation method for gas utilization in gas storage rock pore space |
CN113931620A (en) * | 2020-07-13 | 2022-01-14 | 中国石油天然气股份有限公司 | Method and device for calculating gas suction index model of gas injection development oil reservoir and storage medium |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104500030A (en) * | 2014-11-21 | 2015-04-08 | 中国石油天然气股份有限公司 | Ultra-high pressure gas reservoir production dynamic abnormal data diagnosis and correction method |
CN104989351A (en) * | 2015-06-01 | 2015-10-21 | 四川大学 | Dryness, temperature and pressure coupling predicting method in oil-gas well gas injection process |
-
2016
- 2016-06-02 CN CN201610390275.XA patent/CN107463761A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104500030A (en) * | 2014-11-21 | 2015-04-08 | 中国石油天然气股份有限公司 | Ultra-high pressure gas reservoir production dynamic abnormal data diagnosis and correction method |
CN104989351A (en) * | 2015-06-01 | 2015-10-21 | 四川大学 | Dryness, temperature and pressure coupling predicting method in oil-gas well gas injection process |
Non-Patent Citations (5)
Title |
---|
LUO RUILAN等: "a new method for forecasting dynamic reserve of fractured gas well in low-permeability and tight gas reservoir", 《CPS/SPE INTERNATIONAL OIL & GAS CONFERENCE AND EXHIBITION IN CHINA》 * |
孙仁远等: "煤层气井的井底流压计算方法", 《石油钻采工艺》 * |
孙贺东: "《油气井现代产量递减分析方法及应用》", 31 July 2013, 《石油工业出版社》 * |
孙贺东等: "油气井现代产量递减分析方法及在油气藏动态监测中的应用", 《2014年全国天然气学术年会》 * |
胡俊坤等: "高压超高压气藏天然气等温压缩系数的研究", 《石油天然气学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108804819A (en) * | 2018-06-10 | 2018-11-13 | 西南石油大学 | A kind of low permeability gas reservoirs dynamic holdup evaluation method |
CN112395724A (en) * | 2019-07-30 | 2021-02-23 | 中国石油天然气股份有限公司 | Method and system for predicting hydraulic fracturing stratum effect |
CN111155989A (en) * | 2020-01-14 | 2020-05-15 | 成都理工大学 | Method for predicting drying range of reservoir around injection and production well of gas reservoir type gas storage |
CN113323648A (en) * | 2020-02-12 | 2021-08-31 | 中国石油天然气股份有限公司 | Method and device for determining unobstructed flow of gas well |
CN113567316A (en) * | 2020-04-29 | 2021-10-29 | 中国石油天然气股份有限公司 | Dynamic evaluation method for gas utilization in gas storage rock pore space |
CN113567316B (en) * | 2020-04-29 | 2023-09-26 | 中国石油天然气股份有限公司 | Dynamic assessment method for gas utilization of pore space of gas storage rock |
CN113931620A (en) * | 2020-07-13 | 2022-01-14 | 中国石油天然气股份有限公司 | Method and device for calculating gas suction index model of gas injection development oil reservoir and storage medium |
CN113931620B (en) * | 2020-07-13 | 2023-06-30 | 中国石油天然气股份有限公司 | Method, device and storage medium for calculating gas injection development oil reservoir gas suction index model |
CN113294149A (en) * | 2020-08-05 | 2021-08-24 | 中国石油天然气股份有限公司 | Method and device for determining correlation between bottom hole flowing pressure and injection and production amount of injection and production well |
CN113294149B (en) * | 2020-08-05 | 2023-07-25 | 中国石油天然气股份有限公司 | Method and device for determining association relation between bottom hole flow pressure and injection and production amount of injection and production well |
CN111927444A (en) * | 2020-08-31 | 2020-11-13 | 中国石油集团渤海钻探工程有限公司 | Method for evaluating gas injection capacity of depleted oil-gas reservoir gas storage |
CN113047815A (en) * | 2021-04-02 | 2021-06-29 | 西南石油大学 | Reasonable production allocation method for gas well |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107463761A (en) | gas storage gas injection well control dynamic evaluation method | |
Ranjith et al. | Effective parameters for sand production in unconsolidated formations: An experimental study | |
CN108804819A (en) | A kind of low permeability gas reservoirs dynamic holdup evaluation method | |
CN105740563B (en) | Preferential channel identification method for secondary development of mature oil field | |
CN105840187A (en) | Method for calculating staged fracturing productivity of compact reservoir horizontal well | |
CN106761733A (en) | A kind of horizontal wells in heavy oil reservoir steam soak initial productivity Forecasting Methodology | |
CN107313759B (en) | Hypotonic heavy crude reservoir straight well thermal recovery pressure distribution forecasting method and system | |
CN103590824A (en) | Productivity calculation method for tight gas horizontal wells modified by multi-stage fracturing | |
Chai et al. | Modelling prefabricated vertical drain improved ground in plane strain analysis | |
CN107679293A (en) | One kind is applied to Rockfill Dam material interior mechanics parameter and the mutually preferable method of roller compaction construction parameter | |
Rivière et al. | Experimental and numerical assessment of transient stream–aquifer exchange during disconnection | |
CN104963661B (en) | A kind of gas condensate reservoir production split method and device | |
CN203929557U (en) | A kind of gas bearing shale crack develops and seepage flow characteristics proving installation | |
Liu et al. | Numerical modeling of gas flow in coal using a modified dual-porosity model: a multi-mechanistic approach and finite difference method | |
CN105089566A (en) | Gas well system production allocation method | |
Lou et al. | High-precision nonisothermal transient wellbore drift flow model suitable for the full flow pattern domain and full dip range | |
Wang et al. | Analytical solutions for a quad-linear flow model derived for multistage fractured horizontal wells in tight oil reservoirs | |
CN104989385B (en) | The HTHP oil gas straight well perforating parameter optimization method calculated based on skin factor | |
CN108661631A (en) | A kind of shale gas well yield prediction technique | |
Boulon et al. | Influence of rock joint degradation on hydraulic conductivity | |
CN110321575A (en) | Gas condensate reservoir Dynamic Productivity Calculation prediction technique | |
CN114201932A (en) | Well testing simulation method for tight reservoir fracturing well under complex condition | |
CN108629471A (en) | The method that gas storage gas well interaction note is adopted | |
CN105298488A (en) | Diversion capacity testing method under non-continuous filling mode | |
CN103397874B (en) | The solid sand evaluation experimental device of chemical sand control |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171212 |