[发明专利]岩浆岩体侵入过程的反演方法及装置

说明书

本发明公开了一种岩浆岩体侵入过程的反演方法及装置，其中方法主要包括步骤：(1)加载岩浆岩体侵入剖面矢量数据与岩浆岩体发育时序表，构建岩浆岩体邻接矩阵与岩浆岩体发育时序集合；(2)创建邻接岩浆岩体要素的分块集合；(4)基于分块岩浆岩体邻接数量恢复岩浆岩体；(5)构建岩浆岩体数据集；(6)基于岩浆岩体侵入过程中间帧制作岩浆岩体侵入过程反演动画。本发明正确率更高，效率更高，自动化程度更高，可以直观反演出岩浆岩体的侵入过程。

技术领域

本发明涉及地理信息和地质学领域，尤其涉及一种岩浆岩体侵入过程的反演方法及装置。

背景技术

侵入岩是指当上覆岩层压力减轻时，岩浆在侵入过程中变冷、结晶而形成的岩石。由于岩浆侵入过程中的缓慢冷却，有足够的时间使矿物结晶，因此侵入岩晶粒粗大，具有显晶质结构。侵入岩侵入于早先形成的岩石中时，“最省力”的方式是沿裂隙侵入并使其横截面有较小的周长，主体沿侵入方向延伸，虽形态多样，但多为近圆柱状。

目前，岩浆岩体的侵入过程反演主要是地质工作者依靠专家知识，通过人工绘制侵入过程示意图或手工制作动画，来进行侵入过程的动态表达。相关反演方法，效率低下，投入大，且反演过程的正确性参差不齐。

发明内容

发明目的：本发明针对现有技术存在的问题，提供一种正确率和效率更高的岩浆岩体侵入过程的反演方法及装置。

技术方案：本发明所述的岩浆岩体侵入过程的反演方法包括：

(1)加载岩浆岩体侵入剖面矢量数据与岩浆岩体发育时序表，构建岩浆岩体的邻接关系矩阵与发育时序集合SK；

(2)从集合SK中获取任一岩浆岩体s_a；

(3)根据邻接关系矩阵获取s_a的任一邻接岩浆岩体s_b，并分别获取s_a、s_b中所有岩浆岩体分块，形成分块集合PA、PB；

(4)从分块集合PA中提取出邻接个数大于等于2的分块，对于每一分块，根据与分块集合PB中相邻的分块的左侧和右侧相接线，构建岩浆岩体相接面，作为当前分块的邻接岩浆岩体；

(5)从分块集合PA中提取出邻接个数为1的分块，对于每一分块，根据与分块集合PB中相邻的分块的相接线，构建岩浆岩体相接面，作为当前分块的邻接岩浆岩体；

(6)循环执行步骤(3)-(5)，直至s_a的所有邻接岩浆岩体被遍历，将所有分块的邻接岩浆岩体融合，作为s_a的邻接岩浆岩体数据集GeoRock_a；

(7)循环执行步骤(2)-(6)，直至集合SK中所有元素被遍历，得到所有岩浆岩体的邻接岩浆岩体数据集GeoRock；

(8)基于数据集GeoRock，制作岩浆岩体侵入过程的中间帧，生成侵入过程反演动画。

进一步的，步骤(1)具体包括：

(1-1)加载岩浆岩体侵入剖面矢量数据，得到岩浆岩体集合S＝{s_i|i＝1,2,…,sn}，s_i为第i个岩浆岩体，i为岩浆岩体序号，sn代表岩浆岩体数量；

(1-2)创建大小为sn*sn的岩浆岩体邻接矩阵；

(1-3)根据岩浆岩体的邻接关系，对岩浆岩体邻接矩阵进行赋值，其中，当两个岩浆岩体邻接时，对应岩浆岩体邻接矩阵元素赋值为1；否则，赋值为0；