CN112986247B - 基于矿物含量和分布的岩石脆性分析方法、系统及设备 - Google Patents

基于矿物含量和分布的岩石脆性分析方法、系统及设备 Download PDF

Info

Publication number
CN112986247B
CN112986247B CN202110469932.0A CN202110469932A CN112986247B CN 112986247 B CN112986247 B CN 112986247B CN 202110469932 A CN202110469932 A CN 202110469932A CN 112986247 B CN112986247 B CN 112986247B
Authority
CN
China
Prior art keywords
rock
distribution
analyzed
autocorrelation
brittleness
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
Application number
CN202110469932.0A
Other languages
English (en)
Other versions
CN112986247A (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.)
Institute of Geology and Geophysics of CAS
Original Assignee
Institute of Geology and Geophysics of CAS
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 Institute of Geology and Geophysics of CAS filed Critical Institute of Geology and Geophysics of CAS
Priority to CN202110469932.0A priority Critical patent/CN112986247B/zh
Publication of CN112986247A publication Critical patent/CN112986247A/zh
Application granted granted Critical
Publication of CN112986247B publication Critical patent/CN112986247B/zh
Priority to US17/471,195 priority patent/US11282183B1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84—Systems specially adapted for particular applications
    • 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
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/50—Processing the detected response signal, e.g. electronic circuits specially adapted therefor using auto-correlation techniques or cross-correlation techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/01—Indexing codes associated with the measuring variable
    • G01N2291/011—Velocity or travel time
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/023—Solids
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/023—Solids
    • G01N2291/0232—Glass, ceramics, concrete or stone
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/025—Change of phase or condition
    • G01N2291/0258—Structural degradation, e.g. fatigue of composites, ageing of oils
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/028—Material parameters
    • G01N2291/02827—Elastic parameters, strength or force
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/04—Wave modes and trajectories
    • G01N2291/042—Wave modes
    • G01N2291/0421—Longitudinal waves
    • 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/20—Special algorithmic details
    • G06T2207/20021—Dividing image into blocks, subimages or windows
    • 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

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明属于介质脆性分布表征领域，具体涉及了一种基于矿物含量和分布的岩石脆性分析方法、系统及设备，旨在解决现有岩石脆性分析方法仅考虑内部矿物含量而忽略矿物分布导致分析结果的准确性和精度达不到预期的问题。本发明包括：对岩石图像数字化，并对数字化后的二维数组进行自相关分析；以自相关相关系数变化曲线的纵坐标值下降为1/e时的横坐标值作为自相关长度；构建描述矿物分布的自相关函数；进行岩石样品重构，获得密度及弹性模量分布；进行网格剖分，获得待分析岩石的不同尺度剖分的脆性指数分布。本发明充分考虑到岩石内部颗粒相对位置、孔隙结构变化对岩石脆性特征性质产生的影响，岩石脆性分析结果的精度和准确性更高。

Description

