CN107575219A - A kind of shale gas reservoir formation fracture pressure gradient computational methods - Google Patents
A kind of shale gas reservoir formation fracture pressure gradient computational methods Download PDFInfo
- Publication number
- CN107575219A CN107575219A CN201710840885.XA CN201710840885A CN107575219A CN 107575219 A CN107575219 A CN 107575219A CN 201710840885 A CN201710840885 A CN 201710840885A CN 107575219 A CN107575219 A CN 107575219A
- Authority
- CN
- China
- Prior art keywords
- fracture pressure
- formation
- density
- pressure gradient
- dimensions
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The present invention relates to shale gas reservoir formation fracture pressure gradient computational methods, collect the Logging Curves data of full well section, calculate matrix density ρ ma, determine formation average porosity Φ and fluid density ρ w；Collect shale gas well Reservoir Section depth, vertical depth, sound wave AC, density DEN, neutron CNL borehole log datas；Formation average denstiy is calculated by ρ=(1 Φ) * ρ ma+ Φ * ρ w；Fracture pressure COEFFICIENT K is determined by being fitted the method returned：K=P0+P1*AC+P2*DEN+P3*CNL, shale reservoir formation fracture pressure gradient is calculated according to formula F RAC=K* ρ, then use formula：FP=FRAC*H calculates fracture pressure, last output result.The present invention is closer in the multiple wells of Fuling shale gas field applications pages rock gas well 200, the formation fracture pressure gradient that the formation fracture pressure gradient of calculating obtains with practice of construction, and error is less than 15%.
Description
Technical field
The present invention relates to a kind of method for being calculated using Conventional Logs and determining shale reservoir formation fracture pressure, tool
Body is related to a kind of shale gas reservoir formation fracture pressure gradient computational methods.
Background technology
Shale gas is a kind of new cleaning fuel, and is a kind of important unconventional gas resource.Shale gas is developed
Need with " well factory " pattern drilling level well, need progress multistage massive hydraulic fracture to change to horizontal well horizontal segment shale gas reservoir
Make, cost is high, and risk is big.The precision that shale gas reservoir formation fracture gradient and formation fracture pressure parameter are calculated or predicted is very
It is important, directly affect reservoir reconstruction effect.
Traditional formation fracture pressure gradient calculation method mainly has method of Eaton, Ma Xiusi and Kai Lifa, Christman
Method.These methods are mainly formed in the sixties in last century to during the seventies, and solution is that sandstone, carbonate rock etc. are normal
Advise reservoir formation fracture gradient and formation fracture pressure parameter calculates or forecasting problem.
With the discovery in domestic shale gas field and the expansion of shale gas exploration and development experimental scale, shale gas horizontal well multistage
Pressure break cost is high, risk is big, and the formation fracture pressure gradient, formation fracture pressure computational accuracy to shale gas reservoir require more next
Higher, the problem of traditional formation fracture pressure gradient calculation method is difficult in adapt to production requirement, and calculation error is big, is more and more prominent
Go out.
The content of the invention
The purpose of the present invention is to be directed to the above-mentioned state of the art, it is desirable to provide one kind can utilize Logging Curves data meter
The method for calculating reliable shale gas reservoir formation fracture pressure.
The implementation of the object of the invention is shale gas reservoir formation fracture pressure gradient computational methods, to concretely comprise the following steps：
1) the Logging Curves data of full well section are collected, Logging Curves data includes natural gamma, neutron, close
Degree, sound wave, Resistivity log；Matrix density ρ ma are calculated using Logging Curves data using weighted mean method,
Using full well section Logging Curves data and the core porosity analysis of data of region prospect pit core hole, the average hole in stratum is determined
Porosity Φ and fluid density ρ w；
Matrix density ρ ma dimensions are g/cm3,
Formation average porosity Φ dimensions are percentage,
Fluid density ρ w dimensions are g/cm3；
2) shale gas well Reservoir Section depth, vertical depth, sound wave AC, density DEN, neutron CNL borehole log datas are collected；
Density DEN dimensions are g/cm3,
Sound wave AC dimensions are μ s/m,
Neutron CNL dimensions are percentage,
3) it is average close to calculate stratum by formation rock skeletal density ρ ma, porosity Φ and fluid density ρ w obtained by step 1)
Spend ρ：
ρ=(1- Φ) * ρ ma+ Φ * ρ w；
Formation average denstiy ρ dimensions are g/cm3,
4) it is collected that shale reservoir section sound wave, density DEN, neutron well logging curve data have been pressed with reference to area by step 2)
Well fracture pressure is split, fracture pressure COEFFICIENT K is determined by being fitted the method returned：K=P0+P1*AC+P2*DEN+P3*CNL,
P0 to P3 is regional experience coefficient in formula, and relevant with fracturing section formation lithology, formation lithology is relevant with log；
5) formula is passed through：FRAC=K* ρ calculate shale gas Reservoir Section formation fracture pressure gradient FRAC；
Shale reservoir formation fracture pressure gradient FPAC dimensions are MPa/100m in formula,
Formation average denstiy ρ dimensions are g/cm3,
6) result of calculation is exported.
The present invention in the multiple wells of Fuling shale gas field applications pages rock gas well 200, the formation fracture pressure gradient of calculating with
The formation fracture pressure gradient that practice of construction obtains is closer to, and error is less than 15%.
Brief description of the drawings
Fig. 1 is workflow block diagram of the present invention.
Embodiment
Reference picture 1, of the invention concretely comprises the following steps：
1) the Logging Curves data of full well section are collected, Logging Curves data includes natural gamma, neutron, close
The logs such as degree, sound wave, resistivity；Matrix density ρ is calculated using density log curve data using weighted mean method
Ma, determine that stratum is averaged using the core porosity analysis of data of full well section Logging Curves data and region prospect pit core hole
Porosity Φ and fluid density ρ w；
Matrix density ρ ma dimensions are g/cm3,
Formation average porosity Φ dimensions are percentage,
Fluid density ρ w dimensions are g/cm3。
Wherein Logging Curves data is exclusion bad well and corrupt data data.Fluid density ρ w typically take 1.1
±0.1g/cm3。
2) shale gas well Reservoir Section depth, vertical depth, sound wave AC, density DEN, neutron CNL borehole log datas are collected；
Density DEN dimensions are g/cm3,
Sound wave AC dimensions are μ s/m,
Neutron CNL dimensions are percentage.
3) it is average close to calculate stratum by formation rock skeletal density ρ ma, porosity Φ and fluid density ρ w obtained by step 1)
Spend ρ：
ρ=(1- Φ) * ρ ma+ Φ * ρ w；
Formation average denstiy ρ dimensions are g/cm3,
4) it is collected that shale reservoir section sound wave, density DEN, neutron well logging curve data have been pressed with reference to area by step 2)
Well fracture pressure is split, fracture pressure COEFFICIENT K is determined by being fitted the method returned：K=P0+P1*AC+P2*DEN+P3*CNL,
Fracture pressure coefficient is relevant with fracturing section formation lithology, and formation lithology is relevant with log, by with fractured well
Data combination data fitting method can determine regional empirical coefficient P0-P3, so as to calculate formation fracture pressure coefficient.
5) formula is passed through：FRAC=K* ρ calculate shale gas Reservoir Section formation fracture pressure gradient FRAC；
Shale reservoir formation fracture pressure gradient FPAC dimensions are MPa/100m in formula,
Formation average denstiy ρ dimensions are g/cm3,
6) result of calculation is exported.
The present invention is described in detail with specific embodiment below.(embodiment preferably by above the step of write)
Example one：Certain gas field JY2-3 wells
1) the full well section Logging Curves data of the well (including sound wave, density, neutron curve) are collected, using weighted average
It is 2.67g/cm that method calculates the average skeletal density ρ ma of rock using density log curve data3, it is bent using full well section conventional logging
The core porosity analysis of data of line number evidence and region prospect pit core hole determines that formation average porosity Φ is 0.05, fluid density
ρ w are 1.05g/cm3；
2) it is that 2980.0-4490.0m vertical depths are 2448.0-2521.0m to determine shale gas reservoir well section, and collects shale
Reservoir Section sound wave, density, neutron well logging curve data；
3) ρ=(1- Φ) * ρ ma+ Φ * ρ w are passed through；It is 2.589g/cm to calculate formation average denstiy ρ3；
4) determine that fracture pressure coefficient k (has with regional stratum by fracturing section sound wave, density, neutron well logging curve data
The regional experience coefficient of pass), (AC is sound wave to fracture pressure coefficient formulas, and DEN is density, and CNL is neutron, and P0-P3 is normal
Number)：K=P0+P1*AC+P2*DEN+P3*CNL.By being fitted to the section pressure break data of 260 multiple wells of Fuling Prefecture more than 4000
Analysis, p1=-0.004274867, p2=-0.046441009, p3=0.004972802, p0 are drawn using fitting homing method
=2.019452937.It is 0.977 that the average K values of the fracturing section, which are calculated,.
5) formula is passed through：It is 2.537 that FRAC=K* ρ, which calculate fracturing section formation fracture pressure gradient, mean burst pressure
63.06MPa；
6) fracture pressure gradient and fracture pressure result of calculation are exported according to user's request.
Calculating JY2-3 well 2980.0-4490.0m well section shale gas-bearing formation formation fracture pressures using the present invention is
63.06MPa.Completion pressing crack construction, it is 63.48MPa to survey the average formation fracture pressure of 20 sections of the well, and error is less than 10%, two
Person's result of calculation is very close, open-flow capacity at initial stage day producing natural gas about 116.39 × 104m3, fracturing transformation effect is notable.
From example 1, formation fracture pressure gradient, formation fracture pressure and the measured result that the present invention calculates approach,
Error is less than 10%, meets site operation needs, has stronger application value.
Example two：Certain shale gas field JY13-2 well
JY13-2 wells are a bite long horizontal sections horizontal wells, and gas reservoir horizontal segment is 2630-4078m, and gas reservoir vertical depth is 2330.0-
2412.0m, horizontal segment pass through in gas reservoir.
1) the full well section Logging Curves data of the well (including sound wave, density, neutron curve) are collected, using weighted average
It is 2.67g/cm that method calculates the average skeletal density ρ ma of rock using density log curve data3, it is bent using full well section conventional logging
The core porosity analysis of data of line number evidence and region prospect pit core hole determines that formation average porosity Φ is 0.05, fluid density
ρ w are 1.05g/cm3；
2) it is 2630-4078m, vertical depth 2332.0-2412.0m to determine shale gas reservoir well section, and collects shale storage
Interval sound wave, density, neutron well logging curve data；
3) ρ=(1- Φ) * ρ ma+ Φ * ρ w are passed through；It is 2.589g/cm to calculate formation average denstiy ρ3；
4) determine that fracture pressure coefficient k (has with regional stratum by fracturing section sound wave, density, neutron well logging curve data
The regional experience coefficient of pass), (AC is sound wave to fracture pressure coefficient formulas, and DEN is density, and CNL is neutron, and P0-P3 is normal
Number)：K=P0+P1*AC+P2*DEN+P3*CNL, by being fitted to the section pressure break data of 260 multiple wells of Fuling Prefecture more than 4000
Analysis, p1=-0.004274867, p2=-0.046441009, p3=0.004972802, p0 are drawn using fitting homing method
=2.019452937, it is 1.036 that the average K values of the fracturing section, which are calculated,；
5) formula is passed through：It is 2.137 that FRAC=K* ρ, which calculate fracturing section formation fracture pressure gradient, mean burst pressure
64.2MPa；
6) fracture pressure gradient and fracture pressure result of calculation are exported according to user's request.
It is using the JY13-2 well 2630-4078m horizontal segment shale gas-bearing formations stratum mean burst pressure of the invention calculated
64.2MPa.Completion fracturing work, point 18 sections of constructions, construction actual measurement stratum mean burst pressure is 60.0MPa, and initial stage is without choked flow
Measure day producing natural gas about 111.02 × 104m3, fracturing transformation effect is obvious.
From example 1, formation fracture pressure gradient, formation fracture pressure and the measured result that the present invention calculates approach,
Error is less than 10%, instructs pressing crack construction to work well.
Claims (3)
1. shale gas reservoir formation fracture pressure gradient computational methods, it is characterised in that：Concretely comprise the following steps：
1) the Logging Curves data of full well section are collected, Logging Curves data includes natural gamma, neutron, density, sound
Ripple, Resistivity log；Matrix density ρ ma are calculated using Logging Curves data using weighted mean method, utilized
The core porosity analysis of data of full well section Logging Curves data and region prospect pit core hole, determines formation average porosity
Φ and fluid density ρ w；
Matrix density ρ ma dimensions are g/cm3,
Formation average porosity Φ dimensions are percentage,
Fluid density ρ w dimensions are g/cm3；
2) shale gas well Reservoir Section depth, vertical depth, sound wave AC, density DEN, neutron CNL borehole log datas are collected；
Density DEN dimensions are g/cm3,
Sound wave AC dimensions are μ s/m,
Neutron CNL dimensions are percentage,
3) formation rock skeletal density ρ ma, porosity Φ and fluid density ρ w calculate formation average denstiy ρ obtained by step 1)：
ρ=(1- Φ) * ρ ma+ Φ * ρ w,
Formation average denstiy ρ dimensions are g/cm3；
4) it is collected that shale reservoir section sound wave, density DEN, neutron well logging curve data combine area fractured well by step 2)
Fracture pressure, fracture pressure COEFFICIENT K is determined by being fitted the method returned：K=P0+P1*AC+P2*DEN+P3*CNL,
P0 to P3 is regional experience coefficient in formula, and relevant with fracturing section formation lithology, formation lithology is relevant with log；
5) formula is passed through：FRAC=K* ρ calculate shale gas Reservoir Section formation fracture pressure gradient FRAC；
Shale reservoir formation fracture pressure gradient FPAC dimensions are MPa/100m in formula,
Formation average denstiy ρ dimensions are g/cm3；
6) result of calculation is exported.
2. shale gas reservoir formation fracture pressure gradient computational methods according to claim 1, it is characterised in that：Step 1)
Middle Logging Curves data is exclusion bad well and corrupt data data.
3. shale gas reservoir formation fracture pressure gradient computational methods according to claim 1, it is characterised in that：Fluid is close
Degree ρ w take 1.1 ± 0.1g/cm3。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710840885.XA CN107575219B (en) | 2017-09-15 | 2017-09-15 | Shale gas reservoir stratum fracture pressure gradient calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710840885.XA CN107575219B (en) | 2017-09-15 | 2017-09-15 | Shale gas reservoir stratum fracture pressure gradient calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107575219A true CN107575219A (en) | 2018-01-12 |
| CN107575219B CN107575219B (en) | 2020-08-07 |
Family
ID=61036319
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710840885.XA Active CN107575219B (en) | 2017-09-15 | 2017-09-15 | Shale gas reservoir stratum fracture pressure gradient calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107575219B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108708715A (en) * | 2018-05-02 | 2018-10-26 | 中石化石油工程技术服务有限公司 | The survey logging method of fast prediction shale gas-bearing formation formation fracture pressure gradient |
| CN108801538A (en) * | 2018-06-15 | 2018-11-13 | 中国石油大学(华东) | A method of calculating decline of pressure gradient |
| CN110410069A (en) * | 2019-08-15 | 2019-11-05 | 中石化石油工程技术服务有限公司 | Shale gas horizontal well formation fracture pressure gradient detecting well prediction technique |
| CN110469321A (en) * | 2019-08-05 | 2019-11-19 | 中石化江汉石油工程有限公司 | A kind of logging method of determining stratum fracture pressure gradient |
| CN113550740A (en) * | 2020-04-24 | 2021-10-26 | 中国石油化工股份有限公司 | Method for realizing continuous calculation of longitudinal pressure gradient of single well |
| CN115773106A (en) * | 2021-09-08 | 2023-03-10 | 中国石油化工股份有限公司 | Method for identifying evolution process of overpressure sealing box by using paleo-pressure attenuation gradient |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5282384A (en) * | 1992-10-05 | 1994-02-01 | Baroid Technology, Inc. | Method for calculating sedimentary rock pore pressure |
| CN1588128A (en) * | 2004-10-22 | 2005-03-02 | 石油大学(北京) | Method for predicting collapse pressure and bursting pressure for borehole to be drilled section using earthquake record |
| CN101025084A (en) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | Method for predetecting formation pore pressure under drill-bit while drilling |
| CN102733789A (en) * | 2012-07-06 | 2012-10-17 | 崔彦立 | Staged fracturing construction yield increment method for waterpower in deep thickened oil deposit thick-bedded sandstone storage layer |
| CN104573344A (en) * | 2014-12-25 | 2015-04-29 | 中国海洋石油总公司 | Method for acquiring gas content of shale reservoir through well logging data |
| US20150176387A1 (en) * | 2013-12-20 | 2015-06-25 | Schlumberger Technology Corporation | Perforation strategy |
| CN105257252A (en) * | 2015-06-08 | 2016-01-20 | 中国石油集团川庆钻探工程有限公司 | Method for optimizing shale gas horizontal well clustering perforation well section by using logging information |
| CN106285646A (en) * | 2016-09-09 | 2017-01-04 | 中国石油大学(华东) | Drilling well loss horizon recognition methods based on Multi-information acquisition |
-
2017
- 2017-09-15 CN CN201710840885.XA patent/CN107575219B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5282384A (en) * | 1992-10-05 | 1994-02-01 | Baroid Technology, Inc. | Method for calculating sedimentary rock pore pressure |
| CN1588128A (en) * | 2004-10-22 | 2005-03-02 | 石油大学(北京) | Method for predicting collapse pressure and bursting pressure for borehole to be drilled section using earthquake record |
| CN101025084A (en) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | Method for predetecting formation pore pressure under drill-bit while drilling |
| CN102733789A (en) * | 2012-07-06 | 2012-10-17 | 崔彦立 | Staged fracturing construction yield increment method for waterpower in deep thickened oil deposit thick-bedded sandstone storage layer |
| US20150176387A1 (en) * | 2013-12-20 | 2015-06-25 | Schlumberger Technology Corporation | Perforation strategy |
| CN104573344A (en) * | 2014-12-25 | 2015-04-29 | 中国海洋石油总公司 | Method for acquiring gas content of shale reservoir through well logging data |
| CN105257252A (en) * | 2015-06-08 | 2016-01-20 | 中国石油集团川庆钻探工程有限公司 | Method for optimizing shale gas horizontal well clustering perforation well section by using logging information |
| CN106285646A (en) * | 2016-09-09 | 2017-01-04 | 中国石油大学(华东) | Drilling well loss horizon recognition methods based on Multi-information acquisition |
Non-Patent Citations (1)
| Title |
|---|
| 华北石油管理局干部处: "《华北油田保护油气层技术》", 31 August 1996 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108708715A (en) * | 2018-05-02 | 2018-10-26 | 中石化石油工程技术服务有限公司 | The survey logging method of fast prediction shale gas-bearing formation formation fracture pressure gradient |
| CN108801538A (en) * | 2018-06-15 | 2018-11-13 | 中国石油大学(华东) | A method of calculating decline of pressure gradient |
| CN108801538B (en) * | 2018-06-15 | 2020-10-30 | 中国石油大学(华东) | A method for calculating the pressure decay gradient |
| CN110469321A (en) * | 2019-08-05 | 2019-11-19 | 中石化江汉石油工程有限公司 | A kind of logging method of determining stratum fracture pressure gradient |
| CN110410069A (en) * | 2019-08-15 | 2019-11-05 | 中石化石油工程技术服务有限公司 | Shale gas horizontal well formation fracture pressure gradient detecting well prediction technique |
| CN110410069B (en) * | 2019-08-15 | 2022-08-19 | 中石化石油工程技术服务有限公司 | Shale gas horizontal well stratum fracture pressure gradient logging prediction method |
| CN113550740A (en) * | 2020-04-24 | 2021-10-26 | 中国石油化工股份有限公司 | Method for realizing continuous calculation of longitudinal pressure gradient of single well |
| CN115773106A (en) * | 2021-09-08 | 2023-03-10 | 中国石油化工股份有限公司 | Method for identifying evolution process of overpressure sealing box by using paleo-pressure attenuation gradient |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107575219B (en) | 2020-08-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107575219A (en) | A kind of shale gas reservoir formation fracture pressure gradient computational methods | |
| CN109594968B (en) | Shale gas multistage fracturing horizontal well post-fracturing fracture parameter evaluation method and system | |
| CN103912269B (en) | Method for determining formation fracture pressure gradient logging of shale gas reservoir | |
| Dong et al. | Resource evaluation for shale gas reservoirs | |
| CN107577831B (en) | Method for calculating scale of karst cave of fracture-cavity carbonate oil-gas reservoir | |
| Parsegov et al. | Physics-driven optimization of drained rock volume for multistage fracturing: field example from the Wolfcamp Formation, Midland Basin | |
| CN106869911B (en) | Evaluation method for describing compressibility of shale reservoir | |
| CN111219175B (en) | Optimization method for fracture matching of fractured carbonate rock acid fracturing considering stress sensitivity | |
| CN111927417B (en) | Shale gas sectional fracturing horizontal well group reserve utilization condition evaluation method | |
| Bachman et al. | Examining the effects of stress dependent reservoir permeability on stimulated horizontal Montney gas wells | |
| Zhang et al. | An evaluation method of volume fracturing effects for vertical wells in low permeability reservoirs | |
| CN105893679B (en) | Well Test Data Analysis Method is corrected in low yield horizontal well afterflow | |
| CN108894777A (en) | A kind of separate stratum fracfturing commingling production oil-gas reservoir reservoir properties and characteristic of crack determination method for parameter | |
| CN104695950A (en) | Prediction method for volcanic rock oil reservoir productivity | |
| CN116128084A (en) | Prediction method for volume fracture network control reserves of tight oil reservoir horizontal well | |
| CN110469321A (en) | A kind of logging method of determining stratum fracture pressure gradient | |
| CN114444005B (en) | A method for calculating formation pressure of heterogeneous oil reservoir in a mine | |
| CN113792932B (en) | Shale gas yield prediction method utilizing microseism-damage-seepage relation | |
| CN111090127A (en) | Carbonate rock geothermal reservoir productivity prediction method and system | |
| CN114997083A (en) | Method for calculating chart of natural gas reserves of abnormally high-pressure water-bearing gas reservoir | |
| CN106227995A (en) | Shale gas reservoir single well productivity computational methods | |
| Sun | Characterization and simulation of discrete fracture networks in unconventional shale reservoirs | |
| KR101818098B1 (en) | Method for estimating volume of clay in rocks | |
| CN105221143A (en) | A kind of coal seam permeability computational methods based on mining data | |
| CN112267872B (en) | Method for determining volume of gas well drilling karst cave system by using production dynamic data |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |