CN112464586A - Shale gas well control reserve calculation method - Google Patents
Shale gas well control reserve calculation method Download PDFInfo
- Publication number
- CN112464586A CN112464586A CN202011302455.0A CN202011302455A CN112464586A CN 112464586 A CN112464586 A CN 112464586A CN 202011302455 A CN202011302455 A CN 202011302455A CN 112464586 A CN112464586 A CN 112464586A
- Authority
- CN
- China
- Prior art keywords
- pressure
- well control
- well
- mpa
- reserve
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 43
- 238000004519 manufacturing process Methods 0.000 claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 21
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 18
- 239000000463 material Substances 0.000 claims abstract description 12
- 239000000126 substance Substances 0.000 claims abstract description 4
- 230000006835 compression Effects 0.000 claims description 15
- 238000007906 compression Methods 0.000 claims description 15
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 4
- 239000011435 rock Substances 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 238000013178 mathematical model Methods 0.000 abstract description 2
- 230000035699 permeability Effects 0.000 abstract 1
- 238000011084 recovery Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000000611 regression analysis Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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/08—Thermal analysis or thermal optimisation
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a shale gas well control reserve calculation method, which aims at the characteristics of low porosity and low permeability of a shale gas reservoir, and has the defects that the well closing pressure recovery speed is slow when the formation pressure is obtained, the production is influenced, and the existing shale gas reservoir single well reserve calculation method is difficult to apply. The idea of adopting a material balance equation, firstly assuming single well control reserves and calculating the formation pressure by using single well production data is provided, and the formation pressure can be obtained under the condition of not closing the well. And then, a flowing substance balance equation is adopted, and the control reserves of the single well are obtained through regression by utilizing the stratum pressure and the single well production data obtained through calculation. And finally, comparing the assumed storage quantity value with the calculated storage quantity value, and if the assumed storage quantity value meets the error limit, obtaining the shale gas reservoir single-well control storage quantity and the formation pressure value. In addition, the invention also establishes an automatic calculation mathematical model, converts the trial calculation process into automatic calculation, reduces the workload, reduces the operation difficulty and improves the working efficiency of the reserve calculation.
Description
Technical Field
The invention relates to the field of shale gas well control reserve rating, in particular to a shale gas well control reserve calculation method.
Background
At present, shale gas reservoir reserves commonly used at home and abroad are calculated by a material balance method, an unstable well testing method and the like, and the methods all need an important parameter of average formation pressure. However, the shut-in pressure of the shale gas reservoir which is a low-porosity and low-permeability gas reservoir is recovered slowly, and long-time shut-in is needed to obtain the parameters. This results in very limited application of these methods, which brings about inconvenience for field application, and thus the prior art has yet to be improved and enhanced.
Disclosure of Invention
In view of the defects of the prior art, the invention provides a novel shale gas well reserve control method, which can calculate the single-well control reserve by using production dynamic data without closing a well and measuring pressure, and improves the practicability of the method.
In order to achieve the purpose, the invention adopts the following technical scheme:
the invention provides a shale gas well control reserve calculation method, which comprises the following steps:
s1, acquiring basic data of a single well, wherein the basic data comprises static data: irreducible water saturation, langmuir volume, langmuir pressure, porosity, shale density, gas reservoir original temperature, gas reservoir original pressure; phase state data: pressure and compression factor, pressure and gas viscosity; production data: daily gas production, accumulated gas production and bottom hole flow pressure.
S2, assuming single well to control the reserves, and calculating the average formation pressure by using a material balance equation;
the material balance equation is:
in the formula: p is the formation pressure, MPa; z*To correct the bias factor; z is a deviation factor; gpFor cumulative production, 108m3(ii) a G is single well controlled reserve, 108m3;swcIrreducible water saturation, f; rhobIs the density of rock, g/cm3;VLIs the Langmuir volume, m3/t;pLLane pressure, MPa; phi is the free porosity, f; t is the formation temperature, K; p is a radical ofscPressure under standard conditions, MPa; t isscTemperature under standard conditions, K; zscIs the deviation factor under standard conditions, i is the initial mark.
S3, fitting regression to calculate the single well control reserve according to the average formation pressure calculated in the step S2 by using a flowing substance balance equation and gas well production data;
the flow material balance equation is:
in the formula: μ is the gas viscosity, mPas; m (p) is the pseudo-pressure, MPa2/mpa·s;tca *Simulating time for material balance, d; c. CgIs a free gas compression coefficient, MPa-1;caTo adsorb a gas compression coefficient, MPa-1;ct *Is a comprehensive compression coefficient, MPa-1;qgIs the daily gas production, m3/d;kpssThe slope of the fit for equation (3); m (p)wf) Simulating pressure, MPa, for bottom hole flowing pressure2/mpa·s。
S4, establishing a relative error model according to the function relation between the single-well control reserves assumed in the step S2 and the single-well control reserves calculated in the step S3;
the relative error model is:
lb≤x≤ub (10)
in the formula: lb is the lower limit of the single well control reserve, the default is the current cumulative yield, 108m3(ii) a ub is the upper limit of the single well controlled reserve, and the default is 10 × 108m3(ii) a x is the single well control reserve assumed in step S2, g (x) is the single well control reserve calculated in step S3; f (x) represents a relative error.
And S5, calculating a reasonable single-well control reserve G by adopting a binary search method for the relative error model in the step S4.
The binary search method calculation of the relative error model comprises the following steps:
(3) Judging (ub-lb) orIf it is less than the limit value, and if any value is less than the limit value, the calculation is ended,and (3) the reserve G is reasonably controlled by the single well, otherwise, the step (1) is skipped until the calculation result meets the precision requirement.
Compared with the prior art, the shale gas well control reserve calculation method provided by the invention does not need to carry out well shut-in pressure measurement, and can obtain the single-well control reserve and the average formation pressure only by analyzing the existing production data of the gas well under the condition of not influencing the production plan of the gas well; meanwhile, the method for automatically calculating the control reserves greatly reduces the workload of reserve calculation and realizes the automation of reserve calculation.
Drawings
FIG. 1 is a flow chart of shale gas well reserves calculation provided by the present invention.
FIG. 2 is a diagram of the results of an automatically computed mathematical model provided by the present invention.
FIG. 3 is a flow chart of the automated calculation provided by the present invention.
FIG. 4 is a production dynamic diagram for an X1 well in an example of the present invention.
FIG. 5 is a calculated fit graph of shale gas well reserves of an X1 well in an embodiment of the invention.
Detailed Description
The invention provides a shale gas well control reserve calculation method, and in order to make the purpose, technical scheme and effect of the invention clearer and clearer, the invention is further described in detail below by referring to the attached drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Referring to fig. 1, the shale gas well control reserve calculation method provided by the embodiment of the invention includes the following steps:
s1, acquiring basic data of a single well, wherein the basic data comprises static data: irreducible water saturation, langmuir volume, langmuir pressure, porosity, shale density, gas reservoir original temperature, gas reservoir original pressure; phase state data: pressure and compression factor, pressure and gas viscosity; production data: daily gas production, accumulated gas production and bottom hole flow pressure;
specifically, the method comprises the following steps:
analyzing and obtaining a relational expression between pressure p and a compression factor Z by adopting a regression analysis method and utilizing data between the pressure and the compression factor, wherein Z is f (p);
analyzing a relation between pressure p and gas viscosity mu by using data between the pressure and the gas viscosity by adopting a regression analysis method, wherein mu is g (p);
s2, assuming single well to control the reserves, and calculating the average formation pressure by using a material balance equation;
specifically, the method comprises the following steps:
in the formula: p is the formation pressure, MPa; z*To correct the bias factor; z is a deviation factor; gpFor cumulative production, 108m3(ii) a G is single well controlled reserve, 108m3;swcIrreducible water saturation, f; rhobIs the density of rock, g/cm3;VLIs the Langmuir volume, m3/t;pLLane pressure, MPa; phi is the free porosity, f; t is the formation temperature, K; p is a radical ofscPressure under standard conditions, MPa; t isscTemperature under standard conditions, K; zscIs the deviation factor under standard conditions, i is the initial mark.
S3, fitting regression to calculate the single well control reserve according to the average formation pressure calculated in the step S2 by utilizing a flowing substance balance equation and gas well production data;
specifically, the method comprises the following steps:
Z=f(p) (8)
μ=g(p) (9)
in the formula: μ is the gas viscosity, mPas; m (p) is the pseudo-pressure, MPa2/mpa·s;tca *Simulating time for material balance, d; c. CgIs a free gas compression coefficient, MPa-1;caTo adsorb a gas compression coefficient, MPa-1;ct *Is a comprehensive compression coefficient, MPa-1;qgIs the daily gas production, m3/d;kpssThe slope of the fit for equation (3); m (p)wf) Simulating pressure, MPa, for bottom hole flowing pressure2/mpa·s。
Are respectively provided withIs the y-axis, andfor the x axis, the relation between the two is obtained by regression, and the slope is kpssThen, the formula (3) shows that:
s4, establishing a relative error model according to the function relation between the single-well control reserves assumed in the step S2 and the single-well control reserves calculated in the step S3;
specifically, the method comprises the following steps:
the assumed single-well controlled reserve is defined as x, the single-well controlled reserve calculated from the production data is defined as y, and y is g (x) as shown in equations (1), (3) and (10).
The relative error model is now established as follows:
lb≤x≤ub (12)
in the formula: lb is the lower limit of the single well control reserve, the default is the current cumulative yield, 108m3(ii) a ub is the upper limit of the single well controlled reserve, and the default is 10 × 108m3;
Through a large number of experimental analyses, fig. 2 can be obtained, from which fig. 2 an empirical conclusion can be drawn: f (lb) and f (ub) have opposite signs, and f (x) there must be a zero between lb and ub, which is a reasonable single well control reserve G.
And S5, calculating a reasonable single-well control reserve G by adopting a binary search method for the relative error model in the step S4. The calculation flow is shown in FIG. 3;
specifically, the method comprises the following steps:
(3) Judging (ub-lb) orIf it is less than the limit value, and if any value is less than the limit value, the calculation is ended,and (3) the reserve G is reasonably controlled by the single well, otherwise, the step (1) is skipped until the calculation result meets the precision requirement.
In specific application, collecting static data of a certain shale gas well: irreducible water saturation, langmuir volume, langmuir pressure, porosity, shale density, gas reservoir original temperature, gas reservoir original pressure; phase state data: pressure and compression factor, pressure and gas viscosity; production dynamic data: daily gas production, accumulated gas production and bottom hole flow pressure; the production dynamics data is shown in figure 4.
By adopting the shale gas well reserves calculation method provided by the invention, a shale gas well reserves calculation fitting graph can be obtained, as shown in fig. 5, the graph shows that the method provided by the invention has high fitting degree, and the correlation coefficient reaches 0.99.
By adopting the automatic calculation method provided by the invention, an iterative process of shale gas well reserves calculation can be obtained, and as shown in table 1, the automatic calculation method provided by the invention has high calculation speed, only 13 iterations are needed, the time is consumed for 29.36 seconds, and the relative error can reach 0.003691407% precision.
In conclusion, the shale gas well control reserve calculation method adopted by the embodiment can obtain high precision and high calculation speed, greatly improves the working efficiency, and solves the problem that the production of a gas well is influenced because the well closing back pressure is needed in the reserve calculation process.
It should be understood that equivalents and modifications of the technical solution and inventive concept thereof may occur to those skilled in the art, and all such modifications and alterations should fall within the scope of the appended claims.
Claims (6)
1. A shale gas well control reserve calculation method is characterized by comprising the following steps:
s1, acquiring basic data of a single well;
s2, assuming single well to control the reserves, and calculating the average formation pressure by using a material balance equation;
s3, fitting regression to calculate the single well control reserve according to the average formation pressure calculated in the step S2 by using a flowing substance balance equation and gas well production data;
s4, establishing a relative error model according to the function relation between the single-well control reserves assumed in the step S2 and the single-well control reserves calculated in the step S3;
and S5, calculating reasonable single-well control reserves by adopting a binary search method for the relative error model in the step S4.
2. The shale gas well control reserve calculation method of claim 1, wherein the basic data of the single well in step S1 comprises static data: irreducible water saturation, langmuir volume, langmuir pressure, porosity, shale density, gas reservoir original temperature, gas reservoir original pressure; phase state data: pressure and compression factor, pressure and gas viscosity; production data: daily gas production, accumulated gas production and bottom hole flow pressure.
3. The shale gas well control reserve calculation method of claim 1, wherein the material balance equation in step S2 is:
in the formula: p is the formation pressure, MPa; z*To correct the bias factor; z is a deviation factor; gpFor cumulative production, 108m3(ii) a G is single well controlled reserve, 108m3;swcIrreducible water saturation, f; rhobIs the density of rock, g/cm3;VLIs the Langmuir volume, m3/t;pLLane pressure, MPa; phi is the free porosity, f; t is the formation temperature, K; p is a radical ofscPressure under standard conditions, MPa; t isscTemperature under standard conditions, K; zscIs the deviation factor under standard conditions, i is the initial mark.
4. The shale gas well control reserve calculation method of claim 1, wherein the flowing material balance equation in step S3 is:
in the formula: μ is the gas viscosity, mPas; m (p) is the pseudo-pressure, MPa2/mpa·s;tca *Simulating time for material balance, d; c. CgIs a free gas compression coefficient, MPa-1;caTo adsorb a gas compression coefficient, MPa-1;ct *Is a comprehensive compression coefficient, MPa-1;qgIs the daily gas production, m3/d;kpssThe slope of the fit for equation (3); m (p)wf) Simulating pressure, MPa, for bottom hole flowing pressure2/mpa·s。
5. The shale gas well control reserve calculation method of claim 1, wherein the relative error model in step S4 is:
lb≤x≤ub (10)
in the formula: lb is the lower limit of the single well control reserve, the default is the current cumulative yield, 108m3(ii) a ub is the upper limit of the single well controlled reserve, and the default is 10 × 108m3(ii) a x is the single well control reserve assumed in step S2, g (x) is the single well control reserve calculated in step S3; f (x) represents a relative error.
6. The shale gas well control reserve calculation method of claim 5, wherein the calculation steps of the binary search method for the relative error model are as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011302455.0A CN112464586A (en) | 2020-11-19 | 2020-11-19 | Shale gas well control reserve calculation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011302455.0A CN112464586A (en) | 2020-11-19 | 2020-11-19 | Shale gas well control reserve calculation method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112464586A true CN112464586A (en) | 2021-03-09 |
Family
ID=74837767
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011302455.0A Pending CN112464586A (en) | 2020-11-19 | 2020-11-19 | Shale gas well control reserve calculation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112464586A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113818873A (en) * | 2021-10-20 | 2021-12-21 | 中国石油大学(北京) | Method and device for calculating dynamic reserves of gas reservoir under action of hidden gas-containing area |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104632187A (en) * | 2013-11-14 | 2015-05-20 | 中国石油化工股份有限公司 | Method for determining dynamic reserve volume of water production coal seam gas well |
CN106484933A (en) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | A kind of method and system for determining shale gas well well control dynamic holdup |
CN108612525A (en) * | 2018-04-19 | 2018-10-02 | 重庆科技学院 | A kind of gas reservoir protection Reserve Estimation Method |
CN108804819A (en) * | 2018-06-10 | 2018-11-13 | 西南石油大学 | A kind of low permeability gas reservoirs dynamic holdup evaluation method |
-
2020
- 2020-11-19 CN CN202011302455.0A patent/CN112464586A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104632187A (en) * | 2013-11-14 | 2015-05-20 | 中国石油化工股份有限公司 | Method for determining dynamic reserve volume of water production coal seam gas well |
CN106484933A (en) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | A kind of method and system for determining shale gas well well control dynamic holdup |
CN108612525A (en) * | 2018-04-19 | 2018-10-02 | 重庆科技学院 | A kind of gas reservoir protection Reserve Estimation Method |
CN108804819A (en) * | 2018-06-10 | 2018-11-13 | 西南石油大学 | A kind of low permeability gas reservoirs dynamic holdup evaluation method |
Non-Patent Citations (6)
Title |
---|
刘尧文等: "利用流动物质平衡确定低渗气藏单井控制储量", 《天然气勘探与开发》 * |
刘建仪等: "考虑应力敏感和压裂液影响的页岩气井动态产能评价方法", 《科学技术与工程》 * |
杨龙: "页岩气藏单井动态储量计算方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
王董东等: "页岩气藏可采储量计算新方法", 《大庆石油地质与开发》 * |
蔡振华等: "致密砂岩气藏多段压裂水平气井动态控制储量计算方法", 《长江大学学报(自科版)》 * |
赵坤山: "压降法计算地质储量在低渗气藏中的研究及应用", 《天然气技术与经济》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113818873A (en) * | 2021-10-20 | 2021-12-21 | 中国石油大学(北京) | Method and device for calculating dynamic reserves of gas reservoir under action of hidden gas-containing area |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108612525B (en) | Gas reservoir dynamic reserve calculation method | |
CN109281663B (en) | Method for determining permeability of coal bed gas reservoir | |
CN107506948A (en) | A kind of shale oil and gas composite volume analysis method based on dynamic aerial drainage volume | |
CN106932324B (en) | A kind of method of determining High water cut sandstone oil reservoir reservoir permeability changing rule | |
CN110878690B (en) | Method for determining shale gas well production decreasing curve based on carbon isotope evolution | |
CN111287740A (en) | Method for calculating dynamic reserves of abnormal high-pressure gas reservoir based on real strain | |
CN112464586A (en) | Shale gas well control reserve calculation method | |
CN111502652A (en) | Yield decreasing and production dynamic prediction method for three-hole medium gas reservoir horizontal well | |
CN112554864A (en) | Method for calculating single-well control reserve of water-producing gas well | |
CN114997080A (en) | Method for explaining unstable productivity test data of fracture-cavity gas reservoir | |
CN110295894B (en) | Method for establishing horizontal well productivity prediction model | |
CN113863920A (en) | Method for detecting volume of gas channeling channel | |
CN109710881B (en) | Heterogeneous well region inter-well interference rapid analysis and calculation method | |
CN116738672A (en) | Method for establishing complex biological reef bottom water and gas reservoir numerical simulation model | |
CN116467958A (en) | Salt lake brine numerical model construction and water supplementing and mineral dissolving efficiency calculation method | |
CN114252381B (en) | Method for evaluating pollution degree of fractured reservoir horizontal well drilling fluid after solid-phase pollution | |
CN112989721B (en) | Rapid calculation method for reconstruction volume of compact reservoir volume fracturing horizontal well | |
CN109236274A (en) | A kind of well choosing method for high-water-cut oil-producing well water blockoff | |
CN115060875A (en) | Method for determining hydrate reservoir production pressure interval based on Darcy's law | |
CN111287741B (en) | Rapid calculation method for permeability of compact oil reservoir volume fracturing transformation area | |
CN114658391A (en) | Dynamic pressure regulating method and related device for drainage and production shaft of coal-bed gas well | |
CN111927422A (en) | Shale gas well adsorption output contribution rate calculation method | |
CN115288646B (en) | Connectivity analysis method, device, medium and terminal for fracturing horizontal well | |
CN112464476B (en) | Method for evaluating reservoir residual oil saturation distribution field | |
CN112780266B (en) | New method for calculating inaccessible pore volume |
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 |