CN109901238B - High-stress formation resistivity correction method based on stress difference resistivity experiment - Google Patents
High-stress formation resistivity correction method based on stress difference resistivity experiment Download PDFInfo
- Publication number
- CN109901238B CN109901238B CN201910148795.3A CN201910148795A CN109901238B CN 109901238 B CN109901238 B CN 109901238B CN 201910148795 A CN201910148795 A CN 201910148795A CN 109901238 B CN109901238 B CN 109901238B
- Authority
- CN
- China
- Prior art keywords
- stress
- resistivity
- porosity
- experiment
- formation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a high-stress formation resistivity correction method based on a stress difference resistivity experiment, which comprises the following steps of: obtaining the formation resistivity and the maximum and minimum horizontal principal stress; obtaining the porosity of the rock core, and establishing a porosity calculation model by combining density logging data; selecting rock cores with different physical properties of a target layer to carry out a conventional rock-electricity experiment to obtain a cementation index m value; selecting rock cores with different physical properties of a target layer to carry out stress difference resistivity experiments, obtaining rock core resistivity measurement values under different stress differences, and establishing a stress difference corrected resistivity model; and correcting the resistivity model by utilizing the stress difference, and solving the formation resistivity without stress influence point by point through the porosity, the resistivity, the maximum and minimum horizontal main stress and the cementation index m value. In practical application, the method realizes resistivity abnormity under strong extrusion stress, and obtains the formation resistivity under stress influence elimination through logging information point by point, thereby providing basic parameters for accurate evaluation of the fluid property of the high-stress formation reservoir.
Description
Technical Field
The invention belongs to the field of processing and explaining well logging data of a clastic rock oil-gas reservoir in oil exploration, and particularly relates to a high-stress stratum resistivity correction method based on a stress difference resistivity experiment.
Background
Researchers at home and abroad develop a lot of experimental research works on the resistivity of rocks, and Brace and the like consider that the pore state and the volume expansion of the rocks are the main reasons for the change of the resistivity of the rocks in the process from loading to cracking. The Chengduan and the like research the change condition of the resistivity of the granite sample through repeated stress loading. The stress environment of stratum rock is simulated by utilizing 2-dimensional elastic constraint force, the change of the resistivity of the rock is measured by adopting a quadrupole method, the resistance of the rock is considered to show a descending trend in different stress stages along with the increase of loading stress, and the resistivity is reduced more quickly at the position close to the main fracture of the rock. And the like compares saturated crude oil samples with saturated water samples through simulation of reservoir environments to obtain that the rock resistivity increases along with the increase of pressure and changes in an exponential form of e, but the change gradient is small. The researches such as Zhang Ningsheng and the like find that for the core in the same block, when confining pressure is in the elastic deformation range of the core, the pressure and the resistivity of the core are in a linear relation. The influence of the magnitude and direction of the ground stress on the logging response of the reservoir resistivity of the reservoir-depressed Kerasu structure of the garage vehicle is researched based on logging information, and the results show that the ground stress has an exponential increasing trend along with the increase of the horizontal stress difference, the correlation between the formation resistivity and the horizontal main stress difference is good when the included angle between the horizontal maximum main stress direction and the fracture trend is small, and the correlation between the formation resistivity and the horizontal main stress difference is poor when the included angle between the horizontal maximum main stress direction and the fracture trend is large. The Wangwei and other researches analyze resistivity experiments under triaxial stress, and the core resistivity and the horizontal stress difference are in an exponential relationship. At present, most of research mainly focuses on the relationship between the ground stress and the resistivity, a stress difference correction resistivity model is not established all the time, the fluid property of a resistivity abnormal reservoir under strong extrusion stress cannot be judged, and the oil saturation of the oil field reservoir cannot be accurately evaluated.
Disclosure of Invention
The invention aims to provide a method for correcting the resistivity of a high-stress stratum based on a stress difference resistivity experiment, so as to solve the problems.
In order to achieve the purpose, the invention adopts the following technical scheme:
a high stress stratum resistivity correction method based on a stress difference resistivity experiment comprises the following steps:
Selecting rock cores with different grain diameters of a target layer to carry out a physical property analysis experiment, obtaining the porosity of the rock core, and establishing a porosity calculation model by combining density logging data;
selecting rock cores with different physical properties of a target layer to carry out a conventional rock-electricity experiment to obtain a cementation index m value;
selecting rock cores with different physical properties of a target layer to carry out stress difference resistivity experiments, obtaining rock core resistivity measurement values under different stress differences, and establishing a stress difference corrected resistivity model;
and 5, correcting the resistivity model by utilizing the stress difference, and solving the formation resistivity without stress influence point by point through the porosity, the resistivity, the maximum and minimum horizontal main stress and the cementation index m value.
Further, in the step 2, the core physical property experiment is carried out according to a flow specified by a standard of a core analysis method SY/T5336-2006; the porosity calculation model is as follows:
POR＝a×DEN+c (1)
in the formula: POR is porosity; DEN is density logging; a. and c is a coefficient in the model formula, and is obtained by performing least square fitting on physical property analysis experimental data and density logging data.
Further, in the step 3, the core rock electricity experiment is carried out according to the flow specified in the standard of core analysis method SY/T5385-2007.
Further, in step 4, the stress difference correction resistivity model is as follows:
ΔSH12＝SH1-SH2(4)
in the formula: RD is the formation resistivity for eliminating stress influence; RT is formation measurement resistivity under the influence of stress; Δ SH12Is the maximum minimum level principal stress difference; SH (hydrogen sulfide)1Is the maximum horizontal principal stress; SH (hydrogen sulfide)2Is the minimum horizontal principal stress; POR is porosity; m is a cementation index; n is a stress response index, b is a stress constant factor, and the stress response index and the stress constant factor are obtained by fitting stress difference resistivity experimental data by a least square method.
Further, in step 5, the stress difference is used for correcting a resistivity model, namely formula (3), and the formation resistivity for eliminating the stress influence is obtained point by point through the porosity, the resistivity, the maximum and minimum horizontal principal stress and the cementation index m value.
Compared with the prior art, the invention has the following technical effects:
the invention realizes the resistivity abnormity under strong extrusion stress, obtains the formation resistivity under the stress influence elimination through logging information point by point, provides accurate resistivity parameters for the fluid property evaluation of a high-stress formation reservoir, avoids the fluid property judgment error caused by the stress influence on the resistivity, and provides necessary reservoir parameters for determining the oil saturation of the reservoir in oil field development.
Drawings
FIG. 1 is a flow chart of a high stress formation resistivity correction method based on a stress difference resistivity experiment provided by the invention
FIG. 2 is a diagram of a porosity calculation model provided in an embodiment of the present invention
FIG. 3 is a diagram of a determined cementation index in a rock-electricity experiment provided in an embodiment of the invention
FIG. 4 is a graph of stress response exponent and constant factor determination provided in an embodiment of the present invention
FIG. 5 is a graph illustrating the effect of correcting formation resistivity to eliminate stress effects provided in embodiments of the present invention.
Detailed Description
The following provides a more detailed description of the embodiments of the present invention, with reference to the accompanying drawings.
Referring to fig. 1, a method for correcting resistivity of a high-stress formation based on a stress difference resistivity experiment provided by an embodiment of the invention includes the following steps:
step 101: selecting a well with density, array induction and dipole acoustic logging information in a research area to obtain formation resistivity and maximum and minimum horizontal principal stress;
step 102: selecting rock cores with different grain diameters of a target layer to carry out a physical property analysis experiment, obtaining the porosity of the rock core, and establishing a porosity calculation model by combining density logging data;
step 103: selecting rock cores with different physical properties of a target layer to carry out a conventional rock-electricity experiment to obtain a cementation index m value;
step 104: selecting rock cores with different physical properties of a target layer to carry out stress difference resistivity experiments, obtaining rock core resistivity measurement values under different stress differences, and establishing a stress difference corrected resistivity model;
step 105: and correcting the resistivity model by utilizing the stress difference, and solving the formation resistivity without stress influence point by point through the porosity, the resistivity, the maximum and minimum horizontal main stress and the cementation index m value.
In support of the technical problem to be solved by the present invention, the following steps are performed by further describing the specific implementation of the present embodiment:
the method comprises the steps of firstly, selecting a well with density, array induction and dipole acoustic logging information in a research area, and obtaining formation resistivity RT, maximum horizontal principal stress SH1 and minimum horizontal principal stress SH 2.
Selecting cores with different particle sizes of a target layer to carry out a physical property analysis experiment to obtain the porosity of the core, and combining density logging information to obtain a porosity calculation model (shown in figure 2) through least square fitting:
POR＝-0.4075×DEN+1.0898 (1)
in the formula: POR is porosity; DEN is density logging.
And step three, selecting rock cores with different physical properties of a target layer to carry out a conventional rock-electricity experiment, and obtaining the cementation index m value of 1.685 (shown in figure 3) through least square fitting.
Selecting cores with different physical properties of a target layer to carry out stress difference resistivity experiments to obtain measured values of the core resistivity under different stress differences, and applying a model formula (2) to obtain a stress response index n of 0.0439 and a stress constant factor b of 0.0187 through least square fitting (as shown in figure 4).
And fifthly, correcting the resistivity model by using the stress difference, namely formula (3), and calculating the formation resistivity for eliminating the stress influence point by point according to the porosity, the resistivity, the maximum and minimum horizontal main stress and the cementation index m value.
ΔSH12＝SH1-SH2(4)
In the formula: RD is the formation resistivity for eliminating stress influence; RT is formation measurement resistivity under the influence of stress; Δ SH12Is the maximum minimum level principal stress difference; SH (hydrogen sulfide)1Is the maximum horizontal principal stress; SH (hydrogen sulfide)2Is the minimum horizontal principal stress; POR is porosity.
The method is programmed to realize the processing modularization of the method, and the specific processing result is shown in fig. 5. As can be seen from the figure, by utilizing the method, the formation resistivity influenced by the stress can be effectively eliminated point by point through array induction logging and density logging, and accurate resistivity parameters are provided for the evaluation of the fluid properties of the high-stress formation reservoir.
A well with density, array induction, dipole sonic logging data is known, and the logging techniques known in the art collect density logging data for use in evaluating formation porosity. The acquisition array induction is used for evaluating the formation resistivity, and the dipole acoustic logging is used for evaluating the formation anisotropy, the ground stress condition and the like.
Finally, it should be noted that the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting, and although the present invention has been described in detail with reference to examples, it should be understood by those skilled in the art that modifications or equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, which should be covered by the claims of the present invention.
Claims (3)
1. A high stress stratum resistivity correction method based on a stress difference resistivity experiment is characterized by comprising the following steps:
step 1, selecting a well with density, array induction and dipole acoustic logging information in a research area to obtain formation resistivity RT and maximum horizontal principal stress SH1With minimum horizontal principal stress SH2；
Selecting rock cores with different grain diameters of a target layer to carry out a physical property analysis experiment, obtaining the porosity of the rock core, and establishing a porosity calculation model by combining density logging data;
selecting rock cores with different physical properties of a target layer to carry out a conventional rock-electricity experiment to obtain a cementation index m value;
selecting rock cores with different physical properties of a target layer to carry out stress difference resistivity experiments, obtaining rock core resistivity measurement values under different stress differences, and establishing a stress difference corrected resistivity model;
step 5, correcting the resistivity model by utilizing the stress difference, and solving the formation resistivity without stress influence point by point through the porosity, the resistivity, the maximum and minimum horizontal main stress and the cementation index m value;
in step 4, the stress difference correction resistivity model is as follows:
ΔSH12＝SH1-SH2(4)
in the formula: RD is the formation resistivity for eliminating stress influence; RT is formation measurement resistivity under the influence of stress; Δ SH12Is the maximum minimum level principal stress difference; SH (hydrogen sulfide)1Is the maximum horizontal principal stress; SH (hydrogen sulfide)2Is the minimum horizontal principal stress; POR is porosity; m is a cementation index; n is a stress response index, b is a stress constant factor, and the stress response index and the stress constant factor are obtained by fitting stress difference resistivity experimental data by a least square method;
in the step 5, the resistivity model, namely formula (3), is corrected by using the stress difference, and the formation resistivity without stress influence is obtained point by point through the porosity, the resistivity, the maximum and minimum horizontal principal stress and the cementation index m value.
2. The method for correcting the resistivity of the high-stress stratum based on the stress difference resistivity experiment as claimed in claim 1, wherein in the step 2, the core physical property experiment is performed according to a flow specified in the standard of core analysis method SY/T5336-2006; the porosity calculation model is as follows:
POR＝a×DEN+c (1)
in the formula: POR is porosity; DEN is density logging; a. and c is a coefficient in the model formula, and is obtained by performing least square fitting on physical property analysis experimental data and density logging data.
3. The method for correcting the resistivity of the high-stress stratum based on the stress difference resistivity experiment as claimed in claim 1, wherein in the step 3, the core rock electricity experiment is carried out according to a flow specified in a core analysis method SY/T5385-2007 standard.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910148795.3A CN109901238B (en) | 2019-02-28 | 2019-02-28 | High-stress formation resistivity correction method based on stress difference resistivity experiment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910148795.3A CN109901238B (en) | 2019-02-28 | 2019-02-28 | High-stress formation resistivity correction method based on stress difference resistivity experiment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109901238A CN109901238A (en) | 2019-06-18 |
| CN109901238B true CN109901238B (en) | 2020-09-04 |
Family
ID=66945909
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910148795.3A Active CN109901238B (en) | 2019-02-28 | 2019-02-28 | High-stress formation resistivity correction method based on stress difference resistivity experiment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109901238B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111255446B (en) * | 2020-01-19 | 2022-09-20 | 长江大学 | Resistivity correction method based on stratum simulation |
| CN111350490B (en) * | 2020-03-16 | 2023-11-28 | 中国石油天然气集团有限公司 | Stratum true parameter solving method based on logging data |
| CN114782276B (en) * | 2022-04-29 | 2023-04-11 | 电子科技大学 | Resistivity imaging dislocation correction method based on adaptive gradient projection |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101725344A (en) * | 2008-10-29 | 2010-06-09 | 中国石油天然气集团公司 | Method for determining litho-electric parameters |
| CN102012526A (en) * | 2010-09-09 | 2011-04-13 | 四川德阳西德电器有限公司 | Method for discriminating type of reservoir fluid by using resistivity data |
| CN102323198A (en) * | 2011-06-02 | 2012-01-18 | 中国石油天然气集团公司 | Method and system for correcting rock electricity experiment parameters m and n of core |
| CN103422853A (en) * | 2012-05-24 | 2013-12-04 | 中国石油天然气集团公司 | Method and device for determining water saturation of stratum |
| CN104101905A (en) * | 2013-04-11 | 2014-10-15 | 中国石油天然气集团公司 | Reservoir classification method based on rock electricity parameters |
| CN104453874A (en) * | 2014-10-23 | 2015-03-25 | 中国石油天然气集团公司 | Glutenite reservoir oil saturation calculation method based on nuclear magnetic resonance |
| CN106772674A (en) * | 2016-12-01 | 2017-05-31 | 中国石油天然气股份有限公司 | Method for calculating oil saturation of fractured sandstone reservoir |
| CN106802433A (en) * | 2016-12-28 | 2017-06-06 | 中国石油天然气股份有限公司 | Method and device for correcting formation resistivity |
| CN106842329A (en) * | 2016-12-22 | 2017-06-13 | 中国石油天然气股份有限公司 | Method and device for acquiring oil saturation of reservoir |
| CN107780923A (en) * | 2017-11-01 | 2018-03-09 | 中石化石油工程技术服务有限公司 | A kind of foundation of the water-saturation model based on Shale Correction, emulation mode |
| CN109025982A (en) * | 2018-07-16 | 2018-12-18 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | Classification evaluation method, device and the terminal device of shale gas exploitation interval |
| CN109061752A (en) * | 2018-06-26 | 2018-12-21 | 西南石油大学 | A kind of resistivity curve bearing calibration on the stratum containing grey matter |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9238962B2 (en) * | 2010-12-21 | 2016-01-19 | Schlumberger Technology Corporation | Pore pressure from spectroscopy and sonic data |
| US10473802B2 (en) * | 2017-05-10 | 2019-11-12 | Baker Hughes, A Ge Company, Llc | Downhole energy calibration of gamma ray detectors |
-
2019
- 2019-02-28 CN CN201910148795.3A patent/CN109901238B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101725344A (en) * | 2008-10-29 | 2010-06-09 | 中国石油天然气集团公司 | Method for determining litho-electric parameters |
| CN102012526A (en) * | 2010-09-09 | 2011-04-13 | 四川德阳西德电器有限公司 | Method for discriminating type of reservoir fluid by using resistivity data |
| CN102323198A (en) * | 2011-06-02 | 2012-01-18 | 中国石油天然气集团公司 | Method and system for correcting rock electricity experiment parameters m and n of core |
| CN103422853A (en) * | 2012-05-24 | 2013-12-04 | 中国石油天然气集团公司 | Method and device for determining water saturation of stratum |
| CN104101905A (en) * | 2013-04-11 | 2014-10-15 | 中国石油天然气集团公司 | Reservoir classification method based on rock electricity parameters |
| CN104453874A (en) * | 2014-10-23 | 2015-03-25 | 中国石油天然气集团公司 | Glutenite reservoir oil saturation calculation method based on nuclear magnetic resonance |
| CN106772674A (en) * | 2016-12-01 | 2017-05-31 | 中国石油天然气股份有限公司 | Method for calculating oil saturation of fractured sandstone reservoir |
| CN106842329A (en) * | 2016-12-22 | 2017-06-13 | 中国石油天然气股份有限公司 | Method and device for acquiring oil saturation of reservoir |
| CN106802433A (en) * | 2016-12-28 | 2017-06-06 | 中国石油天然气股份有限公司 | Method and device for correcting formation resistivity |
| CN107780923A (en) * | 2017-11-01 | 2018-03-09 | 中石化石油工程技术服务有限公司 | A kind of foundation of the water-saturation model based on Shale Correction, emulation mode |
| CN109061752A (en) * | 2018-06-26 | 2018-12-21 | 西南石油大学 | A kind of resistivity curve bearing calibration on the stratum containing grey matter |
| CN109025982A (en) * | 2018-07-16 | 2018-12-18 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | Classification evaluation method, device and the terminal device of shale gas exploitation interval |
Non-Patent Citations (2)
| Title |
|---|
| "库车前陆盆地强挤压应力条件下的测井电阻率校正方法";袁龙 等;《石油地球物理勘探》;20181031;第53卷(第5期);正文第1086页左栏第2段-1093页左栏第2段 * |
| "高温高压下地应力对电阻率影响实验研究";王伟 等;《测井技术》;20171231;第41卷(第6期);第642-647页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109901238A (en) | 2019-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2447469B1 (en) | Quantitative calculation method for hydrocarbon saturation of fractured reservoir | |
| CN102434152B (en) | Method for calculating oil saturation of reservoir | |
| CN105114064B (en) | Method for Determining Saturation of Tight Sandstone Reservoir | |
| CN104278991B (en) | Saline Lake Facies hydrocarbon source rock organic carbon and the polynary well logging computational methods of hydrocarbon potential | |
| CN109901238B (en) | High-stress formation resistivity correction method based on stress difference resistivity experiment | |
| CN105386756B (en) | A method of brittle formation porosity is calculated using dependent variable | |
| CN103926617B (en) | Seam hole reservoir body detection method and detection device | |
| CN105588883B (en) | Three-dimensional rock mechanical parameter acquisition method and system | |
| CN111897011B (en) | Reservoir pore characteristic determination method, device and equipment | |
| CN108019204B (en) | Neutron-gamma gas layer identification method based on physical property and formation water mineralization correction | |
| CN102606138B (en) | Method for correcting dielectric constant of electromagnetic wave resistivity logger while drilling by means of phase-amplitude method | |
| CN105447762A (en) | Calculation method for low permeability reservoir flooding information of fluid replacement | |
| CN106897531B (en) | Quantitative evaluation method for permeability of low-permeability limestone reservoir | |
| CN107829731B (en) | Clay alteration volcanic porosity correction method | |
| CN112698399A (en) | Gravel well seismic-logging linkage constraint efficient reservoir quantitative prediction method and system | |
| CN112145165B (en) | Microcrack-pore type reservoir dynamic and static permeability conversion method | |
| CN106285652B (en) | Method for determining shale free gas saturation | |
| CN107861162B (en) | Microelectrode logging data-based natural crack identification method and system | |
| CN108957554A (en) | Seismic inversion method in a kind of geophysical exploration | |
| CN112115592A (en) | Large-inclination/horizontal well resistivity correction method based on rock physics experiment | |
| CN109667576B (en) | High-salinity-formation-factor low-resistance oil layer logging identification method | |
| CN109885927B (en) | Continuous inversion method for stratum radial resistivity | |
| CN104514553A (en) | Recognition method of gas reservoir in lithic sandstone layer and application thereof | |
| CN109655394B (en) | Nuclear magnetic resonance T2 spectrum permeability calculation method under constraint of throat structure parameters | |
| CN112324421A (en) | Method for calculating saturation before and after flooding of low-resistivity heavy oil reservoir |
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 |