岩浆岩体侵入过程的反演方法及装置制造方法及图纸

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

【技术实现步骤摘要】
岩浆岩体侵入过程的反演方法及装置
本专利技术涉及地理信息和地质学领域，尤其涉及一种岩浆岩体侵入过程的反演方法及装置。
技术介绍
侵入岩是指当上覆岩层压力减轻时，岩浆在侵入过程中变冷、结晶而形成的岩石。由于岩浆侵入过程中的缓慢冷却，有足够的时间使矿物结晶，因此侵入岩晶粒粗大，具有显晶质结构。侵入岩侵入于早先形成的岩石中时，“最省力”的方式是沿裂隙侵入并使其横截面有较小的周长，主体沿侵入方向延伸，虽形态多样，但多为近圆柱状。目前，岩浆岩体的侵入过程反演主要是地质工作者依靠专家知识，通过人工绘制侵入过程示意图或手工制作动画，来进行侵入过程的动态表达。相关反演方法，效率低下，投入大，且反演过程的正确性参差不齐。
技术实现思路
专利技术目的：本专利技术针对现有技术存在的问题，提供一种正确率和效率更高的岩浆岩体侵入过程的反演方法及装置。技术方案：本专利技术所述的岩浆岩体侵入过程的反演方法包括：(1)加载岩浆岩体侵入剖面矢量数据与岩浆岩体发育时序表，构建岩浆岩体的邻接关系矩阵与发育时序集合SK；(2)从集合本文档来自技高网...

【技术保护点】
1.一种岩浆岩体侵入过程的反演方法，其特征在于该方法包括：/n(1)加载岩浆岩体侵入剖面矢量数据与岩浆岩体发育时序表，构建岩浆岩体的邻接关系矩阵与发育时序集合SK；/n(2)从集合SK中获取任一岩浆岩体s

【技术特征摘要】
1.一种岩浆岩体侵入过程的反演方法，其特征在于该方法包括：
(1)加载岩浆岩体侵入剖面矢量数据与岩浆岩体发育时序表，构建岩浆岩体的邻接关系矩阵与发育时序集合SK；
(2)从集合SK中获取任一岩浆岩体sa；
(3)根据邻接关系矩阵获取sa的任一邻接岩浆岩体sb，并分别获取sa、sb中所有岩浆岩体分块，形成分块集合PA、PB；
(4)从分块集合PA中提取出邻接个数大于等于2的分块，对于每一分块，根据与分块集合PB中相邻的分块的左侧和右侧相接线，构建岩浆岩体相接面，作为当前分块的邻接岩浆岩体；
(5)从分块集合PA中提取出邻接个数为1的分块，对于每一分块，根据与分块集合PB中相邻的分块的相接线，构建岩浆岩体相接面，作为当前分块的邻接岩浆岩体；
(6)循环执行步骤(3)-(5)，直至sa的所有邻接岩浆岩体被遍历，将所有分块的邻接岩浆岩体融合，作为sa的邻接岩浆岩体数据集GeoRocka；
(7)循环执行步骤(2)-(6)，直至集合SK中所有元素被遍历，得到所有岩浆岩体的邻接岩浆岩体数据集GeoRock；
(8)基于数据集GeoRock，制作岩浆岩体侵入过程的中间帧，生成侵入过程反演动画。


2.根据权利要求1所述的岩浆岩体侵入过程的反演方法，其特征在于：步骤(1)具体包括：
(1-1)加载岩浆岩体侵入剖面矢量数据，得到岩浆岩体集合S＝{si|i＝1,2,…,sn}，si为第i个岩浆岩体，i为岩浆岩体序号，sn代表岩浆岩体数量；
(1-2)创建大小为sn*sn的岩浆岩体邻接矩阵；
(1-3)根据岩浆岩体的邻接关系，对岩浆岩体邻接矩阵进行赋值，其中，当两个岩浆岩体邻接时，对应岩浆岩体邻接矩阵元素赋值为1；否则，赋值为0；
(1-4)加载岩浆岩体发育时序表，从中获取到岩浆岩体的发育时序关系，并按照发育时序从新到老的顺序，对岩浆岩体集合S重新排序，得到岩浆岩体发育时序集合SK＝{sx1,sx2,sx3,。。。}，其中，x1,x2,x3分别为第1,2,3个发育的岩浆岩体的序号，取值[1,sn]。


