1.一种侵入岩体发育时序获取方法，其特征在于该方法包括：
(1)根据侵入岩体的地质剖面矢量面图层，构建岩体集合和岩体邻接矩阵；
(2)根据岩体邻接矩阵从岩体集合中读取任意两个邻接岩体sa、sb，并分别获取这两个岩体的分块集合pa和pb；
(3)基于两侧分布规则或半包围规则进行邻接岩体sa、sb的切割关系判别；
(4)循环执行步骤(2)-(3)，直至完成所有邻接岩体切割关系判别；
(5)根据岩体的切割关系，生成岩体切割关系矩阵；
(6)基于岩体切割关系矩阵，运用归并排序方法，生成岩体的发育时序。
2.根据权利要求1所述的侵入岩体发育时序获取方法，其特征在于：步骤(1)包括：
(1-1)加载侵入岩体的地质剖面矢量面图层，得到所有岩体集合s＝{sn|n＝1,2,…,sn}，sn表示第n个岩体，sn表示岩体数量；
(1-2)创建大小为sn*sn的岩体邻接矩阵；
(1-3)根据地质剖面矢量面图层，判定不同岩体间的邻接关系，当两个岩体邻接时，将对应的岩体邻接矩阵元素值赋值为1；否则，赋值为0。
3.根据权利要求1所述的侵入岩体发育时序获取方法，其特征在于：步骤(2)包括：
(2-1)按照从左向右、从上向下顺序，从岩体邻接矩阵中读取两个尚未读取过的岩体sa、sb；其中a、b表示岩体序号，且a、b∈{1,2,…,sn}；
(2-2)若sa、sb对应的岩体邻接矩阵元素值为1，则执行步骤(2-3)；否则返回执行步骤(2-1)；
(2-3)分别读取岩体sa、sb中所有的分块，形成对应的分块集合pa＝{aα|α＝1,2,…,an}和pb＝{bβ|β＝1,2,…,bn}，其中，aα表示岩体sa第α个分块，bβ表示岩体sb第β个分块，an、bn分别为sa、sb中所有的分块数量。
4.根据权利要求1所述的侵入岩体发育时序获取方法，其特征在于：步骤(3)包括：
(3-1)从分块集合pa中提取出邻接个数大于等于2的分块，存入子集sa＝{ai|i＝1,2,…,ai}，其中，ai表示子集sa中第i个岩体分块，ai为分块数量；
(3-2)从子集sa中读取任意一个元素ai；
(3-3)从分块集合pb中获取与ai邻接的分块，存入子集sb＝{bj|j＝1,2,…,bj}，其中，bj表示子集sb中第j个岩体分块，bj为邻接ai的分块个数；
(3-4)判断子集sb中任意两个分块bj、bj+1与ai是否满足两侧分布规则，若是，则判定sa切割sb，并执行步骤(3-7)，否则执行步骤(3-5)；
(3-5)归并分块集合pb，构建sb的整块岩体rb；
(3-6)判断ai与rb是否满足半包围规则，若是，则判定sa切割sb，并执行步骤(3-7)，否则，判定sa、sb切割关系未知，并执行步骤(3-7)；
(3-7)返回执行步骤(3-2)，直至子集sa中所有元素都被遍历。
5.根据权利要求4所述的侵入岩体发育时序获取方法，其特征在于：步骤(3-4)中两侧分布规则的判断方法的包括：
(3-4-1)获取bj、bj+1的外接矩形的中心点
(3-4-2)根据下式，基于端点pta，pwa计算长边的斜率k：
(3-4-3)根据下式，获取过中心点
(3-4-4)根据下式，计算中心点
(3-4-5)若r≤0，则表示bj、bj+1位于ai两侧，满足两侧分布规则，判定sa切割sb；如果r＞0，表示bj、bj+1位于ai同侧，不满足两侧分布规则。
6.根据权利要求4所述的侵入岩体发育时序获取方法，其特征在于：步骤(3-6)中半包围规则的判断方法的包括：
(3-6-1)获取ai的外接矩形
(3-6-2)从外接矩形
(3-6-3)从外接矩形frb的四个角点中，获取外接矩形frb的最大纵坐标值yrmax、最小纵坐标值yrmin、最大横坐标值xrmax和最小横坐标值xrmin；
(3-6-4)根据下式计算用于判断外接矩形
t＝(yamax-yrmax)(yamin-yrmin)(xamax-xrmax)(xamin-xrmin)
(3-6-5)若t>0，表示外接矩形
(3-6-6)根据下式计算用于判断中心点
(3-6-7)若u>0，则表示中心点
7.根据权利要求1所述的侵入岩体发育时序获取方法，其特征在于：步骤(6)具体包括：
(6-1)基于岩体切割关系矩阵，采用归并排序方法，表示按照从老到新的顺序对岩体进行排序；其中，新老的判断准则为：对任意两个岩体，切割的岩体为新，被切割的岩体为老；
(6-2)按照岩体的排序序号，生成对应岩体的发育时序；
(6-3)将生成的发育时序添加到对应岩体的时序属性timeid内。