CN109657299A - Shale gas reservoir recovery method - Google Patents
Shale gas reservoir recovery method Download PDFInfo
- Publication number
- CN109657299A CN109657299A CN201811449422.1A CN201811449422A CN109657299A CN 109657299 A CN109657299 A CN 109657299A CN 201811449422 A CN201811449422 A CN 201811449422A CN 109657299 A CN109657299 A CN 109657299A
- Authority
- CN
- China
- Prior art keywords
- well
- area
- npv
- production
- npva
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000011084 recovery Methods 0.000 title claims abstract description 10
- 238000004519 manufacturing process Methods 0.000 claims abstract description 48
- 230000001186 cumulative effect Effects 0.000 claims abstract description 18
- 230000008901 benefit Effects 0.000 claims description 24
- 238000012360 testing method Methods 0.000 claims description 15
- 230000015572 biosynthetic process Effects 0.000 claims description 14
- 230000003068 static effect Effects 0.000 claims description 11
- 230000035699 permeability Effects 0.000 claims description 10
- 238000001179 sorption measurement Methods 0.000 claims description 10
- 238000005553 drilling Methods 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 7
- 230000008859 change Effects 0.000 abstract description 2
- 238000005094 computer simulation Methods 0.000 abstract description 2
- 238000011161 development Methods 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 241000669244 Unaspis euonymi Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000000126 substance Substances 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Feedback Control In General (AREA)
- Solid-Sorbent Or Filter-Aiding Compositions (AREA)
Abstract
A kind of shale gas reservoir recovery method, the following steps are included: using staged fracturing horizontal wells shale gas reservoir, gas flows into crack from stratum, flows into horizontal wellbore from crack, and the corresponding relationship of gas output per well, cumulative production and time is calculated using Life cycle Dynamic Simulation Method.Change Area of a well by well control area width xe, specifies the corresponding relationship of different Area of a well Yu individual well cumulative gas production.Individual well cumulative production curve combination net present value (NPV) NPV model, obtains the corresponding relationship curve of different Area of a well Yu productivity effect relationship.According to productivity effect relation curve, best Area of a well is obtained using graphical method.Shale gas reservoir exploitation is carried out according to best Area of a well.
Description
Technical field
The present invention relates to shale gas reservoir development fields, more particularly, to a kind of shale gas reservoir recovery method.
Background technique
Shale gas reservoir is a kind of continuous type gas mineral reserve, and essential characteristic is that distribution area is big, the variation of local geology condition
Not significant, exploitation is usually influenced by regional environment and geological conditions.For this feature, the development of shale gas reservoir
Mode integrally verifies+development technique breakthrough+step-by-step movement implementation.Once development technique makes a breakthrough, has the base of benefit exploitation
This condition, the investment of each well be substantially it is determining, the decision that at this moment operator is faced be how stepping and implementation issue.More
Specifically, the control area of distribution every mouth well is how, so that the reserves on per unit obtain utmostly
On utilization, obtain greatest benefit the problem of.
Therefore, in order to solve many insufficient and defects of the above-mentioned prior art, it is necessary to study a kind of shale gas reservoir exploitation
Method.
Summary of the invention
The present invention is completed in view of at least one above problem, the present invention passes through shale gas well control area and individual well
The corresponding relationship curve of cumulative gas production obtains Area of a well and Production Gain relation curve in conjunction with economic model, thus
Obtain shale gas reservoir recovery method.
Specifically, according to an aspect of the present invention, a kind of shale gas reservoir recovery method is provided, it is characterised in that including following
Step:
A formation fracture dynamic parameter, including Number of Fractures n) are obtained using well testing, productivity testf, fracture length Lf, crack
Flow conductivity Fc, stratum inner region permeability Km1With stratum outskirt permeability Km2；
Stratum static parameter, including original formation pressure P are obtained using well log interpretation, static pressure testi, formation porosityHorizontal well length DfAnd well control area width xeWith length ye；
Gas PVT parameter and youth's Ge Miaoer isothermal adsorption characteristics parameter are obtained using high pressure property and adsorption isotherm experiment
VLAnd PL；
B) the length y of fixed control areae, by changing width xeChange control plane product value, is measured in conjunction with step A)
Parameter determines the corresponding relationship of cumulative gas production and production time under different Area of a well by gas output per well model,
Obtain Area of a well-cumulative gas production relation curve；
C the corresponding relationship of net present value (NPV) NPV and production time under different Area of a well) are obtained by economic model,
In
Wherein, GP, jFor the annual output in jth year, irFor Annual Percentage Rate, j is jth year, and n is the production cycle, and FC is fixed total throwing
Money, CwellFor individual well drilling cost, CfractureFor the main seam pressure break cost of single cluster, nfFor Number of Fractures；
D Area of a well and net present value (NPV) NPV relation curve and Area of a well and average list) are further obtained
Plane accumulates benefit NPVa relation curve；
NPVa=NPV/S
Determine Area of a well and increment benefit D relation curve,
D=dNPV/dS
It determines optimum efficiency point S2, makes to meet
dNPV/dS|S=s2=NPV/S |S=s2Or NPV/S |S=s2=Max
As S=S2, the obtained benefit of Area of a well is maximum, i.e. NPVa=K；
E work area area A) is measured, Area of a well is set as S2, sets N=INT (A/S2), in the area A of work area
N mouthfuls of wells of even setting and Area of a well are S2, shale gas reservoir exploitation are carried out, to the residual area S in the area A of work areaIt is remainingIt gives over to
It exploits in the future, SIt is remaining≥0。
Compared with prior art, the beneficial effects of the present invention are:
The present invention obtains the pass corresponding with the production time net present value (NPV) NPV under different Area of a well by economic model
System, and then obtain Area of a well and net present value (NPV) NPV relation curve and Area of a well and per unit benefit
NPVa relation curve determines the i.e. best Area of a well of optimum efficiency point S2, to provide foundation for shale gas reservoir exploitation.
Detailed description of the invention
Fig. 1 is that a kind of best Area of a well of the shale gas of preferred embodiment determines flow chart according to the present invention.
Fig. 2 is the staged fracturing horizontal well Production development under a kind of Area of a well of preferred embodiment according to the present invention
Physical model.
Fig. 3 is that cumulative production-time under a kind of different Area of a well of preferred embodiment according to the present invention is corresponding
Relational graph.
Area of a well and cumulative production pair under the different production times that Fig. 4 is a kind of preferred embodiment according to the present invention
Answer relational graph.
Fig. 5 is a kind of corresponding single hop crack overall cost of different pressure break length of preferred embodiment according to the present invention.
Fig. 6 is net present value (NPV)-time corresponding pass under a kind of different Area of a well of preferred embodiment according to the present invention
System's figure.
Area of a well is corresponding with net present value (NPV) under the different production times that Fig. 7 is a kind of preferred embodiment according to the present invention
Relational graph.
Fig. 8 is a kind of Area of a well of preferred embodiment and productivity effect graph of relation according to the present invention.
Fig. 9 is a kind of analysis chart of preferred embodiment order well Optimal Control area really according to the present invention.
Figure 10 is that a kind of well spacing of preferred embodiment and Area of a well, well spacing density conversion indicate according to the present invention
Example.
Figure 11 is a kind of typical well underlying parameter datatable example of preferred embodiment according to the present invention.
Figure 12 is a kind of schematic diagram of the opportunity cost of preferred embodiment according to the present invention.
Specific embodiment
With reference to the accompanying drawing, preferred forms of the invention, specific implementation here are described by preferred embodiment
Mode is to explain the present invention in detail, should not be construed as limiting the invention, and is not departing from spirit and reality of the invention
In the case where matter range, various changes and modifications can be made, these should all be included within protection scope of the present invention.
Embodiment 1
Referring to attached drawing 1-12, it is preferable that the present invention provides one kind to provide a kind of shale gas reservoir recovery method, feature
Be the following steps are included:
A formation fracture dynamic parameter, including Number of Fractures n) are obtained using well testing, productivity testf, fracture length Lf, crack
Flow conductivity Fc, stratum inner region permeability Km1With stratum outskirt permeability Km2；
Stratum static parameter, including original formation pressure P are obtained using well log interpretation, static pressure testi, formation porosityHorizontal well length DfAnd well control area width xeWith length ye；
Gas PVT parameter and youth's Ge Miaoer isothermal adsorption characteristics parameter are obtained using high pressure property and adsorption isotherm experiment
VLAnd PL；
B) the length y of fixed control areae, by changing width xeChange control plane product value, is measured in conjunction with step A)
Parameter determines the corresponding relationship of cumulative gas production and production time under different Area of a well by gas output per well model,
Obtain Area of a well-cumulative gas production relation curve；
C the corresponding relationship of net present value (NPV) NPV and production time under different Area of a well) are obtained by economic model,
In
Wherein, GP, jFor the annual output in jth year, irFor Annual Percentage Rate, j is jth year, and n is the production cycle, and FC is fixed total throwing
Money, CwellFor individual well drilling cost, CfractureFor the main seam pressure break cost of single cluster, nfFor Number of Fractures；
D Area of a well and net present value (NPV) NPV relation curve and Area of a well and average list) are further obtained
Plane accumulates benefit NPVa relation curve；
NPVa=NPV/S
Determine Area of a well and increment benefit D relation curve,
D=dNPV/dS
It determines optimum efficiency point S2, makes to meet
dNPV/dS|S=s2=NPV/S |S=s2Or NPV/S |S=s2=Max
As S=S2, the obtained benefit of Area of a well is maximum, i.e. NPVa=K；
E work area area A) is measured, Area of a well is set as S2, sets N=INT (A/S2), in the area A of work area
N mouthfuls of wells of even setting and Area of a well are S2, shale gas reservoir exploitation are carried out, to the residual area S in the area A of work areaIt is remainingIt gives over to
It exploits in the future, SIt is remaining≥0。
Preferably, N=INT (A/S2) is set specifically, being rounded to (A/S2), such as A=100, S2=0.45, (A/S2)
=222.22, then N=INT (A/S2)=222.That is, being uniformly arranged 222 mouthfuls of lists within the scope of 100 square metres of work area
The gas well that well control area is 0.45 square metre, SIt is remaining=0.1 square metre, then remaining 0.1 quadrature is continued to employ.
Preferably, referring to attached drawing 1-12, it is preferable that the present invention also provides one kind to provide a kind of shale gas reservoir exploitation side
Method, it is characterised in that the following steps are included:
A formation fracture dynamic parameter, including Number of Fractures n) are obtained using well testing, productivity testf, fracture length Lf, crack
Flow conductivity Fc, stratum inner region permeability Km1With stratum outskirt permeability Km2；
Stratum static parameter, including original formation pressure P are obtained using well log interpretation, static pressure testi, formation porosityHorizontal well length DfAnd well control area width xeWith length ye；
Gas PVT parameter and youth's Ge Miaoer isothermal adsorption characteristics parameter are obtained using high pressure property and adsorption isotherm experiment
VLAnd PL；
B) the length y of fixed control areae, by changing width xeChange control plane product value, is measured in conjunction with step A)
Parameter determines the corresponding relationship of cumulative gas production and production time under different Area of a well by gas output per well model,
Obtain Area of a well-cumulative gas production relation curve；
C the corresponding relationship of net present value (NPV) NPV and production time under different Area of a well) are obtained by economic model,
In
Wherein, GP, jFor the annual output in jth year, irFor Annual Percentage Rate, j is jth year, and n is the production cycle, and FC is fixed total throwing
Money, CwellFor individual well drilling cost, CfractureFor the main seam pressure break cost of single cluster, nfFor Number of Fractures；
D Area of a well and net present value (NPV) NPV relation curve and Area of a well and average list) are further obtained
Plane accumulates benefit NPVa relation curve；
NPVa=NPV/S
Determine Area of a well and increment benefit D relation curve,
D=dNPV/dS
It determines optimum efficiency point S2, makes to meet
dNPV/dS|S=s2=NPV/S |S=s2Or NPV/S |S=s2=Max
As S=S2, the obtained benefit of Area of a well is maximum, i.e. NPVa=K；When S is below or above S2, machine is defined
Meeting cost P,
P=K-NPVa
E work area area A) is measured, when A/S2 is integer, Area of a well takes S2；
When A/S2 is non-integer, being rounded to A/S2 is N, determines and is evenly arranged corresponding individual well when N mouthfuls of wells in the area A of work area
Corresponding Area of a well S2 "=A/ (N+1) when N+1 mouthfuls of wells is evenly arranged in control area S2 '=A/N, work area area A,
Determine that S takes S2 ' and S2 respectively according to NPVa=NPV/S ' when corresponding NPVa ' and NPVa ", according to P=K-NPVa
Determine the P ' and P " when NPVa takes NPVa ' and NPVa " respectively,
When P ' is less than or equal to P ", Area of a well S2 takes S2 ', and otherwise Area of a well S2 takes S2 " F) setting list
Well control area is S2, carries out shale gas reservoir exploitation.
Preferably, control area is well control Area length and width xeProduct.
Preferably, gas output per well model specifically: according to physical model shown in Fig. 2, it is raw to construct Life cycle gas well
Dynamic mathematical models are produced, the parameter of step A) measurement/underlying parameter was substituted into mathematical model f (), gas well accumulation can be obtained
The corresponding relationship of yield G and time t.
Gp(t)=f (Km1, Km2, Lf, Fc, nf；xe, ye；φm, Pi；PVT parameter；VL, PL)
Advantageously, the present invention obtains net present value (NPV) NPV and production time under different Area of a well by economic model
Corresponding relationship, and then obtain Area of a well and net present value (NPV) NPV relation curve and Area of a well and mean unit face
Product benefit NPVa relation curve, determines the i.e. best Area of a well of optimum efficiency point S2, to mention for shale gas reservoir exploitation
For foundation.
Embodiment 2
Referring to Fig. 1-12, it is preferable that the present invention provides a kind of shale gas reservoir recovery methods, with Sichuan Province China basin shale
Illustrate implementation process for the typical development well of gas producing region mouthful.Specifically, method includes the following steps:
It in step A), is tested and is tested using mining site and tested, obtained the well static geologic parameter accordingly, utilize pressure break
Construction and dynamic monitoring, obtain the corresponding dynamic project parameter of the well, specifically include:
A1 dynamic data) is obtained using well testing, productivity test, relevant stratum, crack dynamic parameter are obtained by analysis,
Including Number of Fractures nf, fracture length Lf, fracture condudtiviy Fc, stratum inner region permeability Km1With stratum outskirt permeability Km2；
A2 well log interpretation, static pressure test) are utilized, stratum static parameter, including original formation pressure P are obtained by analysisi、
Formation porosityHorizontal well length DfAnd the width x of well control areaeWith length ye。
A3 gas PVT parameter, youth's Ge Miaoer isothermal adsorption characteristics parameter) are obtained using high pressure property and adsorption isotherm experiment
VLAnd PL。
The above sound state Data Summary is in Figure 11.
In step B), fixed control area length, by changing control area width xeTo change control plane product value.It is excellent
Selection of land, ratio is set as constant I between fracture length and area width herex=Lf/xe.Control area is set in the present embodiment
xe=100~450m, corresponding fracture length Lf=80~360m, the underlying parameter provided in conjunction with step A) substitute into individual well and produce gas
Amount model is calculated.Fig. 3 reflects the corresponding relationship of cumulative production and production time under different Area of a well, individual well control
Area processed is bigger, and cumulative production is higher.Fig. 4 reflects the pass corresponding with cumulative production of Area of a well under the different production times
System, the production cycle is longer, and Area of a well influences cumulative production more significant.
In step C), economic model is introduced, wherein net present value (NPV) formula is such as:
Here, GP, jFor the annual output in jth year, irFor Annual Percentage Rate, j is jth year, and n is the production cycle, and FC is fixed total throwing
It provides (equipment cost etc.), CwellFor individual well drilling cost, CfractureFor single cluster (main seam) pressure break cost, nfFor Number of Fractures.
In step C1), individual well overall cost (in addition to pressure break) is 40,000,000 yuan, Annual Percentage Rate 10%, gas price 1.2
Member/m3, the period is 20 years.The different exponential increases of pressure break length cost, such as Fig. 5.
In step C2), it is only existing that individual well is obtained according to the individual well cumulative production curve combination formula (1) provided in step B)
Value.Fig. 6 shows the corresponding relationship of net present value (NPV) and production time under different Area of a well, shorter for the production time, single
In the smaller situation of well control area, individual well net present value (NPV) will appear negative value.Fig. 7 shows Area of a well and list under different cycles
Well net present value (NPV) relationship.
In step D), best Area of a well is determined.
In step D1), Area of a well and productivity effect relation curve (curve 1 in Fig. 8) are obtained.
In step D2), Area of a well and per unit benefit (NPVa) relation curve are determined, i.e. scale is imitated
Answer (curve 2 in Fig. 8).
NPVa=NPV/S (2)
In step D3), Area of a well and increment benefit (D) relation curve are determined, i.e. edge effect is (bent in Fig. 8
Line 3).
D=dNPV/dS (3)
In step D4), breakeven point S1 (Fig. 9) is determined, be to meet
NPV|S=s1=0 or NPVa |S=s1=0 (4)
This is Area of a well lowest limit, and individual well Production Gain will not support investment when being lower than this value.
In step D5), determines optimum efficiency point S2 (Fig. 9), make to meet
dNPV/dS|S=s2=NPV/S |S=s2Or NPV/S |S=s2=Max (5)
S2 is determined as best Area of a well.
E Area of a well) is set as S2, carries out shale gas reservoir exploitation.
Preferably, as S=S2, the obtained benefit of Area of a well is maximum, i.e. NPVa=K, when S is below or above S2
When, actual benefit is below this value.Referring to Figure 12, opportunity cost P is determined, define
P=K-NPVa (6)
Then P can measure the size of opportunity cost.Its meaning is the probable lost revenue when S value deviates optimal value.Its
Chance is, if having an opportunity to supply insufficient area or being used to the area for being more than beat new well, can obtain more effects
Benefit.The difference of this better benefit and actual benefit is exactly probable opportunity loss.For continuity oil reservoir, such chance
It is existing for reality, it is probable just to become inevitable.It, can according to P value to the sensibility and P value of S in the asymmetry of S2 point two sides
Effectively to instruct decision process.Preferably, work area area A is measured, when A/S2 is integer, Area of a well is taken
S2；
When A/S2 is non-integer, being rounded to A/S2 is N, determines and is evenly arranged corresponding individual well when N mouthfuls of wells in the area A of work area
Corresponding Area of a well S2 "=A/ (N+1) when N+1 mouthfuls of wells is evenly arranged in control area S2 '=A/N, work area area A,
Determine that S takes S2 ' and S2 respectively according to NPVa=NPV/S ' when corresponding NPVa ' and NPVa ", according to P=K-NPVa
Determine the P ' and P " when NPVa takes NPVa ' and NPVa " respectively,
When P ' is less than or equal to P ", Area of a well S2 takes S2 ', and otherwise Area of a well S2 takes S2 ".
Embodiment 3
Preferably, the present invention provides a kind of shale gas reservoir recovery methods, method includes the following steps:
Step A), according to fig. 2 in physical model, using staged fracturing horizontal wells shale gas reservoir, gas is from stratum
It flows into crack, flow into horizontal wellbore from crack.Gas output per well, cumulative production are calculated using Life cycle Dynamic Simulation Method
With the corresponding relationship of time.
Step B), the present invention changes Area of a well by well control area width xe, specifies different single well-controlled faces
The long-pending corresponding relationship with individual well cumulative gas production.Well control area width (i.e. corresponding well away from) and Area of a well, well spacing density
Transformational relation is summarised in Figure 10.
Step C), individual well cumulative production curve combination net present value (NPV) NPV model, obtains different Area of a well and production is imitated
The corresponding relationship curve of beneficial relationship.The corresponding pressure break cost of difference Area of a well is different in net present value model.
Step D), according to productivity effect relation curve, best Area of a well and/or corresponding is obtained using graphical method
Best techniques parameter.
In conclusion the beneficial effects of the present invention are:
The present invention obtains the pass corresponding with the production time net present value (NPV) NPV under different Area of a well by economic model
System, and then obtain Area of a well and net present value (NPV) NPV relation curve and Area of a well and per unit benefit
NPVa relation curve determines the i.e. best Area of a well of optimum efficiency point S2, to provide foundation for shale gas reservoir exploitation.
The present invention is not limited to above-mentioned specific embodiments.It is understood that not departing from spirit and substance of the present invention model
In the case where enclosing, various changes and modifications can be made, these should all be included within protection scope of the present invention.
Claims (2)
1. a kind of shale gas reservoir recovery method, it is characterised in that the following steps are included:
A formation fracture dynamic parameter, including Number of Fractures n) are obtained using well testing, productivity testf, fracture length Lf, fracture guide
Ability Fc, stratum inner region permeability Km1With stratum outskirt permeability Km2；
Stratum static parameter, including original formation pressure P are obtained using well log interpretation, static pressure testi, formation porosityWater
Horizontal well length DfAnd well control area width xeWith length ye；
Gas PVT parameter and youth's Ge Miaoer isothermal adsorption characteristics parameter V are obtained using high pressure property and adsorption isotherm experimentLWith
PL；
B) the length y of fixed control areae, by changing width xeChange control plane product value, the parameter measured in conjunction with step A)
By gas output per well model, the corresponding relationship of cumulative gas production and production time under different Area of a well is determined, obtain
Area of a well-cumulative gas production relation curve；
C the corresponding relationship of net present value (NPV) NPV and production time under different Area of a well) are obtained by economic model, wherein
Wherein, GP, jFor the annual output in jth year, irFor Annual Percentage Rate, j is jth year, and n is the production cycle, and FC is fixed gross investment,
CwellFor individual well drilling cost, CfractureFor the main seam pressure break cost of single cluster, nfFor Number of Fractures；
D Area of a well and net present value (NPV) NPV relation curve and Area of a well and mean unit face) are further obtained
Product benefit NPVa relation curve；
NPVa=NPV/S
Determine Area of a well and increment benefit D relation curve,
D=dNPV/dS
It determines optimum efficiency point S2, makes to meet
dNPV/dS|S=s2=NPV/S |S=s2Or NPV/S |S=s2=Max
As S=S2, the obtained benefit of Area of a well is maximum, i.e. NPVa=K；
E work area area A) is measured, Area of a well is set as S2, sets N=INT (A/S2), uniformly set in the area A of work area
It sets N mouthfuls of wells and Area of a well is S2, shale gas reservoir exploitation is carried out, to the residual area S in the area A of work areaIt is remainingIt gives over in the future
Exploitation, SIt is remaining≥0。
2. according to the method described in claim 1, it is characterized by: Area of a well is equal to well control Area length and width
Product.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811449422.1A CN109657299B (en) | 2018-11-28 | 2018-11-28 | Shale gas reservoir mining method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811449422.1A CN109657299B (en) | 2018-11-28 | 2018-11-28 | Shale gas reservoir mining method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109657299A true CN109657299A (en) | 2019-04-19 |
| CN109657299B CN109657299B (en) | 2023-01-17 |
Family
ID=66111953
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811449422.1A Active CN109657299B (en) | 2018-11-28 | 2018-11-28 | Shale gas reservoir mining method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109657299B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112241801A (en) * | 2019-07-18 | 2021-01-19 | 中国石油天然气股份有限公司 | Method and device for determining optimal development well spacing of shale gas |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102606129A (en) * | 2012-04-10 | 2012-07-25 | 中国海洋石油总公司 | Method and system for thin interbed oilfield development |
| EP2550422A1 (en) * | 2010-03-25 | 2013-01-30 | Bruce A. Tunget | Pressure controlled well construction and operation systems and methods usable for hydrocarbon operations, storage and solution mining |
| CN103161434A (en) * | 2013-04-07 | 2013-06-19 | 赵万福 | Mining method for low permeability reservoir of shale gas and the like |
| CN104948163A (en) * | 2014-03-24 | 2015-09-30 | 中国石油化工股份有限公司 | Method for measuring shale gas well capacity |
| CN106150477A (en) * | 2015-04-23 | 2016-11-23 | 中国石油化工股份有限公司 | A kind of method determining single well controlled reserves |
-
2018
- 2018-11-28 CN CN201811449422.1A patent/CN109657299B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2550422A1 (en) * | 2010-03-25 | 2013-01-30 | Bruce A. Tunget | Pressure controlled well construction and operation systems and methods usable for hydrocarbon operations, storage and solution mining |
| CN102606129A (en) * | 2012-04-10 | 2012-07-25 | 中国海洋石油总公司 | Method and system for thin interbed oilfield development |
| CN103161434A (en) * | 2013-04-07 | 2013-06-19 | 赵万福 | Mining method for low permeability reservoir of shale gas and the like |
| CN104948163A (en) * | 2014-03-24 | 2015-09-30 | 中国石油化工股份有限公司 | Method for measuring shale gas well capacity |
| CN106150477A (en) * | 2015-04-23 | 2016-11-23 | 中国石油化工股份有限公司 | A kind of method determining single well controlled reserves |
Non-Patent Citations (2)
| Title |
|---|
| WILSON K C,ET AL.: "Optimization of shale gas field development using direct search techniques and reducedphysics models", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 * |
| 郭小哲等: "页岩气压裂水平井组布井参数研究", 《非常规油气》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112241801A (en) * | 2019-07-18 | 2021-01-19 | 中国石油天然气股份有限公司 | Method and device for determining optimal development well spacing of shale gas |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109657299B (en) | 2023-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109441422B (en) | Shale gas well spacing optimization mining method | |
| CN106869911B (en) | Evaluation method for describing compressibility of shale reservoir | |
| CN102041995B (en) | System for monitoring complicated oil deposit flooding conditions | |
| CN109209333B (en) | Shale gas multi-well group efficient mining interval optimization method | |
| CN110704888A (en) | Unconventional oil and gas reservoir encrypted well volume fracturing construction parameter optimization design method | |
| CN103573263B (en) | The horizontal well stream of a kind of tight sandstone reservoir fracturing production is installed with meter method | |
| CN104695950A (en) | Prediction method for volcanic rock oil reservoir productivity | |
| CN109236265B (en) | Method for optimizing tight gas reservoir well pattern | |
| CN103742130A (en) | High permeability band characterization method based on time varying | |
| CN110952976B (en) | Single-well exploitation stable yield potential evaluation method in gas reservoir development mode | |
| CN106408208B (en) | Volume fracturing correctional effect evaluation method | |
| CN111927417A (en) | Shale gas staged fracturing horizontal well group reserve utilization condition evaluation method | |
| Ailin et al. | Technical measures of deliverability enhancement for mature gas fields: A case study of Carboniferous reservoirs in Wubaiti gas field, eastern Sichuan Basin, SW China | |
| CN108952657A (en) | A kind of horizontal well platform fracturing fracture length determining method and device | |
| CN109577967A (en) | Recover of Tight Sandstone Gas Reservoir recovery method | |
| CN114427432A (en) | Method for determining development potential of residual gas in gas reservoir | |
| CN109657299A (en) | Shale gas reservoir recovery method | |
| CN111950111B (en) | Dynamic analysis method suitable for open-bottom carbonate reservoir | |
| CN107704646B (en) | Modeling method after compact reservoir volume modification | |
| CN108875115B (en) | Method for determining rock strength | |
| CN109339745A (en) | Shale gas reservoir recovery method based on optimal Area of a well | |
| CN115860266A (en) | Shale gas/coal bed gas well productivity evaluation method and system and electronic equipment | |
| Olsen et al. | Waterflood performance evaluation in a chalk reservoir with an ensemble of tools | |
| CN106845685B (en) | Parameter optimization method for reducing horizontal well oil testing fracturing operation cost | |
| CN109630103A (en) | Recover of Tight Sandstone Gas Reservoir recovery method based on enrichment region |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Jia Ailin Inventor after: Guo Jianlin Inventor after: Yun Sheng Wei Inventor after: Jia Chengye Inventor after: Wang Junlei Inventor before: Guo Jianlin Inventor before: Jia Ailin Inventor before: Yun Sheng Wei Inventor before: Jia Chengye Inventor before: Wang Junlei |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |