【技术实现步骤摘要】
本专利技术涉及地质解译，具体涉及一种融合地球物理资料与遥感资料或地质踏勘资料的地质构造分析方法及存储介质。
技术介绍
1、查明地质构造对于区域地质调查、矿产勘探和油气勘探都具有十分重要的意义。目前，区域构造信息的提取主要依赖于地质踏勘、遥感解译和地球物理解释。
2、地质踏勘的方法定位准确，是当下地质工作的主要手段，但野外工作时间长、成本高，且受制于自然条件，在植被茂密的地方，工作难以展开。遥感解译可以提供比较准确的线环构造信息，但是无法提取出深部构造信息以及至为重要的控岩构造信息。地球物理，例如航空重力、航空磁测是间接的方法，依赖于地下地质体的物性差异，通过感知物性差异来推断构造，传统上，这种方法依靠线性异常带、串珠状异常、异常错动带等来定位构造，但是由于是间接的方法，因此会产生大量虚假推断。
技术实现思路
1、本专利技术的目的在于提供一种融合多数据源的地质构造分析方法，该分析方法利用重磁资料获取强检测点数量，利用线性构造单元构建缓冲区，再基于缓冲区内强检测点的数量来构建起区域控岩构造骨架，既可以更加精准的定位构造位置，克服地球物理方法多解性的难题，也可以充分利用地球物理方法带来的附属信息，快速建立起构造骨架，划分出重点勘察区域。
2、本专利技术通过下述技术方案实现：
3、融合多数据源的地质构造分析方法，包括以下步骤：
4、s1：获取目标区域的重磁资料和遥感资料或地质踏勘资料；
5、s2：基于所述重磁资料获取至少一个尺度下
6、s3：从所述遥感资料或地质踏勘资料中获取线性构造单元，基于所述线性构造单元建立对应的封闭曲线，所述封闭曲线内的区域作为线性构造单元的缓冲区；
7、s4：累计所述封闭曲线内的各尺度下的强检测点数量，基于所述强检测点数量构建起区域控岩构造骨架。
8、本技术方案中，重磁资料为地球物理原始资料中通过重力法获得的重力资料、以及通过磁法获得的磁测资料。重磁资料通常为二维平面数据，可以表示为tz(x,y,z)，也即观察面位场信息，其中x和y为二维平面上的横坐标和纵坐标，z为重磁测量值的高度。
9、本技术方案中，重磁资料包含的观察面位场信息用于获得目标区域内的、一个或多个尺度下的各个检测点的总模，并可以根据总模的强度将所有检测点分为强检测点和非强检测点。在一个或多个实施例中，当增大尺度后，强检测点数量变化较小、或者几乎不再变化时，即可不再增加尺度，在部分优选的实施例中，尺度的数量为8～10个。在一个或多个实施例中，强检测点和非强检测点可以基于经验设置阈值，大于等于阈值的为强检测点，小于阈值的为非强检测点，强检测点和非强检测点也可以通过自适应双阈值法来将各检测点按照总模强度分为至少两类。在一个实施例中，非强检测点还可以进一步划分为弱检测点和非检测点。
10、本技术方案中，从遥感资料或者从地质踏勘资料中解译出线性构造单元信息为现有技术。基于解译出的一个或多个线性构造单元可以对应地建立起每个线性构造单元的封闭曲线，封闭曲线内的区域则作为该线性构造单元的缓冲区。
11、本技术方案中，需要判断强检测点和封闭曲线的相对位置，从而统计每条封闭曲线内包含的各个尺度下的强检测点的总数量。对于强检测点总数量大于预设值的封闭曲线，则表明该封闭曲线对应的线性构造单元可以同时被重磁位场方法和遥感、地质踏勘观察到，反之，则说明该封闭曲线对应的线性构造单元难以被重磁位场方法观察到。最后，基于能够同时被重磁位场方法和遥感、地质踏勘观察到的各线性构造单元来构建起区域控岩，即磁性岩，主要指岩浆岩、磁性变质岩，来构造骨架。在一个或多个实施例中，可以根据封闭曲线内的强检测点总数来确定对应的线性构造单元的控制等级，并以该控制等级作为线宽来制图，从而能够快速、直观地掌握各线性构造单元能够被多数据源共同标记的程度。
12、本技术方案中，通过重磁资料获得强检测点数量，通过遥感资料或地质踏勘资料获得的线性构造单元构建缓冲区，再基于缓冲区内的强检测点总数来直观地反应出各线性构造单元能够同时被重磁位场方法和遥感、地质踏勘观察到，以及观察到的程度，能够解决重磁资料多解性的问题，更加精准地定位构造位置，同时快速地、直观地划分出重点勘察区域，具有广泛地应用价值。
13、作为本专利技术的一种优选实施方式，所述步骤s2包括以下步骤：
14、s21：确定尺度因子队列，所述尺度因子队列包括至少一个尺度因子；
15、s22：从所述重磁资料中获取观察面位场信息，基于所述尺度因子和观察面位场信息得到当前尺度下处理后的位场信息；
16、s23：根据处理后的位场信息在x方向、y方向上的梯度，计算目标区域内各检测点的总模和梯度方向；
17、s24：判断各检测点的总模是否在所述梯度方向上为极大值，若是则保留该检测点的总模，若否则将该检测点的总模设置为零；
18、s25：将所有总模非零的检测点按照总模的强度划分为强检测点和非强检测点；
19、s26：重复步骤s22至s25，直至获得所述尺度因子队列中所有尺度因子下的强检测点和非强检测点。
20、本技术方案中，在步骤s21中首先确定包括至少一个尺度因子的尺度因子队列，例如，尺度因子队列可以是z＝[z1,z2,z3...zi]。尺度因子队列中的每一个尺度因子均进行步骤s22至s25，以获得所有尺度因子下的强检测点和非强检测点。尺度因子可以看作是相对于观察面位场信息中的重磁测量值高度z在高度方向上的不同延拓。
21、本技术方案中，步骤s22至s25中，对于每一个尺度因子，首先利用尺度因子和观察面位场信息获得当前尺度下处理后的位场信息，之后，根据处理后的位场信息在x方向、y方向上的梯度，可以得到各检测点的总模和梯度方向，在步骤s24中，利用非极大值抑制可以使各检测点沿着梯度方向只保留最大的一个值，从而当连接这些点以后便在空间上形成了边缘，接下来，在步骤s25中，基于总模的强度，通过经验或者计算方法，例如自适应双阈值法能够将总模非零的检测点进行分类，得到强检测点和非强检测点。
22、在一个或多个实施例中，在步骤s24中，可以将梯度方向离散化为5个区间，例如再针对二维总模上的每一个检测点，判断检测点的总模是否在梯度方向上是极大值，如果不是最大值则将该点的梯度总模置于0。
23、本技术方案中，在获得各检测点的总模后，利用非极大值抑制会使得位场数据在边缘方向上建构起一个个薄边，利于后续数据处理，同时，基于总模的强度的划分检测点，能够得到在不同尺度下观察的强检测点，减少了干扰点的影响，进一步提高了重点勘察区域划分的准确性。
24、进一步地，所述步骤s22中，对于任一尺度因子zi，构建光滑函数基于所述光滑函数与观察面位场信息tz(x,y,z)的二维空间褶积得到zi尺度下处理后的位场信息tzi(x,y,z本文档来自技高网...
【技术保护点】
1.融合多数据源的地质构造分析方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的融合多数据源的地质构造分析方法，其特征在于，所述步骤S2包括以下步骤：
3.根据权利要求2所述的融合多数据源的地质构造分析方法，其特征在于，所述步骤S22中，对于任一尺度因子zi，构建光滑函数基于所述光滑函数与观察面位场信息Tz(x,y,z)的二维空间褶积得到zi尺度下处理后的位场信息Tzi(x,y,zi)。
4.根据权利要求3所述的融合多数据源的地质构造分析方法，其特征在于，所述步骤S23中，所述总模fm_zi、梯度方向fa_zi的计算公式为：
5.根据权利要求2所述的融合多数据源的地质构造分析方法，其特征在于，所述步骤S25中，所有总模非零的检测点构成集合P＝{p|p＝fm_zi(x,y)andp≠0}，按照集合P中各检测点的总模强度将集合P划分为n级区间，第k级的区间范围为区间累积概率为其中，max(p)、min(p)分别代表集合P中的总模最大值和总模最小值；
6.根据权利要求5所述的融合多数据源的地质构造分析方法，其特征在于，
7.根据权利要求1所述的融合多数据源的地质构造分析方法，其特征在于，所述步骤S3中，具有s个控制点的线性构造单元对应的封闭曲线L的计算公式为：
8.根据权利要求1所述的融合多数据源的地质构造分析方法，其特征在于，所述步骤S4中，生成一条以强检测点作为起点的水平射线，若所述水平射线与封闭曲线的交点为一个，则所述强检测点位于所述封闭曲线内，若所述水平射线与封闭曲线的交点为两个，则所述强检测点位于所述封闭曲线外。
9.根据权利要求1所述的融合多数据源的地质构造分析方法，其特征在于，所述步骤S4中，所述封闭曲线所对应的线性构造单元的控制等级Ym＝{zi|Cm_zi＝max({Cm_i,i∈所有的检测尺度等级})}，以所述控制等级作为线宽，构建起区域控岩构造骨架，其中，Cm_i为第i个尺度下，所述封闭曲线包含个层级的强检测点数量。
10.一种存储介质，其特征在于，包括存储的计算机程序，其中，在所述计算机程序运行时，控制所述存储介质所在的设备执行如权利要求1～9中任一项所述的融合多数据源的地质构造分析方法。
...【技术特征摘要】
1.融合多数据源的地质构造分析方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的融合多数据源的地质构造分析方法，其特征在于，所述步骤s2包括以下步骤：
3.根据权利要求2所述的融合多数据源的地质构造分析方法，其特征在于，所述步骤s22中，对于任一尺度因子zi，构建光滑函数基于所述光滑函数与观察面位场信息tz(x,y,z)的二维空间褶积得到zi尺度下处理后的位场信息tzi(x,y,zi)。
4.根据权利要求3所述的融合多数据源的地质构造分析方法，其特征在于，所述步骤s23中，所述总模fm_zi、梯度方向fa_zi的计算公式为：
5.根据权利要求2所述的融合多数据源的地质构造分析方法，其特征在于，所述步骤s25中，所有总模非零的检测点构成集合p＝{p|p＝fm_zi(x,y)andp≠0}，按照集合p中各检测点的总模强度将集合p划分为n级区间，第k级的区间范围为区间累积概率为其中，max(p)、min(p)分别代表集合p中的总模最大值和总模最小值；
6.根据权利要求5所述的融合多数据源的地质构造分析方法，其特征在于，遍历所有弱检测点，若所...
【专利技术属性】
技术研发人员：周春波，沈开红，
申请(专利权)人：四川省地球物理调查研究所，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。