3.根据权利要求1所述的岩浆岩体侵入过程的反演方法，其特征在于：步骤(3)具体包括：
(3-1)根据邻接关系矩阵获取sa的所有邻接岩浆岩体，形成邻接岩浆岩体集合NK；
(3-2)从集合NK中读取任一岩浆岩体sb；
(3-3)分别获取岩浆岩体sa、sb中所有的分块，形成对应分块集合PA、PB。


4.根据权利要求1所述的岩浆岩体侵入过程的反演方法，其特征在于：步骤(4)具体包括：
(4-1)从分块集合PA中提取出邻接个数大于等于2的分块，存入子集SA＝{aα|α＝1,2,…,san}，其中，aα为提取的第α个分块，san为提取的分块的数量；
(4-2)从子集SA中获取任一元素aα，从PB中获取与aα邻接的所有分块，形成子集SB＝{bβ|β＝1,2,…,sbn}，其中，bβ为PB中与aα邻接的第β个分块，sbn为与aα邻接的分块的数量；
(4-3)从子集SB中获取任一元素bβ，并获取aα与bβ的相接线Touchline；
(4-4)判断bβ位于aα左侧还是右侧；
(4-5)若bβ位于aα右侧，则将Touchline保存到右侧相接线集合RL中；若bβ位于aα左侧，将Touchline保存到左侧相接线集合LL中；
(4-6)循环执行步骤(4-3)-(4-5)，直至子集SB中所有元素被遍历；
(4-7)从集合RL中获取任一右侧相接线rlk，并从集合LL中寻找到对应的左侧相接线llk；
(4-8)根据相接线rlk和llk构建岩浆岩体相接面，并存入相接面集合SEα；
(4-9)循环执行步骤(4-7)-(4-8)，直至集合RL中所有元素被遍历，得到分块aα的相接面集合SEα；
(4-10)循环执行步骤(4-2)-(4-9)，直至子集SA中所有元素被遍历，得到所有分块的相接面集合SE＝{SEα}。


5.根据权利要求1所述的岩浆岩体侵入过程的反演方法，其特征在于：步骤(4-4)具体包括：
(4-4-1)从子集SB中获取任一元素bβ，并获取bβ的外接矩形的中心点po(xo,yo)以及aα的外接矩形的两条短边线段的端点pa、pb和pc、pd；
(4-4-2)求取pa、pb的中点pe坐标(xe,ye)，以及pc、pd的中点pf坐标(xf,yf)；
(4-4-3)根据下式，计算得到向量与






(4-4-4)采用下式计算bβ与aα的位置关系指标T；



(4-4-5)若T≥0，则判定bβ位于aα右侧，否则判定bβ位于aα左侧。


6.根据权利要求5所述的岩浆岩体侵入过程的反演方法，其特征在于：步骤(4-7)具体包括：
(4-7-1)从集合RL中获取任一元素rlu，获取rlu的起点pst与终点ped；
(4-7-2)求取点pst与点ped的中点pmid，以及点pmid到点pe的直线距离du；
(4-7-3)循环执行步骤(4-7-1)-(4-7-2)，直至得到所有右侧相接线rlu对应的距离du；
(4-7-4)根据距离du的大小，按照从小到大排序对应的相接...

【专利技术属性】
技术研发人员：沈言根，李安波，黄键初，闾国年，
申请(专利权)人：南京师范大学，
类型：发明
国别省市：江苏;32