基于矿物含量和分布的岩石脆性分析方法、系统及设备
技术领域
本发明属于介质脆性分布表征领域，具体涉及了一种基于矿物含量和分布的岩石脆性分析方法、系统及设备。
背景技术
随着油气开发技术研究的逐渐深入，以及非常规油气资源需求的日益增加，对地下介质压裂的需求不断增加。由于岩石压裂过程中，压裂液的注入位置与岩石的脆性特征密切相关，因此需要深入研究岩石的脆性评价标准，为岩石压裂位置的选取提供依据。岩石的脆性特征受岩石中，所含矿物种类、矿物分布，影响较大。因此，研究真实岩石的脆性特征时，需要充分计算矿物类型、矿物分布对岩石弹性性质的影响。然而，常规的岩石脆性分析方法只计算了岩石中，矿物种类对岩石脆性特征的影响，岩石矿物分布差异导致的岩石脆性分布差异，以及差异对岩石脆性特征的影响并未深入研究。
自相关函数是描述物质随机分布状态的有效方法，通过对其功率谱函数的控制，能够表征不同的岩石矿物分布，从而定量的计算矿物分布对岩石各种特征的影响。数字岩芯技术，通过对岩石进行CT扫描，能够显示岩石内部结构。通过数字岩芯技术，可研究岩石内部矿物分布，并且，通过数字岩芯技术与岩石脆性指数表征算法的结合，可反应岩石内部矿物分布对岩石脆性空间分布的影响。在岩石矿物非均匀分布的条件下，研究描述矿物分布特征的自相关函数，并计算脆性指数的分布特征，是研究岩石脆性分布的重要方法，是选择岩石压裂位置的重要依据。
但是，目前岩石脆性的特征的相关分析，只考虑了矿物含量的影响，对矿物分布的影响以及定量描述的有效方法仍然缺失，如何在结合定量表征矿物分布的基础上进行岩石脆性分析，仍是本领域需要研究的问题。
发明内容
为了解决现有技术中的上述问题，即现有岩石脆性分析方法仅考虑内部矿物含量而忽略矿物分布导致分析结果的准确性和精度达不到预期的问题，本发明提供了一种基于矿物含量和分布的岩石脆性分析方法，该方法包括：
步骤S10，进行待分析岩石图像的数字化处理，并进行数字化处理后的二维数组的自相关分析，获得自相关相关系数变化曲线；
步骤S20，将所述自相关相关系数变化曲线的纵坐标值下降为1/e时的横坐标值作为待分析岩石的自相关长度；
步骤S30，根据所述待分析岩石的自相关长度，待分析岩石中矿物以及弹性模量的分布特征，构建描述矿物分布的自相关函数；
步骤S40，根据所述自相关函数进行所述待分析岩石的岩石样品重构，获得待分析岩石的密度及弹性模量分布；
步骤S50，根据所述待分析岩石的密度及弹性模量分布进行重构的岩石样品的网格剖分，获得待分析岩石的不同尺度剖分的脆性指数分布。
在一些优选的实施例中，所述描述矿物分布的自相关函数，其公式表示为：
Figure 100002_DEST_PATH_IMAGE001
其中，
Figure 592709DEST_PATH_IMAGE002
为待分析岩石的自相关长度；
Figure 47217DEST_PATH_IMAGE003
为数据点
Figure 283027DEST_PATH_IMAGE004
与岩芯切片中心
Figure 762550DEST_PATH_IMAGE005
的相对距离，
Figure 809134DEST_PATH_IMAGE006
分别为数据点
Figure 639687DEST_PATH_IMAGE004
的横坐标值和纵坐标值。
在一些优选的实施例中，所述待分析岩石的密度及弹性模量分布，其公式表示为：
Figure 944766DEST_PATH_IMAGE007
其中，
Figure 278796DEST_PATH_IMAGE008
为数据点
Figure 260396DEST_PATH_IMAGE004
的密度，
Figure 578245DEST_PATH_IMAGE009
为待分析岩石的背景密度；
Figure 562381DEST_PATH_IMAGE010
为数据点
Figure 875551DEST_PATH_IMAGE004
的体积模量，
Figure 919730DEST_PATH_IMAGE011
为体积模量背景值；
Figure 334662DEST_PATH_IMAGE012
为数据点
Figure 856910DEST_PATH_IMAGE004
的剪切模量，
Figure 290166DEST_PATH_IMAGE013
为剪切模量背 景值；
Figure 505247DEST_PATH_IMAGE006
分别为数据点
Figure 532108DEST_PATH_IMAGE004
的横坐标值和纵坐标值；
Figure 957581DEST_PATH_IMAGE014
为扰动量，
Figure 386288DEST_PATH_IMAGE015
为岩石非均质性产生的 密度扰动，
Figure 896904DEST_PATH_IMAGE016
为岩石非均质性产生的体积模量扰动，
Figure 145482DEST_PATH_IMAGE017
为岩石非均质性产生的剪切模量扰 动；
Figure 275112DEST_PATH_IMAGE018
代表求平均值操作。
在一些优选的实施例中，所述剪切模量背景值
Figure 433692DEST_PATH_IMAGE013
和所述体积模量背景值
Figure 724996DEST_PATH_IMAGE011
通过新 型孔隙声弹性公式计算获得，其计算方法为：
Figure 585505DEST_PATH_IMAGE019
Figure 518826DEST_PATH_IMAGE020
其中，
Figure DEST_PATH_IMAGE021
和
Figure 764868DEST_PATH_IMAGE022
分别为通过岩石物理实验测量得到的待分析岩石的纵波速度和横波 速度。
在一些优选的实施例中，通过扰动的方差和协方差确定所述扰动量
Figure 227074DEST_PATH_IMAGE014
，所述扰动的 方差和协方差的关系为：
Figure DEST_PATH_IMAGE023
其中，
Figure 43720DEST_PATH_IMAGE024
为扰动的协方差，
Figure 390519DEST_PATH_IMAGE025
为扰动的方差，
Figure 382746DEST_PATH_IMAGE026
为描述矿物分布的自相关 函数。
在一些优选的实施例中，步骤S50中每一个剖分网格，其脆性指数计算方法为：
Figure DEST_PATH_IMAGE027
Figure 140486DEST_PATH_IMAGE028
Figure 850953DEST_PATH_IMAGE029
其中，
Figure 971749DEST_PATH_IMAGE030
、
Figure 84062DEST_PATH_IMAGE031
、
Figure 278283DEST_PATH_IMAGE032
和
Figure 476046DEST_PATH_IMAGE033
分别代表每一个剖分网格的四个节点中的最大杨氏模 量、最小杨氏模量、平均杨氏模量和该网格杨氏模量值的相对差异，
Figure 899068DEST_PATH_IMAGE034
、
Figure 131466DEST_PATH_IMAGE035
、
Figure 371955DEST_PATH_IMAGE036
和
Figure 916069DEST_PATH_IMAGE037
分别 代表每一个剖分网格的四个节点中的最大泊松比、最小泊松比、平均泊松比和该网格泊松 比的相对差异，
Figure 267416DEST_PATH_IMAGE038
代表该节点的脆性指数。
在一些优选的实施例中，步骤S50中每一个剖分网格，其平均杨氏模量
Figure 728222DEST_PATH_IMAGE032
和平均泊 松比
Figure 139611DEST_PATH_IMAGE036
的计算方法为：
Figure 171021DEST_PATH_IMAGE039
Figure 326059DEST_PATH_IMAGE040
其中，
Figure 267470DEST_PATH_IMAGE041
、
Figure 725128DEST_PATH_IMAGE042
、
Figure 384779DEST_PATH_IMAGE043
和
Figure 202562DEST_PATH_IMAGE044
分别为剖分网格的四个节点
Figure 998480DEST_PATH_IMAGE045
、
Figure 751672DEST_PATH_IMAGE046
、
Figure DEST_PATH_IMAGE047
和
Figure 478713DEST_PATH_IMAGE048
的杨氏模量，
Figure DEST_PATH_IMAGE049
、
Figure 834608DEST_PATH_IMAGE050
、
Figure 360399DEST_PATH_IMAGE051
和
Figure 284492DEST_PATH_IMAGE052
分别为剖分网格的四个节点
Figure 653157DEST_PATH_IMAGE045
、
Figure 812743DEST_PATH_IMAGE046
、
Figure 583253DEST_PATH_IMAGE047
和
Figure 848887DEST_PATH_IMAGE048
的泊松比。
本发明的另一方面，提出了一种基于矿物含量和分布的岩石脆性分析系统，该系统包括以下模块：
数字化处理模块，配置为进行待分析岩石图像的数字化处理，获得待分析岩石的二维数组；
自相关分析模块，配置为进行所述二维数组的自相关分析，并将分析获取的自相关相关系数变化曲线的纵坐标值下降为1/e时的横坐标值作为待分析岩石的自相关长度；
自相关函数构建模块，配置为根据所述待分析岩石的自相关长度，待分析岩石中矿物以及弹性模量的分布特征，构建描述矿物分布的自相关函数；
密度及弹性模量分布计算模块，配置为根据所述自相关函数进行所述待分析岩石的岩石样品重构，获得待分析岩石的密度及弹性模量分布；
脆性分析模块，配置为根据所述待分析岩石的密度及弹性模量分布进行重构的岩石样品的网格剖分，获得待分析岩石的不同尺度剖分的脆性指数分布。
本发明的第三方面，提出了一种电子设备，包括：
至少一个处理器；以及
与至少一个所述处理器通信连接的存储器；其中，
所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的基于矿物含量和分布的岩石脆性分析方法。
本发明的第四方面，提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现上述的基于矿物含量和分布的岩石脆性分析方法。
本发明的有益效果：
（1）本发明基于矿物含量和分布的岩石脆性分析方法，通过提取岩石中的矿物分布，可研究岩石中矿物分布对岩石脆性指数的影响，从而提升后续岩石脆性分析的准确性和精度。
（2）本发明基于矿物含量和分布的岩石脆性分析方法，利用数字岩芯技术，通过提取岩石矿物分布结构，得到岩石内部弹性模量扰动分布情况，并以扰动量的形式，叠加在声弹性效应引起的等效弹性模量变化中，建立更加真实的深部储层数值模型，从而进一步提升了后续岩石脆性分析的准确性和精度。
（3）本发明基于矿物含量和分布的岩石脆性分析方法，能够得到更为准确和高精度的岩石脆性指数分布的空间信息，从而能够为岩石压裂位置的选址提供更精确的信息。
附图说明
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述，本申请的其它特征、目的和优点将会变得更明显：
图1是本发明基于矿物含量和分布的岩石脆性分析方法的流程示意图；
图2是本发明基于矿物含量和分布的岩石脆性分析方法一种实施例的岩石图像；
图3是本发明基于矿物含量和分布的岩石脆性分析方法一种实施例的提取和理论计算的自相关系数变化曲线示意图；
图4是本发明基于矿物含量和分布的岩石脆性分析方法一种实施例的密度及弹性模量分布示意图；
图5是本发明基于矿物含量和分布的岩石脆性分析方法一种实施例的网格剖分时每个网格的尺度；
图6是本发明基于矿物含量和分布的岩石脆性分析方法一种实施例的脆性指数分布示意图。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅用于解释相关发明，而非对该发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与有关发明相关的部分。
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
本发明提供一种基于矿物含量和分布的岩石脆性分析方法，本方法充分考虑到岩石内部颗粒相对位置、孔隙结构变化对岩石脆性特征性质产生的影响，可以在给矿物分布的条件下计算岩石脆性指数分布、岩石脆性非均质特征，为岩石的压裂位置选择奠定基础。
本发明的一种基于矿物含量和分布的岩石脆性分析方法，该方法包括：
步骤S10，进行待分析岩石图像的数字化处理，并进行数字化处理后的二维数组的自相关分析，获得自相关相关系数变化曲线；
步骤S20，将所述自相关相关系数变化曲线的纵坐标值下降为1/e时的横坐标值作为待分析岩石的自相关长度；
步骤S30，根据所述待分析岩石的自相关长度，待分析岩石中矿物以及弹性模量的分布特征，构建描述矿物分布的自相关函数；
步骤S40，根据所述自相关函数进行所述待分析岩石的岩石样品重构，获得待分析岩石的密度及弹性模量分布；
步骤S50，根据所述待分析岩石的密度及弹性模量分布进行重构的岩石样品的网格剖分，获得待分析岩石的不同尺度剖分的脆性指数分布。
为了更清晰地对本发明基于矿物含量和分布的岩石脆性分析方法进行说明，下面结合图1对本发明实施例中各步骤展开详述。
本发明第一实施例的基于矿物含量和分布的岩石脆性分析方法，包括步骤S10-步骤S50，各步骤详细描述如下：
步骤S10，进行待分析岩石图像的数字化处理，并进行数字化处理后的二维数组的自相关分析，获得自相关相关系数变化曲线。
如图2所示，为本发明基于矿物含量和分布的岩石脆性分析方法一种实施例的岩石图像，对该岩石图像进行数字化处理，形成二维数组，并进行该二维数组的自相关分析，获得岩石图像的自相关系数变化曲线，0.5mm代表岩石图像的比例尺。
步骤S20，将所述自相关相关系数变化曲线的纵坐标值下降为1/e时的横坐标值作为待分析岩石的自相关长度。
根据自相关系数变化曲线的横坐标和纵坐标值获取岩石的自相关长度，选择曲线纵坐标值下降为1/e（e≈2.71828）时，曲线的横坐标值为岩石的自相关长度。
如图3所示，为本发明基于矿物含量和分布的岩石脆性分析方法一种实施例的提 取和理论计算的自相关系数变化曲线示意图，虚线为根据图2的岩石图像提取的自相关系 数变化曲线，实线为根据真实自相关系数变化曲线拟合得到的理论自相关系数变化曲线， 横坐标和纵坐标的单位均为μm，0.368（1/e≈1/2.71828≈0.368）为曲线纵坐标值下降为1/ e的纵坐标值，
Figure 704847DEST_PATH_IMAGE053
139.4μm代表根据自相关系数变化曲线纵坐标值下降为1/e获取的横坐标 值，即岩石的自相关长度。
步骤S30，根据所述待分析岩石的自相关长度，待分析岩石中矿物以及弹性模量的分布特征，构建描述矿物分布的自相关函数，如式（1）所示：
Figure 684436DEST_PATH_IMAGE001
（1）
其中，
Figure 43873DEST_PATH_IMAGE002
为待分析岩石的自相关长度；
Figure 434403DEST_PATH_IMAGE003
为数据点
Figure 512080DEST_PATH_IMAGE004
与岩芯切片中心
Figure 419993DEST_PATH_IMAGE005
的相对距离，
Figure 276347DEST_PATH_IMAGE006
分别为数据点
Figure 713145DEST_PATH_IMAGE004
的横坐标值和纵坐标值。
步骤S40，根据所述自相关函数进行所述待分析岩石的岩石样品重构，获得待分析岩石的密度及弹性模量分布。
根据步骤S10和S20中提取的矿物以及弹性模量分布，通过空间平均，计算空间背景值，基于步骤S30得到的自相关函数，得到待分析岩石的密度及弹性模量的分布函数，如式（2）所示：
Figure 668331DEST_PATH_IMAGE007
（2）
其中，
Figure 848777DEST_PATH_IMAGE008
为数据点
Figure 58173DEST_PATH_IMAGE004
的密度，
Figure 665872DEST_PATH_IMAGE009
为待分析岩石的背景密度；
Figure 842775DEST_PATH_IMAGE010
为数据点
Figure 561332DEST_PATH_IMAGE004
的体积模量，
Figure 389349DEST_PATH_IMAGE011
为体积模量背景值；
Figure 167949DEST_PATH_IMAGE012
为数据点
Figure 973094DEST_PATH_IMAGE004
的剪切模量，
Figure 619976DEST_PATH_IMAGE013
为剪切模量背 景值；
Figure 928597DEST_PATH_IMAGE006
分别为数据点
Figure 753465DEST_PATH_IMAGE004
的横坐标值和纵坐标值；
Figure 780327DEST_PATH_IMAGE014
为扰动量，
Figure 230900DEST_PATH_IMAGE015
为岩石非均质性产生的 密度扰动，
Figure 394028DEST_PATH_IMAGE016
为岩石非均质性产生的体积模量扰动，
Figure 780010DEST_PATH_IMAGE017
为岩石非均质性产生的剪切模量扰 动；
Figure 682718DEST_PATH_IMAGE018
代表求平均值操作。
如图4所示，为本发明基于矿物含量和分布的岩石脆性分析方法一种实施例的密度及弹性模量分布示意图，对应于图2的岩石图像，横坐标和纵坐标代表待分析岩石的密度及弹性模量分布示意图的宽和高，其单位均为mm。
剪切模量背景值
Figure 546768DEST_PATH_IMAGE013
和所述体积模量背景值
Figure 689037DEST_PATH_IMAGE011
通过新型孔隙声弹性公式计算获得， 其计算方法如式（3）和式（4）所示：
Figure 511499DEST_PATH_IMAGE054
（3）
Figure 122740DEST_PATH_IMAGE020
（4）
其中，
Figure 790482DEST_PATH_IMAGE021
和
Figure 787257DEST_PATH_IMAGE022
分别为通过岩石物理实验测量得到的待分析岩石的纵波速度和横波 速度。
通过扰动的方差和协方差确定所述扰动量
Figure 780621DEST_PATH_IMAGE014
，所述扰动的方差和协方差的关系如 式（5）所示：
Figure 3792DEST_PATH_IMAGE023
（5）
其中，
Figure 583546DEST_PATH_IMAGE024
为扰动的协方差，
Figure 841352DEST_PATH_IMAGE025
为扰动的方差，
Figure 864672DEST_PATH_IMAGE026
为描述矿物分布的自相关 函数。本发明一个实施例中，对于中等分选的砂岩，扰动方差选取35%。
步骤S50，根据所述待分析岩石的密度及弹性模量分布进行重构的岩石样品的网格剖分，获得待分析岩石的不同尺度剖分的脆性指数分布。
每一个剖分网格，其脆性指数计算方法如式（6）-式（8）所示：
Figure 309560DEST_PATH_IMAGE027
（6）
Figure 584683DEST_PATH_IMAGE028
（7）
Figure 572362DEST_PATH_IMAGE029
（8）
其中，
Figure 376370DEST_PATH_IMAGE030
、
Figure 964346DEST_PATH_IMAGE031
、
Figure 512002DEST_PATH_IMAGE032
和
Figure 744400DEST_PATH_IMAGE033
分别代表每一个剖分网格的四个节点中的最大杨氏模 量、最小杨氏模量、平均杨氏模量和该网格杨氏模量值的相对差异，
Figure 96141DEST_PATH_IMAGE034
、
Figure 781200DEST_PATH_IMAGE035
、
Figure 257181DEST_PATH_IMAGE036
和
Figure 344085DEST_PATH_IMAGE037
分别 代表每一个剖分网格的四个节点中的最大泊松比、最小泊松比、平均泊松比和该网格泊松 比的相对差异，
Figure 630841DEST_PATH_IMAGE038
代表该节点的脆性指数。
每一个剖分网格，其平均杨氏模量
Figure 537618DEST_PATH_IMAGE032
和平均泊松比
Figure 692655DEST_PATH_IMAGE036
的计算方法如式（9）和式（10） 所示：
Figure 758700DEST_PATH_IMAGE039
（9）
Figure 75412DEST_PATH_IMAGE040
（10）
其中，
Figure 843386DEST_PATH_IMAGE041
、
Figure 802114DEST_PATH_IMAGE042
、
Figure 722666DEST_PATH_IMAGE043
和
Figure 475858DEST_PATH_IMAGE044
分别为剖分网格的四个节点
Figure 232593DEST_PATH_IMAGE045
、
Figure 729433DEST_PATH_IMAGE046
、
Figure 379857DEST_PATH_IMAGE047
和
Figure 428585DEST_PATH_IMAGE048
的杨氏模量，
Figure 62829DEST_PATH_IMAGE049
、
Figure 474612DEST_PATH_IMAGE050
、
Figure 979542DEST_PATH_IMAGE051
和
Figure 199171DEST_PATH_IMAGE052
分别为剖分网格的四个节点
Figure 320711DEST_PATH_IMAGE045
、
Figure 893775DEST_PATH_IMAGE046
、
Figure 394157DEST_PATH_IMAGE047
和
Figure 925633DEST_PATH_IMAGE048
的泊松比。
如图5所示，为是本发明基于矿物含量和分布的岩石脆性分析方法一种实施例的 网格剖分时每个网格的尺度，网格中点
Figure 127944DEST_PATH_IMAGE055
代表该网格的随性指数
Figure 770278DEST_PATH_IMAGE056
，a 1、a 2、a 3、a 4为剖分网格 角点的脆性指数。
如图6所示，为本发明基于矿物含量和分布的岩石脆性分析方法一种实施例的脆性指数分布示意图，对应于图2的岩石图像，横坐标和纵坐标代表待分析岩石的脆性指数分布示意图的宽和高，其单位均为mm。
上述实施例中虽然将各个步骤按照上述先后次序的方式进行了描述，但是本领域技术人员可以理解，为了实现本实施例的效果，不同的步骤之间不必按照这样的次序执行，其可以同时(并行)执行或以颠倒的次序执行，这些简单的变化都在本发明的保护范围之内。
本发明第二实施例的基于矿物含量和分布的岩石脆性分析系统，该系统包括以下模块：
数字化处理模块，配置为进行待分析岩石图像的数字化处理，获得待分析岩石的二维数组；
自相关分析模块，配置为进行所述二维数组的自相关分析，并将分析获取的自相关相关系数变化曲线的纵坐标值下降为1/e时的横坐标值作为待分析岩石的自相关长度；
自相关函数构建模块，配置为根据所述待分析岩石的自相关长度，待分析岩石中矿物以及弹性模量的分布特征，构建描述矿物分布的自相关函数；
密度及弹性模量分布计算模块，配置为根据所述自相关函数进行所述待分析岩石的岩石样品重构，获得待分析岩石的密度及弹性模量分布；
脆性分析模块，配置为根据所述待分析岩石的密度及弹性模量分布进行重构的岩石样品的网格剖分，获得待分析岩石的不同尺度剖分的脆性指数分布。
所属技术领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统的具体工作过程及有关说明，可以参考前述方法实施例中的对应过程，在此不再赘述。
需要说明的是，上述实施例提供的基于矿物含量和分布的岩石脆性分析系统，仅以上述各功能模块的划分进行举例说明，在实际应用中，可以根据需要而将上述功能分配由不同的功能模块来完成，即将本发明实施例中的模块或者步骤再分解或者组合，例如，上述实施例的模块可以合并为一个模块，也可以进一步拆分成多个子模块，以完成以上描述的全部或者部分功能。对于本发明实施例中涉及的模块、步骤的名称，仅仅是为了区分各个模块或者步骤，不视为对本发明的不当限定。
本发明第三实施例的一种电子设备，包括：
至少一个处理器；以及
与至少一个所述处理器通信连接的存储器；其中，
所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的基于矿物含量和分布的岩石脆性分析方法。
本发明第四实施例的一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现上述的基于矿物含量和分布的岩石脆性分析方法。
所属技术领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的存储装置、处理装置的具体工作过程及有关说明，可以参考前述方法实施例中的对应过程，在此不再赘述。
本领域技术人员应该能够意识到，结合本文中所公开的实施例描述的各示例的模块、方法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，软件模块、方法步骤对应的程序可以置于随机存储器（RAM）、内存、只读存储器（ROM）、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。为了清楚地说明电子硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以电子硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。本领域技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
术语“第一”、 “第二”等是用于区别类似的对象，而不是用于描述或表示特定的顺序或先后次序。
术语“包括”或者任何其它类似用语旨在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备/装置不仅包括那些要素，而且还包括没有明确列出的其它要素，或者还包括这些过程、方法、物品或者设备/装置所固有的要素。
至此，已经结合附图所示的优选实施方式描述了本发明的技术方案，但是，本领域技术人员容易理解的是，本发明的保护范围显然不局限于这些具体实施方式。在不偏离本发明的原理的前提下，本领域技术人员可以对相关技术特征做出等同的更改或替换，这些更改或替换之后的技术方案都将落入本发明的保护范围之内。

Claims (8)

1.一种基于矿物含量和分布的岩石脆性分析方法，其特征在于，该方法包括：
步骤S10，进行待分析岩石图像的数字化处理，并进行数字化处理后的二维数组的自相关分析，获得自相关相关系数变化曲线；
步骤S20，将所述自相关相关系数变化曲线的纵坐标值下降为1/e时的横坐标值作为待分析岩石的自相关长度；
步骤S30，根据所述待分析岩石的自相关长度，待分析岩石中矿物以及弹性模量的分布特征，构建描述矿物分布的自相关函数：
Figure DEST_PATH_IMAGE001
其中，
Figure DEST_PATH_IMAGE002
为待分析岩石的自相关长度；
Figure 496609DEST_PATH_IMAGE003
为数据点
Figure DEST_PATH_IMAGE004
与岩芯切片中心
Figure 820274DEST_PATH_IMAGE005
的相对距离，
Figure DEST_PATH_IMAGE006
分别为数据点
Figure 272115DEST_PATH_IMAGE004
的横坐标值和纵坐标值；
步骤S40，根据所述自相关函数进行所述待分析岩石的岩石样品重构，获得待分析岩石的密度及弹性模量分布：
Figure 932904DEST_PATH_IMAGE007
其中，
Figure DEST_PATH_IMAGE008
为数据点
Figure 278434DEST_PATH_IMAGE004
的密度，
Figure 507421DEST_PATH_IMAGE009
为待分析岩石的背景密度；
Figure DEST_PATH_IMAGE010
为数据点
Figure 305613DEST_PATH_IMAGE004
的体积模量，
Figure 146924DEST_PATH_IMAGE011
为体积模量背景值；
Figure DEST_PATH_IMAGE012
为数据点
Figure 612540DEST_PATH_IMAGE004
的剪切模量，
Figure 746849DEST_PATH_IMAGE013
为剪切模量背景值；
Figure 766758DEST_PATH_IMAGE006
分别为数据点
Figure 769349DEST_PATH_IMAGE004
的横坐标值和纵坐标值；
Figure DEST_PATH_IMAGE014
为岩石非均质性产生的密度扰动，
Figure 964838DEST_PATH_IMAGE015
为岩石非均质性产生的体积模量扰动，
Figure DEST_PATH_IMAGE016
为岩石非均质性产生的剪切模量扰动；
Figure 66786DEST_PATH_IMAGE017
代表求平均值操作；
步骤S50，根据所述待分析岩石的密度及弹性模量分布进行重构的岩石样品的网格剖分，获得待分析岩石的不同尺度剖分的脆性指数分布。
2.根据权利要求1所述的基于矿物含量和分布的岩石脆性分析方法，其特征在于，所述剪切模量背景值
Figure 573991DEST_PATH_IMAGE013
和所述体积模量背景值
Figure 114694DEST_PATH_IMAGE011
通过新型孔隙声弹性公式计算获得，其计算方法为：
Figure DEST_PATH_IMAGE018
Figure 663224DEST_PATH_IMAGE019
其中，
Figure DEST_PATH_IMAGE020
和
Figure 529549DEST_PATH_IMAGE021
分别为通过岩石物理实验测量得到的待分析岩石的纵波速度和横波速度。
3.根据权利要求1所述的基于矿物含量和分布的岩石脆性分析方法，其特征在于，通过扰动的方差和协方差确定扰动量，所述扰动量包括岩石非均质性产生的密度扰动
Figure 524050DEST_PATH_IMAGE014
，岩石非均质性产生的体积模量扰动
Figure 743810DEST_PATH_IMAGE015
，以及岩石非均质性产生的剪切模量扰动
Figure 507367DEST_PATH_IMAGE016
，所述扰动量的方差和协方差的关系为：
Figure DEST_PATH_IMAGE022
其中，
Figure 810172DEST_PATH_IMAGE023
为扰动的协方差，
Figure DEST_PATH_IMAGE024
为扰动的方差，
Figure 432914DEST_PATH_IMAGE025
为描述矿物分布的自相关函数。
4.根据权利要求1所述的基于矿物含量和分布的岩石脆性分析方法，其特征在于，步骤S50中每一个剖分网格，其脆性指数计算方法为：
Figure DEST_PATH_IMAGE026
Figure 721944DEST_PATH_IMAGE027
Figure DEST_PATH_IMAGE028
其中，
Figure 605587DEST_PATH_IMAGE029
、
Figure DEST_PATH_IMAGE030
、
Figure DEST_PATH_IMAGE031
和
Figure DEST_PATH_IMAGE032
分别代表每一个剖分网格的四个节点中的最大杨氏模量、最小杨氏模量、平均杨氏模量和该网格杨氏模量值的相对差异，
Figure DEST_PATH_IMAGE033
、
Figure DEST_PATH_IMAGE034
、
Figure DEST_PATH_IMAGE035
和
Figure DEST_PATH_IMAGE036
分别代表每一个剖分网格的四个节点中的最大泊松比、最小泊松比、平均泊松比和该网格泊松比的相对差异，
Figure 128228DEST_PATH_IMAGE037
代表该节点的脆性指数。
5.根据权利要求4所述的基于矿物含量和分布的岩石脆性分析方法，其特征在于，步骤S50中每一个剖分网格，其平均杨氏模量
Figure 972688DEST_PATH_IMAGE031
和平均泊松比
Figure 924463DEST_PATH_IMAGE035
的计算方法为：
Figure DEST_PATH_IMAGE038
Figure 833251DEST_PATH_IMAGE039
其中，
Figure DEST_PATH_IMAGE040
、
Figure 618804DEST_PATH_IMAGE041
、
Figure DEST_PATH_IMAGE042
和
Figure 950560DEST_PATH_IMAGE043
分别为剖分网格的四个节点
Figure DEST_PATH_IMAGE044
、
Figure 706026DEST_PATH_IMAGE045
、
Figure DEST_PATH_IMAGE046
和
Figure 174048DEST_PATH_IMAGE047
的杨氏模量，
Figure DEST_PATH_IMAGE048
、
Figure 989557DEST_PATH_IMAGE049
、
Figure DEST_PATH_IMAGE050
和
Figure 310073DEST_PATH_IMAGE051
分别为剖分网格的四个节点
Figure 603651DEST_PATH_IMAGE044
、
Figure 50813DEST_PATH_IMAGE045
、
Figure 506065DEST_PATH_IMAGE046
和
Figure 546834DEST_PATH_IMAGE047
的泊松比。
6.一种基于矿物含量和分布的岩石脆性分析系统，其特征在于，该系统包括以下模块：
数字化处理模块，配置为进行待分析岩石图像的数字化处理，获得待分析岩石的二维数组；
自相关分析模块，配置为进行所述二维数组的自相关分析，并将分析获取的自相关相关系数变化曲线的纵坐标值下降为1/e时的横坐标值作为待分析岩石的自相关长度；
自相关函数构建模块，配置为根据所述待分析岩石的自相关长度，待分析岩石中矿物以及弹性模量的分布特征，构建描述矿物分布的自相关函数：
Figure 378523DEST_PATH_IMAGE001
其中，
Figure 211350DEST_PATH_IMAGE002
为待分析岩石的自相关长度；
Figure 712870DEST_PATH_IMAGE003
为数据点
Figure 365568DEST_PATH_IMAGE004
与岩芯切片中心
Figure 949DEST_PATH_IMAGE005
的相对距离，
Figure 157124DEST_PATH_IMAGE006
分别为数据点
Figure 95124DEST_PATH_IMAGE004
的横坐标值和纵坐标值；
密度及弹性模量分布计算模块，配置为根据所述自相关函数进行所述待分析岩石的岩石样品重构，获得待分析岩石的密度及弹性模量分布：
Figure DEST_PATH_IMAGE052
其中，
Figure 235118DEST_PATH_IMAGE008
为数据点
Figure 674190DEST_PATH_IMAGE004
的密度，
Figure 324351DEST_PATH_IMAGE009
为待分析岩石的背景密度；
Figure 823466DEST_PATH_IMAGE010
为数据点
Figure 185177DEST_PATH_IMAGE004
的体积模量，
Figure 37727DEST_PATH_IMAGE011
为体积模量背景值；
Figure 168494DEST_PATH_IMAGE012
为数据点
Figure 572930DEST_PATH_IMAGE004
的剪切模量，
Figure 421938DEST_PATH_IMAGE013
为剪切模量背景值；
Figure 812599DEST_PATH_IMAGE006
分别为数据点
Figure 63451DEST_PATH_IMAGE004
的横坐标值和纵坐标值；
Figure 638789DEST_PATH_IMAGE014
为岩石非均质性产生的密度扰动，
Figure 975093DEST_PATH_IMAGE015
为岩石非均质性产生的体积模量扰动，
Figure 169445DEST_PATH_IMAGE016
为岩石非均质性产生的剪切模量扰动；
Figure 274804DEST_PATH_IMAGE017
代表求平均值操作；
脆性分析模块，配置为根据所述待分析岩石的密度及弹性模量分布进行重构的岩石样品的网格剖分，获得待分析岩石的不同尺度剖分的脆性指数分布。
7.一种电子设备，其特征在于，包括：
至少一个处理器；以及
与至少一个所述处理器通信连接的存储器；其中，
所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现权利要求1-5任一项所述的基于矿物含量和分布的岩石脆性分析方法。
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现权利要求1-5任一项所述的基于矿物含量和分布的岩石脆性分析方法。
CN202110469932.0A 2021-04-29 2021-04-29 基于矿物含量和分布的岩石脆性分析方法、系统及设备 Active CN112986247B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202110469932.0A CN112986247B (zh) 2021-04-29 2021-04-29 基于矿物含量和分布的岩石脆性分析方法、系统及设备
US17/471,195 US11282183B1 (en) 2021-04-29 2021-09-10 Rock brittleness analysis method and system based on mineral content and distribution and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110469932.0A CN112986247B (zh) 2021-04-29 2021-04-29 基于矿物含量和分布的岩石脆性分析方法、系统及设备

