CN107843927A - Shale formation pressure prediction method and device based on well shake joint speed - Google Patents
Shale formation pressure prediction method and device based on well shake joint speed Download PDFInfo
- Publication number
- CN107843927A CN107843927A CN201610835388.6A CN201610835388A CN107843927A CN 107843927 A CN107843927 A CN 107843927A CN 201610835388 A CN201610835388 A CN 201610835388A CN 107843927 A CN107843927 A CN 107843927A
- Authority
- CN
- China
- Prior art keywords
- mrow
- pressure
- interval velocity
- msub
- velocity
- 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 52
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 44
- 238000013528 artificial neural network Methods 0.000 claims description 6
- 230000005251 gamma ray Effects 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000005755 formation reaction Methods 0.000 description 31
- 239000007789 gas Substances 0.000 description 21
- 239000010410 layer Substances 0.000 description 16
- 239000011435 rock Substances 0.000 description 7
- 238000011160 research Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- 238000005056 compaction Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 230000008021 deposition Effects 0.000 description 3
- 239000004575 stone Substances 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000003245 coal Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 206010017076 Fracture Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011005 laboratory method Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000001228 spectrum 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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A kind of shale formation pressure prediction method and device, this method based on well shake joint speed comprise the following steps：Obtain seismic interval velocity；Ask for improved seismic interval velocity；Strata pressure is calculated based on improved seismic interval velocity；Calibrating (base measuring) pressure coefficient is asked for, strata pressure is corrected according to calibrating (base measuring) pressure coefficient.The present invention selects inversion of simulated AC curve interval velocity on the basis of Fillippone pressure prediction empirical equations, and so as to obtain the improvement interval velocity that can more reflect actual formation lithology, using improving, the strata pressure that interval velocity calculates is more reasonable.In addition, the present invention is obtained calibrating (base measuring) pressure coefficient by the statistical fit of interval velocity and pressure correcting coefficient, pressure prediction empirical equation is corrected based on actual work area, actual formation so that the pressure of calculating is more suitable for studying work area, more efficient, reasonable.
Description
Technical field
The present invention relates to unconventional shale gas exploration and development field, more particularly to a kind of page based on well shake joint speed
Rock Formation pressure prediction method and device.
Background technology
Shale gas refers to using methane as main component, contain by principal mode of free state and ADSORPTION STATE in high-carbon mud stone or
Natural gas in person's dead color (black) mud stone.Such clean energy, efficiently, range of deposition is wide in basin, thickness is larger, be after
Another the important Unconventional gas found after coal bed gas.
Global huge (the 456x1012m of shale gas stock number3, 2012), equivalent to compact sandstone gas and the summation of coal bed gas,
It is mainly distributed on North America, Latin America, West Europe, Central and Eastern Europe, the former Soviet Union, the SA on the south middle east and north the Sahara, the Central Asia
With China etc..Although the shale in China and the U.S. is all developed very much, differed greatly in terms of geologic feature, be mainly manifested in page
The geological tectonic conditions of rock Distribution Area, shale deposition type, the shale deposition age, thermal evolution condition, shale maturity period with
And maturity, buried depth, preservation condition etc..
After 21 century, with promotion of the China to energy industry and the increase to petroleum resources demand, shale gas research
Paces start to accelerate.Since shale gas scale probing 2009, petrochina, sinopec and shell China are in COMPONENT IN SOUTH SICHUAN BASIN
Impressive progress is achieved in terms of the Paleozoic marine facies source shale gas exploration of underground, thus China also turns into unique one in addition to north America region
Individual beginning business explores the country of shale gas.Particularly Sinopec Group is in succession on the Wei Jin new fields of Sichuan
Trembled with fear under three folded five sections of Xu jiahes of system under sub- section, Sichuan southwest depression Lve Bao Chang-metal and stone pole structural belt military nine old hole groups,
The group such as Chongqing Peng River, Guizhou Province Bei Dingshan constructions, Fuling Chongqing Jiao masonry dam Longma small stream Duan Jun achieves great shale gas and found, promotes
The business processes of China's shale gas exploitation, are particularly obtained huge in the peak group shale gas of river southeast Jiao's masonry dam area Longma small stream-five
Big breakthrough, scale development is obtained, and complete 5,000,000,000 side's annual capacity construction.
With deepening continuously for shale gas exploration, shale gas and corresponding shale are studied to corresponding geophysical techniques not
Disconnected proposition demand, demand is constantly proposed especially for the pressure prediction techniques related to the actual gas production of shale gas well.At present,
The method for carrying out prediction of formation pressure using seismic technology mainly has four kinds：Be compacted balance method, equivalent depth method (Eaton,
Stone etc.), Fillippone methods, improve Fillippone equations.In terms of shale formation pore pressure prediction, at present
Method based on forecast for seismic data strata pressure is also only limitted on the basis of undercompaction formation theory, and its general principle is just
In normal compacted formation, with the increase of depth, stratum is gradually compacted, and the porosity of formation rock is gradually reduced, and seismic wave exists
Spread speed in rock is gradually accelerated, and in abnormal pressure stratum, the change opposite with normal compaction trend is shown as, hole
Porosity is bigger than the porosity of normal compaction, and rock density is lower than the density value of normal compaction, and seismic wave velocity is than normal compaction
Velocity of wave is small.According to the presence of these abnormal can predicted anomaly high pressures, and the size of its pressure can be estimated.With strata pressure
A kind of method that the development of Forecasting Methodology, Downton and Roure (2010) describe AVAZ while elastic inversion, it can be anti-
Drill to obtain the elastic parameter under isotropic medium background, carry out the forecast analysis of shale fragility and minimal closure pressure, and
Can Simultaneous Inversion obtain describe fracture development density and trend information parameter.Suirez-Rivera etc. (2009) passes through thing
Simulated test is managed, combines isotropism and anisotropic pressure computational methods has calculated formation pore clossing pressure.
According to demand, Shale Pressure Predicting Technique is being updated, for a regional or a set of any side in stratum
Method is accurate or perhaps effective, and this is depended primarily on coincide with actual production data, being capable of effective descriptive study
The technology of area's strata pressure instruction parameter is exactly effective pressure prediction techniques.Along with the research of China's shale gas and exploration and development
Grow in intensity, strata pressure is more accurately and effectively predicted in expectation, is specified pressure condition (reservoir pressure coefficient) and is produced with shale gas well
The relation of tolerance.
The content of the invention
The purpose of the present invention is to propose to a kind of Shale Pressure Forecasting Methodology based on well shake joint speed, to improve pressure
(pressure coefficient) precision of prediction, so as to improve the success rate of shale gas well probing.
An aspect of of the present present invention provides a kind of shale formation pressure prediction method based on well shake joint speed, including following
Step：
Obtain seismic interval velocity；
Ask for improved seismic interval velocity；
Strata pressure is calculated based on the improved seismic interval velocity；
Calibrating (base measuring) pressure coefficient is asked for, the strata pressure is corrected according to the calibrating (base measuring) pressure coefficient.
Preferably, seismic interval velocity is calculated by below equation (6)：
Wherein, VR,nAnd VR,n-1Respectively stratum top interface, the stack velocity of bottom interface,WithRespectively stratum is pushed up
Interface, the root mean sequare velocity of bottom interface, t0,nAnd t0,n-1Respectively stratum top interface, the reflection interval of bottom interface, VnFor earthquake layer
Speed.
Preferably, it is described to ask for improved seismic interval velocity and include：
Based on acoustic logging and Natural Gamma-ray Logging Curves reconstruct pseudo-acoustic curve DT；
Based on the pseudo-acoustic curve DT, pass through the seismic interval velocity of BP neural network algorithm computed improved.
Preferably, multiple wells in work area will be studied and correspond to observed pressure value and meter at the interval velocity and the destination layer of destination layer
The ratio calculated between pressure value carries out index return, asks for the calibrating (base measuring) pressure coefficient k '.
Another aspect of the present invention provides a kind of shale formation pressure prediction device based on well shake joint speed, including：
Interval velocity acquisition module, for obtaining seismic interval velocity；
Interval velocity improves module, for asking for improved seismic interval velocity；
Formation pressure calculation module, for calculating strata pressure based on the improved seismic interval velocity；
Strata pressure correction module, for asking for calibrating (base measuring) pressure coefficient, according to the calibrating (base measuring) pressure coefficient to the stratum
Pressure is corrected.
Preferably, the interval velocity acquisition module calculates seismic interval velocity by below equation (6)：
Wherein, VR,nAnd VR,n-1Respectively stratum top interface, the stack velocity of bottom interface,WithRespectively stratum is pushed up
Interface, the root mean sequare velocity of bottom interface, t0,nAnd t0,n-1Respectively stratum top interface, the reflection interval of bottom interface, VnFor earthquake layer
Speed.
Preferably, it is described to ask for improved seismic interval velocity and include：
Based on acoustic logging and Natural Gamma-ray Logging Curves reconstruct pseudo-acoustic curve DT；
Based on the pseudo-acoustic curve DT, pass through the seismic interval velocity of BP neural network algorithm computed improved.
Preferably, multiple wells in work area will be studied and correspond to observed pressure value and meter at the interval velocity and the destination layer of destination layer
The ratio calculated between pressure value carries out index return, asks for the calibrating (base measuring) pressure coefficient k '.
Compared with prior art, the beneficial effects of the present invention are on the basis of Fillippone pressure prediction empirical equations
On, inversion of simulated AC curve interval velocity is selected, so as to obtain the improvement interval velocity that can more reflect actual formation lithology, using changing
It is more reasonable to enter the strata pressure that interval velocity calculates.In addition, the present invention passes through interval velocity based on actual work area, actual formation
Statistical fit with pressure correcting coefficient obtains calibrating (base measuring) pressure coefficient, and pressure prediction empirical equation is corrected so that calculates
Pressure be more suitable for study work area, it is more efficient, reasonable.
Brief description of the drawings
Disclosure exemplary embodiment is described in more detail in conjunction with the accompanying drawings, the disclosure it is above-mentioned and other
Purpose, feature and advantage will be apparent.
Fig. 1 shows the flow chart of the Shale Pressure Forecasting Methodology based on well shake joint speed according to exemplary embodiment；
Fig. 2 shows what is predicted using the Shale Pressure Forecasting Methodology based on well shake joint speed of exemplary embodiment
JSB constituencies shale formation pressure correcting coefficient crosses well profile figure；
Fig. 3 shows what is predicted using the Shale Pressure Forecasting Methodology based on well shake joint speed of exemplary embodiment
JSB constituencies shale formation pressure correcting coefficient flat distribution map.
Embodiment
Preferred embodiment of the present disclosure is more fully described below with reference to accompanying drawings.Although the disclosure is shown in accompanying drawing
Preferred embodiment, however, it is to be appreciated that may be realized in various forms the disclosure without should be limited by embodiments set forth here
System.On the contrary, these embodiments are provided so that the disclosure is more thorough and complete, and can be complete by the scope of the present disclosure
Ground is communicated to those skilled in the art.
The it is proposed of pressure prediction empirical equation
At present, the Formation pressure prediction method related to seismic data mainly has compacting equilibrium equation method, equivalent depth
Equation, Fillippone methods and improvement Fillippone methods, these four methods are all close with length velocity relation, i.e., speed is
The key of Formation pressure prediction method based on seismic data.The present invention on the basis of Fillippone pressure prediction methods,
The precision of pressure prediction method is improved by improving speed and correction coefficient, to be adapted to the demand of research area's prediction of formation pressure.
Fillippone prediction of formation pressure formula can be expressed as below equation (1)：
Wherein, PpFor strata pressure (MPa)；H is depth (m)；ViFor i-th layer of interval velocity (m/s)；VmaxFor maximum layer speed
Spend (m/s), be rock stratum effecive porosity close to velocity of longitudinal wave when zero；VminIt is rock stratum rigidity for minimum interval velocity (m/s)
Close to velocity of longitudinal wave when zero；K is pressure coefficient, dimensionless；ρAFor formation average denstiy, generally Gardner formula meters are used
Calculate, forecast for seismic data can also be passed through.
Parameter V in formula (1)max、VminIt can be tried to achieve respectively from by below equation (2) and formula (3)：
Vmax=1.4V0+3KT (2)
Vmin=1.4V0+0.5KT (3)
Wherein：
V0=1.4VJ-KT0 (4)
In formula (2)-(5), T0With T be respectively stratum top interface and bottom interface round trip when, VJ0And VJRespectively push up
Interface and the root mean sequare velocity of bottom interface.
The improvement of seismic velocity
Understand that interval velocity has vital work to the pressure prediction based on seismic data based on theory analysis above
With.Therefore, core of the invention be from the angle for improving interval velocity computational accuracy, comprehensive utilization seismic inversion speed,
Shake data, well-log information propose a kind of new, high-precision Formation pressure prediction method.The core of this method is to improve layer speed
The computational accuracy of degree, prediction of formation pressure is carried out using Fillippone formula.In addition, this method combines actual area and stratum
Feature calibrating (base measuring) pressure coefficient k, it is determined that the final pressure prediction formula with Regional suitability.
It is described in detail below with reference to accompanying drawing and is predicted according to the Shale Pressure that joint speed is shaken based on well of exemplary embodiment
Method, it comprises the following steps：
Step 1：Obtain seismic interval velocity
Obtaining the method for spread speed of the seismic wave in subsurface formations mainly includes：Acoustic logging method, VSP methods, ground
Seismic wave method, laboratory method etc..In the exemplary embodiment, the spread speed of seismic wave is obtained using surface seismic method,
Then seismic interval velocity is calculated by DIX formula.
In sub-step 11, using the surface seismic engineering of outdoor, collection receives artificial earthquake signal, passes through seismic signal
The steps such as compiling, first break pickup, normal-moveout spectrum pickup, velocity modeling obtain stack velocity, root mean sequare velocity.The step belongs to existing
Technology, is not the key content of the present invention, therefore is not specifically described herein.
Then, it is horizontal or under conditions of rising and falling less in stratum media in sub-step 12, using DIX formula by sub-step
The stack velocity, the root mean sequare velocity that are obtained in 11 are converted into the interval velocity of seismic interval velocity, i.e. single-layer medium, such as below equation
(6) shown in：
Wherein, VR,nAnd VR,n-1Respectively stratum top interface, the stack velocity of bottom interface,WithRespectively stratum is pushed up
Interface, the root mean sequare velocity of bottom interface, t0,nAnd t0,n-1Respectively stratum top interface, the reflection interval of bottom interface, VnFor layer speed
Degree.
Step 2：Ask for improved seismic interval velocity
Although the method that interval velocity is calculated by DIX formula is simple, in theory, this formula is only applicable to level
The situation of laminar formation and small geophone offset, poor anti jamming capability.In order to improve the seismic interval velocity being converted to using DIX formula
The precision of pressure prediction is carried out, proposes to utilize the curve weight constrained based on well-log information, seismic data in the exemplary embodiment
Structure-onomatopoeia wave parameter inversion method calculates seismic interval velocity.It is realized by following two sub-steps：
Sub-step 21：Based on acoustic logging and Natural Gamma-ray Logging Curves reconstruct pseudo-acoustic curve DT
This method is realized according to the well log curve and activity curve related to stratum characteristic, is intersected by curve,
Reasonable selection characterizes the log parameter of stratum characteristic, the new property parameters of reconstruct reflection stratum characteristic.Sound wave is based in embodiment
Curve and gamma ray curve rebuild the pseudo-acoustic curve of a reflection stratum characteristic, are subsequently used for Application of Logging-constrained Inversion.
Pseudo-acoustic curve DT had both had the low-frequency information of stratum background, and can reflection formation lithology change.Reconstruct pseudo-acoustic curve DT's
Method refers to patent《A kind of pseudo-acoustic curve reconstruct and Sparse Pulse joint inversion method》, Application No.
CN201410213216.6, the content of the patent are all incorporated herein by reference.
Sub-step 22：Based on pseudo-acoustic curve DT, pass through the seismic interval velocity of BP neural network algorithm computed improved
Pseudo-acoustic curve DT had both reflected the seismic velocity information of Different Strata, reflected formation lithology change again.Different Strata
Seismic velocity and pseudo-acoustic curve numerical value into reciprocal relation, gathered by way of artificial earthquake, handle the earthquake number of acquisition
According to the seismic amplitude data that can reflect Different Strata Seismic reflection character, therefore to convert what is obtained by pseudo-acoustic curve DT
Formation velocity is input value (more mouthfuls of drilling wells can provide multigroup formation velocity data in work area), using geological data as constraint,
The improved seismic interval velocity in whole work area is calculated using BP neural network algorithm.
Step 3：Strata pressure is calculated according to improved seismic interval velocity
After improved seismic interval velocity is obtained, formula (1) can be utilized to calculate strata pressure.In initial situation
Under, rule of thumb it can determine pressure coefficient k by formula.For example, under initial situation, can be first in order to not allow k to influence to calculate data
Its assignment 1 is given, then it is modified by step 4 again.
Step 4：Calibrating (base measuring) pressure coefficient is asked for, strata pressure is corrected according to calibrating (base measuring) pressure coefficient
For Fillippone pressure prediction formula, whole district stratum could be improved by being only improved pressure coefficient k precision
Pressure prediction result.Therefore, multiple wells correspond to the interval velocity V of destination layer and observed pressure at the destination layer in research on utilization work area
It is worth calculating ratio (i.e. calibrating (base measuring) pressure coefficient k ') the progress exponential regression analysis between pressure value, goes to seek calibrating (base measuring) pressure coefficient.
In more detail, to study, multiple wells in work area correspond to the interval velocity V of destination layer, observed pressure value is pressed with calculating at the destination layer
Ratio (i.e. calibrating (base measuring) pressure coefficient k ') between force value is variable, exponential relationship therebetween is established, then in step 2
The improved seismic interval velocity of calculating, the observed pressure value that actual drilling test obtains, the pressure value calculated by step 3, are selected
Multi-group data is taken, using the undetermined coefficient in least square method gauge index relation.It can ask for correcting by this method
Pressure coefficient k '.By the calibrating (base measuring) pressure coefficient k ' substitute into formula (1) can calculate correction after strata pressure, the pressure so calculated
Power is more rationally effective.
For example, in one embodiment, draw calibrating (base measuring) pressure coefficient k=0.8e5E-0.5V, wherein v expression interval velocities, 5E-
0.5V represents that 5*10^ { -0.5V }, e represent natural logrithm.
Embodiment
Step 1：It is goal in research stratum to select the peak group shale of river southeast JSB work areas Lower Silurian Series Longma small stream group-five, is surrounded
Target requirement, target of the pressure prediction empirical equation (1) as research is introduced, finds the parameter that can be improved and improve.
Step 2：Initial seismic interval velocity is obtained by Processing Seismic Data.
Step 3：The seismic interval velocity obtained for step 2, it is higher to go out precision using the inversion of simulated AC curve of reconstruct
Improve interval velocity.
Step 4：The improvement interval velocity obtained using step 3, the pressure of shale formation is calculated according to formula (1).
Step 5：In order to improve the applicability of pressure prediction and validity, target is corresponded to by the way that multiple wells in work area will be studied
Observed pressure value and ratio (i.e. calibrating (base measuring) pressure the coefficient k ') progress between pressure value is calculated at the interval velocity of layer and the destination layer
Linear regression, draw calibrating (base measuring) pressure coefficient k ', utilize the calibrating (base measuring) pressure coefficient k ', rule of thumb formula (1) calculates shale formation
Pressure.Fig. 2 and Fig. 3 is shown respectively to be predicted using the Shale Pressure Forecasting Methodology based on well shake joint speed of exemplary embodiment
The JSB constituencies shale formation pressure correcting coefficient gone out crosses well profile figure and flat distribution map.
Above-mentioned technical proposal is a kind of embodiment of the present invention, for those skilled in the art, in this hair
On the basis of bright principle disclosed, it is easy to make various types of improvement or deformation, it is above-mentioned specific to be not limited solely to the present invention
The description of embodiment, therefore description above is simply preferable, and not restrictive meaning.
Claims (8)
1. a kind of shale formation pressure prediction method based on well shake joint speed, comprises the following steps：
Obtain seismic interval velocity；
Ask for improved seismic interval velocity；
Strata pressure is calculated based on the improved seismic interval velocity；
Calibrating (base measuring) pressure coefficient is asked for, the strata pressure is corrected according to the calibrating (base measuring) pressure coefficient.
2. the shale formation pressure prediction method according to claim 1 based on well shake joint speed, wherein, by following
Formula (6) calculates seismic interval velocity：
<mrow>
<msubsup>
<mi>V</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<msubsup>
<mi>V</mi>
<mrow>
<mi>R</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>V</mi>
<mrow>
<mi>R</mi>
<mo>,</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, VR,nAnd VR,n-1Respectively stratum top interface, the stack velocity of bottom interface,WithRespectively stratum top circle
Face, the root mean sequare velocity of bottom interface, t0,nAnd t0,n-1Respectively stratum top interface, the reflection interval of bottom interface, VnFor earthquake layer speed
Degree.
3. the shale formation pressure prediction method according to claim 1 based on well shake joint speed, wherein, it is described to ask for
Improved seismic interval velocity includes：
Based on acoustic logging and Natural Gamma-ray Logging Curves reconstruct pseudo-acoustic curve DT；
Based on the pseudo-acoustic curve DT, pass through the seismic interval velocity of BP neural network algorithm computed improved.
4. the shale formation pressure prediction method according to claim 1 based on well shake joint speed, wherein, work will be studied
Multiple wells correspond to observed pressure value and the ratio calculated between pressure value at the interval velocity and the destination layer of destination layer and carried out in area
Index return, ask for the calibrating (base measuring) pressure coefficient k '.
5. a kind of shale formation pressure prediction device based on well shake joint speed, including：
Interval velocity acquisition module, for obtaining seismic interval velocity；
Interval velocity improves module, for asking for improved seismic interval velocity；
Formation pressure calculation module, for calculating strata pressure based on the improved seismic interval velocity；
Strata pressure correction module, for asking for calibrating (base measuring) pressure coefficient, according to the calibrating (base measuring) pressure coefficient to the strata pressure
It is corrected.
6. the shale formation pressure prediction device according to claim 1 based on well shake joint speed, wherein, the layer speed
Spend acquisition module and seismic interval velocity is calculated by below equation (6)：
<mrow>
<msubsup>
<mi>V</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<msubsup>
<mi>V</mi>
<mrow>
<mi>R</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>V</mi>
<mrow>
<mi>R</mi>
<mo>,</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msub>
<mi>t</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, VR,nAnd VR,n-1Respectively stratum top interface, the stack velocity of bottom interface,WithRespectively stratum top circle
Face, the root mean sequare velocity of bottom interface, t0,nAnd t0,n-1Respectively stratum top interface, the reflection interval of bottom interface, VnFor earthquake layer speed
Degree.
7. the shale formation pressure prediction device according to claim 5 based on well shake joint speed, wherein, it is described to ask for
Improved seismic interval velocity includes：
Based on acoustic logging and Natural Gamma-ray Logging Curves reconstruct pseudo-acoustic curve DT；
Based on the pseudo-acoustic curve DT, pass through the seismic interval velocity of BP neural network algorithm computed improved.
8. the shale formation pressure prediction device according to claim 5 based on well shake joint speed, wherein, work will be studied
Multiple wells correspond to observed pressure value and the ratio calculated between pressure value at the interval velocity and the destination layer of destination layer and carried out in area
Index return, ask for the calibrating (base measuring) pressure coefficient k '.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610835388.6A CN107843927A (en) | 2016-09-20 | 2016-09-20 | Shale formation pressure prediction method and device based on well shake joint speed |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610835388.6A CN107843927A (en) | 2016-09-20 | 2016-09-20 | Shale formation pressure prediction method and device based on well shake joint speed |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107843927A true CN107843927A (en) | 2018-03-27 |
Family
ID=61656738
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610835388.6A Pending CN107843927A (en) | 2016-09-20 | 2016-09-20 | Shale formation pressure prediction method and device based on well shake joint speed |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107843927A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108717202A (en) * | 2018-05-31 | 2018-10-30 | 中国地质大学（武汉） | A kind of shale gas abnormal formation pressure prediction technique considering fluid temperature (F.T.) |
| CN109931054A (en) * | 2018-12-27 | 2019-06-25 | 西南石油大学 | The prediction technique of tight sandstone reservoir pressure |
| CN111060986A (en) * | 2019-10-18 | 2020-04-24 | 中国石油化工股份有限公司 | Formation pressure prediction method and lithologic oil reservoir evaluation method |
| CN111856567A (en) * | 2019-04-30 | 2020-10-30 | 中国石油天然气集团有限公司 | Seismic stack velocity-based compaction coefficient determination method and device |
| CN112100930A (en) * | 2020-11-11 | 2020-12-18 | 中国石油大学(华东) | Formation pore pressure calculation method based on convolutional neural network and Eaton formula |
| CN112241025A (en) * | 2019-07-18 | 2021-01-19 | 中国石油天然气股份有限公司 | Well-seismic combined formation pressure determination method and system |
| CN112906465A (en) * | 2021-01-15 | 2021-06-04 | 阳泉煤业（集团）股份有限公司 | Coal measure stratum acoustic curve reconstruction method and system based on stratum factors |
| CN117434599A (en) * | 2023-08-08 | 2024-01-23 | 浙江大学 | Method for predicting formation pressure based on seismic data |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102129084A (en) * | 2010-12-17 | 2011-07-20 | 中国石油天然气股份有限公司 | Method and device for acquiring seismic velocity in thin reservoir layer through well control |
| CN103792581A (en) * | 2013-12-31 | 2014-05-14 | 刘致水 | Shale gas reservoir crustal stress logging prediction method based on rock physics model |
| CN104267429A (en) * | 2014-09-30 | 2015-01-07 | 中国石油天然气股份有限公司 | Method and device for determining formation pressure |
| CN104483703A (en) * | 2014-12-31 | 2015-04-01 | 中国石油天然气股份有限公司 | Well-to-seismic integration paleo-formation pressure prediction method |
| CN105089652A (en) * | 2014-05-20 | 2015-11-25 | 中国石油化工股份有限公司 | Pseudo-acoustic curve rebuilding and sparse pulse joint inversion method |
| CN105549082A (en) * | 2014-10-29 | 2016-05-04 | 中国石油天然气股份有限公司 | Establishing method and system of three-dimensional geomechanical field of extra-deep carbonate reservoir |
-
2016
- 2016-09-20 CN CN201610835388.6A patent/CN107843927A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102129084A (en) * | 2010-12-17 | 2011-07-20 | 中国石油天然气股份有限公司 | Method and device for acquiring seismic velocity in thin reservoir layer through well control |
| CN103792581A (en) * | 2013-12-31 | 2014-05-14 | 刘致水 | Shale gas reservoir crustal stress logging prediction method based on rock physics model |
| CN105089652A (en) * | 2014-05-20 | 2015-11-25 | 中国石油化工股份有限公司 | Pseudo-acoustic curve rebuilding and sparse pulse joint inversion method |
| CN104267429A (en) * | 2014-09-30 | 2015-01-07 | 中国石油天然气股份有限公司 | Method and device for determining formation pressure |
| CN105549082A (en) * | 2014-10-29 | 2016-05-04 | 中国石油天然气股份有限公司 | Establishing method and system of three-dimensional geomechanical field of extra-deep carbonate reservoir |
| CN104483703A (en) * | 2014-12-31 | 2015-04-01 | 中国石油天然气股份有限公司 | Well-to-seismic integration paleo-formation pressure prediction method |
Non-Patent Citations (2)
| Title |
|---|
| 朱超等: "声波曲线重构技术在火山岩储层预测中的应用", 《石油天然气学报（江汉石油学院学报）》 * |
| 罗胜元等: "渤南洼陷钻前地层压力地震预测", 《石油地球物理勘探》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108717202A (en) * | 2018-05-31 | 2018-10-30 | 中国地质大学（武汉） | A kind of shale gas abnormal formation pressure prediction technique considering fluid temperature (F.T.) |
| CN108717202B (en) * | 2018-05-31 | 2020-01-17 | 中国地质大学（武汉） | Shale gas abnormal formation pressure prediction method considering fluid temperature |
| CN109931054A (en) * | 2018-12-27 | 2019-06-25 | 西南石油大学 | The prediction technique of tight sandstone reservoir pressure |
| CN111856567B (en) * | 2019-04-30 | 2023-09-26 | 中国石油天然气集团有限公司 | Compaction coefficient determination method and device based on earthquake superposition speed |
| CN111856567A (en) * | 2019-04-30 | 2020-10-30 | 中国石油天然气集团有限公司 | Seismic stack velocity-based compaction coefficient determination method and device |
| CN112241025A (en) * | 2019-07-18 | 2021-01-19 | 中国石油天然气股份有限公司 | Well-seismic combined formation pressure determination method and system |
| CN112241025B (en) * | 2019-07-18 | 2023-11-28 | 中国石油天然气股份有限公司 | Well-seismic joint formation pressure determination method and system |
| CN111060986A (en) * | 2019-10-18 | 2020-04-24 | 中国石油化工股份有限公司 | Formation pressure prediction method and lithologic oil reservoir evaluation method |
| CN112100930A (en) * | 2020-11-11 | 2020-12-18 | 中国石油大学(华东) | Formation pore pressure calculation method based on convolutional neural network and Eaton formula |
| CN112100930B (en) * | 2020-11-11 | 2021-02-02 | 中国石油大学(华东) | Formation pore pressure calculation method based on convolutional neural network and Eaton formula |
| CN112906465A (en) * | 2021-01-15 | 2021-06-04 | 阳泉煤业（集团）股份有限公司 | Coal measure stratum acoustic curve reconstruction method and system based on stratum factors |
| CN112906465B (en) * | 2021-01-15 | 2023-12-22 | 阳泉煤业（集团）股份有限公司 | Coal measure stratum acoustic curve reconstruction method and system based on stratum factors |
| CN117434599A (en) * | 2023-08-08 | 2024-01-23 | 浙江大学 | Method for predicting formation pressure based on seismic data |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107843927A (en) | Shale formation pressure prediction method and device based on well shake joint speed | |
| CN103258091B (en) | Unconventionaloil pool hides the method and device that net horizontal section three-dimensional mechanical models for rock mass is set up | |
| CN101872024B (en) | Method for carrying out well design by using time-lapse seismic | |
| CN104453873A (en) | Shale oil-gas economical and effective layer section evaluation method | |
| CN103576195B (en) | A kind of pressure-dependent fissuted medium S-Wave Velocity Predicted Method | |
| CN105277982B (en) | A kind of mud shale total content of organic carbon earthquake prediction method | |
| CN105221141A (en) | A kind of mud shale brittleness index Forecasting Methodology | |
| CN112746837A (en) | Shale oil reservoir exploration data acquisition system and method based on distributed optical fiber sensing | |
| CN101887132A (en) | Method for quantificationally predicting sandstone reservoir fluid saturation by combining well and seism | |
| CN104345345A (en) | Method for forecasting TOC (Total Organic Carbon) content of shale reservoir stratum | |
| CN105242307B (en) | Carbonate complex seismic reservoir porosity acquisition methods and device | |
| CN104181585A (en) | Shear wave estimation method and system in geophysical exploration | |
| CN104749617A (en) | Multi-scale fractured reservoir forward model establishing method | |
| CN104155693A (en) | Angle gather seismic response numerical computation method of reservoir fluid fluidity | |
| CN107817520A (en) | The pressure coefficient Forecasting Methodology and system of marine facies mud shale stratum | |
| KR101415196B1 (en) | Method for selecting fracking intervals of horizontal drilling zone in sweet spot range using resistivity and density logging data in shale play | |
| CN109509111A (en) | The prediction technique and system of prospect pit strata pressure | |
| CN104407381A (en) | Method for acquiring stratum rock brittleness property through radial change of borehole elastic wave velocity | |
| CN110058323A (en) | A kind of tight sand formation brittleness index calculation method | |
| CN214576965U (en) | Shale oil reservoir exploration data acquisition system based on distributed optical fiber sensing | |
| CN107728205B (en) | A kind of Formation pressure prediction method | |
| CN103777245A (en) | Method of quantitative evaluation on oil and gas accumulation conditions based on seismic data | |
| KR101415197B1 (en) | Method for selecting fracking intervals of horizontal drilling zone in sweet spot range using resistivity and neutron logging data in shale play | |
| CN102288996B (en) | Method for quantificationally predicting gas containing probability through prestack automatic fault indication (AFI) inversion | |
| KR20180014294A (en) | Method for estimating veloscity of S wave from other logging 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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180327 |
|
| RJ01 | Rejection of invention patent application after publication |