CN109425904A - A kind of Pore Pressure on Carbonate Rock Formation acquisition methods - Google Patents
A kind of Pore Pressure on Carbonate Rock Formation acquisition methods Download PDFInfo
- Publication number
- CN109425904A CN109425904A CN201710761717.1A CN201710761717A CN109425904A CN 109425904 A CN109425904 A CN 109425904A CN 201710761717 A CN201710761717 A CN 201710761717A CN 109425904 A CN109425904 A CN 109425904A
- Authority
- CN
- China
- Prior art keywords
- characteristic
- rock
- formation
- bulk modulus
- indicate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- 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
- G01V2210/6248—Pore pressure
Abstract
A kind of Pore Pressure on Carbonate Rock Formation acquisition methods, this method determines the formation pore pressure characteristic in carbonate formation to be analyzed at each depth point using the Pore Pressure on Carbonate Rock Formation computation model of building, it establishes formation pore pressure model using pore-fluid acoustic speed, it can more accurately calculate formation pore pressure, it for Drilling Design, selects and determines suitable safe drilling fluid density scientific basis is provided, to effectively prevent the generation of underground complex accident during drilling construction.
Description
Technical field
The present invention relates to oil exploration technology fields, specifically, being related to a kind of Pore Pressure on Carbonate Rock Formation acquisition
Method.
Background technique
Formation pore pressure, especially abnormal pore pressure bring many difficulties to the exploration and exploitation of oil-gas reservoir.At present
For Clastic Stratum of Country Rocks pore pressure, formation mechenism, model, method and prediction are more mature, can accurately calculate broken
Consider rock formation pore pressure, and precision with higher to be worth doing.
And for carbonate formation, at present engineering and theoretical educational circles using it is more be by means of Clastic Stratum of Country Rocks hole
The computation model and method of pressure, calculate Pore Pressure on Carbonate Rock Formation on this basis.But carbonate formation abnormal hole
The formation mechenism of gap pressure is different from Clastic Stratum of Country Rocks, such as computation model simply by Clastic Stratum of Country Rocks pore pressure and side
Method calculates Pore Pressure on Carbonate Rock Formation, it will causes biggish error, or even at abnormal high pressure by normal pressure stratigraphic interpretation
Layer.Therefore calculating Pore Pressure on Carbonate Rock Formation using the computation model of Clastic Stratum of Country Rocks pore pressure and method at present is not
Suitably, it is necessary to form a kind of model and method that can effectively calculate Pore Pressure on Carbonate Rock Formation.
Summary of the invention
To solve the above problems, the present invention provides a kind of Pore Pressure on Carbonate Rock Formation acquisition methods, the method
Include:
Step 1: the rock forming mineral component characteristics data in carbonate formation to be analyzed at each depth point are obtained, according to
The rock forming mineral component characteristics data determine the rock solid volume in the carbonate formation to be analyzed at each depth point
Elasticity modulus characteristic;
Step 2: the velocity of longitudinal wave characteristic of the carbonate formation to be analyzed is obtained, according to the velocity of longitudinal wave
Characteristic calculates the rock matrix bulk modulus characteristic in the carbonate formation to be analyzed at each depth point;
Step 3: according to the saturated rock volume bullet in the carbonate formation to be analyzed got at each depth point
Property modulus characteristic, porosity characteristics data, in conjunction with the rock solid bulk modulus characteristic and rock bone
Frame bulk modulus characteristic calculates the pore-fluid volume bullet in the carbonate formation to be analyzed at each depth point
Property modulus characteristic;
Step 4: the overburden pressure characteristic in carbonate formation to be analyzed at each depth point is obtained, according to
The overburden pressure characteristic, pore-fluid bulk modulus characteristic, saturated rock bulk modulus are special
It is special to levy data, rock solid bulk modulus characteristic, rock matrix bulk modulus characteristic and porosity
Sign data determine the formation pore pressure characteristic in the carbonate formation to be analyzed at each depth point.
According to one embodiment of present invention, in the step 4, the carbon to be analyzed is calculated according to following expression
Formation pore pressure characteristic at the depth point Zhong Ge of Carbonate Rocks stratum:
Wherein, PpIndicate formation pore pressure characteristic, KfIndicate pore-fluid bulk modulus characteristic, Kfr
Indicate rock matrix bulk modulus characteristic, φ indicates porosity, KsIndicate rock solid bulk modulus feature
Data, P0Indicate overburden pressure characteristic.
According to one embodiment of present invention, in said step 1, using VRH mean value method come according to the rock mine
Object component characteristics data determine that the rock solid bulk modulus in the carbonate formation to be analyzed at each depth point is special
Levy data.
According to one embodiment of present invention, in the step 2, velocity of longitudinal wave characteristic is decomposed, is obtained
The skeleton velocity of longitudinal wave characteristic of formation rock in the carbonate formation to be analyzed at each depth point, according to the bone
Frame velocity of longitudinal wave characteristic calculates the rock matrix volume elasticity mould in the carbonate formation to be analyzed at each depth point
Measure feature data.
According to one embodiment of present invention, in the step 2, it is calculated at each depth point according to following expression
Rock matrix bulk modulus characteristic:
Wherein, KfrIndicate rock matrix bulk modulus characteristic, φ indicates porosity, ρsIndicate that rock solid is close
Spend characteristic, νp,frIndicate the skeleton velocity of longitudinal wave characteristic of formation rock, νsIndicate the shear wave velocity feature of formation rock
Data.
According to one embodiment of present invention, in the step 3, according to the carbonate rock to be analyzed got
Velocity of longitudinal wave characteristic, shear wave velocity characteristic and density feature data at the depth point Zhong Ge of stratum determine each depth
Saturated rock bulk modulus characteristic at point.
According to one embodiment of present invention, the saturated rock volume elasticity at each depth point is calculated according to following expression
Modulus characteristic:
Wherein, KsatIndicate saturated rock bulk modulus characteristic, ρ indicates density feature data, νpIndicate longitudinal wave
Velocity characteristic data, νsIndicate shear wave velocity characteristic.
According to one embodiment of present invention, the pore-fluid volume elasticity at each depth point is calculated according to following expression
Modulus characteristic:
Wherein, KsatIndicate saturated rock bulk modulus characteristic, KfIndicate that pore-fluid bulk modulus is special
Levy data, KfrIndicate rock matrix bulk modulus characteristic, φ indicates porosity, KsIndicate rock solid volume elasticity
Modulus characteristic.
According to one embodiment of present invention, it in the step 4, obtains each in the carbonate formation to be analyzed
Density feature data, depth characteristic data at depth point, it is true according to the density feature data and depth characteristic data difference
Overburden pressure characteristic at fixed each depth point.
Method provided by the present invention is determined to be analyzed using the Pore Pressure on Carbonate Rock Formation computation model of building
Formation pore pressure characteristic in carbonate formation at each depth point.Wherein, constructed carbonate formation hole
Calculation of pressure model considers the feature of carbonate rock rock physics structure, analyzes carbonate rock petrophysical parameter to stratum
The influence of pore pressure.This method is from the angle of poroelasticity theory of mechanics, the Pore Pressure on Carbonate Rock Formation of foundation
Computation model is overburden pressure, porosity, pore-fluid bulk modulus, rock matrix bulk modulus and rock
Function between solid volume elasticity modulus.
This method establishes formation pore pressure model using pore-fluid acoustic speed, more accurately can calculate ground
Layer pore pressure, for Drilling Design, selects and determines suitable safe drilling fluid density and provides scientific basis, to effectively prevent
The generation of underground complex accident during drilling construction.
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
It obtains it is clear that understand through the implementation of the invention.The objectives and other advantages of the invention can be by specification, right
Specifically noted structure is achieved and obtained in claim and attached drawing.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is required attached drawing in technical description to do simple introduction:
Fig. 1 is the implementation process signal of Pore Pressure on Carbonate Rock Formation acquisition methods according to an embodiment of the invention
Figure;
Fig. 2 is the schematic diagram of rock solid bulk modulus characteristic series according to an embodiment of the invention;
Fig. 3 is the schematic diagram of rock solid density feature data series according to an embodiment of the invention;
Fig. 4 is the schematic diagram of acoustic speed characteristic series according to an embodiment of the invention;
Fig. 5 is the schematic diagram of porosity characteristics data series according to an embodiment of the invention;
Fig. 6 is the schematic diagram of rock density characteristic series according to an embodiment of the invention;
Fig. 7 is the schematic diagram of rock matrix bulk modulus characteristic series according to an embodiment of the invention;
Fig. 8 is the schematic diagram of saturated rock bulk modulus characteristic series according to an embodiment of the invention;
Fig. 9 is the schematic diagram of pore-fluid bulk modulus characteristic series according to an embodiment of the invention;
Figure 10 is the schematic diagram of overburden pressure characteristic series according to an embodiment of the invention;
Figure 11 is the formation pore pressure of depth segment locating for carbonate formation to be analyzed according to an embodiment of the invention
The schematic diagram of power characteristic series.
Specific embodiment
Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings and examples, how to apply to the present invention whereby
Technological means solves technical problem, and the realization process for reaching technical effect can fully understand and implement.It needs to illustrate
As long as not constituting conflict, each feature in each embodiment and each embodiment in the present invention can be combined with each other,
It is within the scope of the present invention to be formed by technical solution.
Meanwhile in the following description, for illustrative purposes and numerous specific details are set forth, to provide to of the invention real
Apply the thorough understanding of example.It will be apparent, however, to one skilled in the art, that the present invention can not have to tool here
Body details or described ad hoc fashion are implemented.
In addition, step shown in the flowchart of the accompanying drawings can be in the department of computer science of such as a group of computer-executable instructions
It is executed in system, although also, logical order is shown in flow charts, and it in some cases, can be to be different from herein
Sequence execute shown or described step.
The existing method for calculating formation pore pressure focuses primarily upon Accurate Prediction Clastic Stratum of Country Rocks pore pressure, and clast
Rock formation pore pressure sedimentation mechanism is different from carbonate rock, cannot directly use the prediction technique of Clastic Stratum of Country Rocks pore pressure
In calculating Pore Pressure on Carbonate Rock Formation.
In addition, it is existing about calculate carbonate porosity pressure method or by establish effective stress, skeleton and
Relationship between pore-fluid coefficient of bulk compressibility and velocity of longitudinal wave, or velocity of longitudinal wave is decomposed using wavelet transformation to mention
The relational model that pore-fluid velocity of longitudinal wave is taken to establish between formation pore pressure and pore-fluid velocity of longitudinal wave.In these methods
Pore Pressure on Carbonate Rock Formation computation model be all empirical relation model, without theoretical foundation.
For the problems of in the prior art, the present invention provides a kind of new Pore Pressure on Carbonate Rock Formation to obtain
Method is taken, this method is based on poroelasticity theory of mechanics, is asked using the relationship between petrophysical parameter and formation pore pressure
Pore Pressure on Carbonate Rock Formation is taken, with universality and theoretical basis.
Fig. 1 shows the implementation process signal of Pore Pressure on Carbonate Rock Formation acquisition methods provided by the present embodiment
Figure.
As shown in Figure 1, this method obtains each in carbonate formation to be analyzed in step s101 first in the present embodiment
Rock forming mineral component characteristics data at depth point.Specifically, in the present embodiment, this method is advantageously according to the geology got
Data determines depth segment locating for carbonate formation to be analyzed.Certainly, in other embodiments of the invention, this method is also
Depth segment locating for carbonate formation to be analyzed can be determined by other related datas such as geology, well logging.Determine to
After analyzing depth segment locating for carbonate formation, this method can also determine stratum required for calculating formation pore pressure
Performance data in the present embodiment, calculates formation characteristics data required for formation pore pressure and preferably includes: depth of stratum number
According to, acoustic speed data and response data of porosity etc..
After determining depth segment locating for carbonate formation to be analyzed, this method can be in step s101 according to scene
Log data analyzes these landwaste to acquire the landwaste that well head returns out, to obtain carbonate formation to be analyzed
In rock forming mineral component characteristics data series at each depth pointAnd rock solid is close
Spend characteristic seriesWherein, h(i)Indicate the depth of the i-th depth point,Indicate jth kind mine at the i-th depth point
The volume fraction of object component, n indicate the sum of mineral constituent at the i-th depth point,Indicate the rock solid at the i-th depth point
Density feature data.
Specifically, in the present embodiment, for the landwaste at collected each depth point, this method is preferably by X diffractometer
It determines the rock forming mineral component of taken landwaste, and weighs to obtain the density of its rock solid using electronic balance scale, to obtain
Rock forming mineral component characteristics data series in carbonate formation to be analyzed at each depth point
And rock solid density feature data series
Certainly, in other embodiments of the invention, this method can also use other rational methods in step s101
To determine the rock forming mineral component characteristics data and rock solid density in carbonate formation to be analyzed at each depth point
Characteristic, the invention is not limited thereto.
As shown in Figure 1, the rock forming mineral component characteristics number in obtaining carbonate formation to be analyzed at each depth point
According to rear, this method can be in step s 102 according to the rock forming mineral group dtex at depth point each in carbonate formation to be analyzed
Data are levied to determine the rock solid bulk modulus characteristic in carbonate formation to be analyzed at each depth point.
Specifically, in the present embodiment, this method is preferably by VRH mean value method and comes according to carbonate formation to be analyzed
In rock forming mineral component characteristics data at each depth point determine the rock in carbonate formation to be analyzed at each depth point
Solid volume elasticity modulus characteristic.
This method calculates the rock solid bulk modulus characteristic at each depth point according to following expression:
Wherein, KsIndicate rock solid bulk modulus characteristic, KjAnd XjRespectively indicate jth kind mineral constituent
Bulk modulus and volume fraction.
I.e. for the rock solid bulk modulus characteristic at the i-th depth pointFor, exist:
Wherein,Indicate the bulk modulus of jth kind mineral constituent at the i-th depth point.
In the present embodiment, the rock of this method depth segment locating for available carbonate formation to be analyzed in step s 102
Stone solid volume elasticity modulus characteristic seriesAnd rock solid density feature data seriesIts
In, Fig. 2 shows the schematic diagram of rock solid bulk modulus characteristic series in the present embodiment, Fig. 3 shows this reality
Apply the schematic diagram of rock solid density feature data series in example.
In step s 103, this method can obtain the velocity of longitudinal wave characteristic series of carbonate formation to be analyzedIn the present embodiment, this method can also obtain carbonic acid to be analyzed required for subsequent calculating process in step s 103
The shear wave velocity characteristic series on rock salt stratumDensity feature data series { h(i),ρ(i)And porosity spy
Levy data series { h(i),φ(i)}.Wherein,ρ(i)And φ(i)Respectively indicate the velocity of longitudinal wave feature at the i-th depth point
Data, shear wave velocity characteristic, density feature data and porosity characteristics data.Fig. 4 shows sound wave in the present embodiment
The schematic diagram of the characteristic series of speed, Fig. 5 show the schematic diagram of the present embodiment porosity characteristic series, Fig. 6
Show the schematic diagram of rock density characteristic series in the present embodiment.
Then, this method can calculate each depth according to the velocity of longitudinal wave characteristic at each depth point in step S104
Rock matrix bulk modulus characteristic at point.Specifically, in the present embodiment, this method in step S104 preferably
The longitudinal wave characteristic at each depth point is decomposed using EMD method, to extract carbonate formation rock to be analyzed
The skeleton velocity of longitudinal wave characteristic series of stoneWherein,Indicate the skeleton of the formation rock at the i-th depth point
Velocity of longitudinal wave characteristic.
In the skeleton velocity of longitudinal wave characteristic series for obtaining formation rockAfterwards, this method can be vertical according to skeleton
Wave velocity characteristic seriesTo calculate the rock matrix volume in carbonate formation to be analyzed at each depth point
Elasticity modulus characteristic, and obtain the rock matrix bulk modulus feature of depth segment locating for carbonate formation to be analyzed
Data seriesWherein, Fig. 7 shows showing for the series of rock matrix bulk modulus characteristic in the present embodiment
It is intended to.
Specifically, in the present embodiment, this method calculates rock matrix bulk modulus characteristic according to following expression
According to:
Wherein, KfrIndicate rock matrix bulk modulus characteristic, φ indicates porosity, ρsIndicate that rock solid is close
Spend characteristic, νp,frIndicate the skeleton velocity of longitudinal wave characteristic of formation rock, νsIndicate the shear wave velocity feature of formation rock
Data.
I.e. for the rock matrix bulk modulus characteristic at the i-th depth pointFor, exist:
Wherein,φ(i)、WithRespectively indicate the rock matrix bulk modulus feature at the i-th depth point
Data, porosity characteristics data, the shear wave velocity characteristic of the skeleton velocity of longitudinal wave characteristic of formation rock and formation rock
According to.
As shown in Figure 1, this method can obtain each depth in carbonate formation to be analyzed in step s105 in the present embodiment
Saturation delay bulk modulus characteristic at degree point, and according to the saturation delay bulk modulus characteristic got
It is special in conjunction with rock solid bulk modulus characteristic and rock matrix bulk modulus according to, porosity characteristics data
Data are levied, the pore-fluid bulk modulus characteristic in carbonate formation to be analyzed at each depth point is calculated.
In the present embodiment, this method is advantageously according to vertical at each depth point in the carbonate formation to be analyzed got
Wave velocity characteristic, shear wave velocity characteristic and density feature data determine the saturated rock volume bullet at each depth point
Property modulus characteristic, and obtain the saturated rock bulk modulus characteristic of depth segment locating for carbonate formation to be analyzed
According to series
Specifically, in the present embodiment, this method calculates the saturated rock at each depth point advantageously according to following expression
Bulk modulus characteristic:
Wherein, KsatIndicate saturated rock bulk modulus characteristic, ρ indicates density feature data, νpIndicate longitudinal wave
Velocity characteristic data, νsIndicate shear wave velocity characteristic.
I.e. for the saturated rock bulk modulus characteristic at the i-th depth pointFor, exist:
Wherein, ρ(i)、WithIndicate density feature data, velocity of longitudinal wave characteristic and the shear wave at the i-th depth point
Velocity characteristic data.Fig. 8 shows the schematic diagram that rock volume elasticity modulus characteristic series is saturated in the present embodiment.
In step s105, this method is preferably by Biot-Gassmann equation to calculate the hole at each depth point
Bulk Elastic Modulus of Fluid characteristic.That is, in the presence of:
Wherein, KfIndicate pore-fluid bulk modulus characteristic.
For the pore-fluid bulk modulus characteristic at the i-th depth pointFor, exist:
Fig. 9 shows pore-fluid bulk modulus characteristic series in the present embodimentSchematic diagram.
As shown in Figure 1, this method can be also obtained in carbonate formation to be analyzed in step s 106 at each depth point
Overburden pressure characteristic, and in step s 107 according to overburden pressure characteristic, pore-fluid volume elasticity
Modulus characteristic, saturated rock bulk modulus, rock solid bulk modulus characteristic, rock matrix volume bullet
Property modulus characteristic and porosity characteristics data determine the formation pore in carbonate formation to be analyzed at each depth point
Pressure characterization data.
Specifically, in the present embodiment, this method is preferably by accessed in step s 103 in step s 106
Density feature data series { h(i),ρ(i)Calculate the overburden pressure characteristic series { h at each depth point(i),P0 (i)}.Specifically, in the present embodiment, this method calculates the overburden pressure feature at each depth point according to following expression
Data P0:
P0=0.00981 ρ h (9)
For the overburden pressure characteristic P at the i-th depth point0 (i)For, exist:
P0 (i)=0.00981 ρ(i)h(i) (10)
Figure 10 shows overburden pressure characteristic series { h in the present embodiment(i),P0 (i)Schematic diagram.
In the present embodiment, this method can be in step s 107 by overburden pressure characteristic obtained in the above process
According to, pore-fluid bulk modulus characteristic, saturated rock bulk modulus, rock solid bulk modulus feature
Data, rock matrix bulk modulus characteristic and porosity characteristics data substitute into default carbonate formation Pore Pressure
In power computation model, to finally obtain the formation pore pressure characteristic in carbonate formation to be analyzed at each depth point
According to.
In the present embodiment, this method is used in step s 107 to preset Pore Pressure on Carbonate Rock Formation computation model
It is to be obtained based on poroelasticity theory of mechanics and in conjunction with carbonate rock rock matrix feature construction, it preferably can be using such as
Lower expression formula is indicated:
Wherein, PpIndicate formation pore pressure characteristic.
For the formation pore pressure characteristic at the i-th depth pointFor, it can be calculated using following expression
It obtains:
Figure 11 shows the formation pore pressure characteristic of depth segment locating for carbonate formation to be analyzed in the present embodiment
According to seriesSchematic diagram.
As can be seen that method provided by the present embodiment utilizes the carbonate formation Pore Pressure constructed from foregoing description
Power computation model determines the formation pore pressure characteristic in carbonate formation to be analyzed at each depth point.Wherein, institute
The Pore Pressure on Carbonate Rock Formation computation model of building considers the feature of carbonate rock rock physics structure, analyzes carbonic acid
Influence of the rock salt petrophysical parameter to formation pore pressure.This method is established from the angle of poroelasticity theory of mechanics
Pore Pressure on Carbonate Rock Formation computation model be overburden pressure, porosity, pore-fluid bulk modulus, rock
Function between skeleton bulk modulus and rock solid bulk modulus.
This method establishes formation pore pressure model using pore-fluid acoustic speed, more accurately can calculate ground
Layer pore pressure, for Drilling Design, selects and determines suitable safe drilling fluid density and provides scientific basis, to effectively prevent
The generation of underground complex accident during drilling construction.
It should be understood that disclosed embodiment of this invention is not limited to specific structure disclosed herein or processing step
Suddenly, the equivalent substitute for these features that those of ordinary skill in the related art are understood should be extended to.It should also be understood that
It is that term as used herein is used only for the purpose of describing specific embodiments, and is not intended to limit.
" one embodiment " or " embodiment " mentioned in specification means the special characteristic described in conjunction with the embodiments, structure
Or characteristic is included at least one embodiment of the present invention.Therefore, the phrase " reality that specification various places throughout occurs
Apply example " or " embodiment " the same embodiment might not be referred both to.
Although above-mentioned example is used to illustrate principle of the present invention in one or more application, for the technology of this field
For personnel, without departing from the principles and ideas of the present invention, hence it is evident that can in form, the details of usage and implementation
It is upper that various modifications may be made and does not have to make the creative labor.Therefore, the present invention is defined by the appended claims.
Claims (9)
1. a kind of Pore Pressure on Carbonate Rock Formation acquisition methods, which is characterized in that the described method includes:
Step 1: the rock forming mineral component characteristics data in carbonate formation to be analyzed at each depth point are obtained, according to described
Rock forming mineral component characteristics data determine the rock solid volume elasticity in the carbonate formation to be analyzed at each depth point
Modulus characteristic;
Step 2: the velocity of longitudinal wave characteristic of the carbonate formation to be analyzed is obtained, according to the velocity of longitudinal wave feature
Data calculate the rock matrix bulk modulus characteristic in the carbonate formation to be analyzed at each depth point;
Step 3: according to the saturated rock volume elasticity mould in the carbonate formation to be analyzed got at each depth point
Measure feature data, porosity characteristics data, in conjunction with the rock solid bulk modulus characteristic and rock matrix body
Product elasticity modulus characteristic, calculates the pore-fluid volume elasticity mould in the carbonate formation to be analyzed at each depth point
Measure feature data;
Step 4: the overburden pressure characteristic in carbonate formation to be analyzed at each depth point is obtained, according to described
Overburden pressure characteristic, pore-fluid bulk modulus characteristic, saturated rock bulk modulus characteristic
According to, rock solid bulk modulus characteristic, rock matrix bulk modulus characteristic and porosity characteristics number
According to the formation pore pressure characteristic at each depth point in the determination carbonate formation to be analyzed.
2. the method as described in claim 1, which is characterized in that in the step 4, calculated according to following expression described in
Formation pore pressure characteristic in carbonate formation to be analyzed at each depth point:
Wherein, PpIndicate formation pore pressure characteristic, KfIndicate pore-fluid bulk modulus characteristic, KfrIt indicates
Rock matrix bulk modulus characteristic, φ indicate porosity, KsIndicate rock solid bulk modulus characteristic,
P0Indicate overburden pressure characteristic.
3. method according to claim 1 or 2, which is characterized in that in said step 1, using VRH mean value method come root
The rock solid body in the carbonate formation to be analyzed at each depth point is determined according to the rock forming mineral component characteristics data
Product elasticity modulus characteristic.
4. method according to any one of claims 1 to 3, which is characterized in that in the step 2, to velocity of longitudinal wave spy
Sign data are decomposed, and the skeleton velocity of longitudinal wave of the formation rock in the carbonate formation to be analyzed at each depth point is obtained
Characteristic calculates in the carbonate formation to be analyzed at each depth point according to the skeleton velocity of longitudinal wave characteristic
Rock matrix bulk modulus characteristic.
5. method as claimed in claim 4, which is characterized in that in the step 2, calculate each depth according to following expression
Rock matrix bulk modulus characteristic at degree point:
Wherein, KfrIndicate rock matrix bulk modulus characteristic, φ indicates porosity, ρsIndicate that rock solid density is special
Levy data, νp,frIndicate the skeleton velocity of longitudinal wave characteristic of formation rock, νsIndicate the shear wave velocity characteristic of formation rock
According to.
6. such as method according to any one of claims 1 to 5, which is characterized in that in the step 3, according to what is got
Velocity of longitudinal wave characteristic, shear wave velocity characteristic and density in the carbonate formation to be analyzed at each depth point
Characteristic determines the saturated rock bulk modulus characteristic at each depth point.
7. method as claimed in claim 6, which is characterized in that calculate the saturated rock at each depth point according to following expression
Bulk modulus characteristic:
Wherein, KsatIndicate saturated rock bulk modulus characteristic, ρ indicates density feature data, νpIndicate velocity of longitudinal wave
Characteristic, νsIndicate shear wave velocity characteristic.
8. such as method according to any one of claims 1 to 7, which is characterized in that calculate each depth point according to following expression
The pore-fluid bulk modulus characteristic at place:
Wherein, KsatIndicate saturated rock bulk modulus characteristic, KfIndicate pore-fluid bulk modulus characteristic
According to KfrIndicate rock matrix bulk modulus characteristic, φ indicates porosity, KsIndicate rock solid bulk modulus
Characteristic.
9. such as method according to any one of claims 1 to 8, which is characterized in that in the step 4, obtain described wait divide
Analyse density feature data, the depth characteristic data in carbonate formation at each depth point, according to the density feature data and
Depth characteristic data determine the overburden pressure characteristic at each depth point respectively.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710761717.1A CN109425904B (en) | 2017-08-30 | 2017-08-30 | Method for acquiring pore pressure of carbonate rock stratum |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710761717.1A CN109425904B (en) | 2017-08-30 | 2017-08-30 | Method for acquiring pore pressure of carbonate rock stratum |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109425904A true CN109425904A (en) | 2019-03-05 |
CN109425904B CN109425904B (en) | 2020-07-31 |
Family
ID=65504022
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710761717.1A Active CN109425904B (en) | 2017-08-30 | 2017-08-30 | Method for acquiring pore pressure of carbonate rock stratum |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109425904B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101963056A (en) * | 2010-08-19 | 2011-02-02 | 中国石油大学(北京) | Method for predicting carbonate formation pore pressure by using log information |
CN105484739A (en) * | 2015-11-26 | 2016-04-13 | 中国科学院武汉岩土力学研究所 | Carbonate rock formation pore pressure testing method and device |
WO2017074884A1 (en) * | 2015-10-28 | 2017-05-04 | Schlumberger Technology Corporation | Formation evaluation |
-
2017
- 2017-08-30 CN CN201710761717.1A patent/CN109425904B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101963056A (en) * | 2010-08-19 | 2011-02-02 | 中国石油大学(北京) | Method for predicting carbonate formation pore pressure by using log information |
WO2017074884A1 (en) * | 2015-10-28 | 2017-05-04 | Schlumberger Technology Corporation | Formation evaluation |
CN105484739A (en) * | 2015-11-26 | 2016-04-13 | 中国科学院武汉岩土力学研究所 | Carbonate rock formation pore pressure testing method and device |
Non-Patent Citations (2)
Title |
---|
余夫等: "异常高压地层的纵波速度响应特征分析", 《石油钻探技术》 * |
路保平等: "基于流体声速的碳酸盐岩地层孔隙压力求取方法", 《石油钻探技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109425904B (en) | 2020-07-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106855636B (en) | The prototype geological model Seismic forward method appeared based on carbonate reservoir | |
CN107290782B (en) | Reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method | |
CN101414013B (en) | Method for determining underground fluid by seismic data | |
US8898046B2 (en) | Method to improve reservoir simulation and recovery from fractured reservoirs | |
CN107203005A (en) | A kind of method that quantification calculates crack characterising parameter | |
CN102788994A (en) | Reservoir fracture determining method | |
CN104181585A (en) | Shear wave estimation method and system in geophysical exploration | |
CN106772673B (en) | A kind of shale reservoir crustal stress Forecasting Model System | |
CN109655903B (en) | Shale layer transverse wave velocity prediction method and system | |
CN106291682B (en) | A kind of poststack sound impedance inversion method based on base method for tracing | |
CN103914620B (en) | Method for computing distribution of opening spaces of fractures of fault rupture zones | |
CN104749617A (en) | Multi-scale fractured reservoir forward model establishing method | |
CN104570127B (en) | A kind of method of utilization seimic wave velocity Simultaneous Inversion porosity and shale content | |
CN104316958B (en) | A kind of Coherent processing method for identifying different scale stratum breaking | |
CN105089615A (en) | Log data historical retrogression treatment method based on oil reservoir model | |
CN105588883B (en) | Three-dimensional rock mechanics parameters acquisition methods and system | |
Shi et al. | Brittleness index prediction in shale gas reservoirs based on efficient network models | |
CN107728205B (en) | A kind of Formation pressure prediction method | |
CN104570065B (en) | Method for quantitatively inverting porosity by using seismic wave impedance | |
CN103713327A (en) | Shale gas reservoir stratum minimum closing pressure evaluation method based on logging and seismic data | |
CN110095811B (en) | Method and device for constructing and processing velocity model of paste rock stratum | |
CN106501872B (en) | A kind of calculation method and device of fracture reservoir ground stress characteristics | |
CN106353807A (en) | Fracture identification method and device | |
CN103628866B (en) | Method for obtaining stratum pressure coefficients of stratum | |
CN106842316B (en) | Crack determines method and apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |