CN112525803A - Sandstone porosity effective stress coefficient calculation method based on clay shell model - Google Patents

Sandstone porosity effective stress coefficient calculation method based on clay shell model Download PDF

Info

Publication number
CN112525803A
CN112525803A CN202110181206.9A CN202110181206A CN112525803A CN 112525803 A CN112525803 A CN 112525803A CN 202110181206 A CN202110181206 A CN 202110181206A CN 112525803 A CN112525803 A CN 112525803A
Authority
CN
China
Prior art keywords
porosity
clay
pore
effective stress
stress coefficient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110181206.9A
Other languages
Chinese (zh)
Inventor
李溢龙
李晓平
谭晓华
孟展
卢虹林
刘蜀东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN202110181206.9A priority Critical patent/CN112525803A/en
Publication of CN112525803A publication Critical patent/CN112525803A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Abstract

The invention relates to a clay shell model-based sandstone porosity effective stress coefficient calculation method, belonging to the field of rock mechanics effective stress coefficient calculation; the method solves the problems that the existing method for calculating the effective porosity model by using a theoretical model is few, the precision is low and the like; the technical scheme is as follows: providing a clay shell-shaped model, obtaining a clay shell-shaped porosity effective stress coefficient expansion formula by a chain method according to a porosity effective stress coefficient definition formula and in combination with a clay shell-shaped physical model porosity definition formula, and obtaining a clay shell-shaped porosity effective stress coefficient equation; and substituting the porosity, the clay content, the Poisson ratio and the rigidity ratio into a clay shell-shaped porosity effective stress coefficient model, and calculating the porosity effective stress coefficient. Based on the analysis result of the scanning electron microscope, the clay shell model is provided, the clay shell porosity effective stress coefficient model is deduced, the porosity effective stress coefficient is accurately calculated, the calculation effect accuracy is high, and the popularization is strong.

Description

Sandstone porosity effective stress coefficient calculation method based on clay shell model
Technical Field
The invention relates to a clay shell model-based sandstone porosity effective stress coefficient calculation method, and belongs to the field of rock mechanics effective stress coefficient calculation.
Background
All oil and gas reservoir rocks belong to porous media, and the most difference of the porous media from general solid materials is that the porous media have pores and the pores contain fluid, and the pore fluid generates pore fluid pressure. Therefore, the concept of the effective stress law is often used in the research of the mechanical properties of the oil and gas reservoir rock and the physical properties related to the mechanical properties such as stress sensitivity. The research of the effective stress law is the basis for understanding the sensitivity degree of the reservoir stress from the mechanism, and a series of models representing the relationship between the porosity and the permeability and the effective stress are built by many scholars based on the effective stress law. The research on the effective stress rule is essentially the research on the effective stress coefficient, and is very important for determining the effective stress coefficient of porosity.
However, the existing research on the effective stress coefficient of porosity has few documents, and the patent number is CN201710403114.4, namely a method for determining the effective stress coefficient of rock, the method comprises the steps of drawing a relation curve of unloading bias voltage increment and axial/circumferential strain increment under different confining pressures, and performing data fitting on the curve to obtain the effective elastic modulus and Poisson ratio related to the bias voltage level; the axial effective stress coefficient and the annular effective stress coefficient of the rock sample are calculated through a formula, the effective stress coefficient obtained by the method is relatively reliable, but the test is required to be completed, and the cost is high. CN201410320028.3 'a method for determining effective stress coefficient of rock based on pore compression experiment' fits a fitting function of effective stress and pore compression coefficient, measures a pore pressure constant value of a saturated rock core and a pore pressure value after increasing confining pressure, determines the reduction of the pore volume of the rock core, and finally determines the effective stress coefficient value of the rock according to a calculation formula of the effective stress born by a rock framework and the pore compression coefficient of the rock. In 2017, panke in the thesis of master graduate, "research on porosity stress sensitivity of clay-containing tight sandstone", imagines the distribution forms of four kinds of clay, and supposes that not only the clay exists in the pores, but also a part of the clay exists outside the pores, and is directly acted by confining pressure, so that a porosity effective stress calculation model is deduced, but the model is not very consistent with the actual situation, the calculated result is greatly different from the actual real result, and the accuracy of the calculation result is in doubt.
Generally, most of the existing methods for calculating and analyzing the effective stress coefficient of porosity use physical experiments to calculate and measure, and few methods for calculating by using theoretical models have questionable calculation accuracy, so that more accurate calculation methods are needed.
Disclosure of Invention
The invention aims to: in order to solve the problems of few methods for calculating the effective porosity model by using a theoretical model, low precision and the like, the invention provides a clay shell model based on the analysis result of a scanning electron microscope, deduces the clay shell porosity effective stress coefficient model, accurately calculates the effective porosity stress coefficient, has high calculation effect accuracy and strong popularization.
In order to achieve the purpose, the invention provides a clay shell model-based sandstone porosity effective stress coefficient calculation method, which comprises the following steps:
s100, adhering chlorite and illite to the inner wall of a pore in a shell form based on a scanning electron microscope analysis result, providing a clay shell model, and determining the porosity, the clay content, the Poisson ratio and the rigidity ratio of a rock sample containing clay minerals;
s200, obtaining the clay shell porosity effective stress coefficient expansion by a chain method according to the porosity effective stress coefficient definition formula and combining the clay shell physical model porosity definition formula, the main steps are,
s201, determining a porosity effective stress coefficient definition formula and a clay shell physical model porosity definition formula, wherein the porosity effective stress coefficient definition formula is
Figure DEST_PATH_IMAGE001
In the formula
Figure 100002_DEST_PATH_IMAGE002
Is a porosity effective stress coefficient, and has no dimensional quantity;φporosity, dimensionless;P pis pore pressure in MPa;P cis confining pressure in MPa; the porosity of the clay shell-shaped physical model is defined as
Figure DEST_PATH_IMAGE003
In the formula:
Figure 100002_DEST_PATH_IMAGE004
is the pore inner radius, and the unit is mum; b is the outer radius of the pores in μm;
s202, obtaining a clay shell porosity effective stress coefficient expansion form by a chain method, wherein the clay shell porosity effective stress coefficient expansion form is
Figure DEST_PATH_IMAGE005
S300, solving a partial derivative relational expression of the displacement of the inner wall surface and the outer wall surface of the pore to the pore pressure and the confining pressure and a partial derivative relational expression of the porosity to the inner wall surface and the outer wall surface of the pore, mainly comprising the following steps,
s301, calculating a partial derivative relational expression of the displacement of the inner wall surface of the pore to the pore pressure, wherein the partial derivative relational expression of the displacement of the inner wall surface of the pore to the pore pressure is
Figure 100002_DEST_PATH_IMAGE006
In the formula
Figure DEST_PATH_IMAGE007
Is the partial derivative of the displacement of the inner wall surface of the pore to the pore pressure;v cis the Poisson ratio of clay minerals without dimension;μ cthe Lame coefficient of clay mineral is dimensionless;A 1pB 1pis a transition coefficient;
s302, calculating a partial derivative relation formula of the displacement of the inner wall surface of the pore to the confining pressure, wherein the partial derivative relation formula of the displacement of the inner wall surface of the pore to the confining pressure is
Figure 100002_DEST_PATH_IMAGE008
In the formula
Figure DEST_PATH_IMAGE009
The partial derivative of the displacement of the inner wall surface of the pore to the confining pressure;A 1cB 1cis a transition coefficient;
s303, solving a partial derivative relational expression of the displacement of the outer wall surface of the pore to the pore pressure, wherein the partial derivative relational expression of the displacement of the outer wall surface of the pore to the pore pressure is
Figure 100002_DEST_PATH_IMAGE010
In the formula
Figure DEST_PATH_IMAGE011
The partial derivative of the displacement of the outer wall surface of the pore to the pore pressure is obtained; b is the outer radius of the pores in μm;v rthe Poisson ratio of the rock skeleton is a dimensionless quantity;A 2pB 2pis a transition coefficient;
s304, solving a partial derivative relational expression of the displacement of the outer wall surface of the pore to the confining pressure, wherein the partial derivative relational expression of the displacement of the outer wall surface of the pore to the pore pressure is
Figure 100002_DEST_PATH_IMAGE012
In the formula
Figure DEST_PATH_IMAGE013
The partial derivative of the displacement of the outer wall surface of the pore to the confining pressure is obtained;μ rthe Lame coefficient of the rock skeleton is free of dimensional quantity;A 2cB 2cis a transition coefficient;
s305, calculating partial derivative relational expressions of the porosity on the inner wall surface and the outer wall surface of the pore, wherein the partial derivative relational expressions of the porosity on the inner wall surface and the outer wall surface of the pore are respectively
Figure 100002_DEST_PATH_IMAGE014
S400, defining a shell-shaped clay mineral content relational expression and a rigidity ratio of a rock framework to clay to obtain an effective stress coefficient equation of the clay shell-shaped porosity;
s401, defining the relation of the content of the shell-shaped clay minerals as
Figure DEST_PATH_IMAGE015
Defining the rigidity ratio of the rock skeleton to the clay as in the formulaF cThe content of the shell-shaped clay minerals is zero dimensional quantity; c is the radius of the inner wall surface of the rock framework, and the unit is mum; the rigidity ratio of the stone framework to the clay is free of dimensional quantity;
s402, substituting the 5 partial derivative relational expressions obtained in the step S300 into an effective stress coefficient expansion of clay shell porosity to obtain
Figure 100002_DEST_PATH_IMAGE018
S403, substituting the rigidity ratio of the stone skeleton to the clay and the shell-shaped clay mineral content relational expression into the step S402 formula to obtain the effective stress coefficient model of the clay shell-shaped porosity, wherein the effective stress coefficient model of the clay shell-shaped porosity is
Figure 100002_DEST_PATH_IMAGE019
In the formula
Figure 100002_DEST_PATH_IMAGE020
ωIs a constant;
s500, substituting the porosity, the clay content, the Poisson ratio and the rigidity ratio into a clay shell-shaped porosity effective stress coefficient model, and calculating the porosity effective stress coefficient.
In the method for calculating the effective stress coefficient of sandstone porosity based on the clay shell model, when the clay mineral content value of the rock sample is 0 and the rigidity ratio value is 1, the effective stress coefficient of porosity calculated by using the clay shell porosity effective stress coefficient model of the step S403 is equal to 1.
Compared with the prior art, the invention has the following beneficial effects: (1) the effective stress coefficient of the porosity is calculated by using the model, so that the complicated process of physical experiment testing is reduced; (2) the model is consistent with the reality, and the calculation result is more accurate; (3) the popularization is strong.
Drawings
In the drawings:
FIG. 1 is a technical scheme of the method.
Fig. 2 is a schematic diagram of a clay shell model.
Detailed Description
The present invention will be further described with reference to the following embodiments and drawings.
The invention provides a clay shell model-based sandstone porosity effective stress coefficient calculation method, and FIG. 1 is a technical route diagram of the method, and the method comprises the following steps:
s100, based on the analysis result of a scanning electron microscope, adhering chlorite and illite to the inner wall of a pore in a shell form, and providing a clay shell model, wherein the clay shell model is as shown in figure 2, the porosity of a rock sample containing clay minerals is measured by a pore infiltration instrument, the Poisson ratio is measured by Katahara in 1949, the Poisson ratios of the kaolinite, the illite and the chlorite are respectively 0.26, 0.25 and 0.27, and the shear modulus of the clay minerals is uniformly 0.25 multiplied by 104 MPa; the comprehensive elastic modulus and Poisson's ratio of the rock sample are measured by a longitudinal and transverse wave experiment;
s200, obtaining the clay shell porosity effective stress coefficient expansion by a chain method according to the porosity effective stress coefficient definition formula and combining the clay shell physical model porosity definition formula, the main steps are,
s201, determining a porosity effective stress coefficient definition formula and a clay shell physical model porosity definition formula, wherein the porosity effective stress coefficient definition formula is
Figure 267647DEST_PATH_IMAGE001
In the formula
Figure 818714DEST_PATH_IMAGE002
Is a porosity effective stress coefficient, and has no dimensional quantity;φporosity, dimensionless;P pis pore pressure in MPa;P cis confining pressure in MPa; the above-mentionedThe porosity of the clay shell-like physical model is defined as
Figure 352463DEST_PATH_IMAGE003
In the formula:
Figure 174051DEST_PATH_IMAGE004
is the pore inner radius, and the unit is mum; b is the outer radius of the pores in μm;
s202, obtaining a clay shell porosity effective stress coefficient expansion form by a chain method, wherein the clay shell porosity effective stress coefficient expansion form is
Figure 400633DEST_PATH_IMAGE005
S300, solving a partial derivative relational expression of the displacement of the inner wall surface and the outer wall surface of the pore to the pore pressure and the confining pressure and a partial derivative relational expression of the porosity to the inner wall surface and the outer wall surface of the pore, mainly comprising the following steps,
s301, calculating a partial derivative relational expression of the displacement of the inner wall surface of the pore to the pore pressure, wherein the partial derivative relational expression of the displacement of the inner wall surface of the pore to the pore pressure is
Figure 704576DEST_PATH_IMAGE006
In the formula
Figure 776437DEST_PATH_IMAGE007
Is the partial derivative of the displacement of the inner wall surface of the pore to the pore pressure;v cis the Poisson ratio of clay minerals without dimension;μ cthe Lame coefficient of clay mineral is dimensionless;A 1pB 1pis a transition coefficient;
s302, calculating a partial derivative relation formula of the displacement of the inner wall surface of the pore to the confining pressure, wherein the partial derivative relation formula of the displacement of the inner wall surface of the pore to the confining pressure is
Figure 482225DEST_PATH_IMAGE008
In the formula
Figure 867592DEST_PATH_IMAGE009
The partial derivative of the displacement of the inner wall surface of the pore to the confining pressure;A 1cB 1cis a transition coefficient;
s303, solving a partial derivative relational expression of the displacement of the outer wall surface of the pore to the pore pressure, wherein the partial derivative relational expression of the displacement of the outer wall surface of the pore to the pore pressure is
Figure 658831DEST_PATH_IMAGE010
In the formula
Figure 534383DEST_PATH_IMAGE011
The partial derivative of the displacement of the outer wall surface of the pore to the pore pressure is obtained; b is the outer radius of the pores in μm;v rthe Poisson ratio of the rock skeleton is a dimensionless quantity;A 2pB 2pis a transition coefficient;
s304, solving a partial derivative relational expression of the displacement of the outer wall surface of the pore to the confining pressure, wherein the partial derivative relational expression of the displacement of the outer wall surface of the pore to the pore pressure is
Figure 829098DEST_PATH_IMAGE012
In the formula
Figure 397483DEST_PATH_IMAGE013
The partial derivative of the displacement of the outer wall surface of the pore to the confining pressure is obtained;μ rthe Lame coefficient of the rock skeleton is free of dimensional quantity;A 2cB 2cis a transition coefficient;
s305, calculating partial derivative relational expressions of the porosity on the inner wall surface and the outer wall surface of the pore, wherein the partial derivative relational expressions of the porosity on the inner wall surface and the outer wall surface of the pore are respectively
Figure 646324DEST_PATH_IMAGE014
S400, defining a shell-shaped clay mineral content relational expression and a rigidity ratio of a rock framework to clay to obtain an effective stress coefficient equation of the clay shell-shaped porosity;
s401, defining the relation of the content of the shell-shaped clay minerals as
Figure 325567DEST_PATH_IMAGE015
Defining the rigidity ratio of the rock skeleton to the clay as in the formulaF cThe content of the shell-shaped clay minerals is zero dimensional quantity; c is the radius of the inner wall surface of the rock framework, and the unit is mum; the rigidity ratio of the stone framework to the clay is free of dimensional quantity;
s402, substituting the 5 partial derivative relational expressions obtained in the step S300 into an effective stress coefficient expansion of clay shell porosity to obtain
Figure 278107DEST_PATH_IMAGE018
S403, substituting the rigidity ratio of the stone skeleton to the clay and the shell-shaped clay mineral content relational expression into the step S402 formula to obtain the effective stress coefficient model of the clay shell-shaped porosity, wherein the effective stress coefficient model of the clay shell-shaped porosity is
Figure DEST_PATH_IMAGE021
In the formula
Figure 557779DEST_PATH_IMAGE020
ωIs a constant;
s500, substituting the porosity, the clay content, the Poisson ratio and the rigidity ratio into a clay shell-shaped porosity effective stress coefficient model, and calculating the porosity effective stress coefficient.
Further, when the clay mineral content of the rock sample is taken as 0 and the rigidity ratio is taken as 1, the effective stress coefficient of porosity calculated by using the clay shell-shaped porosity effective stress coefficient model of the step S403 is equal to 1.
The method comprises the following steps:
the first verification method comprises the following steps: for homogeneous single-component rock samples, the effective stress coefficient of porosity is 1, and the degree of influence of internal pressure and confining pressure on the porosity is the same. One of the methods for verifying whether the theoretical model is reasonable is to take an extreme value for the parameter and analyze whether the calculation result under the extreme condition is correct. For conventional sandstone samples, the main source of the effective stress coefficient of porosity not equal to 1 is the presence of clay minerals, in other words, the effective stress coefficient of porosity of a one-component homogeneous rock sample should be equal to 1. To change the model to a single component model, there are two approaches: (1) taking the clay mineral content as 0; (2) the stiffness ratio was taken to be 1. Both of them are satisfied to explain the reliability of the theoretical model.
By respectively adopting the methods researched by the invention and the research on the porosity stress sensitivity of clay-containing tight sandstone, extreme values are taken to calculate the effective stress coefficient of porosity, and the calculation results are shown in table 1. None of the results of the calculation of the velcade papers satisfies this phenomenon, and the gap is large.
TABLE 1 theoretical verification of effective stress coefficient of porosity for shell model
Figure DEST_PATH_IMAGE022
And a second verification method:
filling modes of chlorite and illite in clay minerals conform to distribution characteristics of a shell model, a porosity effective stress coefficient determination experiment is designed and developed, porosity values under different confining pressure and pore fluid pressure combinations are determined, regression fitting is carried out on confining pressure, internal pressure and porosity data obtained by the experiment by adopting a binary linear regression method, and the relation formula is
Figure DEST_PATH_IMAGE023
Wherein:a 0a 1a 2-fitting coefficients, dimensionless; -confining pressure, MPa; internal pressure, MPa.
Combining the definition of the effective stress coefficient of the porosity, obtaining the effective stress coefficient expression of the porosity as follows:
Figure DEST_PATH_IMAGE024
according to the above method, three blocks can be calculated by experimentThe effective stress coefficient of the porosity of the rock sample is compared with the calculation result of a theoretical model, such as the result shown in table 2.
TABLE 2 Shell model porosity effective stress coefficient calculation
Figure DEST_PATH_IMAGE025
The comparison result shows that the absolute error between the theoretical calculation result and the measured value is small, the absolute error is less than 1 percent, and the precision is high, so that the reliability of the theoretical model is demonstrated.
Compared with the prior art, the invention has the following beneficial effects: (1) the effective stress coefficient of the porosity is calculated by using the model, so that the complicated process of physical experiment testing is reduced; (2) the model is consistent with the reality, and the calculation result is more accurate; (3) the popularization is strong.
Finally, it should be noted that: although the present invention has been described in detail with reference to the above embodiments, it should be understood by those skilled in the art that: modifications and equivalents may be made thereto without departing from the spirit and scope of the invention and it is intended to cover in the claims the invention as defined in the appended claims.

