CN106570269A - 一种基于数字图像的页岩破裂模式分析方法 - Google Patents
一种基于数字图像的页岩破裂模式分析方法 Download PDFInfo
- Publication number
- CN106570269A CN106570269A CN201610971453.8A CN201610971453A CN106570269A CN 106570269 A CN106570269 A CN 106570269A CN 201610971453 A CN201610971453 A CN 201610971453A CN 106570269 A CN106570269 A CN 106570269A
- Authority
- CN
- China
- Prior art keywords
- shale
- shales
- heterogeneity
- digital
- digital picture
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30184—Infrastructure
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Silicates, Zeolites, And Molecular Sieves (AREA)
Abstract
本发明提供一种基于数字图像的页岩破裂模式分析方法,在细观尺度上利用数字图像处理技术表征页岩中由于方解石矿物的形状、大小及分布对页岩造成的非均匀性,并映射到有限元网格中,结合基于细观结构的岩石破裂过程分析系统RFPA2D‑DIP建立能准确反映页岩细观结构的数值模型。利用该模型进行页岩单轴压缩数值试验,研究在不同方向加载下页岩细观结构对其抗压强度和破坏模式的影响。以解决现有针对页岩抗压强度及破坏模式的研究,很少涉及页岩细观尺度下自身非均匀性引起的应力分布非均匀性和由此诱发的局部破坏问题,难以得到页岩脆性破坏机制的本质的问题。本发明属于岩石力学领域。
Description
技术领域
本发明涉及一种页岩气勘探开发中基于数字图像的页岩破裂模式的新方法,属于岩石力学领域。
背景技术
针对页岩的抗压强度及破坏模式国内外学者做了比较多的研究,特别是针对页岩在不同围压、不同倾角下的强度、变形与破坏模式等有了很细致的分析,但这些研究都是主要基于页岩的宏观尺度破坏,很少涉及页岩细观尺度下自身非均匀性引起的应力分布非均匀性和由此诱发的局部破坏问题,因此,难以得到页岩脆性破坏机制的本质。
发明内容
本发明的目的在于:提供一种基于数字图像的页岩破裂模式分析方法,以解决现有针对页岩抗压强度及破坏模式的研究,很少涉及页岩细观尺度下自身非均匀性引起的应力分布非均匀性和由此诱发的局部破坏问题,难以得到页岩脆性破坏机制的本质的问题。
为解决上述问题,提供一种基于数字图像的页岩破裂模式分析方法,其步骤如下:
第一步,用数码相机对页岩进行拍照,获取页岩的表面图像;
第二步,数字图像处理中采用阈值分割或边缘检测的方法对具有不同特征的图像区域进行划分,根据不同的图像特征,选择合适的颜色空间进行图像分割;
第三步,采用数字图像处理技术表征页岩中方解石矿物与页岩基质的空间分布;
第四步,在方解石与基质内部采用统计的方法来描述其非均匀性,建立与有限元模型网格数据之间的关系;
第五步,建立页岩细-微观尺度耦合分析的页岩损伤数值模型;
第六步,利用RFPA2D的破坏过程分析优势,模拟页岩在单轴载荷作用下的破坏过程,研究页岩抗压强度及破坏形态与不同方位角的变化。
前述分析方法的第二步中,对于页岩的图像,选择HSI颜色空间进行图像分割。
前述分析方法的第四步中,假定页岩与方解石内部的细观单元的弹性模量和强度服从Weibull分布,
式中a代表细观单元的力学性质参数;a0是细观单元力学性质的平均值;m是分布函数的形状参数,其物理意义反映材料的均质性,定义为材料的均匀性系数;f(a)为细观单元力学性质a的统计分布密度。
前述分析方法中,在获取页岩数字图像过程中,利用50mmx50mm的正方形框架为工具,然后固定正方形的中心,逆时针方向每隔15°获取1幅页岩数字图像,一共获取7幅图像,分别为方位角α=0°、15°、30°、45°、60°、75°、90°的数字图像,研究不同方位角下页岩的强度变形特征。
页岩作为一种非均质材料,其宏观力学响应是其内部各细观介质力学响应的总体外在表现,内部细观介质几何分布控制着材料的破裂演化过程。
与现有技术相比,本发明通过数字图像处理技术表征了页岩中常见矿物方解石和页岩基质的空间分布,并且在方解石与基质内部采用统计的方法来描述其非均匀性,然后结合RFPA岩石破裂过程分析系统软件,模拟了页岩在单轴载荷作用下的破坏过程,研究了页岩抗压强度及破裂模式与不同方位角的变化,并分析了页岩声发射空间演化分布规律和分形维数特征与不同方位角的关系。这不仅对海相页岩岩石力学理论体系的建立具有重要的科学价值,而且对南方地区寒武系页岩气的勘探开发具有重要的指导作用和实践意义。
附图说明
图1是实施例的页岩表面图;
图2是图1中AA’扫描线上I值变化曲线图;
图3是实施例阈值分割后的图像;
图4是实施例页岩数字图像转化后的数值模型与网格;
图5是实施例模型加载示意图;
图6是实施例不同方位角下页岩数字图像;
图7是实施例不同方位角下页岩抗压强度和弹性模量变化趋势图;
图8是实施例不同方位角下页岩破裂模式图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将对本发明作进一步地详细描述,
实施例:
本实施例提供一种基于数字图像的页岩破裂模式分析方法,步骤如下:
第一步,用数码相机对页岩进行拍照,获取页岩的表面图像,在图1中较亮的部分为方解石,较暗的部分为页岩。
第二步,数字图像处理中采用阈值分割或边缘检测的方法对具有不同特征的图像区域进行划分,根据不同的图像特征,选择合适的颜色空间进行图像分割。对于页岩的图像,多为灰度图像,可以选择HSI颜色空间进行图像分割。参照图1中AB线穿过的方解石脉与页岩的位置与图2中曲线的变化情况可以看出,页岩的亮度(I值)在145以下,方解石的亮度在145以上波动,通过反复试验的方法,最终选择I=145作为区域分割的阈值进行图像分割根据反复调试,选择I值为145进行分割。
第三步,采用数字图像处理技术表征页岩中方解石矿物与页岩基质的空间分布。图3为分割后的结果,黄色表示方解石,浅蓝色表示页岩。由图3可以看出,经过数字图像处理后图像比较准确的反映出方解石脉在页岩中的形状和空间分布。
第四步,在方解石与基质内部采用统计的方法来描述其非均匀性,建立与有限元模型网格数据之间的关系。假定页岩与方解石内部的细观单元的弹性模量和强度服从Weibull分布,
式中这里a代表细观单元的力学性质参数(弹性模量、强度等);a0是细观单元力学性质的平均值;m是分布函数的形状参数,其物理意义反映材料的均质性,定义为材料的均匀性系数;f(a)为细观单元力学性质a的统计分布密度。通过数字图像处理后的图像转化成的数值模型如图4所示。图中颜色愈亮表示弹性模量愈大
第五步,建立页岩细-微观尺度耦合分析的页岩损伤数值模型。图5是页岩单轴加载示意图。
第六步,利用RFPA2D的破坏过程分析优势,模拟了页岩在单轴载荷作用下的破坏过程,研究了页岩抗压强度及破坏形态与不同方位角的变化。为了研究页岩中细观结构的影响,在获取页岩数字图像过程中,利用50mmx50mm的正方形框架为工具,然后固定正方形的中心,逆时针方向每隔15°获取1幅页岩数字图像,一共获取7幅图像,分别为方位角α=0°、15°、30°、45°、60°、75°、90°的数字图像,如图6所示,方位角α为方解石脉与水平方向的夹角。试验中,保证页岩试件的细观结构基本相同来研究不同方位角下页岩的强度变形特征。图7为页岩在不同方位角下的抗压强度以及弹性模量变化趋势。图8为不同方位角下页岩破裂模式图。
Claims (4)
1.一种基于数字图像的页岩破裂模式分析方法,其特征在于,步骤如下:
第一步,用数码相机对页岩进行拍照,获取页岩的表面图像;
第二步,数字图像处理中采用阈值分割或边缘检测的方法对具有不同特征的图像区域进行划分,根据不同的图像特征,选择合适的颜色空间进行图像分割;
第三步,采用数字图像处理技术表征页岩中方解石矿物与页岩基质的空间分布;
第四步,在方解石与基质内部采用统计的方法来描述其非均匀性,建立与有限元模型网格数据之间的关系;
第五步,建立页岩细-微观尺度耦合分析的页岩损伤数值模型;
第六步,利用RFPA2D的破坏过程分析优势,模拟页岩在单轴载荷作用下的破坏过程,研究页岩抗压强度及破坏形态与不同方位角的变化。
2.根据权利要求1所述一种基于数字图像的页岩破裂模式分析方法,其特征在于:第二步中,对于页岩的图像,选择HSI颜色空间进行图像分割。
3.根据权利要求1所述一种基于数字图像的页岩破裂模式分析方法,其特征在于:第四步中,假定页岩与方解石内部的细观单元的弹性模量和强度服从Weibull分布,
式中a代表细观单元的力学性质参数;a0是细观单元力学性质的平均值;m是分布函数的形状参数,其物理意义反映材料的均质性,定义为材料的均匀性系数;f(a)为细观单元力学性质a的统计分布密度。
4.根据权利要求1所述一种基于数字图像的页岩破裂模式分析方法,其特征在于:在获取页岩数字图像过程中,利用50mmx50mm的正方形框架为工具,然后固定正方形的中心,逆时针方向每隔15°获取1幅页岩数字图像,一共获取7幅图像,分别为方位角α=0°、15°、30°、45°、60°、75°、90°的数字图像,研究不同方位角下页岩的强度变形特征。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610971453.8A CN106570269A (zh) | 2016-10-27 | 2016-10-27 | 一种基于数字图像的页岩破裂模式分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610971453.8A CN106570269A (zh) | 2016-10-27 | 2016-10-27 | 一种基于数字图像的页岩破裂模式分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106570269A true CN106570269A (zh) | 2017-04-19 |
Family
ID=58541625
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610971453.8A Pending CN106570269A (zh) | 2016-10-27 | 2016-10-27 | 一种基于数字图像的页岩破裂模式分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106570269A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107101888A (zh) * | 2017-05-27 | 2017-08-29 | 青岛理工大学 | 一种岩石压缩破坏预测方法 |
CN107133399A (zh) * | 2017-05-02 | 2017-09-05 | 浙江省交通规划设计研究院 | 抗压强度评估方法及装置 |
CN107192600A (zh) * | 2017-05-18 | 2017-09-22 | 西安科技大学 | 一种分析工程扰动影响下煤岩强度劣化耦合特征的方法 |
CN110765572A (zh) * | 2019-09-12 | 2020-02-07 | 中国科学院武汉岩土力学研究所 | 一种杏仁状玄武岩单三轴试验连续非连续数值模拟方法 |
CN112150430A (zh) * | 2020-09-21 | 2020-12-29 | 中国矿业大学(北京) | 一种利用岩石细观结构数字图像的数值分析方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110082678A1 (en) * | 2009-10-01 | 2011-04-07 | Algive Lionnel | Method of optimizing the injection of a reactive fluid into a porous medium |
US20130090902A1 (en) * | 2010-06-28 | 2013-04-11 | Yao Yao | Method and System for Modeling Fractures in Ductile Rock |
CN104615802A (zh) * | 2014-12-09 | 2015-05-13 | 中国石油集团川庆钻探工程有限公司 | 利用岩屑元素分析数据计算页岩脆性矿物含量的方法 |
CN104749619A (zh) * | 2013-12-26 | 2015-07-01 | 中国石油化工股份有限公司 | 基于方位杨氏模量的地层各向异性优势方向预测方法 |
CN105716947A (zh) * | 2016-01-26 | 2016-06-29 | 河海大学 | 一种岩石的细观蠕变破坏过程的测试方法 |
-
2016
- 2016-10-27 CN CN201610971453.8A patent/CN106570269A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110082678A1 (en) * | 2009-10-01 | 2011-04-07 | Algive Lionnel | Method of optimizing the injection of a reactive fluid into a porous medium |
US20130090902A1 (en) * | 2010-06-28 | 2013-04-11 | Yao Yao | Method and System for Modeling Fractures in Ductile Rock |
CN104749619A (zh) * | 2013-12-26 | 2015-07-01 | 中国石油化工股份有限公司 | 基于方位杨氏模量的地层各向异性优势方向预测方法 |
CN104615802A (zh) * | 2014-12-09 | 2015-05-13 | 中国石油集团川庆钻探工程有限公司 | 利用岩屑元素分析数据计算页岩脆性矿物含量的方法 |
CN105716947A (zh) * | 2016-01-26 | 2016-06-29 | 河海大学 | 一种岩石的细观蠕变破坏过程的测试方法 |
Non-Patent Citations (2)
Title |
---|
于庆磊 等: "《基于数字图像的混凝土破坏过程的数值模拟》", 《工程力学》 * |
于庆磊 等: "《岩石细观结构对其变形强度影响的数值分析》", 《岩石力学》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107133399A (zh) * | 2017-05-02 | 2017-09-05 | 浙江省交通规划设计研究院 | 抗压强度评估方法及装置 |
CN107133399B (zh) * | 2017-05-02 | 2020-08-14 | 浙江省交通规划设计研究院 | 抗压强度评估方法及装置 |
CN107192600A (zh) * | 2017-05-18 | 2017-09-22 | 西安科技大学 | 一种分析工程扰动影响下煤岩强度劣化耦合特征的方法 |
CN107192600B (zh) * | 2017-05-18 | 2019-08-16 | 西安科技大学 | 一种分析工程扰动影响下煤岩强度劣化耦合特征的方法 |
CN107101888A (zh) * | 2017-05-27 | 2017-08-29 | 青岛理工大学 | 一种岩石压缩破坏预测方法 |
CN110765572A (zh) * | 2019-09-12 | 2020-02-07 | 中国科学院武汉岩土力学研究所 | 一种杏仁状玄武岩单三轴试验连续非连续数值模拟方法 |
CN112150430A (zh) * | 2020-09-21 | 2020-12-29 | 中国矿业大学(北京) | 一种利用岩石细观结构数字图像的数值分析方法 |
CN112150430B (zh) * | 2020-09-21 | 2023-10-10 | 中国矿业大学(北京) | 一种利用岩石细观结构数字图像的数值分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106570269A (zh) | 一种基于数字图像的页岩破裂模式分析方法 | |
Li | Research progress on evaluation methods and factors influencing shale brittleness: A review | |
Zhang et al. | Cracking processes in rock-like material containing a single flaw under uniaxial compression: a numerical study based on parallel bonded-particle model approach | |
Wu et al. | Numerical simulation and fractal analysis of mesoscopic scale failure in shale using digital images | |
Semnani et al. | Quantifying the heterogeneity of shale through statistical combination of imaging across scales | |
Karakul et al. | Empirical models and numerical analysis for assessing strength anisotropy based on block punch index and uniaxial compression tests | |
Zhao et al. | Study on failure characteristic of rock‐like materials with an open‐hole under uniaxial compression | |
WO2021022465A1 (zh) | 一种基于数字图像技术的粗糙面裂隙生成方法及试验系统 | |
Pan et al. | Study on the effect of micro-geometric heterogeneity on mechanical properties of brittle rock using a grain-based discrete element method coupling with the cohesive zone model | |
Yang et al. | Numerical study on failure behavior of brittle rock specimen containing pre-existing combined flaws under different confining pressure | |
Qiu et al. | Dynamic failure of a phyllite with a low degree of metamorphism under impact Brazilian test | |
Miller et al. | Understanding local deformation in metallic polycrystals using high energy X-rays and finite elements | |
Xiao et al. | Texture synthesis: A novel method for generating digital models with heterogeneous diversity of rock materials and its CGM verification | |
CN103076347B (zh) | 基于原位x射线断层照相的脆性材料力学损伤的测量方法 | |
Zhao et al. | Establishment of numerical cracking constitutive models using 3D reconstruction and X-ray CT images of geomaterials | |
Aligholi et al. | Estimating engineering properties of igneous rocks using semi-automatic petrographic analysis | |
Marinelli et al. | Parameter calibration for high-porosity sandstones deformed in the compaction banding regime | |
Fan et al. | A real-time visual investigation on microscopic progressive fatigue deterioration of granite under cyclic loading | |
Zhou et al. | An innovative acousto-optic-sensing-based triaxial testing system for rocks | |
Cheng et al. | Free desiccation shrinkage process in clayey soils | |
Tian et al. | Quantification of rock heterogeneity and application in predicting rock mechanical properties | |
Liu et al. | Elastic characteristics of digital cores from Longmaxi shale using lattice spring models | |
Chajed et al. | Acoustic Emission (AE) Based Damage Quantification and Its Relation with AE-Based Micromechanical Coupled Damage Plasticity Model for Intact Rocks | |
Wu et al. | Study on a new inversion method for non-uniform distribution of rock material parameters | |
Wu et al. | Damage evolution characteristics of 3D reconstructed bedding-containing shale based on CT technology and digital image processing |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170419 |
|
RJ01 | Rejection of invention patent application after publication |