Publications (2)

Publication Number Publication Date
CN112986247A CN112986247A (zh) 2021-06-18
CN112986247B true CN112986247B (zh) 2021-08-03

Family

ID=76336523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110469932.0A Active CN112986247B (zh) 2021-04-29 2021-04-29 基于矿物含量和分布的岩石脆性分析方法、系统及设备

Country Status (2)

Country Link
US (1) US11282183B1 (zh)
CN (1) CN112986247B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112986247B (zh) * 2021-04-29 2021-08-03 中国科学院地质与地球物理研究所 基于矿物含量和分布的岩石脆性分析方法、系统及设备
CN115619848B (zh) * 2022-10-13 2023-08-11 水利部交通运输部国家能源局南京水利科学研究院 一种评价骨料中碱硅酸反应活性矿物分布特征的方法
CN115357994B (zh) * 2022-10-20 2023-03-17 中国地质大学（北京） 软岩隧道围岩参数空间随机场建模方法、装置及设备
CN116106510B (zh) * 2023-02-28 2024-05-10 山东科技大学 一种基于超声波测试的岩石冻融后脆性程度评价方法
CN116183875B (zh) * 2023-04-21 2023-07-07 煤炭科学研究总院有限公司 岩体波速的确定方法、装置和电子设备
CN116165054B (zh) * 2023-04-21 2023-07-14 煤炭科学研究总院有限公司 岩体力学参数获取方法、装置和电子设备
CN116468724B (zh) * 2023-06-08 2023-08-29 四川亿欣新材料有限公司 一种基于光学传感技术的碳酸钙含量测试方法
CN117370718B (zh) * 2023-10-18 2024-03-29 中国地质大学（北京） 页岩脆性指数的获取方法、装置、终端及存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502971A (zh) * 2014-12-22 2015-04-08 中国石油天然气集团公司 页岩储层总有机碳含量和脆性空间展布的检测方法
CN107290219A (zh) * 2017-07-20 2017-10-24 长江大学 储气库盖层岩石的综合脆性评价方法
WO2018175404A1 (en) * 2017-03-20 2018-09-27 Saudi Arabian Oil Company Determining rock properties
CN109490965A (zh) * 2018-10-15 2019-03-19 长江大学 一种定量评价地层非均匀性的方法及装置
CN109782348A (zh) * 2019-01-22 2019-05-21 太原理工大学 基于多孔介质模型的储层流体与脆性同时反演的识别方法
WO2019117857A1 (en) * 2017-12-11 2019-06-20 Halliburton Energy Services, Inc. Measuring mechanical properties of rock cuttings
CN110851937A (zh) * 2018-07-25 2020-02-28 中国石油化工股份有限公司 页岩储层脆性指数计算方法及系统

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8542674B2 (en) * 2008-12-02 2013-09-24 International Business Machines Corporation System and method to initiate a presence driven peer to peer communications session on non-IMS and IMS networks
US8612156B2 (en) * 2010-03-05 2013-12-17 Vialogy Llc Active noise injection computations for improved predictability in oil and gas reservoir discovery and characterization
US20130013209A1 (en) * 2010-03-11 2013-01-10 Yaping Zhu Predicting anisotropic source rock properties from well data
WO2013148632A1 (en) * 2012-03-29 2013-10-03 Ingrain, Inc. A method and system for estimating properties of porous media such as fine pore or tight rocks
RU2015104162A (ru) * 2012-08-10 2016-09-27 Ингрейн, Инк. Способ повышения точности значений свойств горной породы, полученных из цифровых изображений
CN105445440B (zh) * 2015-11-19 2017-06-16 西南石油大学 一种基于岩屑微观特征的页岩可压裂性评价方法
US10816440B2 (en) * 2017-02-20 2020-10-27 Conocophillips Company Rock mechanical properties from drill cuttings
CN109828031B (zh) * 2019-02-15 2021-07-23 西南石油大学 岩石脆性评价方法及装置
CN111398273B (zh) * 2019-10-09 2021-07-06 天津大学 一种岩石几何-力学参数获取方法及全息扫描系统
CN112986247B (zh) * 2021-04-29 2021-08-03 中国科学院地质与地球物理研究所 基于矿物含量和分布的岩石脆性分析方法、系统及设备

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502971A (zh) * 2014-12-22 2015-04-08 中国石油天然气集团公司 页岩储层总有机碳含量和脆性空间展布的检测方法
WO2018175404A1 (en) * 2017-03-20 2018-09-27 Saudi Arabian Oil Company Determining rock properties
CN107290219A (zh) * 2017-07-20 2017-10-24 长江大学 储气库盖层岩石的综合脆性评价方法
WO2019117857A1 (en) * 2017-12-11 2019-06-20 Halliburton Energy Services, Inc. Measuring mechanical properties of rock cuttings
CN110851937A (zh) * 2018-07-25 2020-02-28 中国石油化工股份有限公司 页岩储层脆性指数计算方法及系统
CN109490965A (zh) * 2018-10-15 2019-03-19 长江大学 一种定量评价地层非均匀性的方法及装置
CN109782348A (zh) * 2019-01-22 2019-05-21 太原理工大学 基于多孔介质模型的储层流体与脆性同时反演的识别方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Evaluation of methods for determining rock brittleness under compression;Wen Tao等;《Journal of Natural Gas Science and Engineering》;20200418;全文 *
基于岩石物理与矿物组成的页岩脆性评价新方法;秦晓艳等;《天然气地球科学》;20161031;第27卷(第10期);全文 *

