本申请实施例提供一种岩石力学参数获取方法及其装置、计算机设备、存储介质，所述岩石力学参数获取方法包括：获取微岩心的三维图像和矿物相成分；根据所述三维图像和所述矿物相成分构建三维孔隙矿物结构模型；对所述三维孔隙矿物结构模型进行切片，得到二维孔隙矿物结构模型；根据所述二维孔隙矿物结构模型构建数字岩心几何模型；根据有限单元法和离散单元法对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的岩石力学参数。本申请实施例采用微岩心样本，通过有限单元法和离散单元法相结合对微岩心样本建立的数字岩心几何模型进行数值模拟，能够准确快速得到岩石力学参数。快速得到岩石力学参数。快速得到岩石力学参数。
【技术实现步骤摘要】
岩石力学参数获取方法及其装置、计算机设备、存储介质
[0001]本专利技术涉及岩石物理学及石油测井
，尤其涉及一种岩石力学参数获取方法及其装置、计算机设备、存储介质。
技术介绍
[0002]相关技术中，采用石油钻井技术获取储层岩样，通过岩样的单轴压缩、三轴压缩试验进行岩样受压破裂过程分析，得到静态岩石力学参数。但是通过石油钻井技术获取的岩样具有成本高、数量少、质量差、周期长等限制，且一些钻井难以获得完整的岩心。
技术实现思路
[0003]本申请实施例的主要目的在于提出一种岩石力学参数获取方法及其装置、计算机设备、存储介质，能够避免石油钻井技术获取岩样的限制，快速获取岩石力学参数。
[0004]为实现上述目的，本申请实施例的第一方面提出了一种岩石力学参数获取方法，所述方法包括：
[0005]获取微岩心的三维图像和矿物相成分；
[0006]根据所述三维图像和所述矿物相成分构建三维孔隙矿物结构模型；
[0007]对所述三维孔隙矿物结构模型进行切片，得到二维孔隙矿物结构模型；
[0008]根据所述二维孔隙矿物结构模型构建数字岩心几何模型；
[0009]根据有限单元法和离散单元法对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的岩石力学参数。
[0010]在一些实施例，所述根据有限单元法和离散单元法对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的岩石力学参数，包括：
[0011]根据有限单元法对所述数字岩心几何模型进行网格划分，得到单元集合；
[0012]对所述单元集合中的网格单元设置网格单元矿物力学参数；
[0013]构建边界约束条件；
[0014]将所述网格单元矿物力学参数作为所述离散单元法的输入参数，根据所述输入参数和所述边界约束条件对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的岩石力学参数。
[0015]在一些实施例，所述对所述单元集合中的网格单元设置网格单元矿物力学参数，包括：
[0016]获取所述数字岩心几何模型中像素点的坐标和像素值；
[0017]根据所述坐标与所述单元集合的网格单元进行匹配，将像素值作为对应的所述网格单元的标记；
[0018]根据所述标记对所述单元集合进行分组，得到分组网格单元；
[0019]对所述分组网格单元设置所述网格单元矿物力学参数。
[0020]在一些实施例，所述网格单元矿物力学参数包括单元边界接触参数，所述将所述
网格单元矿物力学参数作为所述离散单元法的输入参数，根据所述输入参数和所述边界约束条件对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的岩石力学参数，包括：
[0021]将所述单元边界接触参数作为所述离散单元法的输入参数，根据所述输入参数和所述边界约束条件对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的应力数据和应变数据；
[0022]根据所述应力数据和所述应变数据得到所述数字岩心几何模型的岩石力学参数。
[0023]在一些实施例，在所述根据有限单元法和离散单元法对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的岩石力学参数之后，所述岩石力学参数获取方法还包括：
[0024]获取样本岩石力学参数；
[0025]根据所述样本岩石力学参数对所述数字岩心几何模型的岩石力学参数进行参数标定。
[0026]在一些实施例，所述根据所述二维孔隙矿物结构模型构建数字岩心几何模型，包括：
[0027]获取所述二维孔隙矿物结构模型的像素点密度和像素点尺寸；
[0028]根据所述像素点密度和所述像素点尺寸构建所述数字岩心几何模型。
[0029]在一些实施例，所述获取微岩心的三维图像和矿物相成分，包括：
[0030]获取计算机断层扫描装置采集的所述微岩心的三维图像；
[0031]获取X射线衍射装置采集的所述微岩心的矿物相成分。
[0032]本申请实施例的第二方面提出了一种岩石力学参数获取装置，所述装置包括：
[0033]获取模块，用于获取微岩心的三维图像和矿物相成分；
[0034]第一模型构建模块，用于根据所述三维图像和所述矿物相成分构建三维孔隙矿物结构模型；
[0035]切片模块，用于对所述三维孔隙矿物结构模型进行切片，得到二维孔隙矿物结构模型；
[0036]第二模型构建模块，用于根据所述二维孔隙矿物结构模型构建数字岩心几何模型；
[0037]数值模拟模块，用于根据有限单元法和离散单元法对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的岩石力学参数。
[0038]本申请实施例的第三方面提出了一种计算机设备，所述计算机设备包括存储器和处理器，其中，所述存储器中存储有程序，所述程序被所述处理器执行时所述处理器用于执行如本申请第一方面实施例任一项所述的方法。
[0039]本申请实施例的第四方面提出了一种存储介质，该存储介质为计算机可读存储介质，所述存储介质存储有计算机可执行指令，所述计算机可执行指令用于使计算机执行如本申请第一方面实施例任一项所述的方法。
[0040]本申请实施例提出的岩石力学参数获取方法及其装置、计算机设备、存储介质，通过获取微岩心的三维图像和矿物相成分，根据三维图像和矿物相成分构建三维孔隙矿物结构模型，对三维孔隙矿物结构模型进行切片，得到二维孔隙矿物结构模型，根据二维孔隙矿
物结构模型构建数字岩心几何模型，根据有限单元法和离散单元法对数字岩心几何模型的受压破裂过程进行数值模拟，得到数字岩心几何模型的岩石力学参数，本申请实施例通过有限单元法和离散单元法相结合的数值模拟方法模拟数字岩心几何模型的受压破裂过程，能够快速得到岩石力学参数，避免岩样获取、岩心不完整等方面的限制。
附图说明
[0041]图1是本申请实施例提供的岩石力学参数获取方法的第一流程图；
[0042]图2是图1中步骤S150的具体方法的流程图；
[0043]图3是图2中步骤S220的具体方法的流程图；
[0044]图4是图2中步骤S240的具体方法的流程图；
[0045]图5是本申请实施例提供的岩石力学参数获取方法的第二流程图；
[0046]图6是图1中步骤S140的具体方法的流程图；
[0047]图7是本申请实施例提供的岩石力学参数获取装置的模块结构图。
具体实施方式
[0048]为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本申请，并不用于限定本申请。
[0049]需要说明的是，虽然在装置示意图中进行了功能模块划分，在流程图中示出了逻辑顺序，但是在某些情况下，可以本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.岩石力学参数获取方法，其特征在于，所述方法包括：获取微岩心的三维图像和矿物相成分；根据所述三维图像和所述矿物相成分构建三维孔隙矿物结构模型；对所述三维孔隙矿物结构模型进行切片，得到二维孔隙矿物结构模型；根据所述二维孔隙矿物结构模型构建数字岩心几何模型；根据有限单元法和离散单元法对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的岩石力学参数。2.根据权利要求1所述的岩石力学参数获取方法，其特征在于，所述根据有限单元法和离散单元法对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的岩石力学参数，包括：根据有限单元法对所述数字岩心几何模型进行网格划分，得到单元集合；对所述单元集合中的网格单元设置网格单元矿物力学参数；构建边界约束条件；将所述网格单元矿物力学参数作为所述离散单元法的输入参数，根据所述输入参数和所述边界约束条件对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的岩石力学参数。3.根据权利要求2所述的岩石力学参数获取方法，其特征在于，所述对所述单元集合中的网格单元设置网格单元矿物力学参数，包括：获取所述数字岩心几何模型中像素点的坐标和像素值；根据所述坐标与所述单元集合的网格单元进行匹配，将像素值作为对应的所述网格单元的标记；根据所述标记对所述单元集合进行分组，得到分组网格单元；对所述分组网格单元设置所述网格单元矿物力学参数。4.根据权利要求2所述的岩石力学参数获取方法，其特征在于，所述网格单元矿物力学参数包括单元边界接触参数，所述将所述网格单元矿物力学参数作为所述离散单元法的输入参数，根据所述输入参数和所述边界约束条件对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的岩石力学参数，包括：将所述单元边界接触参数作为所述离散单元法的输入参数，根据所述输入参数和所述边界约束条件对所述数字岩心几何模型的受压破裂过程进行数值模拟，得到所述数字岩心几何模型的应力数据和应变...
【专利技术属性】
技术研发人员：龙威，王冠群，李炜，蔡坤鹏，
申请(专利权)人：深圳清华大学研究院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。