CN109916716B - Method for quantitatively evaluating pressurization of pore fluid by pressure stress - Google Patents
Method for quantitatively evaluating pressurization of pore fluid by pressure stress Download PDFInfo
- Publication number
- CN109916716B CN109916716B CN201910141702.4A CN201910141702A CN109916716B CN 109916716 B CN109916716 B CN 109916716B CN 201910141702 A CN201910141702 A CN 201910141702A CN 109916716 B CN109916716 B CN 109916716B
- Authority
- CN
- China
- Prior art keywords
- pore fluid
- stratum
- maximum deformation
- maximum
- fluid pressure
- 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.)
- Expired - Fee Related
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 81
- 239000011148 porous material Substances 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000009933 burial Methods 0.000 claims abstract description 18
- 230000009467 reduction Effects 0.000 claims abstract description 17
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 5
- 239000011435 rock Substances 0.000 claims description 26
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000006731 degradation reaction Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 230000003628 erosive effect Effects 0.000 claims 3
- 238000002679 ablation Methods 0.000 claims 1
- 230000008961 swelling Effects 0.000 abstract 2
- 238000005260 corrosion Methods 0.000 abstract 1
- 230000007797 corrosion Effects 0.000 abstract 1
- 230000037303 wrinkles Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
The invention provides a method for quantitatively evaluating pressurization of pore fluid by pressure stress, which comprises the following steps: calculating the maximum burial depth period pore fluid pressure PmbRecovering the buried depth of stratum in maximum deformation period (the stratum deforms maximally but does not break), recovering the buried depth of stratum in maximum deformation period after the maximum deformation period, forming fault and starting to slide along the fault plane, continuously lifting the stratum, obtaining the vertical fault distance of the current fault, recovering the buried depth of stratum in maximum deformation period, calculating the pressure reduction amount △ P of pore fluid caused by fold, swelling and corrosion of stratum, obtaining the pressure P of pore fluid in maximum deformation periodmfThe method comprises the steps of quantitatively calculating the pressurization P of the pore fluid by the pressure stress, and subtracting the pore fluid pressure in the maximum burial depth period from the sum of the pore fluid pressure in the maximum deformation period and the pore fluid pressure reduction caused by the formation wrinkle swelling and denudation to obtain the pressurization of the pore fluid by the pressure stress.
Description
Technical Field
The invention belongs to the technical field of oil-gas exploration fluid pressure prediction, and particularly relates to a method for quantitatively evaluating pressurization of pore fluid by pressure stress.
Technical Field
The overpressure is closely related to oil and gas, the phenomenon that overpressure is generated by tectonic pressure stress generally exists in the world, and quantitative calculation of the pressure stress has important significance for pore fluid pressurization for finding oil and gas reservoirs.
The existing method for quantitatively calculating the pressurization of the pressure stress fluid mainly comprises the following steps: a method for quantitatively calculating the maximum value of the pressure stress pressurization based on a totally closed system plane stress field (Terzaghi, 1943; Wang et al, 2005); a method for quantitatively calculating pressure stress pressurization based on a semi-closed system plane stress field (Wang et al, 2005); methods for quantitatively calculating compressive stress pressurization based on three-dimensional stress fields (Luo, 2004; Obradors-Prats et al, 2016,2017). The following problems still remain: (1) the plane stress field cannot describe the actual stress state of the rock, and the pressure stress pressurization cannot be accurately evaluated; (2) in a semi-closed system, it is difficult to obtain an accurate partition coefficient due to anisotropy of the sealing ability of the capping layer and the lack of related data; (3) the practicability of the numerical simulation of the pressure stress pressurization in the three-dimensional stress field is poor, and a quantitative calculation model capable of directly calculating the pressure stress pressurization is lacked at present.
Disclosure of Invention
The invention provides a method for quantitatively evaluating pressurization of a pore fluid by pressure stress, which solves the problem of quantitatively calculating the pressurization of the pore fluid by the pressure stress.
The technical scheme of the invention is as follows: a method for quantitatively evaluating pressurization of a pore fluid by a compressive stress, comprising the steps of:
(1) calculating the maximum burial depth period pore fluid pressure PmbThe method specifically comprises the following steps: calculating the maximum burial depth stage pore fluid pressure by using an equivalent depth method, and correcting the pore fluid pressure by using the maximum burial depth stage fluid pressure measured by the regional fluid inclusion, so as to finally obtain the more accurate maximum burial depth stage pore fluid pressure;
(2) recovering the buried depth of the stratum in the maximum deformation period (the stratum deforms maximally but is not broken), breaking the stratum after the maximum deformation period to form a fault and start sliding along the fault surface, continuously lifting the stratum, solving the vertical fault distance of the current fault, and recovering the buried depth of the stratum in the maximum deformation period;
(3) calculating the reduction △ P of pore fluid pressure caused by the collapse and the denudation of the stratum, the overlying load reduction, the stratum temperature reduction, the rock pore rebound and the fluid expansion, the pore fluid pressure reduction, and calculating the reduction of the pore fluid pressure by using a formula:
△ h for denudation thickness, △ T for temperature change during build lift, CpIs the coefficient of pore resilience, ρrDensity of the rock to be degraded αfCoefficient of thermal expansion of pore fluid, βfIs the pore fluid compression coefficient, v is the Poisson's ratio, g is the gravitational acceleration;
(4) obtaining the maximum deformation period pore fluid pressure PmfAfter the maximum deformation period, the rock stratum is broken and begins to slide along the fault plane, and the rock stratum is weakly compacted, so that the rock stratum density is similar to that of the existing rock stratum; obtaining the overburden load in the maximum deformation period by recovering the buried depth of the rock stratum in the maximum deformation period; the stress can be divided into normal stress and shear stress, the normal stress causes the rock stratum to generate volume deformation, the shear stress causes the rock stratum to generate shape deformation, and the maximum deformation period horizontal maximum effective stress sigmaHHorizontal minimum effective stress sigmahAnd vertical effective stress sigmavThe functional relationship of the three is as follows: sigmav＝σh/ν-σHV is Poisson's ratio, the horizontal maximum effective stress and the horizontal minimum effective stress of the target layer in the maximum deformation period can be obtained by using an acoustic emission experiment, the vertical effective stress of the maximum deformation period is further obtained, and the pore fluid pressure in the maximum deformation period is obtained by subtracting the vertical effective stress from the overlying load in the maximum deformation period;
(5) quantitatively calculating the pressurizing P of the pore fluid by the pressure stress and the pore fluid pressure P in the maximum deformation periodmfThe maximum burial depth pore fluid pressure P subtracted from the sum of the pore fluid pressure reductions △ P caused by formation fold swell degradationmbObtaining the pressurization of the pore fluid by the pressure stress, namely the formula:
P＝Pmf+Δp-Pmb。
the invention has the advantages that: the invention provides a quantitative evaluation model of pressurizing pore fluid by pressure stress under a three-dimensional stress field semi-closed system, which solves the defects of the conventional quantitative calculation pressure stress pressurizing method, has the advantages of higher feasibility and practicability, better accordance with geological practice and the like, and solves the problem of pressurizing pore fluid by quantitatively calculating pressure stress.
Drawings
FIG. 1 is a schematic flow diagram of the present invention.
Detailed Description
The present invention will now be described more fully hereinafter with reference to the accompanying drawings, in which a person skilled in the art can, without any creative effort, fully implement the present invention.
The specific implementation mode of the invention is as follows: as shown in fig. 1, a method for quantitatively evaluating the pressurization of a pore fluid by a compressive stress comprises the following steps:
(1) calculating the maximum burial depth period pore fluid pressure PmbThe method specifically comprises the following steps: calculating the maximum burial depth stage pore fluid pressure by using an equivalent depth method, and correcting the pore fluid pressure by using the maximum burial depth stage fluid pressure measured by the regional fluid inclusion, so as to finally obtain the more accurate maximum burial depth stage pore fluid pressure;
(2) recovering the buried depth of the stratum in the maximum deformation period (the stratum deforms maximally but is not broken), breaking the stratum after the maximum deformation period to form a fault and start sliding along the fault surface, continuously lifting the stratum, solving the vertical fault distance of the current fault, and recovering the buried depth of the stratum in the maximum deformation period;
(3) calculating the reduction △ P of pore fluid pressure caused by the collapse and the denudation of the stratum, the overlying load reduction, the stratum temperature reduction, the rock pore rebound and the fluid expansion, the pore fluid pressure reduction, and calculating the reduction of the pore fluid pressure by using a formula:
△ h for denudation thickness, △ T for temperature change during build lift, CpIs the coefficient of pore resilience, ρrDensity of the rock to be degraded αfCoefficient of thermal expansion of pore fluid, βfIs the pore fluid compression coefficient, v is the Poisson's ratio, g is the gravitational acceleration;
(4) obtaining the maximum deformation period pore fluid pressure PmfAfter the maximum deformation period, the rock stratum is broken and begins to slide along the fault plane, and the rock stratum is weakly compacted, so that the rock stratum density is similar to that of the existing rock stratum; obtaining the overburden load in the maximum deformation period by recovering the buried depth of the rock stratum in the maximum deformation period; stress separableThe normal stress causes the rock stratum to generate volume deformation, the shear stress causes the rock stratum to generate shape deformation, and the maximum deformation period horizontal maximum effective stress sigma is normal stress and shear stressHHorizontal minimum effective stress sigmahAnd vertical effective stress sigmavThe functional relationship of the three is as follows: sigmav＝σh/ν-σHV is Poisson's ratio, the horizontal maximum effective stress and the horizontal minimum effective stress of the target layer in the maximum deformation period can be obtained by using an acoustic emission experiment, the vertical effective stress of the maximum deformation period is further obtained, and the pore fluid pressure in the maximum deformation period is obtained by subtracting the vertical effective stress from the overlying load in the maximum deformation period;
(5) quantitatively calculating the pressurizing P of the pore fluid by the pressure stress and the pore fluid pressure P in the maximum deformation periodmfThe maximum burial depth pore fluid pressure P subtracted from the sum of the pore fluid pressure reductions △ P caused by formation fold swell degradationmbObtaining the pressurization of the pore fluid by the pressure stress, namely the formula:
P＝Pmf+Δp-Pmb。
while the preferred embodiments of the invention have been described, it is to be understood that the invention is not limited to the precise embodiments described, and that equipment and structures not described in detail are understood to be practiced as commonly known in the art; any simple modification, equivalent change and modification of the above embodiments according to the technical essence of the present invention by those skilled in the art can be made without departing from the technical scope of the present invention, and still fall within the protection scope of the technical solution of the present invention.
Claims (3)
1. A method for quantitatively evaluating pressurization of a pore fluid by a compressive stress, comprising the steps of:
(1) calculating the maximum burial depth period pore fluid pressure PmbThe method specifically comprises the following steps: calculating the maximum burial depth stage pore fluid pressure by using an equivalent depth method, and correcting the pore fluid pressure by using the maximum burial depth stage fluid pressure measured by the regional fluid inclusion, so as to finally obtain the more accurate maximum burial depth stage pore fluid pressure;
(2) recovering the buried depth of the stratum in the maximum deformation period, breaking the stratum after the maximum deformation period to form a fault and start sliding along the fault surface, continuously lifting the stratum, solving the vertical fault distance of the current fault and recovering the buried depth of the stratum in the maximum deformation period;
(3) the formula is used for calculating the reduction △ P of the pore fluid pressure caused by the collapse and the erosion of the stratum fold bulge, the stratum fold bulge and the erosion, the overburden load is reduced, the stratum temperature is reduced to cause the rebound of rock pores and the expansion of fluid, the pore fluid pressure is reduced, and the formula for calculating the reduction △ P of the pore fluid pressure caused by the collapse and the erosion of the stratum fold bulge is as follows:
wherein △ h is ablation thickness, △ T is temperature change in the structure lifting process, CpIs the coefficient of pore resilience, ρrDensity of the rock to be degraded αfCoefficient of thermal expansion of pore fluid, βfIs the pore fluid compression coefficient, v is the Poisson's ratio, g is the gravitational acceleration;
(4) obtaining the maximum deformation period pore fluid pressure PmfAfter the maximum deformation period, the rock stratum is broken and begins to slide along the fault plane, the rock stratum is weakly compacted, so that the density of the rock stratum is similar to that of the existing rock stratum, the overburden load in the maximum deformation period is obtained by recovering the burial depth of the rock stratum in the maximum deformation period, and the pore fluid pressure in the maximum deformation period is obtained by subtracting the vertical effective stress in the maximum deformation period from the overburden load in the maximum deformation period;
(5) quantitatively calculating the pressurizing P of the pore fluid by the pressure stress and the pore fluid pressure P in the maximum deformation periodmfThe maximum burial depth pore fluid pressure P subtracted from the sum of the pore fluid pressure reductions △ P caused by formation fold swell degradationmbThe pressurization of the pore fluid by the pressure stress is obtained, and the formula is as follows:
P＝Pmf+ΔP-Pmb。
2. the method for quantitatively evaluating the pressurization of pore fluid by compressive stress according to claim 1, wherein the maximum deformation period in the step (2) is a period in which the formation is deformed maximally but not fractured.
3. The method for quantitatively evaluating the pressurization of pore fluid by compressive stress according to claim 2, wherein the vertical effective stress in the maximum deformation period is σvAnd maximum effective stress sigma of maximum deformation period levelHHorizontal minimum effective stress sigmahThe functional relationship of (A) is as follows: sigmav＝σh/ν-σHAnd v is Poisson's ratio, the horizontal maximum effective stress and the horizontal minimum effective stress of the target layer in the maximum deformation period can be obtained by utilizing an acoustic emission experiment, and then the vertical effective stress of the maximum deformation period is sigmav。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910141702.4A CN109916716B (en) | 2019-02-26 | 2019-02-26 | Method for quantitatively evaluating pressurization of pore fluid by pressure stress |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910141702.4A CN109916716B (en) | 2019-02-26 | 2019-02-26 | Method for quantitatively evaluating pressurization of pore fluid by pressure stress |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109916716A CN109916716A (en) | 2019-06-21 |
| CN109916716B true CN109916716B (en) | 2020-04-17 |
Family
ID=66962230
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910141702.4A Expired - Fee Related CN109916716B (en) | 2019-02-26 | 2019-02-26 | Method for quantitatively evaluating pressurization of pore fluid by pressure stress |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109916716B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105468886A (en) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | Method for calculating formation pressure based on petrophysical parameters |
| CN106368686A (en) * | 2015-07-17 | 2017-02-01 | 中国石油化工股份有限公司 | Method for computing stratum pressure based on rock Poisson's ratio |
| CN107179393A (en) * | 2017-06-19 | 2017-09-19 | 中国科学院武汉岩土力学研究所 | A kind of pore water pressure force test method of the fine and close rock soil medium of low water content |
| WO2018111257A1 (en) * | 2016-12-14 | 2018-06-21 | Halliburton Energy Services, Inc. | Hydraulic fracturing methods and systems using gas mixture |
-
2019
- 2019-02-26 CN CN201910141702.4A patent/CN109916716B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105468886A (en) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | Method for calculating formation pressure based on petrophysical parameters |
| CN106368686A (en) * | 2015-07-17 | 2017-02-01 | 中国石油化工股份有限公司 | Method for computing stratum pressure based on rock Poisson's ratio |
| WO2018111257A1 (en) * | 2016-12-14 | 2018-06-21 | Halliburton Energy Services, Inc. | Hydraulic fracturing methods and systems using gas mixture |
| CN107179393A (en) * | 2017-06-19 | 2017-09-19 | 中国科学院武汉岩土力学研究所 | A kind of pore water pressure force test method of the fine and close rock soil medium of low water content |
Non-Patent Citations (2)
| Title |
|---|
| 库车坳陷构造挤压增压定量评价新方法;张凤奇等;《中国石油大学学报( 自然科学版)》;20110831;1-7 * |
| 异常孔隙流体压力的成因及其贡献探讨;付广等;《海相油气地质》;19991231;46-50 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109916716A (en) | 2019-06-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20200386081A1 (en) | Optimized design method for temporary blocking agent to promote uniform expansion of fractures produced by fracturing in horizontal wells | |
| CN109324345B (en) | Method for recovering porosity of rock in oil and gas accumulation period of superimposed basin | |
| CN109162701B (en) | Coal seam open hole well fracture pressure prediction method | |
| CN112698392A (en) | Multilayer system reservoir three-dimensional drawing and oil-gas spatial distribution and constant volume method | |
| CN110348031B (en) | Numerical simulation method for horizontal well fracturing near-wellbore fracture distortion form | |
| CN111157346B (en) | Analysis system and method for design and correction of water curtain system of underground water-sealed oil storage | |
| CN113790044A (en) | Method for evaluating pore pressure of multi-mechanism composite overpressure formation | |
| CN106289964A (en) | A kind of mud shale Gu comprcssive strength determines method | |
| CN117291059A (en) | Method for calculating recoverable reserves and water invasion parameters of high-temperature high-pressure water-bearing reservoir | |
| CN109916716B (en) | Method for quantitatively evaluating pressurization of pore fluid by pressure stress | |
| CN113919201A (en) | Multi-scale expansion grid self-adaption method for hydraulic fracturing fracture | |
| Park | Geomechanical Simulation of Bayou Choctaw Strategic Petroleum Reserve-Model Calibration. | |
| US10794179B2 (en) | Poromechanical impact on yield behavior in unconventional reservoirs | |
| CN110850470A (en) | Method for calculating ground stress difference coefficient by using seismic data | |
| CN106842329B (en) | The acquisition methods and device of reservoir oil saturation | |
| CN107478768B (en) | A kind of mud shale Gu shearing strength determines method | |
| CN113944459B (en) | Method and device for predicting resource amount of shale in-situ conversion recoverable oil gas | |
| CN107016219B (en) | Early warning method and system for carbonate reservoir drilling emptying | |
| Isotton et al. | On the possible contribution of clayey inter-layers to delayed land subsidence above producing aquifers | |
| CN112651190B (en) | Method for representing dual-medium seepage flow rate through fracture density | |
| CN114169204B (en) | Sand control opportunity determination method for offshore oil and gas field development and production | |
| CN113255037B (en) | Method for estimating floating amount of double-mode shield tunnel segment in upper-soft and lower-hard stratum | |
| Gambolati et al. | Numerical modeling of natural land subsidence over sedimentary basins undergoing large compaction | |
| Park | Assessing Potential Risk of Loss in Integrity of the Sidewall of a Cavern Located Near the Edge of the Salt Dome | |
| EP3658750B1 (en) | Drained reservoir volume diagnostics from mandel-cryer pressure signal |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200417 |