一种岩石有限应变测量方法、系统、电子设备及存储介质技术方案

本发明专利技术公开一种岩石有限应变测量方法、系统、电子设备及存储介质，涉及数据处理技术领域，包括：获取岩石图像，标记应变标志体的中心点，创建点集，绘制Fry点云图；定义等面积的同心内椭圆和外椭圆环的密度比，所述密度比为处于外椭圆环内的中心点数与处于内椭圆内的中心点数的比；在Fry点云图中绘制空心椭圆，设定空心椭圆参数的取值范围，在取值范围内遍历使密度比最大的空心椭圆参数组合，从而得到最佳应变椭圆；根据最佳应变椭圆分析得到岩石有限应变测量结果，由此实现对岩石有限应变的测量分析。分析。分析。

【技术实现步骤摘要】
一种岩石有限应变测量方法、系统、电子设备及存储介质


[0001]本专利技术涉及数据处理
，特别是涉及一种岩石有限应变测量方法、系统、电子设备及存储介质。

技术介绍

[0002]岩石有限应变测量是现代构造地质学研究的重要手段，其利用岩石内应变标志体的形态和分布来确定岩石的应变状态，构成岩石变形定量研究的基础，反映岩石变形前和变形后几何形态差异的实体称为应变标志体，常用的如化石、鲕粒等。
[0003]目前有较多应变测量方法被提出，其中一类关注应变标志体单体或集合体变形前后的形态或方向变化，如反向轮法、莫尔圆图解法等，另一类关注应变标志体的相对位置关系，如心对心法等。Fry法将心对心法发展为一种更为简便的测量方法，其中，需要说明的是，Fry法是英国地质学家Fry N于1979年提出的一种岩石有限应变测量方法，系本领域公知常识，中文曾译名为弗莱法，但后续文献多仍以英文表述。该方法通过平移统计应变标志体中心点，形成一个空心椭圆来代表应变椭圆。
[0004]Fry法利用透明纸进行手工制图，操作过程简单，但经多次重复，耗时长、效率低，且该方法操作过程中可能出现误差的地方主要有两处，即标志点的选择和空心椭圆的确定；前者主要由于样品表面不平整或亚颗粒的存在所造成，但基于大量的统计，这种误差的影响是微乎其微的；而空心椭圆作为Fry法最后确定应变的直接依据，通常由人为划定，主观性较强、误差大。
[0005]为了解决这一问题，有研究提出基于计算程序编程语言的自动实现方法，但是对于一些设计年代较早的语言，对可视化支持较差，即使能解决自动化的问题，最终仍是由人工确定应变椭圆，仍存在主观性强的问题；而且自动识别方法存在一个关键挑战
‑‑
伪影，即一个与真实空心椭圆存在交角的、范围较小的应变标志点空缺区，伪影的存在会对结果造成误差。

技术实现思路

