CN108133086A - Water Fractured Gas Wells fracture half-length's inversion method is produced in a kind of stress sensitive reservoir - Google Patents
Water Fractured Gas Wells fracture half-length's inversion method is produced in a kind of stress sensitive reservoir Download PDFInfo
- Publication number
- CN108133086A CN108133086A CN201711304191.0A CN201711304191A CN108133086A CN 108133086 A CN108133086 A CN 108133086A CN 201711304191 A CN201711304191 A CN 201711304191A CN 108133086 A CN108133086 A CN 108133086A
- Authority
- CN
- China
- Prior art keywords
- water
- gas
- pressure
- phase
- stress sensitive
- 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.)
- Withdrawn
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 166
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000035699 permeability Effects 0.000 claims abstract description 85
- 238000005325 percolation Methods 0.000 claims abstract description 43
- 238000004519 manufacturing process Methods 0.000 claims abstract description 24
- 230000006872 improvement Effects 0.000 claims abstract description 12
- 239000012530 fluid Substances 0.000 claims abstract description 11
- 230000000704 physical effect Effects 0.000 claims abstract description 7
- 230000000694 effects Effects 0.000 claims abstract description 6
- 238000012417 linear regression Methods 0.000 claims abstract description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 34
- 230000010354 integration Effects 0.000 claims description 12
- 238000000205 computational method Methods 0.000 claims description 7
- 239000011148 porous material Substances 0.000 claims description 7
- 238000013459 approach Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 5
- 238000005452 bending Methods 0.000 claims description 4
- 239000011435 rock Substances 0.000 claims description 3
- 238000009738 saturating Methods 0.000 claims description 3
- 239000002352 surface water Substances 0.000 claims description 3
- 230000008595 infiltration Effects 0.000 claims description 2
- 238000001764 infiltration Methods 0.000 claims description 2
- 239000007788 liquid Substances 0.000 claims description 2
- 230000006835 compression Effects 0.000 claims 1
- 238000007906 compression Methods 0.000 claims 1
- 238000012360 testing method Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000012466 permeate Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 210000000867 larynx Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Theoretical Computer Science (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The present invention relates to water Fractured Gas Wells fracture half-length's inversion method is produced in a kind of stress sensitive reservoir, include the following steps:(1) gas well daily creation data, reservoir properties, physical properties of fluids, stress sensitive data are compiled;(2) the equivalent relative permeability of air water under the conditions of different pressures is calculated;(3) it calculates and considers that parameter is intended in the improvement of stress sensitive reservoir air water two-phase complicated percolation feature;(4) it draws and utilizes time graph identification stratum linear flow under improved radical sign;(5) data point in layer linear flow stage carries out linear regression over the ground, based on obtained straight slope, is calculated using inverse model, realizes half long inverting of fracture.The complicated percolations feature such as phase percolation curve stress sensitive, absolute permeability stress sensitive, slippage effect when present invention can consider air water two phase fluid flow in stress sensitive reservoir, the error caused by Traditional calculating methods is eliminated, is evaluated after can be widely used for the pressure of stress sensitive reservoir production water Fractured Gas Wells.
Description
Technical field
Hydraulic pressure, which is produced, the present invention relates to assessment technique field after hydraulic fracturing pressure, in especially a kind of stress sensitive reservoir splits gas
Well fracture half-length's inversion method.
Background technology
In stress sensitive reservoir production process, since reservoir stress condition constantly changes, reservoir pore throat character can be sent out therewith
It is raw to change.And for the water-producing gas well in stress sensitive reservoir, pore throat character can lead to Gas-water phases with the variation of stress
The complicated percolations such as absolute permeability stress sensitive, phase percolation curve stress sensitive, dynamic slippage effect spy is showed in flow event
Sign.On the other hand, the permeability of stress sensitive reservoir is usually relatively low, needs that by fracturing reform industrial gas could be obtained.And
The fracture half-length that pressure break is formed be evaluate fracturing effect important indicator, the formulation to subsequent capability forecasting and development plan
It has a very important significance.The presence of Gas-water phases complicated percolation feature in stress sensitive reservoir so that production water fractured well is split
There are larger differences with conventional gas well for half long inverting of seam.
But presently, there are some fracture half-length's inversion methods, do not consider stress sensitive reservoir Gas-water phases
Complicated percolation feature, particularly ignores the stress sensitive of air-water phase percolation curve, obtained fracture half-length's value and actual value it
Between usually there are certain difference, therefore be not suitable in stress sensitive reservoir produce water Fractured Gas Wells fracture half-length's inverting.
Invention content
The technical problem to be solved by the present invention is to:In order to overcome deficiency in the prior art, the present invention provides a kind of stress
In sensitive reservoir produce water Fractured Gas Wells fracture half-length's inversion method, this method considered the gas such as phase percolation curve stress sensitive-
Water two phase fluid flow feature, so as to improve the accuracy of crack inversion result.
The technical solution adopted by the present invention to solve the technical problems is:Water Fractured Gas Wells are produced in a kind of stress sensitive reservoir
Fracture half-length's inversion method has following steps:
(1), the daily creation data of gas well, reservoir properties, physical properties of fluids, air-water mutually ooze data, the receipts of stress sensitive data
Collection arranges;
(2), consider phase percolation curve stress sensitive, absolute permeability stress sensitive and dynamic slippage effect, calculate not
With the equivalent air-water relative permeability of Gas-water phases complicated percolation feature under pressure condition;
(3), it using equivalent air-water phase percolation curve, with reference to gas well liquid loading data, according to the following formula, calculates and considers stress sensitive
The improved pseudopressure of reservoir Gas-water phases complicated percolation feature and pseudotime;
Wherein, ψtwoTo consider the improved pseudopressure MPa/cp of stress sensitive reservoir Gas-water phases complicated percolation feature;
ttwoTo consider the improved pseudotime d of stress sensitive reservoir Gas-water phases complicated percolation feature;P is pressure MPa;paFor reference
Pressure MPa;T is actual time d;taFor reference time d;ρgFor gas phase lower density kg/m3;ρwGround lower density kg/m for water3;
ρwscGround standard density kg/m for water3;For the gas phase corresponding to average pressure lower density kg/m3;For mean pressure
Water phase ground standard density kg/m corresponding to power3;krgEFor the equivalent relative permeability of gas phase;krwEFor the equivalent opposite infiltration of water phase
Rate;For the equivalent relative permeability of gas phase corresponding to average pressure;It is equivalent opposite for the water phase corresponding to average pressure
Permeability;μgFor gas viscosity cp;μwViscosity cp for water;μgiFor the gas viscosity cp under original formation pressure;It is average
Gas viscosity corresponding to pressure, cp;Viscosity cp for the water corresponding to average pressure;Ct-twoiFor under original formation pressure
Air water two-phase system compressibility MPa-1;For the air water two-phase system compressibility MPa corresponding to average pressure-1;
(4), it based on the improved pseudopressure obtained and pseudotime is calculated, draws and utilizes time graph under improved radical sign
Identify stratum linear flow;
(5), linear regression is carried out to the scatterplot on curve, according to obtained straight slope is returned, utilizes inverse model
It is calculated, so as to fulfill half long inverting of fracture.
In above-mentioned steps (1), the necessary data compiled is needed to include:Include water yield qw, gas production qg, water-gas ratio
fw, casing pressure pcOr flowing bottomhole pressure (FBHP) pwfGas well Field Production Data inside;Hole under the conditions of reservoir thickness h and original formation pressure
Porosity φi, absolute permeability ki, gas slip factor biWait reservoir physical parameters;Fluid viscosity μ under the conditions of different pressuresg、μw
And compressed coefficient Cg、CwEtc. physical properties of fluids parameter;Phase percolation curve endpoint value (S corresponding to original formation pressurewci、Sgri、
krgendi、krwendi), the phase percolation curves parameter such as the non-isometrical coefficient lambda of reservoir pore throat and bending coefficient η;Reservoir absolute permeability, hole
Stress sensitive coefficients (α, γ, C, D, E, F) of degree and phase percolation curve endpoint value etc.;
In step (2), the computation model of equivalent air-water relative permeability is:
Wherein,AndCalculation formula be respectively:
Wherein, krgE-p1<Sw>For a certain pressure (p1) and saturation degree (Sw) under the equivalent relative permeability of gas phase;krwE-p1<
Sw>For a certain pressure (p1) and saturation degree (Sw) under the equivalent relative permeability of water phase;krg-p1<Sw>For a certain pressure (p1) and it is full
With degree (Sw) under gas phase relative permeability;krw-p1<Sw>For a certain pressure (p1) and saturation degree (Sw) under water phase opposite permeate
Rate;biFor the slip factor corresponding to original formation pressure;B is slip factor regression coefficient;α is permeability stress sensitive coefficients
MPa-1;krgendiEndpoint value is mutually oozed for the gas phase under original formation pressure;SwFor water saturation;SwciFor under original formation pressure
Irreducible water saturation;SgriFor the residual gas saturation under original formation pressure;λ is hollow billet profile exponent;η is bent for hollow billet
Coefficient;krwendiEndpoint value is mutually oozed for the water phase under original formation pressure;C mutually oozes endpoint value stress sensitive coefficients, unit for gas phase
For MPa-1;D mutually oozes endpoint value stress sensitive coefficients, unit MPa for water phase-1;E is irreducible water saturation stress sensitive coefficients,
Unit is MPa-1;F be residual gas saturation stress sensitive coefficients, unit MPa-1;P is to specify strata pressure, unit MPa;
piFor original formation pressure, unit MPa;
In step (3), the computational methods of pseudopressure are:
(a), using each time point t corresponding to gas production qgAnd water yield qw, calculate corresponding water yield of each time point with
Gas production ratio qw/qg;
(b), with time t1For, first pressure spot p is selected in limit of integration1, based on gas phase and water phase physical parameter
With the relational expression of pressure, related property parameter value (μ corresponding to the pressure value is calculatedg、Bg、μw、Bw);
(c), using following formula, the pressure p is calculated1Under the conditions of water phase and the equivalent relative permeability of gas phase ratio krwE/krgE;
(d), using the equivalent phase percolation curve under the conditions of the different pressures established in step (2), each pressure is further determined that
The equivalent relative permeability of water phase and the equivalent relative permeability ratio (k of gas phase under the conditions of powerrwE/krgE) with the variation of water saturation
Feature;
(e), with reference to step (c) and step (d) as a result, determining the corresponding water saturation value of the pressure spot, and utilize
The saturation degree calculates pressure spot p1The corresponding equivalent relative permeability k of gas phasergEAnd the equivalent relative permeability k of water phaserwE;
(f), next pressure spot p is chosen in limit of integration2, analyzed using similar approach, obtain pressure spot p2It is right
The equivalent relative permeability of gas phase and the equivalent relative permeability ... of water phase answered and so on, if the pressure chosen in limit of integration
Force quantity is enough, you can establishes pressure value and gas phase and the equivalent relative permeability k of water phasergEOr krwERelational expression, so as to
The numerical integration realized to relevant parameter is substituted into following formula, and obtains changing suitable for stress sensitive reservoir Gas-water phases seepage flow
Into pseudopressure;
Wherein, qgFor surface gas rate, m3/d;qwFor surface water rate, m3/d;BgFor gas volume factor;BwFor water
Volume factor.
In step (3), the computational methods of pseudotime are:
(a), time point t is selected in limit of integration1, it is true by pressure propagation range formula or stream material equilibria method
Mean reservoir pressure in constant-pressure coverage or reservoirAnd then determine the corresponding related air-water physical property of the average pressure
Parameter value (μg、Bg、μw、Bw);
(c), based on the equivalent phase percolation curve in step (2), the equivalent of water phase and gas phase under the conditions of the average pressure is obtained
Relative permeability ratioWith the variation characteristic of water saturation;
(d), it with reference to step (b) and (c) as a result, average staturation value corresponding to determining average pressure, and then calculates
Average pressureCorresponding average air equivalent relative permeabilityAnd average water equivalent relative permeability
(e), next time point t is chosen in limit of integration2, analyzed using similar approach, obtain t2Moment corresponds to
Gas phase average equivalent relative permeability and water phase average equivalent relative permeability ..., if being selected in limit of integration
The time point quantity taken is enough, you can settling time and average gas phase and the equivalent relative permeability of water phasePass
It is formula, so as to substitute into following formula the numerical integration realized to relevant parameter, and obtains and be suitable for stress sensitive reservoir Gas-water phases
The improvement pseudotime of seepage flow;
Wherein, t is actual time, d;For mean reservoir pressure, MPa;It is equivalent for the gas phase corresponding to average pressure
Relative permeability;For the equivalent relative permeability of gas phase corresponding to average pressure.
In step (3), the computational methods of Gas-water phases seepage flow system compressibility in the reservoir in pseudotime calculation formula
For:
Wherein, Ct-twoFor air water two-phase system compressibility, MPa-1;SgFor gas phase saturation;CgFor gas compressibility factor,
MPa-1;CwFor the compressed coefficient of water, MPa-1;CpFor pore space compressibility of rock, MPa-1。
In step (4), when drawing time graph under improved radical sign, when gas well is determines air content production, abscissa isAnd when gas well gas production is non-constant, abscissa isAnd ordinate isBecome gas production feelings
Abscissa under conditionIt is calculated using following formula.
Wherein, ψtwoiFor the improvement pseudopressure corresponding to original formation pressure;ψtwowfFor the improvement corresponding to flowing bottomhole pressure (FBHP)
Pseudopressure;qtwoFor air-water total output;tsl-twoFor the improved material balance pseudotime corresponding to the stratum linear flow stage, use
Determine production conditions in varying flow rate condition is converted into, unit is day (d);N is the total quantity of data point.
Wherein, mtwo-sBy the straight line that time graph overlying strata linear flow number of stages strong point returns under improved radical sign
Slope;xfFor fracture half-length m;φ is porosity;H is gas-bearing net pay m;kiIt is oozed for original under the conditions of original formation pressure
Saturating rate mD.
Pass through above step, you can realize the inverting to producing water Fractured Gas Wells fracture half-length in stress sensitive reservoir.
The beneficial effects of the invention are as follows:The present invention compared with current inversion method, has due to taking more than technical method
Advantages below and advantageous effect:
1. this method has considered the complicated percolations feature such as phase percolation curve stress sensitive, therefore relative to current main-stream
Fracture half-length's inversion method, is more nearly the practical seepage flow situation of reservoir, and obtained fracture half-length is more accurate;
2. this method employs semi-analytic method, easier relative to numerical analysis method, speed is asked for faster;
3. this method is analyzed based on daily creation data, solve traditional well testing inversion method and excessively rely on now
The drawbacks of field well testing test, without carrying out closing well test, significantly reduce the cost of fracture half-length's inverting.
Description of the drawings
The present invention is further described with reference to the accompanying drawings and examples.
Fig. 1 is flow chart of the present invention.
Fig. 2 is the aqueous vapor when flowing bottomhole pressure (FBHP) data in the embodiment of the present invention.
Fig. 3 is the equivalent phase percolation curve in the embodiment of the present invention.
Fig. 4 is time graph under improved radical sign in the embodiment of the present invention.
Specific embodiment
In conjunction with the accompanying drawings, the present invention is further explained in detail.These attached drawings are simplified schematic diagram, only with
Illustration illustrates the basic structure of the present invention, therefore it only shows composition related to the present invention.
It is the stream that water Fractured Gas Wells fracture half-length's inversion method is produced in a kind of stress sensitive reservoir of the present invention as shown in Figure 1
Cheng Tu, its step are as follows:
In step (1), the necessary data compiled includes:Water yield qw, gas production qg, water-gas ratio fw, flowing bottomhole pressure (FBHP)
pwfGas well Field Production Data Deng including;Porosity φ under the conditions of reservoir thickness h and original formation pressurei, Absolute permeation
Rate ki, gas slip factor biWait reservoir physical parameters;Fluid viscosity μ under the conditions of different pressuresg、μwAnd compressed coefficient Cg、Cw
Etc. physical properties of fluids parameter;Phase percolation curve endpoint value (S corresponding to original formation pressurewci、Sgri、krgendi、krwendi), reservoir hole
The phase percolation curves parameter such as the non-isometrical coefficient lambda of larynx and bending coefficient η;Reservoir absolute permeability, porosity and phase percolation curve endpoint value
Stress sensitive coefficients (α, γ, C, D, E, F) etc..Wherein, the when flowing bottomhole pressure (FBHP) data such as Fig. 2 institutes of the aqueous vapor in the present embodiment
Show, the other underlying parameters in the part being collected into are as shown in table 1.
Table 1
In step (2), equivalent relative permeability is calculated based on following formula, and draw corresponding equivalent phase percolation curve, as a result
As shown in Figure 3:
Wherein,AndCalculation formula be respectively:
Wherein, krgE-p1<Sw>For a certain pressure (p1) and saturation degree (Sw) under the equivalent relative permeability of gas phase;krwE-p1<
Sw>For a certain pressure (p1) and saturation degree (Sw) under the equivalent relative permeability of water phase;krg-p1<Sw>For a certain pressure (p1) and it is full
With degree (Sw) under gas phase relative permeability;krw-p1<Sw>For a certain pressure (p1) and saturation degree (Sw) under water phase opposite permeate
Rate;biFor the slip factor corresponding to original formation pressure, MPa;B is slip factor regression coefficient;α is permeability stress sensitive
Coefficient, MPa-1;krgendiEndpoint value is mutually oozed for the gas phase under original formation pressure;SwFor water saturation;SwciFor prime stratum
Irreducible water saturation under pressure;SgriFor the residual gas saturation under original formation pressure;λ is hollow billet profile exponent;η is hair
Pipe bending coefficient;krwendiEndpoint value is mutually oozed for the water phase under original formation pressure;C mutually oozes endpoint value stress sensitive system for gas phase
Number;D mutually oozes endpoint value stress sensitive coefficients for water phase;E is irreducible water saturation stress sensitive coefficients;F is residual gas saturation
Stress sensitive coefficients;P is to specify strata pressure, unit MPa;piFor original formation pressure, unit MPa;
In step (3), the improvement pseudopressure for considering Gas-water phases complicated percolation feature is calculated according to below scheme:
(a), using each time point t corresponding to gas production qgAnd water yield qw, calculate corresponding water yield of each time point with
Gas production ratio qw/qg;
(b), with time t1For, first pressure spot p is selected in limit of integration1, based on gas phase and water phase physical parameter
With the relational expression of pressure, related property parameter value (μ corresponding to the pressure value is calculatedg、Bg、μw、Bw);
(c), using following formula, the pressure p is calculated1Under the conditions of water phase and the equivalent relative permeability of gas phase ratio krwE/krgE;
(d), using the equivalent phase percolation curve under the conditions of the different pressures established in step (2), each pressure is further determined that
The equivalent relative permeability of water phase and the equivalent relative permeability ratio (k of gas phase under the conditions of powerrwE/krgE) with the variation of water saturation
Feature;
(e), with reference to step (c) and step (d) as a result, determining the corresponding water saturation value of the pressure spot, and utilize
The saturation degree calculates pressure spot p1The corresponding equivalent relative permeability k of gas phasergEAnd the equivalent relative permeability k of water phaserwE;
(f), next pressure spot p is chosen in limit of integration2, analyzed using similar approach, obtain pressure spot p2It is right
The equivalent relative permeability of gas phase and the equivalent relative permeability ... of water phase answered and so on, if the pressure chosen in limit of integration
Force quantity is enough, you can establishes pressure value and gas phase and the equivalent relative permeability k of water phasergEOr krwERelational expression, so as to
The numerical integration realized to relevant parameter is substituted into following formula, and obtains changing suitable for stress sensitive reservoir Gas-water phases seepage flow
Into pseudopressure,
Wherein, ψtwoTo consider the improved pseudopressure MPa/cp of stress sensitive reservoir Gas-water phases complicated percolation feature;p
For pressure MPa;paFor reference pressure MPa;ρgFor gas phase lower density kg/m3;ρwGround lower density kg/m for water3;ρwscFor water
Ground standard density kg/m3;krgEFor the equivalent relative permeability of gas phase;krwEFor the equivalent relative permeability of water phase;μgFor gas viscosity
cp;μwViscosity cp for water;qgFor surface gas rate m3/d;qwFor surface water rate m3/d;BgFor gas volume factor;BwFor water
Volume factor.
In step (3), the computational methods of pseudotime are:
(a), time point t is selected in limit of integration1, it is true by pressure propagation range formula or stream material equilibria method
Mean reservoir pressure in constant-pressure coverage or reservoirAnd then determine the corresponding related air-water physical property of the average pressure
Parameter value (μg、Bg、μw、Bw);
(c), based on the equivalent phase percolation curve in step (2), the equivalent of water phase and gas phase under the conditions of the average pressure is obtained
Relative permeability ratioWith the variation characteristic of water saturation;
(d), it with reference to step (b) and (c) as a result, average staturation value corresponding to determining average pressure, and then calculates
Average pressureCorresponding average air equivalent relative permeabilityAnd average water equivalent relative permeability
(e), next time point t is chosen in limit of integration2, analyzed using similar approach, obtain t2Moment corresponds to
Gas phase average equivalent relative permeability and water phase average equivalent relative permeability ..., if being selected in limit of integration
The time point quantity taken is enough, you can settling time and average gas phase and the equivalent relative permeability of water phasePass
It is formula, so as to substitute into following formula the numerical integration realized to relevant parameter, and obtains and be suitable for stress sensitive reservoir Gas-water phases
The improvement pseudotime of seepage flow,
Wherein,T is actual time d;For
Mean reservoir pressure MPa;For the equivalent relative permeability of gas phase corresponding to average pressure;For corresponding to average pressure
The equivalent relative permeability of gas phase;Ct-twoFor air water two-phase system compressibility MPa-1;SgFor gas phase saturation;CgFor gas pressure
Contracting coefficient MPa-1;CwCompressed coefficient MPa for water-1;CpFor pore space compressibility of rock MPa-1, ttwoTo consider stress sensitive storage
The improved pseudotime d of layer Gas-water phases complicated percolation feature;taFor reference time d;For the gas phase corresponding to average pressure
Ground lower density kg/m3;For the water phase ground standard density kg/m corresponding to average pressure3;For corresponding to average pressure
The equivalent relative permeability of gas phase;For the equivalent relative permeability of water phase corresponding to average pressure;μgFor gas viscosity cp;
μwViscosity cp for water;For the gas viscosity cp corresponding to average pressure;Viscosity for the water corresponding to average pressure
cp;Ct-twoiFor the air water two-phase system compressibility MPa under original formation pressure-1;For the gas corresponding to average pressure
Water two-phase system compressibility MPa-1。
In step (4), due to the present embodiment be fixed output quota tolerance production, using abscissa asAnd ordinate isTime graph under improved radical sign is drawn, as shown in Figure 4.
Wherein, ψtwoiFor the improvement pseudopressure corresponding to original formation pressure;ψtwowfFor the improvement corresponding to flowing bottomhole pressure (FBHP)
Pseudopressure;qtwoFor air-water total output.
In step (5), based on the obtained straight slope m of recurrencetwo-s, utilize the computation model of fracture half-lengthThe fracture half-length that production water Fractured Gas Wells in stress sensitive reservoir are calculated is 101.02m.
Wherein, mtwo-sBy the straight line that time graph overlying strata linear flow number of stages strong point returns under improved radical sign
Slope;xfFor fracture half-length m;φ is porosity;H is gas-bearing net pay m;kiIt is oozed for original under the conditions of original formation pressure
Saturating rate mD.
Using above-mentioned desirable embodiment according to the present invention as enlightenment, by above-mentioned description, relevant staff is complete
Various changes and amendments can be carried out without departing from the scope of the technological thought of the present invention' entirely.The technology of this invention
Property range is not limited to the content on specification, it is necessary to determine its technical scope according to right.
Claims (8)
1. water Fractured Gas Wells fracture half-length's inversion method is produced in a kind of stress sensitive reservoir, which is characterized in that include the following steps:
(1), it is whole mutually to ooze data, the collection of stress sensitive data for the daily creation data of gas well, reservoir properties, physical properties of fluids, air-water
Reason;
(2), consider phase percolation curve stress sensitive, absolute permeability stress sensitive and dynamic slippage effect, calculate different pressures
The equivalent air-water relative permeability of Gas-water phases complicated percolation feature under the conditions of power;
(3), it using equivalent air-water phase percolation curve, with reference to gas well liquid loading data, according to the following formula, calculates and considers stress sensitive reservoir
The improved pseudopressure of Gas-water phases complicated percolation feature and pseudotime;
Wherein, ψtwoTo consider the improved pseudopressure MPa/cp of stress sensitive reservoir Gas-water phases complicated percolation feature;ttwoFor
Consider the improved pseudotime d of stress sensitive reservoir Gas-water phases complicated percolation feature;P is pressure MPa;paFor reference pressure
MPa;ρgFor gas phase lower density kg/m3;T is actual time d;taFor reference time d;ρwGround lower density kg/m for water3;ρwsc
Ground standard density kg/m for water3;For the gas phase corresponding to average pressure lower density kg/m3;For average pressure institute
Corresponding water phase ground standard density kg/m3;krgEFor the equivalent relative permeability of gas phase;krwEFor the equivalent relative permeability of water phase;For the equivalent relative permeability of gas phase corresponding to average pressure;It opposite oozes for the water phase corresponding to average pressure is equivalent
Saturating rate;μgFor gas viscosity cp;μwViscosity cp for water;μgiFor the gas viscosity cp under original formation pressure;For mean pressure
Gas viscosity cp corresponding to power;Viscosity cp for the water corresponding to average pressure;Ct-twoiFor under original formation pressure
Air water two-phase system compressibility MPa-1;For the air water two-phase system compressibility MPa corresponding to average pressure-1;
(4), it based on the improved pseudopressure obtained and pseudotime is calculated, draws and is identified using time graph under improved radical sign
Stratum linear flow;
(5), the data point in layer linear flow stage carries out linear regression over the ground, according to returning obtained straight slope, using anti-
It drills model to be calculated, so as to fulfill half long inverting of fracture.
2. producing water Fractured Gas Wells fracture half-length's inversion method in stress sensitive reservoir according to claim 1, feature exists
In:The necessary data compiled is needed to include in step (1):Include water yield qw, gas production qg, water-gas ratio fw, casing pressure pc
Or flowing bottomhole pressure (FBHP) pwfGas well Field Production Data inside;Porosity φ under the conditions of reservoir thickness h and original formation pressurei、
Absolute permeability ki, gas slip factor biReservoir physical parameter;Fluid viscosity μ under the conditions of different pressuresg、μwAnd compression system
Number Cg、CwPhysical properties of fluids parameter;Phase percolation curve endpoint value (S corresponding to original formation pressurewci、Sgri、krgendi、krwendi), storage
The non-isometrical coefficient lambda of layer pore throat and bending coefficient η phase percolation curve parameters;Reservoir absolute permeability, porosity and phase percolation curve endpoint
The stress sensitive coefficients (α, γ, C, D, E, F) of value.
3. producing water Fractured Gas Wells fracture half-length's inversion method in stress sensitive reservoir according to claim 1, feature exists
In:The computation model of equivalent air-water relative permeability is in step (2):
Wherein,AndCalculation formula be respectively:
Wherein, krgE-p1<Sw>For a certain pressure (p1) and saturation degree (Sw) under the equivalent relative permeability of gas phase;krwE-p1<Sw>For
A certain pressure (p1) and saturation degree (Sw) under the equivalent relative permeability of water phase;krg-p1<Sw>For a certain pressure (p1) and saturation degree
(Sw) under gas phase relative permeability;krw-p1<Sw>For a certain pressure (p1) and saturation degree (Sw) under water phase relative permeability;bi
For the slip factor MPa corresponding to original formation pressure;B is slip factor regression coefficient;α is permeability stress sensitive coefficients
MPa-1;krgendiEndpoint value is mutually oozed for the gas phase under original formation pressure;SwFor water saturation;SwciFor under original formation pressure
Irreducible water saturation;SgriFor the residual gas saturation under original formation pressure;λ is hollow billet profile exponent;η is bent for hollow billet
Coefficient;krwendiEndpoint value is mutually oozed for the water phase under original formation pressure;C mutually oozes endpoint value stress sensitive coefficients, unit for gas phase
For MPa-1;D mutually oozes endpoint value stress sensitive coefficients, unit MPa for water phase-1;E is irreducible water saturation stress sensitive coefficients,
Unit is MPa-1;F be residual gas saturation stress sensitive coefficients, unit MPa-1;P is to specify strata pressure, unit MPa;
piFor original formation pressure, unit MPa.
4. producing water Fractured Gas Wells fracture half-length's inversion method in stress sensitive reservoir according to claim 1, feature exists
In:The computational methods of pseudopressure are in step (3):
(a), using each time point t corresponding to gas production qgAnd water yield qw, calculate corresponding water yield of each time point and aerogenesis
Measure ratio qw/qg;
(b), with time t1For, first pressure spot p is selected in limit of integration1, based on gas phase and water phase physical parameter and pressure
The relational expression of power calculates related property parameter value (μ corresponding to the pressure valueg、Bg、μw、Bw);
(c), using following formula, the pressure p is calculated1Under the conditions of water phase and the equivalent relative permeability of gas phase ratio krwE/krgE
(d), using the equivalent phase percolation curve under the conditions of the different pressures established in step (2), each pressure strip is further determined that
The equivalent relative permeability of water phase and the equivalent relative permeability ratio (k of gas phase under partrwE/krgE) special with the variation of water saturation
Sign;
(e), with reference to step (c) and step (d) as a result, determining the corresponding water saturation value of the pressure spot, and full using this
And degree, calculate pressure spot p1The corresponding equivalent relative permeability k of gas phasergEAnd the equivalent relative permeability k of water phaserwE;
(f), next pressure spot p is chosen in limit of integration2, analyzed using similar approach, obtain pressure spot p2It is corresponding
The equivalent relative permeability of gas phase and the equivalent relative permeability ... of water phase and so on, if the pressure spot chosen in limit of integration
Quantity is enough, you can establishes pressure value and gas phase and the equivalent relative permeability k of water phasergEOr krwERelational expression, so as to substitute into
The numerical integration to relevant parameter is realized in following formula, and obtains the improvement plan suitable for stress sensitive reservoir Gas-water phases seepage flow
Pressure;
Wherein, qgFor surface gas rate, m3/d;qwFor surface water rate, m3/d;BgFor gas volume factor;BwVolume for water
Coefficient.
5. producing water Fractured Gas Wells fracture half-length's inversion method in stress sensitive reservoir according to claim 1, feature exists
In:The computational methods of pseudotime are in step (3):
(a), time point t is selected in limit of integration1, determine to press by pressure propagation range formula or stream material equilibria method
Mean reservoir pressure in power coverage or reservoirAnd then determine the corresponding related air-water physical parameter of the average pressure
It is worth (μg、Bg、μw、Bw);
(b), using following formula, the equivalent relative permeability ratio of water phase and gas phase under the conditions of average pressure is calculated
(c), based on the equivalent phase percolation curve in step (2), the equivalent opposite of water phase and gas phase under the conditions of the average pressure is obtained
Permeability ratioWith the variation characteristic of water saturation;
(d), it with reference to step (b) and (c) as a result, average staturation value corresponding to determining average pressure, and then calculates average
PressureCorresponding average air equivalent relative permeabilityAnd average water equivalent relative permeability
(e), next time point t is chosen in limit of integration2, analyzed using similar approach, obtain t2Moment corresponding gas
The equivalent relative permeability of phase average and water phase average equivalent relative permeability ... and so on, if chosen in limit of integration
Time point quantity is enough, you can settling time and average gas phase and the equivalent relative permeability of water phaseRelationship
Formula so as to substitute into following formula the numerical integration realized to relevant parameter, and is obtained and is oozed suitable for stress sensitive reservoir Gas-water phases
The improvement pseudotime of stream.
Wherein, t is actual time, d;For mean reservoir pressure, MPa;It is equivalent opposite for the gas phase corresponding to average pressure
Permeability;For the equivalent relative permeability of gas phase corresponding to average pressure.
6. producing water Fractured Gas Wells fracture half-length's inversion method in stress sensitive reservoir according to claim 5, feature exists
In:The computational methods of Gas-water phases seepage flow system compressibility are in reservoir in step (3) in pseudotime calculation formula:
Wherein, Ct-twoFor air water two-phase system compressibility MPa-1;SgFor gas phase saturation;CgFor gas compressibility factor MPa-1;
CwCompressed coefficient MPa for water-1;CpFor pore space compressibility of rock MPa-1。
7. producing water Fractured Gas Wells fracture half-length's inversion method in stress sensitive reservoir according to claim 1, feature exists
In:When time graph under improved radical sign is drawn in step (4), when gas well is determines air content production, abscissa is
And when gas well gas production is non-constant, abscissa isAnd ordinate isBecome in the case of gas production
AbscissaIt is calculated using following formula,
Wherein, ψtwoiFor the improvement pseudopressure corresponding to original formation pressure;ψtwowfIntend pressure for the improvement corresponding to flowing bottomhole pressure (FBHP)
Power;qtwoFor air-water total output;tsl-twoFor the improved material balance pseudotime corresponding to the stratum linear flow stage, for will
Varying flow rate condition, which is converted into, determines production conditions, and unit is day;N is the total quantity of data point.
8. water Fractured Gas Wells fracture half-length's inversion method is produced in stress sensitive reservoir according to claim 1, it is characterised in that:
In step (5), based on the obtained slope m of straight line of recurrencetwo-s, the computation model of fracture half-length is
Wherein, mtwo-sIt is oblique by time graph overlying strata linear flow number of stages strong point returns under improved radical sign straight line
Rate;xfFor fracture half-length m;φ is porosity;H is gas-bearing net pay m;kiFor the original infiltration under the conditions of original formation pressure
Rate mD.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711304191.0A CN108133086A (en) | 2017-12-11 | 2017-12-11 | Water Fractured Gas Wells fracture half-length's inversion method is produced in a kind of stress sensitive reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711304191.0A CN108133086A (en) | 2017-12-11 | 2017-12-11 | Water Fractured Gas Wells fracture half-length's inversion method is produced in a kind of stress sensitive reservoir |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108133086A true CN108133086A (en) | 2018-06-08 |
Family
ID=62389100
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711304191.0A Withdrawn CN108133086A (en) | 2017-12-11 | 2017-12-11 | Water Fractured Gas Wells fracture half-length's inversion method is produced in a kind of stress sensitive reservoir |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108133086A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108830020A (en) * | 2018-07-12 | 2018-11-16 | 西南石油大学 | A method of the micro- Fracturing Technology crack extension of simulation offshore oilfield based on heat flow piercement theory |
CN111027211A (en) * | 2019-12-10 | 2020-04-17 | 中国石油大学(华东) | Dense reservoir fluid simulation method considering pore throat limitation mechanism and application of dense reservoir fluid simulation method in gas injection development simulation |
CN111608647A (en) * | 2020-04-14 | 2020-09-01 | 中国石油化工股份有限公司 | Method for predicting formation pressure of water injection and polymer injection well region |
CN111691869A (en) * | 2019-03-15 | 2020-09-22 | 中国石油化工股份有限公司 | Method for predicting yield change trend of high-water-content compact gas reservoir |
CN112855127A (en) * | 2019-11-28 | 2021-05-28 | 北京国双科技有限公司 | Gas well accumulated liquid identification method and device |
CN113672851A (en) * | 2020-05-14 | 2021-11-19 | 中国石油化工股份有限公司 | Gas reservoir water production gas well parameter prediction method and device, electronic equipment and medium |
CN116579263A (en) * | 2023-05-17 | 2023-08-11 | 中国石油大学(北京) | Comprehensive analysis method based on oil and gas well drainage dynamic data |
-
2017
- 2017-12-11 CN CN201711304191.0A patent/CN108133086A/en not_active Withdrawn
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108830020A (en) * | 2018-07-12 | 2018-11-16 | 西南石油大学 | A method of the micro- Fracturing Technology crack extension of simulation offshore oilfield based on heat flow piercement theory |
CN108830020B (en) * | 2018-07-12 | 2019-05-31 | 西南石油大学 | A method of the micro- Fracturing Technology crack extension of simulation offshore oilfield |
CN111691869A (en) * | 2019-03-15 | 2020-09-22 | 中国石油化工股份有限公司 | Method for predicting yield change trend of high-water-content compact gas reservoir |
CN112855127A (en) * | 2019-11-28 | 2021-05-28 | 北京国双科技有限公司 | Gas well accumulated liquid identification method and device |
CN111027211A (en) * | 2019-12-10 | 2020-04-17 | 中国石油大学(华东) | Dense reservoir fluid simulation method considering pore throat limitation mechanism and application of dense reservoir fluid simulation method in gas injection development simulation |
CN111608647A (en) * | 2020-04-14 | 2020-09-01 | 中国石油化工股份有限公司 | Method for predicting formation pressure of water injection and polymer injection well region |
CN113672851A (en) * | 2020-05-14 | 2021-11-19 | 中国石油化工股份有限公司 | Gas reservoir water production gas well parameter prediction method and device, electronic equipment and medium |
CN116579263A (en) * | 2023-05-17 | 2023-08-11 | 中国石油大学(北京) | Comprehensive analysis method based on oil and gas well drainage dynamic data |
CN116579263B (en) * | 2023-05-17 | 2024-01-30 | 中国石油大学(北京) | Comprehensive analysis method based on oil and gas well drainage dynamic data |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108133086A (en) | Water Fractured Gas Wells fracture half-length's inversion method is produced in a kind of stress sensitive reservoir | |
CN109488276B (en) | Method for predicting shale gas yield of water-producing shale gas well transformed through hydraulic fracturing | |
CN104018829B (en) | A kind of method that utilization coal bed gas well creation data measures air water phase percolation curve | |
CN111353205B (en) | Method for calculating formation pressure and dynamic productivity of water-producing gas well of tight gas reservoir | |
CN104834807B (en) | A kind of stress sensitive reservoir relative permeability computational methods based on fractal theory | |
CN108133087A (en) | A kind of Gas-water phases seepage stress sensitivity reservoir original permeability inversion method | |
CN103912269B (en) | Method for determining formation fracture pressure gradient logging of shale gas reservoir | |
CN107462936B (en) | Utilize the method for pressure monitoring Data Inversion low permeability reservoir non-Darcy percolation law | |
CN105893679B (en) | Well Test Data Analysis Method is corrected in low yield horizontal well afterflow | |
CN109488282B (en) | Method for determining lower limit of movable storage capacity physical property | |
CN114818533B (en) | Shale oil and gas reservoir fracture parameter determination method and device based on drainage and production data | |
CN108959767A (en) | A kind of narrow river channel type gas reservoir difference well type condensate injury method for numerical simulation | |
CN105822298A (en) | Method for acquiring absolute open flow of shale gas layer based on gas productivity index | |
CN110321575A (en) | Gas condensate reservoir Dynamic Productivity Calculation prediction technique | |
CN104500016A (en) | Novel method for analyzing reservoir physical properties by using pressure drop section data in shale gas reservoir fracturing construction | |
CN107169684B (en) | Development dynamic calculation method under constant liquid volume production condition of multilayer commingled production reservoir | |
CN116006132A (en) | Water intrusion calculation method taking influence of condensation and water evolution in gas into consideration in water-flooding gas reservoir | |
CN110173249A (en) | A kind of straight well separate stratum fracfturing seam based on interlayer equilibrium displacement is long to determine method | |
CN108153944A (en) | A kind of stress sensitive reservoir water-producing gas well natural gas single well controlled reserves determine method | |
CN108132204A (en) | Unconventional reservoir air-water phase percolation curve computational methods in a kind of varied stress field | |
CN108664677B (en) | Oil and gas well production data analysis method | |
CN110424947B (en) | Method for measuring size of micro annular gap of first interface of well cementation | |
CN107247816B (en) | Method for judging farthest construction radius of cold recovery chemical viscosity reduction for heavy oil reservoir | |
CN111963149A (en) | Post-fracturing formation pressure calculation method considering stagnant ground fluid volume pressurization | |
CN111287741B (en) | Rapid calculation method for permeability of compact oil reservoir volume fracturing transformation area |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180608 |