CN106324680A - Stratum rupture pressure prediction method - Google Patents
Stratum rupture pressure prediction method Download PDFInfo
- Publication number
- CN106324680A CN106324680A CN201610682924.3A CN201610682924A CN106324680A CN 106324680 A CN106324680 A CN 106324680A CN 201610682924 A CN201610682924 A CN 201610682924A CN 106324680 A CN106324680 A CN 106324680A
- Authority
- CN
- China
- Prior art keywords
- stratum
- pressure
- individual well
- rock
- formation
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 108
- 239000011435 rock Substances 0.000 claims abstract description 58
- 239000011148 porous material Substances 0.000 claims abstract description 37
- 238000012937 correction Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 12
- 239000012530 fluid Substances 0.000 claims description 4
- 239000000758 substrate Substances 0.000 claims description 2
- 238000005755 formation reaction Methods 0.000 description 84
- 238000005553 drilling Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 5
- 238000001615 p wave Methods 0.000 description 5
- 238000005056 compaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000008398 formation water Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 241000219000 Populus Species 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013211 curve analysis Methods 0.000 description 1
- 235000021185 dessert Nutrition 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention provides a stratum rupture pressure prediction method. The stratum rupture pressure prediction method comprises steps of obtaining a transverse wave curve of a stratum at a single well through earthquake between single wells in a work area; obtaining an elastic parameter of the single well according to the transverse wave curve; obtaining overlaying rock pressure, the stratum formation pore pressure and rock tension resistance strength of the stratum at the single well; obtaining a tectonic stress coefficient of the single well through a fracture micro-earthquake; obtaining the stratum rupture pressure of the stratum at the single well according to the tectonic stress coefficient, the overlaying rock stratum pressure, the stratum pore pressure and the rock tension resistance strength; performing pre-stack elastic parameter inversion on the work area to obtain the elastic parameter of the work area by using the elastic parameter of the stratum of the single well as a constraint condition; and fitting the stratum rupture pressure of the stratum at the single well and the elastic parameter of the work area to obtain the stratum rupture pressure of the work area. The stratum rupture pressure prediction method can improve the prediction accuracy of the stratum rupture pressure.
Description
Technical field
The application relates to oil-gas field development technical field, especially relates to a kind of formation fracture pressure prediction method.
Background technology
The stratum exposed in well, can make formation fracture, this pressure limit value claim bearing when pressure reaches a certain limit
For formation fracture pressure.Formation fracture pressure is one of underlying parameter of drilling well and fracturing design, predicts exactly
Formation fracture pressure is for preventing drilling failure, reducing drilling cost and guarantee that pressure break success is significant.
In conventional reservoir petroleum drilling engineering, can be that balance pressure reduction drilling well carries at drilling phase prediction formation fracture pressure
For strata pressure foundation, in order to reasonably select wellbore mud density, blowout prevention leakproof, reduce the generation of drilling engineering accident, carry
High drilling efficiency, shortening drilling period, reduction drilling cost.In unconventionaloil pool is explored, unconventionaloil pool maximum production
It is decided by geology dessert and engineering factor, when only both reach optimum optimization, guarantee shale gas maximum production.And ground
Layer fracture pressure is the important parameter in engineering factor, and prediction formation fracture pressure is conducive to reasonably optimizing horizontal well pressure exactly
Split conceptual design, promote unconventional yield.
Research shows, formation fracture pressure is with rock mechanics parameters, formation pore pressure, overburden pressure, level should
Power is relevant, additionally, shear wave information plays most important effect in formation fracture pressure prediction.Formation fracture pressure the most both at home and abroad
Power prediction is mainly based upon well logging, carries out choosing suitable model calculating formation fracture by directly or indirectly asking for shear wave information
Pressure.And tradition logging method accident shear wave, and usually utilize empirical equation or model to carry out indirect predictions, the most this
The precision of the shear wave information that mode obtains is the highest, thus the final prediction easily affecting formation fracture pressure, the most easily affect ground
The precision of layer fracture pressure.
Summary of the invention
The purpose of the embodiment of the present application is to provide a kind of formation fracture pressure prediction method, to improve formation fracture pressure
The precision of prediction.
For reaching above-mentioned purpose, the embodiment of the present application provides a kind of formation fracture pressure prediction method, including following step
Rapid:
The shear wave curve on stratum at described individual well is obtained by the crosshole seismic of individual well in work area；
According to the elastic parameter on stratum at individual well described in described shear wave curve acquisition；
Obtain the overburden pressure on stratum, formation pore pressure and tensile strength of rock at described individual well；
Obtaining the pressure break curve on stratum at described individual well by pressure break microseism, obtaining the structure on stratum at described individual well should
Force coefficient；
According to described tectonic stress coefficient, described overburden pressure, described formation pore pressure and described rock anti-tensile
Intensity, obtains the formation fracture pressure on stratum at described individual well；
At described individual well, the elastic parameter on stratum is as constraints, and described work area is carried out prestack elastic parameter inversion,
Obtain the elastic parameter in described work area；
The formation fracture pressure on stratum and the elastic parameter in described work area at individual well described in matching, it is thus achieved that the ground in described work area
Layer fracture pressure.
The formation fracture pressure prediction method of the embodiment of the present application, described by work area individual well crosshole seismic obtain institute
State the shear wave curve on stratum at individual well, including:
Obtain vertical seismic profiling (VSP) VSP curve and the interval transit time curve of the crosshole seismic of individual well in work area, and according to institute
State interval transit time curve described in VSP curvature correction；
According to the interval transit time curve after correction, obtain the shear wave curve on stratum at described individual well.
The formation fracture pressure prediction method of the embodiment of the present application, described according to individual well described in described shear wave curve acquisition at
The elastic parameter on stratum, including:
Described elastic parameter includes: p-wave impedance, S-wave impedance, stratum Poisson's ratio and stratum Young's modulus.
The formation fracture pressure prediction method of the embodiment of the present application, described by ground at the pressure break microseism described individual well of acquisition
The pressure break curve of layer, obtains the tectonic stress coefficient on stratum at described individual well, including:
The pressure break curve on stratum at described individual well is obtained by pressure break microseism, and according to described pressure break curve and described
At individual well, the actual measurement stratum Poisson's ratio on stratum, actual measurement overburden pressure, actual measurement formation pore pressure and actual measurement rock anti-tensile are strong
Degree, obtains the tectonic stress coefficient on stratum at described individual well.
The formation fracture pressure prediction method of the embodiment of the present application, described according to described tectonic stress coefficient, described overlying
Rock pressure, described formation pore pressure and described tensile strength of rock, obtain the formation fracture pressure on stratum at described individual well,
Realize especially by below equation:
Pf=η1*Pp+η2*(Po-Pp)+St
In formula, η is tectonic stress coefficient, PfFor formation fracture pressure, PpFor formation pore pressure, PoFor overlying rock pressure
Power, StFor tensile strength of rock.
The formation fracture pressure prediction method of the embodiment of the present application, the overlying rock pressure on stratum at the described individual well of described acquisition
Power, including:
According to described elastic parameter and formula Po=H [(1-φ) ρma+φρf], at the described individual well of acquisition, stratum is upper
Rock stratum pressure；
In formula, PoFor overlying rock pressure, H is stratum depth of stratum, and φ is overlying rock average pore；ρmaIt is upper overlying strata
The average density of layer substrate；ρfIt it is the average fluid density in overlying rock hole.
The formation fracture pressure prediction method of the embodiment of the present application, the formation pore pressure on stratum at the described individual well of described acquisition
Power, including:
The formation pore pressure on stratum at described individual well is obtained according to Fillippone formula.
The formation fracture pressure prediction method of the embodiment of the present application, the formation pore pressure on stratum at the described individual well of described acquisition
Power, including:
The formation pore pressure on stratum at described individual well is obtained according to Eaton formula.
The formation fracture pressure prediction method of the embodiment of the present application, the acquisition Rock Anti on stratum at the described individual well of described acquisition
Zhang Qiangdu, including:
According to formula St=(0.0045E+0.0035EVsh)/12, the acquisition rock anti-tensile obtaining stratum at described individual well is strong
Degree；
In formula, StFor tensile strength of rock, E is stratum Young's modulus, VshFor shale content.
The formation fracture pressure prediction method of the embodiment of the present application, according to the interval transit time curve after correction, obtains described
The shear wave curve on stratum at individual well, including:
According to the petrophysical model on stratum at the interval transit time curve after correction and described individual well, it is thus achieved that at described individual well
The shear wave curve on stratum.
The embodiment of the present application obtains the shear wave curve on stratum at individual well based on work area crosshole seismic, secondly according to shear wave curve
Obtain the elastic parameter on stratum at described individual well；Then obtain the overburden pressure on stratum at individual well, formation pore pressure and
Tensile strength of rock；Then the tectonic stress coefficient on stratum at individual well is obtained by pressure break microseism；Secondly according to tectonic stress
Coefficient, overburden pressure, formation pore pressure and described tensile strength of rock, obtain the formation fracture pressure on stratum at individual well
Power；Then as constraints, work area is carried out prestack elastic parameter inversion, it is thus achieved that work area with the elastic parameter on stratum at individual well
Elastic parameter；The elastic parameter in the formation fracture pressure on stratum and work area at last matching individual well, it is thus achieved that the stratum in described work area
Fracture pressure.For than the shear wave curve preparation method in conventional formation fracture pressure prediction method, due to the application
Embodiment obtains the shear wave curve shear wave information the most accurately on stratum at individual well based on work area crosshole seismic, therefore, by this
The precision of the formation fracture pressure that application embodiment dopes is improved.
Accompanying drawing explanation
Accompanying drawing described herein is used for providing being further appreciated by the embodiment of the present application, constitutes the embodiment of the present application
A part, is not intended that the restriction to the embodiment of the present application.In the accompanying drawings:
Fig. 1 is the formation fracture pressure prediction method flow diagram of the embodiment of the present application.
Fig. 2 is that p-wave impedance and the formation fracture pressure of the application one embodiment crosses matching schematic diagram；
Fig. 3 is that S-wave impedance and the formation fracture pressure of the application one embodiment crosses matching schematic diagram；
Fig. 4 is formation fracture pressure exhibition portion, the work area schematic diagram that the application one embodiment obtains.
Detailed description of the invention
For making the purpose of the embodiment of the present application, technical scheme and advantage clearer, below in conjunction with embodiment and attached
Figure, is described in further details the embodiment of the present application.Here, the schematic description and description of the embodiment of the present application is used for
Explain the embodiment of the present application, but be not intended as the restriction to the embodiment of the present application.
Below in conjunction with the accompanying drawings, the detailed description of the invention of the embodiment of the present application is described in further detail.
With reference to shown in Fig. 1, the embodiment of the present application provides a kind of formation fracture pressure prediction method and comprises the following steps:
Step S101, obtain the shear wave curve on stratum at described individual well by the crosshole seismic of individual well in work area.
It is in the embodiment of the present application, described by the shear wave curve on stratum at the crosshole seismic acquisition individual well of individual well in work area,
Specifically may include that
1) (Vertical Seismic Profile, vertical seismic cuts open, to obtain the VSP of the crosshole seismic of individual well in work area
Face) curve and interval transit time curve；
2), according to interval transit time curve described in described VSP curvature correction, in order to obtain one closer to seismic band
Interval transit time curve, thus follow-up acquisition shear wave information the most accurately.
First timing analyzes well-log information, unified with nature gamma (GR), natural potential (SP), hole diameter (CAL) well logging song
Obvious problematic well section found out by line, calculates sand shale baseline, analyzes log distribution situation, finds out abnormal section, convenient
Followed by correction, it is then determined that correlation model, carry out parameter in conjunction with rock core information and real logging data and arrange and quality prison
Control, utilizes model inversion completion log, the difference before and after contrast primitive curve analysis reconstruct, then carries out parameter adjustment again
Secondary reconstruct, is repeated, until primitive curve and reconstruct curve are at key hole section zero difference.
3), according to the interval transit time curve after correction, the shear wave curve on stratum at described individual well is obtained.Concrete, can root
According to the petrophysical model on stratum at the interval transit time curve after correction and described individual well, it is thus achieved that the shear wave on stratum at described individual well
Curve.At wherein said individual well, the petrophysical model on stratum can be in conjunction with rock core information, real logging data and Greenberg-
Castagna formula obtains, and the equation gives the petrophysical model of multi mineral, saturation and formation water.
Step S102, according to the elastic parameter on stratum at individual well described in described shear wave curve acquisition.
In the embodiment of the present application, described elastic parameter includes: p-wave impedance, S-wave impedance, stratum Poisson's ratio and stratum poplar
Family name's modulus etc..These parameters can corresponding be obtained by conventionally calculation formula, such as:
P-wave impedance Ip: Ip=Vp*ρ；
S-wave impedance Is: Is=Vs*ρ；
Stratum Poisson ratioσ: σ=[(Vp/Vs)2-2]/[2*(Vp/Vs)2-2]；
Young's modulus E:
Shaliness index IGR:
Shale content Vsh:
In above-mentioned formula, VPFor longitudinal wave velocity, VsFor transverse wave speed, ρ is density of earth formations, GR, GRmin、GRmaxIt is respectively stratum
Natural gamma value, sand layers natural gamma value and shale layer natural gamma value, GCUR is Xi Erqi index, and general old stratum is 2,
New stratum is 3.7.
Step S103, obtain the overburden pressure on stratum, formation pore pressure and tensile strength of rock at described individual well.
In the embodiment of the present application, at the described individual well of described acquisition, the overburden pressure on stratum, specifically may include that
According to described elastic parameter and formula Po=H [(1-φ) ρma+φρf], at the described individual well of acquisition, stratum is upper
Rock stratum pressure；
In formula, PoFor overlying rock pressure, H is depth of stratum, and φ is overlying rock average pore；ρmaIt it is overlying rock base
The average density of matter；ρfIt it is the average fluid density in overlying rock hole.
In the application one embodiment, at the described individual well of described acquisition, the overburden pressure on stratum, is specifically as follows root
The formation pore pressure on stratum at described individual well is obtained according to Fillippone formula.Wherein, Fillippone formula is as follows:
In formula: PFillipponeFor the formation pore pressure of prediction, MPa；PoFor overlying formation pressure, MPa；VminFirm for rock
Property close to formation velocity when zero, is similar to fluid velocity when porosity reaches the upper limit, m/s；VmaxConnect for rock porosity
Velocity of longitudinal wave when being bordering on zero, is similar to matrix velocity, m/s；VinstFor formation velocity, for m/s.VmaxAnd VminTypically can lead to
Cross root mean sequare velocity approximation to set up.
In another embodiment of the application, at the described individual well of described acquisition, the overburden pressure on stratum, is specifically as follows
The formation pore pressure on stratum at described individual well is obtained according to Eaton formula.Wherein, Eaton formula is as follows:
Wherein, PoFor overlying rock pressure, σ is stratum Poisson's ratio, VnormalFor normal compaction speed, m/s；VinstFor stratum speed
Degree, as formation velocity Vinst=VnormalTime, stratum normal compaction, effective pressure meets σ=Po-Pw, PwFor hydrostatic pressure, MPa,
PwFormula P can be passed throughw=0.098 ρwH obtains, ρwFor formation water relative density, h is water-column；And n is Eaton parameter, its
Numerical value is relevant with target area, can be obtained by individual well observed pressure.
In the embodiment of the present application, the tensile strength of rock on stratum at the described individual well of described acquisition, specifically include:
According to formula St=(0.0045E+0.0035EVsh)/12, the acquisition rock anti-tensile obtaining stratum at described individual well is strong
Degree；In formula, StFor tensile strength of rock, E is stratum Young's modulus, VshFor shale content.
Step S104, obtain the pressure break curve on stratum at described individual well by pressure break microseism, obtain ground at described individual well
The tectonic stress coefficient of layer.
In the embodiment of the present application, specifically can obtain the pressure break curve on stratum at described individual well by pressure break microseism, from
Pressure break curve can read the formation fracture pressure of each each fracturing section, therefore according to ground at described pressure break curve and described individual well
The actual measurement stratum Poisson's ratio of layer, actual measurement overburden pressure, actual measurement formation pore pressure and actual measurement tensile strength of rock (specifically come
Say, parameter will be surveyed and substitute into equation below), the tectonic stress coefficient obtaining stratum at described individual well (i.e. obtains η1And η2), it is public
Formula is:
Pf=η1*Pp+η2*(Po-Pp)+St
In formula, η is tectonic stress coefficient, PfFor formation fracture pressure, PpFor formation pore pressure, PoFor overlying rock pressure
Power, StFor tensile strength of rock.
Step S105, according to described tectonic stress coefficient, described overburden pressure, described formation pore pressure and described
Tensile strength of rock, obtains the formation fracture pressure on stratum at described individual well.Described according to described tectonic stress coefficient, described on
Rock stratum pressure, described formation pore pressure and described tensile strength of rock, obtain the formation fracture pressure on stratum at described individual well
Power, realizes especially by below equation:
Pf=η1*Pp+η2*(Po-Pp)+St
In formula, η is tectonic stress coefficient, PfFor formation fracture pressure, PpFor formation pore pressure, PoFor overlying rock pressure
Power, StFor tensile strength of rock.
Step S106, with the elastic parameter on stratum at described individual well as constraints, described work area is carried out prestack elastic
Parametric inversion, it is thus achieved that the elastic parameter in described work area.
In the embodiment of the present application, to carry out prestack elastic parameter anti-for available Zuo Pulizi (Zoeppritz) approximate formula
Drill.
The formation fracture pressure on stratum and the elastic parameter in described work area at individual well described in step S107, matching, it is thus achieved that institute
State the formation fracture pressure in work area.In an example gonosome embodiment, p-wave impedance and formation fracture pressure are through the matching that crosses
After fitting result, can as in figure 2 it is shown, wherein, abscissa is p-wave impedance, and unit is (m/s*g/cc), and vertical coordinate is for rupturing
Pressure, unit is Mpa.In an example gonosome embodiment, S-wave impedance and formation fracture pressure plan after the matching that crosses
Close result, can as it is shown on figure 3, wherein, abscissa is S-wave impedance, and unit is (m/s*g/cc), and vertical coordinate is fracture pressure, single
Position is Mpa.Finally can obtain the planar distribution of formation fracture pressure (such as Fig. 4 institute according to the formation fracture pressure curve simulated
Show).Ellipse is two different horizontal well groups in the diagram, and the formation fracture pressure utilizing the embodiment of the present application to dope shows
There is the biggest difference in these two groups horizontal well group formation fracture pressures, H2 well group formation fracture pressure is less than the horizontal well group of YS108, this
It coincide preferably with later stage scene pressing crack construction, it was demonstrated that the reliability of the embodiment of the present application.
The embodiment of the present application obtains the shear wave curve on stratum at individual well based on work area crosshole seismic, secondly according to shear wave curve
Obtain the elastic parameter on stratum at described individual well；Then obtain the overburden pressure on stratum at individual well, formation pore pressure and
Tensile strength of rock；Then the tectonic stress coefficient on stratum at individual well is obtained by pressure break microseism；Secondly according to tectonic stress
Coefficient, overburden pressure, formation pore pressure and described tensile strength of rock, obtain the formation fracture pressure on stratum at individual well
Power；Then as constraints, work area is carried out prestack elastic parameter inversion, it is thus achieved that work area with the elastic parameter on stratum at individual well
Elastic parameter；The elastic parameter in the formation fracture pressure on stratum and work area at last matching individual well, it is thus achieved that the stratum in described work area
Fracture pressure.For than the shear wave curve preparation method in conventional formation fracture pressure prediction method, due to the application
Embodiment obtains the shear wave curve shear wave information the most accurately on stratum at individual well based on work area crosshole seismic, therefore, by this
The precision of the formation fracture pressure that application embodiment dopes is improved.
Particular embodiments described above, has been carried out the most in detail purpose, technical scheme and the beneficial effect of the application
Describe in detail bright, be it should be understood that the specific embodiment that the foregoing is only the embodiment of the present application, be not used to limit this Shen
Protection domain please, all within spirit herein and principle, any modification, equivalent substitution and improvement etc. done, all should wrap
Within being contained in the protection domain of the application.
Claims (10)
1. a formation fracture pressure prediction method, it is characterised in that comprise the following steps:
The shear wave curve on stratum at described individual well is obtained by the crosshole seismic of individual well in work area；
According to the elastic parameter on stratum at individual well described in described shear wave curve acquisition；
Obtain the overburden pressure on stratum, formation pore pressure and tensile strength of rock at described individual well；
The tectonic stress coefficient on stratum at described individual well is obtained by pressure break microseism；
According to described tectonic stress coefficient, described overburden pressure, described formation pore pressure and described tensile strength of rock,
Obtain the formation fracture pressure on stratum at described individual well；
At described individual well, the elastic parameter on stratum is as constraints, and described work area is carried out prestack elastic parameter inversion, it is thus achieved that
The elastic parameter in described work area；
The formation fracture pressure on stratum and the elastic parameter in described work area at individual well described in matching, it is thus achieved that the stratum in described work area is broken
Split pressure.
Formation fracture pressure prediction method the most according to claim 1, it is characterised in that described by individual well in work area
Crosshole seismic obtains the shear wave curve on stratum at described individual well, including:
Obtain vertical seismic profiling (VSP) VSP curve and the interval transit time curve of the crosshole seismic of individual well in work area, and according to described VSP
Interval transit time curve described in curvature correction；
According to the interval transit time curve after correction, obtain the shear wave curve on stratum at described individual well.
Formation fracture pressure prediction method the most according to claim 1, it is characterised in that described elastic parameter includes: vertical
Natural impedance, S-wave impedance, stratum Poisson's ratio and stratum Young's modulus.
Formation fracture pressure prediction method the most according to claim 1, it is characterised in that described obtained by pressure break microseism
Take the tectonic stress coefficient on stratum at described individual well, including:
The pressure break curve on stratum at described individual well is obtained by pressure break microseism, and according to described pressure break curve and described individual well
The actual measurement stratum Poisson's ratio on stratum, place, actual measurement overburden pressure, actual measurement formation pore pressure and actual measurement tensile strength of rock, obtain
Take the tectonic stress coefficient on stratum at described individual well.
Formation fracture pressure prediction method the most according to claim 1, it is characterised in that described according to described tectonic stress
Coefficient, described overburden pressure, described formation pore pressure and described tensile strength of rock, obtain stratum at described individual well
Formation fracture pressure, realizes especially by below equation:
Pf=η1*Pp+η2*(Po-Pp)+St
In formula, η is tectonic stress coefficient, PfFor formation fracture pressure, PpFor formation pore pressure, PoFor overburden pressure, St
For tensile strength of rock.
Formation fracture pressure prediction method the most according to claim 1, it is characterised in that ground at the described individual well of described acquisition
The overburden pressure of layer, including:
According to described elastic parameter and formula Po=H [(1-φ) ρma+φρf], obtain the upper overlying strata on stratum at described individual well
Stressor layer；
In formula, PoFor overlying rock pressure, H is depth of stratum, and φ is overlying rock average pore；ρmaIt is overlying rock substrate
Average density；ρfIt it is the average fluid density in overlying rock hole.
Formation fracture pressure prediction method the most according to claim 1, it is characterised in that ground at the described individual well of described acquisition
The formation pore pressure of layer, including:
The formation pore pressure on stratum at described individual well is obtained according to Fillippone formula.
Formation fracture pressure prediction method the most according to claim 1, it is characterised in that ground at the described individual well of described acquisition
The formation pore pressure of layer, including:
The formation pore pressure on stratum at described individual well is obtained according to Eaton formula.
Formation fracture pressure prediction method the most according to claim 1, it is characterised in that ground at the described individual well of described acquisition
The acquisition tensile strength of rock of layer, including:
According to formula St=(0.0045E+0.0035EVsh)/12, obtain the acquisition tensile strength of rock on stratum at described individual well；
In formula, StFor tensile strength of rock, E is stratum Young's modulus, VshFor shale content.
Formation fracture pressure prediction method the most according to claim 2, it is characterised in that according to correction after sound wave time
Difference curve, obtains the shear wave curve on stratum at described individual well, including:
According to the petrophysical model on stratum at the interval transit time curve after correction and described individual well, it is thus achieved that stratum at described individual well
Shear wave curve.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610682924.3A CN106324680A (en) | 2016-08-18 | 2016-08-18 | Stratum rupture pressure prediction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610682924.3A CN106324680A (en) | 2016-08-18 | 2016-08-18 | Stratum rupture pressure prediction method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106324680A true CN106324680A (en) | 2017-01-11 |
Family
ID=57743826
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610682924.3A Pending CN106324680A (en) | 2016-08-18 | 2016-08-18 | Stratum rupture pressure prediction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106324680A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107817520A (en) * | 2017-09-20 | 2018-03-20 | 中国石油化工股份有限公司 | The pressure coefficient Forecasting Methodology and system of marine facies mud shale stratum |
| CN108205158A (en) * | 2018-01-26 | 2018-06-26 | 成都理工大学 | A kind of formation pore pressure Forecasting Methodology and system based on index constraint |
| CN108798662A (en) * | 2018-06-23 | 2018-11-13 | 东北石油大学 | Coal seam directional well borehole wall Fracturing Pressure Prediction method |
| CN110761764A (en) * | 2019-10-11 | 2020-02-07 | 中国石油天然气股份有限公司 | Liquid carbon dioxide fracturing method |
| CN113777668A (en) * | 2020-06-10 | 2021-12-10 | 中国石油化工股份有限公司 | Geostress calculation method and device for tight gas reservoir of sand-shale interbed |
| CN115618658A (en) * | 2022-12-16 | 2023-01-17 | 新疆石油管理局有限公司 | Stratum fracture pressure dynamic prediction method and device based on elastic-plastic constitutive |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103837893A (en) * | 2014-03-10 | 2014-06-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Acoustic logging curve correction method based on VSP data |
| CN104863577A (en) * | 2015-04-09 | 2015-08-26 | 中国石油大学(北京) | Method for forecasting formation pore pressure by utilizing propagation time of seismic longitudinal waves |
| CN104989360A (en) * | 2015-07-10 | 2015-10-21 | 中国石油天然气股份有限公司 | Characterization method of low-permeability reservoir dynamic fracture |
-
2016
- 2016-08-18 CN CN201610682924.3A patent/CN106324680A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103837893A (en) * | 2014-03-10 | 2014-06-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Acoustic logging curve correction method based on VSP data |
| CN104863577A (en) * | 2015-04-09 | 2015-08-26 | 中国石油大学(北京) | Method for forecasting formation pore pressure by utilizing propagation time of seismic longitudinal waves |
| CN104989360A (en) * | 2015-07-10 | 2015-10-21 | 中国石油天然气股份有限公司 | Characterization method of low-permeability reservoir dynamic fracture |
Non-Patent Citations (5)
| Title |
|---|
| 伍友佳: "《油藏地质学（第二版）》", 29 February 2004, 石油工业出版社 * |
| 刘红岐 等: "《现场测井与工程应用基础》", 30 April 2013, 石油工业出版社 * |
| 夏庆龙 等: "《稠油油田储层精细描述技术》", 31 January 2010, 石油工业出版社 * |
| 孙宁: "《钻井手册（第二版）》", 30 September 2013, 石油工业出版社 * |
| 李敏 等: ""岩石力学参数试验与地层破裂压力预测研究"", 《石油钻采工艺》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107817520A (en) * | 2017-09-20 | 2018-03-20 | 中国石油化工股份有限公司 | The pressure coefficient Forecasting Methodology and system of marine facies mud shale stratum |
| CN107817520B (en) * | 2017-09-20 | 2020-05-05 | 中国石油化工股份有限公司 | Method and system for predicting pressure coefficient of marine facies shale stratum |
| CN108205158A (en) * | 2018-01-26 | 2018-06-26 | 成都理工大学 | A kind of formation pore pressure Forecasting Methodology and system based on index constraint |
| CN108798662A (en) * | 2018-06-23 | 2018-11-13 | 东北石油大学 | Coal seam directional well borehole wall Fracturing Pressure Prediction method |
| CN110761764A (en) * | 2019-10-11 | 2020-02-07 | 中国石油天然气股份有限公司 | Liquid carbon dioxide fracturing method |
| CN110761764B (en) * | 2019-10-11 | 2022-02-01 | 中国石油天然气股份有限公司 | Liquid carbon dioxide fracturing method |
| CN113777668A (en) * | 2020-06-10 | 2021-12-10 | 中国石油化工股份有限公司 | Geostress calculation method and device for tight gas reservoir of sand-shale interbed |
| CN115618658A (en) * | 2022-12-16 | 2023-01-17 | 新疆石油管理局有限公司 | Stratum fracture pressure dynamic prediction method and device based on elastic-plastic constitutive |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106324680A (en) | Stratum rupture pressure prediction method | |
| Häring et al. | Characterisation of the Basel 1 enhanced geothermal system | |
| CN107766662B (en) | Shale gas horizontal well test sectional evaluation method | |
| US10810331B2 (en) | System for predicting induced seismicity potential resulting from injection of fluids in naturally fractured reservoirs | |
| Patterson et al. | Integration of microseismic data, completion data, and production data to characterize fracture geometry in the Permian Basin | |
| Konstantinovskaya et al. | Strike-slip reactivation of a high-angle normal fault induced by increase of reservoir pore pressure: insight from 3D coupled reservoir-geomechanical modeling | |
| Olsen et al. | Quantifying proppant transport for complex fractures in unconventional formations | |
| Kellogg | Geology and petroleum of the Mancos B formation Douglas Creek Arch area Colorado and Utah | |
| MacKay et al. | Integration of outcrop, subsurface, and microseismic interpretation for rock-mass characterization: An example from the Duvernay Formation, Western Canada | |
| Patterson | Characterizing complex fracture geometry through data integration in the Permian basin | |
| Prankada et al. | Analysis of wellbore stability by pore pressure prediction using seismic velocity | |
| Major et al. | Identifying fracture orientation in a mature carbonate platform reservoir | |
| Ren et al. | Quantitative characterization of interlayer fractures in carbonate rocks based on finite element numerical simulation | |
| Weingarten | On the interaction between fluids and earthquakes in both natural and induced seismicity | |
| Povstyanova et al. | Stress Testing with Improved Wireline Formation Tester in Low Permeability Unconventional Formation | |
| Peza et al. | 3-D Integrated Workflow for Understanding the Fracture Interference and Its Impact into the Gas Production of the Woodford Shale | |
| Li et al. | Hydraulic fracturing-induced microseismicity controlled by rock brittleness and natural fractures in Tongren, Guizhou, China | |
| Xu et al. | Geomechanical characteristics and stimulation of Dibei deep tight sandstone reservoirs in the Kuqa DeprEssion of Tarim Basin | |
| Grazulis | Analysis of stress and geomechanical properties in the Niobrara Formation of Wattenberg Field, Colorado, USA | |
| Zeng et al. | A new method for assessing stage-based hydraulic fracturing quality in Wolfcamp formation of Permian Basin | |
| Perbawa et al. | Integration of Seismic Inversion, Pore Pressure Prediction, and TOC Prediction in Preliminary Study of Shale Gas Exploration | |
| Saadatnia et al. | Geomechanical modelling and cap-rock integrity of one of the southwest Iranian giant carbonate oil field | |
| Clemons et al. | Seismic attributes: Exploiting seismic data to understand heterogeneous reservoir performance in the Eagle Ford Shale, south Texas, USA | |
| Gang | An integrated approach to characterize hydraulic fracturing-induced seismicity in shale reservoirs | |
| Ganaeva | The Sakhaling Offshore Field 4D Geomechanical Modeling for Fracturing Planning |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170111 |
|
| RJ01 | Rejection of invention patent application after publication |