The Forecasting Methodology of formation pore pressure between a kind of salt
Technical field
The present invention relates to shale oil exploration and development technology, particularly relate to the Forecasting Methodology of formation pore pressure between a kind of salt.
Background technology
Research shows that shale oil reservoir between salt " dessert district " is main by strata pressure control, is laminated the most in detail between research salt
Power has highly important realistic meaning.
Prediction of formation pressure can be divided into two classes: one is to utilize well-log information to carry out prediction of formation pressure, and two is to utilize earthquake
Data does prediction of formation pressure.
The first kind utilizes well-log information predictably stressor layer to conclude 3 kinds of methods: A1 traditional method, is characterized in setting up pressure
Real Trendline, then estimates strata pressure according to the degree of measured value deviation Trendline;Or set up measured value and ground lamination
Empirical relation between power, and estimate strata pressure with this.The method is applicable to sand-mudstone formation.The traditional method that A2 improves, changes
Enter and be a little mainly around how setting up what normal compaction trend line was carried out.A3 comprehensive parameters method, its advantage is to need not set up
Normal compaction trend line, it was predicted that consider the impact of shale content in model, utilizes porosity, shale content curve, adds
The control result of lithology is more identical with geological conditions.The method is applicable to determine the abnormal stratum Pore Pressure of the mechanism of being differently formed
Power.Generally speaking, the result precision that the first kind obtains is higher, effect preferable, but such method can only be applied after drilling well terminates,
And it is unpredictable to the strata pressure below shaft bottom.
Equations of The Second Kind utilizes seismic data to do prediction of formation pressure and may be summarized to be 3 kinds of methods: B1 classical way, also referred to as
Fillippone method, its advantage is to be independent of normal compaction trend line, calculates simple, applicable.The method can be directly by earthquake
Interval velocity calculates strata pressure, and the area that, area few for drilling well is bigger carries out pressure prediction, is the method for unique feasible.B2
The classical way improved, improvement is how to improve the computational accuracy of seismic interval velocity.B3 many seismic properties method, its feature
It is to combine multiple seismic properties, by considering that genesis of pressure sets the different weight coefficient factors and is analyzed.Equations of The Second Kind method
On the basis of being all built upon the data such as drilling well, well logging, geology, earthquake are carried out comprehensive study, data is the most comprehensive, research is the deepest
Enter, it was predicted that the degree of reliability of result is the highest.
At present, most of prediction of formation pressure is all to carry out from the angle of mechanics, and few people are from geological angle
Degree macroscopically holds abnormal pressure Forming Mechanism, and selects and set up corresponding Forecasting Methodology.
Summary of the invention
The technical problem to be solved in the present invention is for defect of the prior art, it is provided that formation pore pressure between a kind of salt
The Forecasting Methodology of power.
The technical solution adopted for the present invention to solve the technical problems is: the prediction side of formation pore pressure between a kind of salt
Method, comprises the following steps:
1) choosing the interior different configuration position of stratum depression and survey between salt has interval transit time, density curve, physical properties of rock to test
Data and the drilling well of strata pressure, set up drilling well log data and measured formation pressure data storehouse;Described physical properties of rock is surveyed
Examination data includes porosity and shale content data;
2) measured formation pressure data database data is utilized, the theoretical and formula P according to Terzaghi effective stressp=Po-Pe,
Calculate the vertical effective stress of stress test well section;
P in formulaeFor vertical effective stress, P0For overburden pressure, PpFor formation pore pressure;
3) add up the log value of relevant pressure testing well section, make log value relevant to vertical effective stress
Property plate, determine vertical effective stress and porosity, shale content, interval transit time and the dependency of density, and set up mathematical modulo
Type;
In formula, VpFor SVEL, unit km/s;
For calculating porosity, numerical range 0-1;
VshFor shale content, numerical range 0-1;
PeFor vertical effective stress, unit K Bar, convert 1KBar=100Mpa;
P0For overlying formation pressure, unit K Bar;
PpFor pore pressure, unit K Bar;
A0, A1, A2, A3, A4For regional experience model parameter;
4) Multiple non-linear is utilized to determine mathematical model parameter;
5) relevant between ratio and the interval velocity of correspondent section of the strata pressure that statistics makes Measured formation pressure and prediction
Property plate, determines that both have exponential relationship, and sets up corresponding mathematical model;
PpActual measurement=C (Vp) * PpPrediction
In formula, C (Vp) is correction coefficient;
6) according to dependency plate between observed pressure and ratio and the interval velocity of correspondent section, curvilinear regression matching is utilized to obtain
Correction coefficient to model;
C (Vp)=0.9548e(0.00001*Vp), wherein Vp is SVEL;
7) obtain final mask, and use the strata pressure Pp of model prediction stepout wellPrediction:
The beneficial effect comprise that: between the salt that the present invention proposes, stratum calculates formation pore pressure Forecasting Methodology,
Consider the impact of shale content, utilized porosity, shale content curve, add the control knot of lithology and structural factor
Fruit is more identical with geological conditions.Practical is strong, it was predicted that accuracy rate high.
Accompanying drawing explanation
Below in conjunction with drawings and Examples, the invention will be further described, in accompanying drawing:
Fig. 1 is the workflow diagram of the embodiment of the present invention;
Fig. 2 is vertical effective stress and the porosity dependency plate of the embodiment of the present invention;
Fig. 3 is vertical effective stress and the shale content of the embodiment of the present invention;
Fig. 4 is vertical effective stress and the interval transit time dependency plate of the embodiment of the present invention;
Fig. 5 is vertical effective stress and the sound wave density dependence plate of the embodiment of the present invention;
Fig. 6 is pressure ratio and the interval velocity graph of a relation of the embodiment of the present invention.
Detailed description of the invention
In order to make the purpose of the present invention, technical scheme and advantage clearer, below in conjunction with embodiment, to the present invention
It is further elaborated.Should be appreciated that specific embodiment described herein, only in order to explain the present invention, is not used to limit
Determine the present invention.
As it is shown in figure 1, the Forecasting Methodology of formation pore pressure between a kind of salt, comprise the following steps:
1) choose the interior different configuration position of Yan Jian stratum, Qianjiang depression and survey has interval transit time, density curve and physical properties of rock
Test data (porosity, shale content etc.) and the drilling well of strata pressure, be laminated with setting up drilling well log data and actual measurement
Force data storehouse;
2) measured formation pressure data database data is utilized, the theoretical and formula P according to Terzaghi (1943) effective stressp=
Po-Pe, calculate the vertical effective stress of stress test well section;P in formulaeFor vertical effective stress, P0For overlying rock pressure
Power, PpFor formation pore pressure;
3) utilize drilling well log data statistics relevant pressure testing well section log value, make log value with
The dependency plate of vertical effective stress, is shown in accompanying drawing 2 to 5, determines vertical effective stress and porosity, shale content, sound wave
The time difference and density have preferable dependency, and founding mathematical models;
In formula, Vp is SVEL, unit km/s;
For calculating porosity, numerical range 0-1;
VshFor shale content, numerical range 0-1;
PeFor vertical effective stress, unit K Bar, convert 1KBar=100Mpa;
P0For overlying formation pressure, unit K Bar;
PpFor pore pressure, unit K Bar;
A0, A1, A2, A3, A4For regional experience model parameter;
4) Multiple non-linear is utilized to determine model parameter;
5) relevant between ratio and the interval velocity of correspondent section of the strata pressure that statistics makes Measured formation pressure and prediction
Property plate, is shown in Fig. 6, determines that both have exponential relationship, and set up corresponding mathematical model;
Pp actual measurement=C (Vp) * Pp predicts;
In formula, C (Vp) is correction coefficient;
6) according to dependency plate between observed pressure and ratio and the interval velocity of correspondent section, curvilinear regression matching is utilized to obtain
Correction coefficient to model;
C (Vp)=0.9548e(0.00001*Vp);
7) final mask is formed, and with the strata pressure of this model prediction stepout well:
The application of the present invention:
Embodiment 1:
The present invention applies at Qianjiang Depression Wangchang area Wang Xi 10-3 well.Mainly comprise the steps:
The actual measurement GR curve of Wang Xi 10-3 well, interval transit time, density is used to calculate observed pressure section shale respectively
Content Vsh=0.56, speed Vp=5.27km/s, porosityAdd up the seismic interval velocity v of this interval
=3324m/s.Pressure prediction model is utilized to calculate pore pressure
Pressure prediction result calculates relative error with observed pressure value (33.25Mpa), and relative error is 0.72%
Embodiment 2:
The present invention applies at Qianjiang Depression extensively China area wide 25 oblique-2 wells.Mainly comprise the steps:
Use wide 25 tiltedly-2 actual measurement GR curve, interval transit time, density calculate observed pressure section shale respectively and contain
Amount Vsh=0.3, speed Vp=3.43km/s, porosityAdd up the seismic interval velocity v=of this interval
3109.4m/s.Pressure prediction model is utilized to calculate pore pressure
Pressure prediction result calculates relative error with observed pressure value (31.93Mpa), and relative error is 1.91%.
From example it can be seen that the accuracy rate of the Forecasting Methodology prediction of the present invention is high.
It should be appreciated that for those of ordinary skills, can be improved according to the above description or be converted,
And all these modifications and variations all should belong to the protection domain of claims of the present invention.