Claims (2)

1. The sandstone porosity effective stress coefficient calculation method based on the clay shell model is characterized by comprising the following steps of:
s100, determining the porosity, the clay content, the Poisson ratio and the rigidity ratio of a clay mineral-containing rock sample;
s200, obtaining the clay shell porosity effective stress coefficient expansion by a chain method according to the porosity effective stress coefficient definition formula and combining the clay shell physical model porosity definition formula, the main steps are,
s201, determining a porosity effective stress coefficient definition formula and a clay shell physical model porosity definition formula, wherein the porosity effective stress coefficient definition formula is
Figure 92098DEST_PATH_IMAGE001
In the formula
Figure DEST_PATH_IMAGE002
Is a porosity effective stress coefficient, and has no dimensional quantity;φporosity, dimensionless;P pis pore pressure in MPa;P cis confining pressure in MPa; the porosity of the clay shell-shaped physical model is defined as
Figure 150184DEST_PATH_IMAGE003
In the formula:
Figure DEST_PATH_IMAGE004
is the pore inner radius, and the unit is mum; b is the outer radius of the pores in μm;
s202, obtaining a clay shell porosity effective stress coefficient expansion form by a chain method, wherein the clay shell porosity effective stress coefficient expansion form is
Figure 832052DEST_PATH_IMAGE005
S300, solving a partial derivative relational expression of the displacement of the inner wall surface and the outer wall surface of the pore to the pore pressure and the confining pressure and a partial derivative relational expression of the porosity to the inner wall surface and the outer wall surface of the pore, mainly comprising the following steps,
s301, calculating a partial derivative relational expression of the displacement of the inner wall surface of the pore to the pore pressure, wherein the partial derivative relational expression of the displacement of the inner wall surface of the pore to the pore pressure is
Figure DEST_PATH_IMAGE006
In the formula
Figure 633785DEST_PATH_IMAGE007
Is the partial derivative of the displacement of the inner wall surface of the pore to the pore pressure;v cis the Poisson ratio of clay minerals without dimension;μ cis the Lame coefficient of clay mineral and has no dimensionAn amount;A 1pB 1pis a transition coefficient;
s302, calculating a partial derivative relation formula of the displacement of the inner wall surface of the pore to the confining pressure, wherein the partial derivative relation formula of the displacement of the inner wall surface of the pore to the confining pressure is
Figure DEST_PATH_IMAGE008
In the formula
Figure 854682DEST_PATH_IMAGE009
The partial derivative of the displacement of the inner wall surface of the pore to the confining pressure;A 1cB 1cis a transition coefficient;
s303, solving a partial derivative relational expression of the displacement of the outer wall surface of the pore to the pore pressure, wherein the partial derivative relational expression of the displacement of the outer wall surface of the pore to the pore pressure is
Figure DEST_PATH_IMAGE010
In the formula
Figure 847784DEST_PATH_IMAGE011
The partial derivative of the displacement of the outer wall surface of the pore to the pore pressure is obtained; b is the outer radius of the pores in μm;v rthe Poisson ratio of the rock skeleton is a dimensionless quantity;A 2pB 2pis a transition coefficient;
s304, solving a partial derivative relational expression of the displacement of the outer wall surface of the pore to the confining pressure, wherein the partial derivative relational expression of the displacement of the outer wall surface of the pore to the pore pressure is
Figure DEST_PATH_IMAGE012
In the formula
Figure 521342DEST_PATH_IMAGE013
The partial derivative of the displacement of the outer wall surface of the pore to the confining pressure is obtained;μ rthe Lame coefficient of the rock skeleton is free of dimensional quantity;A 2cB 2cis a transition coefficient;
s305, calculating the porosity to the inner wall of the poreThe partial derivative relation between the surface and the outer wall surface of the pore, and the partial derivative relation between the porosity and the inner wall surface of the pore and the outer wall surface of the pore are respectively
Figure DEST_PATH_IMAGE014
S400, defining a shell-shaped clay mineral content relation and a rigidity ratio of a rock framework to clay to obtain a clay shell-shaped porosity effective stress coefficient model;
s401, defining the relation of the content of the shell-shaped clay minerals as
Figure 657925DEST_PATH_IMAGE015
Defining the rigidity ratio of the rock skeleton to the clay as in the formulaF cThe content of the shell-shaped clay minerals is zero dimensional quantity; c is the radius of the inner wall surface of the rock framework, and the unit is mum; the rigidity ratio of the stone framework to the clay is free of dimensional quantity;
s402, substituting the 5 partial derivative relational expressions obtained in the step S300 into an effective stress coefficient expansion of clay shell porosity to obtain
Figure DEST_PATH_IMAGE018
S403, substituting the rigidity ratio of the stone skeleton to the clay and the shell-shaped clay mineral content relational expression into the step S402 formula to obtain the effective stress coefficient model of the clay shell-shaped porosity, wherein the effective stress coefficient model of the clay shell-shaped porosity is
Figure DEST_PATH_IMAGE019
In the formula
Figure DEST_PATH_IMAGE020
ωIs a constant;
s500, substituting the porosity, the clay content, the Poisson ratio and the rigidity ratio into a clay shell-shaped porosity effective stress coefficient model, and calculating the porosity effective stress coefficient.
2. The clay shell model-based sandstone porosity effective stress coefficient calculation method according to claim 1, wherein the method comprises the following steps: and when the clay mineral content of the rock sample is taken as 0 and the rigidity ratio is taken as 1, calculating by using the clay shell-shaped porosity effective stress coefficient model obtained in the step S403 to obtain the porosity effective stress coefficient equal to 1.
CN202110181206.9A 2021-02-10 2021-02-10 Sandstone porosity effective stress coefficient calculation method based on clay shell model Pending CN112525803A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110181206.9A CN112525803A (en) 2021-02-10 2021-02-10 Sandstone porosity effective stress coefficient calculation method based on clay shell model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110181206.9A CN112525803A (en) 2021-02-10 2021-02-10 Sandstone porosity effective stress coefficient calculation method based on clay shell model

Publications (1)

Publication Number Publication Date
CN112525803A true CN112525803A (en) 2021-03-19

Family

ID=74975658

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110181206.9A Pending CN112525803A (en) 2021-02-10 2021-02-10 Sandstone porosity effective stress coefficient calculation method based on clay shell model

Country Status (1)

Country Link
CN (1) CN112525803A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112800632A (en) * 2021-04-01 2021-05-14 西南石油大学 Sandstone reservoir permeability effective stress coefficient calculation method based on granular model

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4782501A (en) * 1987-09-04 1988-11-01 Mobil Oil Corporation Method for determining pore volume compressibility of a porous material
CN102519999A (en) * 2011-11-11 2012-06-27 中国石油大学(北京) Nuclear magnetic resonance analyzer and nuclear magnetic resonance measuring method
US20140122034A1 (en) * 2011-12-09 2014-05-01 Jonathan M. Hanson Drill bit body rubbing simulation
CN104089823A (en) * 2014-07-07 2014-10-08 中国石油大学(北京) Method for determining effective rock stress coefficient based on pore compression experiment
CN104133050A (en) * 2014-07-23 2014-11-05 中国科学院武汉岩土力学研究所 Method for simultaneously testing effective stress coefficient and porosity of porous rock under stress conditions
CN106368686A (en) * 2015-07-17 2017-02-01 中国石油化工股份有限公司 Method for computing stratum pressure based on rock Poisson's ratio
CN106483056A (en) * 2016-09-28 2017-03-08 西南石油大学 A kind of shale porosity measurement method based on longitudinal wave velocity and measurement apparatus
CN106841001A (en) * 2017-01-17 2017-06-13 西南石油大学 A kind of tight sand porosity based on reservoir quality Analysis The Main Control Factor, Permeability Prediction method
CN107829731A (en) * 2017-11-06 2018-03-23 陈国军 A kind of volcanic rock porosity correction method of clay alteration
US20180357347A1 (en) * 2013-02-04 2018-12-13 Comsol Ab Apparatus and method for defining coupled systems on spatial dimensions and extra dimensions
CN109241588A (en) * 2018-08-21 2019-01-18 北京大学 A kind of analogy method of the monolete extension based on quasi-continuous geomechanics model
CN110427732A (en) * 2019-09-11 2019-11-08 河海大学 A kind of bilayer heterogeneous soil model of a layered structure construction method
CN110702585A (en) * 2019-11-05 2020-01-17 重庆科技学院 Rock compression coefficient calculation method
CN110909487A (en) * 2019-12-06 2020-03-24 西南石油大学 Method for establishing fine anisotropic shale rock physical model
CN111460601A (en) * 2019-12-06 2020-07-28 西南石油大学 Orthotropic formation ground stress prediction method based on rock physics modeling
CN111999227A (en) * 2020-08-19 2020-11-27 中国石油大学(北京) Stress sensitivity testing method and device

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4782501A (en) * 1987-09-04 1988-11-01 Mobil Oil Corporation Method for determining pore volume compressibility of a porous material
CN102519999A (en) * 2011-11-11 2012-06-27 中国石油大学(北京) Nuclear magnetic resonance analyzer and nuclear magnetic resonance measuring method
US20140122034A1 (en) * 2011-12-09 2014-05-01 Jonathan M. Hanson Drill bit body rubbing simulation
US20180357347A1 (en) * 2013-02-04 2018-12-13 Comsol Ab Apparatus and method for defining coupled systems on spatial dimensions and extra dimensions
CN104089823A (en) * 2014-07-07 2014-10-08 中国石油大学(北京) Method for determining effective rock stress coefficient based on pore compression experiment
CN104133050A (en) * 2014-07-23 2014-11-05 中国科学院武汉岩土力学研究所 Method for simultaneously testing effective stress coefficient and porosity of porous rock under stress conditions
CN106368686A (en) * 2015-07-17 2017-02-01 中国石油化工股份有限公司 Method for computing stratum pressure based on rock Poisson's ratio
CN106483056A (en) * 2016-09-28 2017-03-08 西南石油大学 A kind of shale porosity measurement method based on longitudinal wave velocity and measurement apparatus
CN106841001A (en) * 2017-01-17 2017-06-13 西南石油大学 A kind of tight sand porosity based on reservoir quality Analysis The Main Control Factor, Permeability Prediction method
CN107829731A (en) * 2017-11-06 2018-03-23 陈国军 A kind of volcanic rock porosity correction method of clay alteration
CN109241588A (en) * 2018-08-21 2019-01-18 北京大学 A kind of analogy method of the monolete extension based on quasi-continuous geomechanics model
CN110427732A (en) * 2019-09-11 2019-11-08 河海大学 A kind of bilayer heterogeneous soil model of a layered structure construction method
CN110702585A (en) * 2019-11-05 2020-01-17 重庆科技学院 Rock compression coefficient calculation method
CN110909487A (en) * 2019-12-06 2020-03-24 西南石油大学 Method for establishing fine anisotropic shale rock physical model
CN111460601A (en) * 2019-12-06 2020-07-28 西南石油大学 Orthotropic formation ground stress prediction method based on rock physics modeling
CN111999227A (en) * 2020-08-19 2020-11-27 中国石油大学(北京) Stress sensitivity testing method and device

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
STANISLAV GLUBOKOVSKIKH 等: ""Effect of micro‐inhomogeneity on the effective stress coefficients and undrained bulk modulus of a poroelastic medium: a double spherical shell model"", 《GEOPHYSICAL PROSPECTING》 *
张睿: ""页岩储层有效应力及应力敏感机理研究"", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112800632A (en) * 2021-04-01 2021-05-14 西南石油大学 Sandstone reservoir permeability effective stress coefficient calculation method based on granular model

