[发明专利]一种断层发育时序的获取方法

说明书

本发明公开了一种断层发育时序的获取方法，具体包括步骤：首先，基于断层分段间的邻接关系，计算断层分段切割关系矩阵；其次，基于断层分段切割关系，计算断层切割关系矩阵；然后，基于断层切割关系矩阵，计算断层发育时序；最后，计算断层分段发育时序编码。本发明实现了基于切割律快速获取断层发育时序，与现有技术相比，可以有效提高判别效率和判别质量，对于复杂断层发育区的断层发育过程表达与时空模拟具有重要的研究与应用价值。

技术领域

本发明涉及地理信息技术和地质学领域，尤其涉及一种断层发育时序的获取方法。

背景技术

断层是岩层或岩体顺破裂面发生明显位移的构造，在地壳中广泛发育。地壳断块沿断层的突然运动是地震发生的主要原因，且在断层带上由于岩石破碎和易于风化侵蚀，而沿断层线常常发育沟谷、泉或湖泊。区域断层发育时序的准确解析，是断层构造时空演化特征解析与时空表达的前提。区域断层发育时序代表了研究区域内断层发育的先后顺序，对于区域构造恢复与时空模拟，具有重要研究意义。此外，对于重大工程规划、地质灾害防治、地下空间开发等诸多应用也具有重要的促进作用。

发明内容

发明目的：本发明针对现有技术存在的问题，提供一种断层发育时序的获取方法，可以有效获取断层的发育时序，用于区域构造恢复与时空模拟和其他应用中，有效提高判别效率和判别质量。

技术方案：本发明所述的断层发育时序的获取方法包括：

(1)读取断层线矢量数据到断层分段集合F，并获取断层分段的走向集合D和所属断层编号集合FID；

(2)根据所述断层分段集合F基于断层分段间的邻接关系，生成断层分段切割关系矩阵A；

(3)基于所述断层分段切割关系矩阵A和所述所属断层编号集合FID，生成表征断层切割关系的断层切割关系矩阵B；

(4)基于断层切割关系矩阵B，获取断层的发育时序，形成断层发育时序序列L；

(5)基于断层发育时序序列L，生成断层分段的发育时序编码。

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

(1-1)读取断层线矢量数据到断层分段集合F＝{f_i|i＝1,2,…,fn}；其中，f_i表示第i个断层分段，fn表示断层分段数量；

(1-2)获取集合F中所有断层分段的走向属性，存入走向集合D＝{d_i|i＝1,2,…,fn}；其中，d_i为断层分段f_i的走向；

(1-3)获取集合F中所有断层分段所属断层编号，存入所属断层编号集合FID＝{fid_i|i＝1,2,…,fn}；其中，fid_i为断层分段f_i所属断层编号。

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

(2-1)基于断层分段数量，创建一个断层分段切割关系矩阵A＝{a_i,j|i＝1,2,…,fn,j＝1,2,…,fn}，其中，a_i,j表示断层分段f_i、f_j的切割关系，其初值设为0，fn表示断层分段数量；

(2-2)从断层分段集合F中读取任意两条断层分段f_i、f_j，若f_i、f_j相交，则计算f_i、f_j交点p_ij，并执行步骤(2-3)；否则执行步骤(2-4)；