Also Published As

Publication number Publication date
US11282183B1 (en) 2022-03-22
CN112986247A (zh) 2021-06-18

Similar Documents

Publication Publication Date Title
CN112986247B (zh) 基于矿物含量和分布的岩石脆性分析方法、系统及设备
Goda et al. New scaling relationships of earthquake source parameters for stochastic tsunami simulation
Salinas et al. Comparative assessment of predictions in ungauged basins–Part 2: Flood and low flow studies
Quintana-Seguí et al. Validation of a new SAFRAN-based gridded precipitation product for Spain and comparisons to Spain02 and ERA-Interim
Kepner et al. An automated cluster finder: the adaptive matched filter
Garnero et al. Comparisons between different interpolation techniques
CN102147479B (zh) 一种储层空间物性参数的建模方法
CN103282908A (zh) 用于表征储藏层评估不确定性的系统和方法
Gyasi‐Agyei Copula‐based daily rainfall disaggregation model
Feng et al. A training image evaluation and selection method based on minimum data event distance for multiple-point geostatistics
Deshmukh et al. Physio‐climatic controls on vulnerability of watersheds to climate and land use change across the US
CN102576370B (zh) 空隙度分析的系统和方法
Yanai et al. Quantifying uncertainty in forest nutrient budgets
CN110244363B (zh) 一种用于预测缝洞型储集层资源量的方法
Fathian et al. Temporal trends in precipitation using spatial techniques in GIS over Urmia Lake Basin, Iran
Yamamoto et al. A comparison of different survey methods for assessing gap parameters in old-growth forests
Goslee Behavior of vegetation sampling methods in the presence of spatial autocorrelation
Lane et al. A large-sample investigation into uncertain climate change impacts on high flows across Great Britain
Wechsler Digital Elevation Model (DEM) uncertainty: evaluation and effect on topographic parameters
Borkowski et al. Geostatistical modelling as an assessment tool of soil pollution based on deposition from atmospheric air
CN116842416A (zh) 一种适用珊瑚礁砂基础管桩侧摩阻力的计算方法
Xuelei et al. Pedodiversity analysis in Hainan Island
Erdede et al. Examining the interpolation methods used in forming the digital elevation models
Eris et al. Streamflow map of the eastern black sea region, turkey
Crema Aoristic approaches and voxel models for spatial analysis

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