[0006]为了解决上述问题，本专利技术提出了一种岩石有限应变测量方法、系统、电子设备及存储介质，以密度比为优化目标，在取值范围内遍历使密度比最大的椭圆参数，从而得到最佳应变椭圆，以此实现对岩石有限应变的测量分析。
[0007]为了实现上述目的，本专利技术采用如下技术方案：第一方面，本专利技术提供一种岩石有限应变测量方法，包括：获取岩石图像，并标记应变标志体的中心点，以此创建点集，根据点集中各中心点的坐标，绘制Fry点云图；定义等面积的同心内椭圆和外椭圆环的密度比，所述密度比为处于外椭圆环内的中心点数与处于内椭圆内的中心点数的比；在Fry点云图中绘制空心椭圆，设定空心椭圆参数的取值范围，在取值范围内遍历
使密度比最大的空心椭圆参数组合，从而得到最佳应变椭圆；根据最佳应变椭圆分析得到岩石有限应变测量结果。
[0008]作为可选择的实施方式，创建点集的过程包括：对应变标志体的中心点形成原始点集，将原始点集中所有中心点移动设定距离，记录中心点移动后的坐标，从而得到新的点集。
[0009]作为可选择的实施方式，设定距离为，每次移动后，第n个中心点移动至点（0,0），并记录其余中心点移动后的坐标；分别是第n个中心点的横坐标和纵坐标。
[0010]作为可选择的实施方式，点集为，其中，为点集中第m个中心点的坐标。
[0011]作为可选择的实施方式，根据Fry点云图设定空心椭圆参数的取值范围，使密度比最大的空心椭圆参数组合为最佳应变椭圆的参数，由此得到最佳应变椭圆。
[0012]作为可选择的实施方式，所述空心椭圆参数包括半长轴、半短轴和旋转角。
[0013]作为可选择的实施方式，根据最佳应变椭圆的参数确定包括应变方向和应变值的岩石有限应变测量结果；其中，旋转角代表应变方向，半短轴与半长轴的比值代表应变值。
[0014]第二方面，本专利技术提供一种岩石有限应变测量系统，包括：图绘制模块，被配置为获取岩石图像，并标记应变标志体的中心点，以此创建点集，根据点集中各中心点的坐标，绘制Fry点云图；密度比定义模块，被配置为定义等面积的同心内椭圆和外椭圆环的密度比，所述密度比为处于外椭圆环内的中心点数与处于内椭圆内的中心点数的比；拟合模块，被配置为在Fry点云图中绘制空心椭圆，设定空心椭圆参数的取值范围，在取值范围内遍历使密度比最大的空心椭圆参数组合，从而得到最佳应变椭圆；分析模块，被配置为根据最佳应变椭圆分析得到岩石有限应变测量结果。
[0015]第三方面，本专利技术提供一种电子设备，包括存储器和处理器以及存储在存储器上并在处理器上运行的计算机指令，所述计算机指令被处理器运行时，完成第一方面所述的方法。
[0016]第四方面，本专利技术提供一种计算机可读存储介质，用于存储计算机指令，所述计算机指令被处理器执行时，完成第一方面所述的方法。
[0017]与现有技术相比，本专利技术的有益效果为：本专利技术提出一种岩石有限应变测量方法，标记应变标志体的中心点后创建点要素类，以此绘制Fry点云图，然后设定应变椭圆半长轴、半短轴和旋转角的取值范围，在取值范围内遍历使密度比最大的椭圆参数，从而得到最佳应变椭圆，以此实现对岩石有限应变的测量分析，以密度比为优化目标，自动校验椭圆参数并得到最佳参数，解决现有方法中主观性强的问题。
[0018]本专利技术附加方面的优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本专利技术的实践了解到。
附图说明
[0019]构成本专利技术的一部分的说明书附图用来提供对本专利技术的进一步理解，本专利技术的示意性实施例及其说明用于解释本专利技术，并不构成对本专利技术的不当限定。
[0020]图1（a）为Fry法原理中构成的一系列圆的示意图；图1（b）为Fry法原理图中岩石均匀变形后圆变为椭圆的示意图；图2为本专利技术实施例1提供的岩石有限应变测量方法流程图；图3为点密度法原理示意图；图4为本专利技术实施例1提供的砂岩图像示意图；图5为本专利技术实施例1提供的基于砂岩图像得到的最佳空心椭圆示意图；图6为本专利技术实施例1提供的砾岩图像示意图；图7为本专利技术实施例1提供的基于砾岩图像得到的最佳空心椭圆示意图。
具体实施方式
[0021]下面结合附图与实施例对本专利技术做进一步说明。
[0022]应该指出，以下详细说明都是示例性的，旨在对本专利技术提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本专利技术所属
的普通技术人员通常理解的相同含义。
[0023]需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本专利技术的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
[0024]在不冲突的情况下，本专利技术中的实施例及实施例中的特征可以相互组合。
[0025]实本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种岩石有限应变测量方法，其特征在于，包括：获取岩石图像，并标记应变标志体的中心点，以此创建点集，根据点集中各中心点的坐标，绘制Fry点云图；定义等面积的同心内椭圆和外椭圆环的密度比，所述密度比为处于外椭圆环内的中心点数与处于内椭圆内的中心点数的比；在Fry点云图中绘制空心椭圆，设定空心椭圆参数的取值范围，在取值范围内遍历使密度比最大的空心椭圆参数组合，从而得到最佳应变椭圆；根据最佳应变椭圆分析得到岩石有限应变测量结果。2.如权利要求1所述的一种岩石有限应变测量方法，其特征在于，创建点集的过程包括：对应变标志体的中心点形成原始点集，将原始点集中所有中心点移动设定距离，记录中心点移动后的坐标，从而得到新的点集。3.如权利要求2所述的一种岩石有限应变测量方法，其特征在于，设定距离为，每次移动后，第n个中心点移动至点（0,0），并记录其余中心点移动后的坐标，分别是第n个中心点的横坐标和纵坐标。4.如权利要求3所述的一种岩石有限应变测量方法，其特征在于，点集为，其中，为点集中第m个中心点的坐标。5.如权利要求1所述的一种岩石有限应变测量方法，其特征在于，根据Fry点云图设定空心椭圆参数的取值范围，使密度比最大的空心椭圆参数组合为最佳应变椭圆的参数，由此得到最佳应变椭圆。6.如权利要求1所述的一种岩石有限应变...

【专利技术属性】
技术研发人员：常杰，倪永进，杨斌，宋香锁，梁吉坡，仵康林，
申请(专利权)人：山东省地质科学研究院，
类型：发明
国别省市：