CN108982818A - 一种描述外载作用岩体微细观破裂与微裂隙形成的方法 - Google Patents
一种描述外载作用岩体微细观破裂与微裂隙形成的方法 Download PDFInfo
- Publication number
- CN108982818A CN108982818A CN201811006299.6A CN201811006299A CN108982818A CN 108982818 A CN108982818 A CN 108982818A CN 201811006299 A CN201811006299 A CN 201811006299A CN 108982818 A CN108982818 A CN 108982818A
- Authority
- CN
- China
- Prior art keywords
- model
- fractional order
- micro
- rock
- kelvin
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Food Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Remote Sensing (AREA)
- General Health & Medical Sciences (AREA)
- Geology (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Medicinal Chemistry (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Adornments (AREA)
Abstract
本发明涉及的是一种描述外载作用岩体微细观破裂与微裂隙形成的方法,具体为:一、从细化角度研究岩石破裂问题,在岩石微观颗粒之间建立黏结模型,将颗粒简化为圆形颗粒,以两颗粒圆心为端点的建立等效微梁模型;二、建立基于分数阶的黏壶元件的本构关系;三、建立分数阶kelvin模型,由Abel黏壶元件与弹簧元件并联组成,四、建立由弹簧元件与黏壶元件串联的Maxwll模型;五、得到分数阶Burgers模型本构关系。本发明主要用于解决现阶段岩石复杂特性无法用宏观本构精确描述,以及利用分数阶建立的细观本构模型。
Description
技术领域
本发明涉及石油与天然气工程中的页岩、致密砂岩、煤岩等非常规油气储层的微观裂缝扩展领域,具体涉及一种描述外载作用岩体微细观破裂与微裂隙形成的方法。
背景技术
岩石内部存在大量规模不等的构造及复杂多变的构造(如孔洞缝的多样组合和特定岩性组合等)相比宏观尺度对储层岩石的均匀性、连续性、各向同性等认识而言,其在微观尺度上则显示出非均匀、不连续及各向异性的特征。由于岩石的结构特性,其内部存在大量的微观和宏观裂隙,岩石在外载作用下会依次经过裂隙压实阶段,线弹性阶段,微裂隙扩展阶段,大变形和累积破坏阶段,以及残余应力阶段。岩石受力经过微裂缝局部成核并持续发展,最后微裂缝迅速扩展、成核、逐步形成宏观破坏。岩石内部众多裂纹端部的应力集中是引起岩石开裂,进而发生解体破坏的原因。大量研究表明,岩石的宏观变形与破坏是受岩石内部的细观缺陷控制的。岩石裂纹萌生、扩展、贯通及相互作用的内部损伤演化过程和规律及其与应力之间的对应关系,从而应用于工程中指导生产。岩石的变形与破坏主要与岩石和宏观裂隙有关,而宏观裂隙是有众多的微裂隙的集合体。因此受载荷作用下的岩石微观裂隙的起裂、扩展过程,对非常规油气藏的开采过程中有重要的意义,尤其是页岩压裂形成缝网的过程中有重要意义。
岩石受外界载荷、变形、温度等条件的影响,会呈现与时间有关的变形、流动和破坏的规律性。所以建立合理的本构关系成为描述裂缝扩展的关键环节。到现在为止,岩石一直在套用传统连续介质模型,还没有它特有的本构模型。而传统连续介质力学方法直接将本构关系建立在宏观连续层次基础之上,而不是细观结构基础之上,这使得传统连续介质模型本身存在一定的局限性。岩石不同于其他金属材料,具有非常复杂的细观结构。在细观尺度上,岩石由矿物颗粒组成,颗粒之间相互胶结,形成界面网络结构。这种细观结构对岩石宏观力学属性具有决定性作用。
目前,关于岩石变形的研究仍然是将岩石假设为理想弹性体,通过泊松比来描述两个方向的变形。岩石的体积变形特性与理想弹性材料有很大区别,其最大的特点是体积变形的剪缩剪胀特性,如果只通过泊松比来反映这些复杂的现象就似乎不太合理。分数阶微积分的阶数是连续变化的,能够很好地体现自然界的连续性。基于分数阶微积分建立的物理模型,能够良好描述具有记忆及时间依赖性的物理现象,且具有参数少、适用性广泛等优点。分数阶微积分在力学中应用的较早成果始于黏弹性材料的本构模型研究,也是目前最活跃的领域之一,这其中包括黏弹性阻尼器研究、黏弹性流体研究和岩土流变模型研究等。因此,如何用分数阶描述岩石微观本构关系是亟待解决的技术问题。
发明内容
本发明的目的是提供一种描述外载作用岩体微细观破裂与微裂隙形成的方法,这种描述外载作用岩体微细观破裂与微裂隙形成的方法用于解决现阶段岩石复杂特性无法用宏观本构精确描述的问题。
本发明解决其技术问题所采用的技术方案是:这种描述外载作用岩体微细观破裂与微裂隙形成的方法:
一、从细化角度研究岩石破裂问题,在岩石微观颗粒之间建立黏结模型,将颗粒简化为圆形颗粒,以两颗粒圆心为端点的建立等效微梁模型:模型长 L=RA+RB,高L=2×(RA+RB)/2=RA+RB,宽度或厚度为t;
二、建立基于分数阶的黏壶元件的本构关系:
式中:η为黏性系数,k为常数,Γ(n+1)为n+1的阶乘,ε为应变,σ为应力,d是微分算子,D(·)为Riemann-Liouville型分数阶微积分算子;
三、建立分数阶kelvin模型,由Abel黏壶元件与弹簧元件并联组成,步骤如下:
σk1(t)=Ekεk1(t) (4)
分数阶kelvin模型为:
σk(t)=σk1+σk2=Ekεk1(t)+kDnεk2(t)=(Ek+kDn)εk(t) (6)
式中:Ek为kelvin体的弹性模量,k为常数;
四、建立由弹簧元件与黏壶元件串联的Maxwll模型,步骤如下:
σm1(t)=Emεm1(t) (7)
σm2(t)=ηmDεm2(t) (8)
式中:Em为kelvin体的弹性模量,k为常数;
五、将Maxwell模型与分数阶Kelvin模型串联,应变相加,应力相同,上式两边分别相加,当元件串联时,由ε=ε1+ε2,σ1=σ2,得到σ=σm=σk,ε=εm+εk,得到分数阶Burgers模型本构关系为:
式中:η为黏性系数,k为常数,ε为应变,σ为应力,Ek为kelvin体的弹性模量,Em为Maxwll体的弹性模量。
本发明具有以下有益效果:
1.本发明提供了一种基于分数阶建立微观颗粒之间黏结的本构关系,使用分数阶元件构建本构模型更加准确的描述岩石损伤、蠕变等性质,解决了岩石宏观本构关系描述岩石破裂过程中的不足。
2.本发明提供了一种基于分数阶方法描述裂缝扩展的方法,充分利用分数阶特性描述岩石破裂过程;提供了基于分数阶的描述的从微观角度描述了岩石微观裂缝起裂、扩展、成核提供了新方法。
附图说明
图1是本发明颗粒黏结模型示意图;
图2是颗粒受力时的模型示意图;
图3是颗粒之间接触黏结的本构关系模型示意图。
具体实施方式
下面结合附图对本发明作进一步的说明:
本发明自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的,解释本发明。
这种描述外载作用岩体微细观破裂与微裂隙形成的方法:岩石微观颗粒之间形成的黏结,定义黏结模型包括了颗粒的物理力学特性,颗粒黏结材料***共同描述岩石类材料的细观力学行为。将岩石颗粒简化为圆形颗粒,颗粒之间形成的黏结建立本构关系。
在建立颗粒之间的本构模型时,用分数阶黏壶元件替换黏壶元件。分数阶导数模型和整数阶导数模型的本质区别在于,对于时间而言,整数结导数所表征的是一个物理或力学过程某时刻的变化或某种性质,而对于分数阶导数所表征的性质则与整个发展历史有关。
根据经典固体和流体力学,理想固体的应力应变关系满足胡克定律σ(t)-ε(t),理想流体满足牛顿粘性定律σ(t)-d1ε(t)/dt1,如果将σ(t)-ε(t)改写为σ(t)-d0ε(t)/dt0,则有充分理由认为介于理想固体和理想流体之间的岩体材料采用分数阶微分形式表示的应力-应变关系式为:
步骤一:如图1所示建立微观颗粒之间的黏结模型。对于接触黏结模型,以颗粒的A,B圆心为端点的建立等效微梁模型。长L=RA+RB,高 L=2×(RA+RB)/2=RA+RB,宽(厚度)为t。
接触的黏结键受外力压缩,接触黏结受到外力挤压发生形变(如图2),基于分数阶建立其本构模型(如图3),其本构模型是由分数阶的Kelvin模型与 Maxwll模型串联组成的。
步骤二:建立基于分数阶的kelvin模型,先建立黏性元件的分数阶形式(Abel 元件)的本构关系如(1、2)式:
式中:η为黏性系数,k为常数,Γ(n+1)为n+1的阶乘,ε为应变,σ为应力。
步骤三:分数阶kelvin模型由Abel黏壶元件与弹簧元件并联组成。元件并联时应力相加,应变相同。对于分数阶的kelvin模型本构关系如(6)式:
σk=σk1+σk2,εk=εk1=εk2 (3)
σk1(t)=Ekεk1(t) (4)
σk(t)=σk1+σk2=Ekεk1(t)+kDnεk2(t)=(Ek+kDn)εk(t) (6)
步骤四:建立由弹簧元件与黏壶元件串联的Maxwll模型,串联的弹簧(应力σm1(t)、应变εm1(t))和定常黏壶(应力σm2(t)、应变εm2(t))本构关系如(7、8)式,在等式两边同作用算子(Ek+kDn),如(9、10)式。
σm1(t)=Emεm1(t) (7)
σm2(t)=ηmDεm2(t) (8)
步骤五:将Maxwell模型与分数阶Kelvin模型串联。应变相加,应力相同,上式两边分别相加,当元件串联时,由ε=ε1+ε2,σ1=σ2。可得σ=σm=σk,ε=εm+εk,得到分数阶Burgers模型本构关系如(12)式。
整理得:
式中:η为黏性系数,k为常数,ε为应变,σ为应力,Ek为kelvin体的弹性模量,Em为Maxwll体的弹性模量。
本发明主要用于解决现阶段岩石复杂特性无法用宏观本构精确描述,以及利用分数阶建立的细观本构模型。为解决岩石细观破裂提供了可靠的依据,对研究岩石细观结构有重要意义。
Claims (1)
1.一种描述外载作用岩体微细观破裂与微裂隙形成的方法,其特征在于包括如下步骤:
一、从细化角度研究岩石破裂问题,在岩石微观颗粒之间建立黏结模型,将颗粒简化为圆形颗粒,以两颗粒圆心为端点的建立等效微梁模型:模型长L=RA+RB,高L=2×(RA+RB)/2=RA+RB,宽度或厚度为t;
二、建立基于分数阶的黏壶元件的本构关系:
式中:η为黏性系数,k为常数,Γ(n+1)为n+1的阶乘,ε为应变,σ为应力,d为微分符号,D(·)为Riemann-Liouville型分数阶微积分算子;
三、建立分数阶kelvin模型,由Abel黏壶元件与弹簧元件并联组成,步骤如下:
σk1(t)=Ekεk1(t) (4)
分数阶kelvin模型为:
σk(t)=σk1+σk2=Ekεk1(t)+kDnεk2(t)=(Ek+kDn)εk(t) (6)
式中:Ek为kelvin体的弹性模量,k为常数;
四、建立由弹簧元件与黏壶元件串联的Maxwll模型,步骤如下:
σm1(t)=Emεm1(t) (7)
σm2(t)=ηmDεm2(t) (8)
式中:Em为kelvin体的弹性模量,k为常数;
五、将Maxwell模型与分数阶Kelvin模型串联,应变相加,应力相同,上式两边分别相加,当元件串联时,由ε=ε1+ε2,σ1=σ2,得到σ=σm=σk,ε=εm+εk,得到分数阶Burgers模型本构关系为:
式中:η为黏性系数,k为常数,ε为应变,σ为应力,Ek为kelvin体的弹性模量,Em为Maxwll体的弹性模量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811006299.6A CN108982818B (zh) | 2018-08-30 | 2018-08-30 | 一种描述外载作用岩体微细观破裂与微裂隙形成的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811006299.6A CN108982818B (zh) | 2018-08-30 | 2018-08-30 | 一种描述外载作用岩体微细观破裂与微裂隙形成的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108982818A true CN108982818A (zh) | 2018-12-11 |
CN108982818B CN108982818B (zh) | 2021-04-06 |
Family
ID=64548177
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811006299.6A Expired - Fee Related CN108982818B (zh) | 2018-08-30 | 2018-08-30 | 一种描述外载作用岩体微细观破裂与微裂隙形成的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108982818B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116611269A (zh) * | 2023-07-19 | 2023-08-18 | 中国矿业大学(北京) | 弹黏性元件模型构建方法、求解方法、设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7374659B1 (en) * | 2004-06-22 | 2008-05-20 | Asphalt Technology, Llc. | Methods and systems for modifying asphalts |
CN103942387A (zh) * | 2014-04-16 | 2014-07-23 | 四川大学 | 一种基于变分数阶导数建立岩石蠕变本构模型的新方法 |
CN107290215A (zh) * | 2017-06-23 | 2017-10-24 | 中国矿业大学 | 一种用于涂层织物膜材黏弹性本构行为的预测方法 |
CN107843722A (zh) * | 2017-07-12 | 2018-03-27 | 南通大学 | 一种基于随机有限元的纤维沥青混凝土细观行为分析方法 |
-
2018
- 2018-08-30 CN CN201811006299.6A patent/CN108982818B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7374659B1 (en) * | 2004-06-22 | 2008-05-20 | Asphalt Technology, Llc. | Methods and systems for modifying asphalts |
CN103942387A (zh) * | 2014-04-16 | 2014-07-23 | 四川大学 | 一种基于变分数阶导数建立岩石蠕变本构模型的新方法 |
CN107290215A (zh) * | 2017-06-23 | 2017-10-24 | 中国矿业大学 | 一种用于涂层织物膜材黏弹性本构行为的预测方法 |
CN107843722A (zh) * | 2017-07-12 | 2018-03-27 | 南通大学 | 一种基于随机有限元的纤维沥青混凝土细观行为分析方法 |
Non-Patent Citations (6)
Title |
---|
YAO ZHAOMING等: "Particle Swarm Fractional Order Burgers Model of Accumulated Pore Pressure Subject to Constant Cyclic Loading of Principal Stress Axis", 《2ND INTERNATIONAL CONFERENCE ON ARCHITECTURAL ENGINEERING AND NEW MATERIALS》 * |
何明明等: "基于分数阶微积分岩石的动态变形行为研究", 《岩土工程学报》 * |
刘利: "软岩变形本构特性研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
尹应梅等: "基于动态蠕变实验的沥青混合料本构关系研究", 《功能材料》 * |
浦少云等: "低动应力下岩石分数阶Burgers本构模型", 《长江科学院院报》 * |
王志方等: "采用分数阶导数描述胶凝原油的流变模型", 《中国石油大学学报( 自然科学版) 》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116611269A (zh) * | 2023-07-19 | 2023-08-18 | 中国矿业大学(北京) | 弹黏性元件模型构建方法、求解方法、设备及存储介质 |
CN116611269B (zh) * | 2023-07-19 | 2023-10-13 | 中国矿业大学(北京) | 弹黏性元件模型构建方法、求解方法、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108982818B (zh) | 2021-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shan et al. | Study on the triaxial unloading creep mechanical properties and damage constitutive model of red sandstone containing a single ice-filled flaw | |
Lee et al. | Borehole breakouts induced in arkosic sandstones and a discrete element analysis | |
Rahmati et al. | Review of sand production prediction models | |
CN109374867B (zh) | 一种基于离散元的砂砾岩水力压裂模拟方法 | |
Klerck et al. | Discrete fracture in quasi-brittle materials under compressive and tensile stress states | |
CN110705165A (zh) | 一种构建岩石材料弹塑性-损伤耦合力学本构模型的方法 | |
Lin et al. | Anisotropic characteristic of irregular columnar-jointed rock mass based on physical model test | |
Perkins et al. | Stability and failure of spherical cavities in unconsolidated sand and weakly consolidated rock | |
Raisianzadeh et al. | 2D simulation of breakage of angular particles using combined DEM and XFEM | |
Yang et al. | Model I cohesive zone models of different rank coals | |
CN109357921A (zh) | 一种参数可控的缝洞油***造岩心制作方法 | |
CN108982818A (zh) | 一种描述外载作用岩体微细观破裂与微裂隙形成的方法 | |
González et al. | Pulse fracture simulation in shale rock reservoirs: DEM and FEM–DEM approaches | |
Ma et al. | Influence of pore pressure on tensile fracture growth in rocks: a new explanation based on numerical testing | |
Chazallon et al. | A constitutive model coupling elastoplasticity and damage for cohesive‐frictional materials | |
Wang et al. | Poroelastic versus poroplastic modeling of hydraulic fracturing | |
Wise et al. | Numerical Analysis of Wellbore Integrity and Cement Sheath Debonding for Wells in the Eugene Island OPD, Gulf of Mexico | |
Crawford et al. | Constitutive modeling of deformation and permeability: Relationships between critical state and micromechanics | |
Knauth et al. | Modeling of the mechanical and hydraulic behavior of salt rocks on a discontinuum-mechanical basis | |
Gao et al. | Parallel finite element simulations of 3D hydraulic fracture propagation using a coupled hydro-mechanical interface element | |
Mandl | Modelling incipient tectonic faulting in the brittle crust of the earth | |
He et al. | Failure mode of cement sheath in salt cavern gas storge wellbore based on coupling plasticity and damage evolution | |
Wang et al. | Burgers creep model used for describing and predicting the creep behaviour of a rock under uniaxial and triaxial compression test conditions | |
Liu et al. | Dynamic mechanical properties of jointed soft rock samples subjected to cyclic triaxial loading: a FEM-DEM-based study | |
Palencia Yrausquin et al. | Laboratory study and numerical simulation of mechanism of failure on hydraulic fracturing in unconsolidated sandstones |
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 |
Granted publication date: 20210406 Termination date: 20210830 |
|
CF01 | Termination of patent right due to non-payment of annual fee |