Similar Documents

Publication Publication Date Title
CN201242522Y (en) Apparatus for measuring low-leakage and ultra-low leakage rock core coke porosity
US8224629B2 (en) Method of modelling a saturation dependant property in a sample
CN110702585B (en) Rock compression coefficient calculation method
CN110348590A (en) A kind of shale adsorbed gas content calculation method
CN112525803A (en) Sandstone porosity effective stress coefficient calculation method based on clay shell model
CN102323199B (en) Porosity determination experiment design method and design system thereof
CN106354955A (en) Sliding bearing rigidity recognition method based on mill vibration mode parameters
CN106897531A (en) A kind of method for quantitatively evaluating of hyposmosis limestone reservoir permeability
CN110700820A (en) Dessert classification method for compact oil reservoir in northern part of Songliao basin
CN112214935A (en) Method and system for calculating fluid phase state in porous medium
CN106840790A (en) CO is tested based on tubule consolidated model long2The method and system of crude oil MMP
CN110344822B (en) Determination method of water saturation of calcareous compact reservoir
CN106092820A (en) A kind of mud shale gas absorption tolerance bearing calibration based on density and device
CN111666665A (en) Method for improving accuracy of dynamic and static Young modulus conversion model of carbonatite
CN111751259A (en) Method and device for determining effective porosity of irregular rock sample
CN113624847B (en) Method for establishing prediction model of shale hydration damage coefficient and prediction method
CN114441401A (en) Determination method and device for overburden porosity, computer equipment and storage medium
CN111693427B (en) Analysis method for mobility of oil and gas reservoir fluid
CN106840791A (en) Back-up sand tubule consolidated model preparation method long and system
CN108090264B (en) Consolidation model test similarity criterion establishing method
CN114077784A (en) Determination method and device of overburden porosity, storage medium and computer equipment
Olek et al. Determination of quasi-filtration phase of consolidation based on experimental and theoretical course of the uniaxial deformation and distribution of pore pressure
CN114441405B (en) Quantitative evaluation method for secondary hole increasing amplitude based on compaction and cementation hole decreasing trend
CN108956423A (en) Determine rock two dimension porosity-three-dimensional apertures porosity quantitative relationship method
CN114492235B (en) Stratum pore fluid pressure prediction method suitable for carbonate rock stratum

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20210319