CN105352787B - 一种测定岩性非衰减蠕变阶段在不同应力下的失稳时间的方法 - Google Patents
一种测定岩性非衰减蠕变阶段在不同应力下的失稳时间的方法 Download PDFInfo
- Publication number
- CN105352787B CN105352787B CN201510801486.3A CN201510801486A CN105352787B CN 105352787 B CN105352787 B CN 105352787B CN 201510801486 A CN201510801486 A CN 201510801486A CN 105352787 B CN105352787 B CN 105352787B
- Authority
- CN
- China
- Prior art keywords
- msub
- rock
- creep
- mfrac
- mrow
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种测定岩性非衰减蠕变阶段在不同应力下失稳时间的方法，包括如下步骤：步骤a，对研究区块岩石进行取样；步骤b，对所取岩样进行蠕变力学实验，步骤c，将测定的参数代入蠕变方程；步骤d，根据蠕变方程中的应力、应变及时间关系可以推断出对应应力条件下岩石的失稳时间。本发明能很好的描述岩石的蠕变特性并测定岩性非衰减蠕变阶段在不同应力下的失稳时间，对于研究油气井井壁及其他岩石工程稳定具有重要意义。
Description
技术领域
属于石油天然气勘探开发及其他岩石工程领域，具体涉及一种从能量转换角度描述岩性蠕变的方法，通过判别分析，建立能够有效描述岩石第三蠕变阶段的本构模型，解决了常规蠕变模型难以解释、描述岩石蠕变第三阶段的问题，
背景技术
流变作为岩石的一项重要力学性质与许多工程质量息息相关，在油气钻井工程及其他岩土工程中，由岩石流变导致工程失稳所造成的施工安全及经济损失不可小觑。
国内外学者对岩石的流变特性做了大量的理论及实验研究。毛东风等从理论上对均匀地应力场下及非均匀地应力场下蠕变引起的岩压外载及因远近含水量不同而引起的变化进行了分析；范翔宇等对煤岩储气层岩石蠕变特性与本构模型进行了研究，建立了饱依丁–汤姆逊体与宾汉姆体串联模型；杨逾等建立了西原加速模型；何峰等对破碎煤岩峰后渗流-蠕变进行了实验研究；李金和等通过单轴压缩蠕变实验研究了大理岩蠕变模型参数；陈绍杰等采用LS-DYNA数值方法分析了深部煤柱的蠕变支撑效应；周长冰等建立了高温三轴应力条件下的煤岩蠕变模型；赵斌等对不同应力水平条件下煤岩蠕变破坏特征进行了研究；刘钦等对灰质页岩蠕变特性进行了试验研究；杨红伟等对三轴压缩条件下岩石孔隙水压力分级加载蠕变试验进行了蠕变特性及模型研究；赵洪宝等对处于不同载荷水平作用下的红砂岩试样剪切蠕变力学特性进行了系统的研究，并基于试验结果尝试解释了一些边坡失稳发生时的现象；黄达等开展相同初始高应力状态条件下恒轴压分级卸围压三轴卸荷蠕变试验。
以上理论和方法基本是从力学的角度研究了岩石的蠕变特性，为后人研究岩石蠕变特性打下了一定的基础。由于从不同角度研究岩石的蠕变特性，取得的效果会有所差别。
岩石蠕变通常分为两种形式，即衰减蠕变和非衰减蠕变。在岩石所受荷载小于岩石的长期强度时，岩石发生衰减蠕变；在岩石所受荷载大于长期强度时，岩石发生非衰减蠕变。
岩石的工程失稳很多情况是由非衰减蠕变引起的。现在从能量转换的角度分析岩石非衰减蠕变的三个阶段：(1)减速蠕变阶段，该阶段主要以弹性形变为主。其原因是岩石要破裂，首先要破坏各基元的价键。而岩石基元粒子之间是通过各种化学键相结合的，这些价键的结合力是吸引力和排斥力共同作用的结果。在此阶段，外力功主要转变为弹性势能(由于岩石先天性的缺陷，因此该阶段并非表现理想弹性体的特点)。(2)稳定蠕变阶段，主要以塑性变形为主。其原因是岩石弹性势能达到一定程度，价键开始断裂，岩石的不同组元内开始大量产生各种缺陷，使得岩石所受的外力功主要转化为塑性变形能。使该阶段主要表现为塑性变形的特点，使组元结构发生不可逆的畸变。(3)加速蠕变阶段，该阶段发展迅速，并且在短时间内发生破坏。主要是因为岩石在外力功的作用下，岩石内部积蓄的塑性变形能已经达到极限，岩石产生宏观裂纹并迅速扩展。在此阶段，外力功主要转化为动能，导致该阶段应变迅速增大，并使岩石快速破裂。
在油气钻采及其他岩石工程开挖过程中，随着时间的推移，井壁及其他岩石工程失稳事故时有发生，很大一部分是由岩石蠕变失稳造成的。传统的岩石蠕变模型对于描述蠕变的全过程效果不理想，主要存在以下问题：
(1)无法有效描述岩石蠕变的第三阶段(即岩石的加速蠕变阶段)。
(2)只能解释岩石蠕变的现象而不能解释岩石发生蠕变的内在原因及无法解释岩石的蠕变为什么会分为三个阶段。
西原模型是由胡克体、开尔文体和理想粘塑性体串联而成(见图1)，能够较全面反映岩石的弹-粘弹-粘塑性特征。
模型中，当岩石所受荷载(σ)小于岩石的长期强度(σS)时，只有胡克体和开尔文体发生形变；当岩石所受荷载大于岩石长期强度时，除了胡克体及开尔文体发生形变外，理想粘塑性体也会发生形变。
图1中E1、E2分别为胡克体和开尔文体中的弹性模量；该模型对应的蠕变方程为：
(1)当时，蠕变方程：
(2)当时，蠕变方程：
式中：ε为岩石应变(％)；σ0为岩石所受荷载(MPa)；t为时间(h)，η1、η2分别为开尔文体和理想粘塑性体中的黏性系数。
从上述理论分析可知，在岩石的蠕变模型中，西原模型能够很好的描述大部分岩石的衰减蠕变过程，但是对于描述非衰减蠕变具有局限性。其主要原因是西原模型的元件组成只能解释出岩石在外力功的情况下转变为弹性势能，塑性势能两个阶段，对于第三阶段，即外力功转换为动能阶段无法解释。
发明内容
本发明目的是提供一种测定岩性非衰减蠕变阶段在不同应力下的失稳时间的方法，解决常规蠕变模型难以解释、描述岩石蠕变第三阶段的问题。
本发明的技术方案为：一种测定岩性非衰减蠕变阶段在不同应力下的失稳时间的方法，其特征在于，包括如下步骤：
步骤a:对研究区块岩石进行取样；
步骤b:对所取岩样进行蠕变力学实验；
步骤c:将测定的参数代入蠕变方程；蠕变方程为：
其中，ε为岩石应变，σ0为岩石所受荷载，σS为岩石的长期强度，E1为胡克体的弹性模量，E2为开尔文体的弹性模量，η1为开尔文体的黏性系数，η2为理想粘塑性体的黏性系数，t岩石蠕变时间，ρ为岩样的密度，L为岩样的高度，λ为功能转换率；
步骤d:根据蠕变方程中的应力、应变及时间关系推断出对应应力条件下岩石的失稳时间。
本发明与现有技术相比具有如下优点：
(1)本发明能够解释岩石蠕变三个阶段的功能转换关系。
(2)本发明能很好的测定岩性非衰减蠕变阶段在不同应力下的失稳时间，对于研究油气井井壁及其他岩石工程稳定具有重要意义。
附图说明
图1为西原模型图；
图2为本发明改进的西原模型图；
图3为8MPa煤岩蠕变实验实测值与理论值比较图；
图4为10MPa煤岩蠕变实验实测值与理论值比较图；
图5为12MPa煤岩蠕变实验实测值与理论值比较；
图6为38.83MPa二辉橄榄岩蠕变实验实测值与理论值比较图。
具体实施方式
本发明的蠕变方程的推导过程如下：
取岩石试件在单位时间内的微元应变为△ε，应变速率为△υ，根据功能关系推导出如下公式：
对公式(1)进行积分可得：
公式(1)、(2)中σ为岩样所受荷载；S为岩样所受荷载面积；L为岩样的高度；λ为功能转换率；M所受荷载岩样的质量；ρ为岩样的密度；t为时间；σ·S·L·Δε为荷载所做的功；ε为应变；为岩样变形破坏时的动能；为岩样变形破坏时的加速度。
因此，在第三蠕变阶段岩石的能量转换关系应遵循公式(2)。
将以上公式(2)代入原有模型(以西原模型为例)本构方程得到新的蠕变模型(改进的西原模型)，其本构方程为：
蠕变方程为：
①当σ≦σs时，只有广义开尔文体发生形变
②当σ﹥σs时，除了广义开尔文体发生形变外，还有理想粘塑性体及绕组也要发生形变。
采用WDW-200微机控制电子万能实验机对制作好的岩石试件进行蠕变力学实验，分析出不同轴向荷载条件下岩石的应变与时间关系。
(1)将岩石试件安放在加载台上，调节横梁，使上承压板与试件上端接触；
(2)以0.03mm/min的加载速率施加轴向荷载至设定的荷载值；
(3)保持轴向荷载不变，进行岩石蠕变实验，并记录实验过程中岩石试件的轴向应变与时间。
实施例1
为了验证两种模型对煤岩蠕变特性的适用性，本次实验采用山西晋城地区的煤岩进行煤岩蠕变力学实验研究。运用matlab软件对实验数据进行拟合，得到各拟合参数见表1及拟合结果见图3、图4和图5。
表1:煤岩蠕变实验拟合参数表
从表1中可以看出：在相同的荷载条件下，改进的西原模型与实际结果的拟合程度要高于西原模型。
从图3、图4和图5可以看出，改进西原模型所描述的曲线与实际结果更加吻合，尤其能更直观的描述出煤岩加速蠕变阶段的曲线特性。
实施例2
本实验岩石试件取自金川集团股份有限公司某矿区的橄榄岩试样。运用matlab软件对该岩石试件蠕变实验结果进行拟合，得到各拟合参数见表2及拟合结果见图6。
表2：38.83MPa橄榄岩蠕变实验拟合参数表
从表2可以可看出：在相同荷载条件下，改进的西原模型对于实测值的拟合程度高于西原模型。
图6直观的反映了两种模型对于研究橄榄岩蠕变特性所取得的效果。从图中不难看出，改进后的西原模型能更好的表现出橄榄岩的加速蠕变阶段特性，而西原模型很难表现出该阶段的特点。
以上所述实施例仅表达了本申请的具体实施方式，其描述较为具体和详细，但并不能因此而理解为对本申请保护范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请技术方案构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。
Claims (1)
1.一种测定岩性非衰减蠕变阶段在不同应力下的失稳时间的方法，其特征在于，包括如下步骤：
步骤a:对研究区块岩石进行取样；
步骤b:对所取岩石进行蠕变力学实验；
步骤c:将测定的参数代入蠕变方程；蠕变方程为：
<mrow>
<mi>&epsiv;</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>&sigma;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>&sigma;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&eta;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&eta;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>t</mi>
<mo>+</mo>
<mfrac>
<mi>&lambda;</mi>
<mrow>
<mn>2</mn>
<msup>
<mi>&rho;L</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mi>s</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msup>
<mi>t</mi>
<mn>3</mn>
</msup>
</mrow>
其中，ε为岩石应变，σ0为岩石所受荷载，σs为岩石的长期强度，E1为胡克体的弹性模量，E2为开尔文体的弹性模量，η1为开尔文体的黏性系数，η2为理想粘塑性体的黏性系数，t为岩石蠕变时间，ρ为岩石的密度，L为岩石的高度，λ为功能转换率；
步骤d:根据蠕变方程中的岩石所受荷载、应变及时间关系推断出对应应力条件下岩石的失稳时间。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510801486.3A CN105352787B (zh) | 2015-11-19 | 2015-11-19 | 一种测定岩性非衰减蠕变阶段在不同应力下的失稳时间的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510801486.3A CN105352787B (zh) | 2015-11-19 | 2015-11-19 | 一种测定岩性非衰减蠕变阶段在不同应力下的失稳时间的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105352787A CN105352787A (zh) | 2016-02-24 |
| CN105352787B true CN105352787B (zh) | 2017-11-14 |
Family
ID=55328794
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510801486.3A Active CN105352787B (zh) | 2015-11-19 | 2015-11-19 | 一种测定岩性非衰减蠕变阶段在不同应力下的失稳时间的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105352787B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106885733A (zh) * | 2017-03-15 | 2017-06-23 | 中国电建集团华东勘测设计研究院有限公司 | 一种由快速蠕变试验直接确定岩体长期强度的方法 |
| CN107748111B (zh) * | 2017-10-13 | 2020-08-04 | 华北水利水电大学 | 一种岩体结构面长期抗剪强度的确定方法 |
| CN109187199B (zh) * | 2018-09-18 | 2021-01-22 | 中国石油大学(华东) | 单轴压缩下加锚岩体蠕变特性的黏弹性理论分析方法 |
| CN110610062B (zh) * | 2019-10-21 | 2021-09-17 | 安徽理工大学 | 一种长期静荷载作用下的岩石弹性模量计算方法 |
| CN114398805B (zh) * | 2022-03-25 | 2022-07-08 | 四川省水利水电勘测设计研究院有限公司 | 一种水岩耦合作用下裂隙岩石蠕变模型构建方法及系统 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104280299A (zh) * | 2014-10-11 | 2015-01-14 | 中国石油大学(华东) | 基于细观组构模拟的岩石蠕变性质预测方法 |
| CN104833593A (zh) * | 2015-01-23 | 2015-08-12 | 河海大学 | 应用于软岩洞室预警系统的围岩蠕变破坏预警方法 |
-
2015
- 2015-11-19 CN CN201510801486.3A patent/CN105352787B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105352787A (zh) | 2016-02-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105352787B (zh) | 一种测定岩性非衰减蠕变阶段在不同应力下的失稳时间的方法 | |
| Lang et al. | Creep experiment and rheological model of deep saturated rock | |
| Yang et al. | Numerical study on failure behavior of brittle rock specimen containing pre-existing combined flaws under different confining pressure | |
| CN102221504B (zh) | 一种测定破裂后岩石三轴强度的加载路径方法 | |
| Zhong et al. | Macro-fracture mode and micro-fracture mechanism of shale | |
| CN108827774B (zh) | 煤岩脆性评价方法 | |
| Zhang et al. | Dynamic behavior of frozen soil under uniaxial strain and stress conditions | |
| CN105021458B (zh) | 一种含油泥页岩杨氏模量的定量评价方法 | |
| Zhou et al. | Mechanical behavior of red sandstone under cyclic point loading | |
| CN204165873U (zh) | 用于测试采动过程卸荷煤岩体渗透率的系统 | |
| Liu et al. | Numerical simulation on effects of embedded crack on rock fragmentation by a tunnel boring machine cutter | |
| Li et al. | Experimental study of compaction creep model of broken rock | |
| Zhou et al. | Degradation mechanism of rock under impact loadings by integrated investigation on crack and damage development | |
| ZHAO et al. | Study on the shear strength of deep reconstituted soils | |
| Wang et al. | Damage evolution model of cemented tailing backfill based on acoustic emission energy | |
| Zhu et al. | Nonlinear constitutive model and discrete-element-method modeling of synthetic Methane Hydrate sand | |
| Li et al. | Deformation and failure mechanism of surrounding rock considering the influence of tensile yield on mechanical properties of rock mass | |
| Zhang et al. | Experimental study on the effect of confining pressure and interval time on the fatigue mechanical properties of mudstone | |
| Shao-Liang | Characterization and Application: Tri-parameter Strength Criterion for Rocks Based on Ultimate Stress Ratio | |
| He et al. | Experimental Study on Influence of Unloading Rate on Creep Characteristics of Marble under High Stress | |
| Lee et al. | Evaluation of cementation effect of sand using cone resistance | |
| Moravej et al. | Thermomechanical Behaviour of Rock Salt | |
| Li et al. | An experimental study of loading rate effects on mechanical characteristics of cemented paste backfill | |
| Zhou et al. | Research on the Dimensional Radius Control Parameters of Crack Plastic Zone Based on the Coal Cutting Cavity Crushing Model | |
| Yao et al. | Monotonic and cyclic analysis of granular soils |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |