基于三维激光扫描技术的地质露头裂缝识别方法与流程

1.本发明涉及油气地质勘探及工程测量技术领域，尤其涉及一种基于三维激光扫描技术的地质露头裂缝识别方法。


背景技术：

2.随着全球性油气勘测和工程测量技术水平的提高，裂缝研究已经逐渐成为一个重要的勘测新领域。在油气勘探领域，裂缝既是油气储集的空间，也是油气运移的通道，在油气藏的勘探开发中有着非常重要的地位。但是由于裂缝的成因复杂，控制和影响的因素众多，形成和发育随机，且油气目的层常位于地下深部，裂缝的产状和规模很难直接测量，研究难度较高。通过选取与油气目的层的地质情况相似的典型地质露头剖面进行裂缝研究，可为地下目的层的裂缝研究提供可靠依据。
3.野外地质露头直接测量裂缝信息准确且直接，是进行裂缝研究的最佳选择之一，但是往往耗时、危险，甚至很多地方无法依靠人工完成。因此，寻找有效、安全而又准确的提取裂缝相关信息的新方法非常必要。


技术实现要素：

4.为了克服上述现有技术的不足，本发明的目的在于提供一种基于三维激光扫描技术的地质露头裂缝识别方法，利用三维激光扫描技术获取地质露头剖面的点云数据，在点云数据上对裂缝位置、产状和规模进行精准识别，能够为油气储层的裂缝识别与预测提供科学指导，具有重要的理论价值和广阔的应用前景，具有快速性、高精度、智能化水平高的优点。
5.为了实现以上技术目的，本发明提供以下技术方案。
6.一种基于三维激光扫描技术的地质露头裂缝识别方法，具体步骤如下：
7.一、用三维激光扫描仪获取地质露头剖面三维点云数据；
8.二、对步骤一获取的地质露头剖面三维点云数据进行滤波去噪；
9.三、对步骤二去噪处理后的三维点云数据进行二维转换；
10.四、对步骤三得到的二维点云数据进行裂缝轮廓的识别；
11.五、对步骤四识别出的裂缝进行包括产状、长度、宽度和面密度的关键参数提取。
12.所述步骤一的具体方法为：
13.采用三维激光扫描仪对地质露头剖面进行扫描，获取地质露头剖面表面的三维点云数据；以三维激光扫描仪为坐标中心建立三维笛卡尔坐标系，并标定点云数据的笛卡尔坐标；地质露头剖面表面的三维点云数据的坐标标定为p(xi，yi，zi)，i＝1,2,3......n，n为点云数据中点的数量；x、y、z为地质露头剖面表面点的三维坐标，i代表第i个点。
14.所述步骤二的具体方法为：
15.利用点云处理软件对步骤一获取的三维点云数据中偏离露头主体的孤立点和离散点直接手动删除。
16.所述点云处理软件包括cyclone，realworks或pointtools。
17.所述步骤三的具体方法为：
18.将步骤二滤波去噪后的三维点云数据经投影变化后，转为二维数据，具体方法如下：
19.第1步：以三维激光扫描仪为坐标中心，建立三维柱坐标系，并建立三维柱坐标系与步骤一中的三维笛卡尔坐标系的转换关系；将笛卡尔空间坐标系中的三维点云对应到柱坐标系中；
20.第2步：将三维柱坐标系的柱面展开，构建二维柱面坐标系，并建立三维柱坐标系与二维柱面坐标系的转换关系；将柱坐标系中的三维点云数据转换为二维数据。
21.所述步骤四的具体方法为：
22.利用delaunay三角剖分改进的alpha shapes模型对步骤三得到的二维点云数据进行裂缝轮廓识别，具体实现步骤如下：
23.第1步：建立步骤三中点集的delaunay三角网；
24.第2步：对delaunay三角网的各边依次进行是否符合alpha shape模型判断，当delaunay三角网中某条边的长度大于da，则该边不符合alpha shapes要求，删除该边；
25.第3步：对其余边依次进行是否符合alpha shape模型判断，求出每条边以da为直径的两个外接圆圆心位置；建立该边对应的检索格网，对检索格网内的任意点到两个外接圆圆心的距离进行判断，识别出符合alpha shapes要求的边，边长≤da即为符合，直至所有边均判断完毕；
26.第4步：所有符合alpha shapes要求的边，构成的轮廓线即为裂缝的边缘线。
27.所述步骤四第3步中da值的选取计算公式为：
[0028][0029]
其中，w是所选取点云的宽度，h是所选取点云的长度，np是点云中点的数量，3是通过试验获取的经验值。
[0030]
所述步骤五的具体方法为：
[0031]
第1步：提取裂缝产状
[0032]
用最小二乘法对平面进行拟合，空间平面的方程用ax by cz 1＝0来表示，由平面方程可以确定裂缝面的法向量为n(a，b，c)，当法向量的z坐标小于0，则给法向量乘以
‑
1将其z坐标变为正值，即保证法向量是垂直于裂缝面向上；根据法向量的三个坐标判断裂缝面的方位，排除多点共线、平面水平以及平面垂直这三种特殊情况，确定裂缝面为倾斜平面，即可以计算倾向和倾角之后计算裂缝面的空间方位，倾向为法向量在水平面上的投影向量与y轴正方向之间的夹角，倾角为法向量与法向量在水平面上的投影向量之间的夹角；
[0033]
第2步：分段测量提取裂缝的长度和宽度
[0034]
先根据裂缝骨架即裂缝内全部最大内切圆的圆心集合与裂缝上下边界点的距离，将整个裂缝分若干段，再针对每段裂缝，确定该段裂缝的最小外接矩形，再根据矩形边界，完成裂缝长度和宽度的计算；
[0035]
第3步:裂缝面密度的计算。
[0036]
面密度指单位面积内裂缝的累积长度，在计算出裂缝的长度后，可由下式阶段裂
缝面密度：
[0037][0038]
其中f为面密度，单位l/m，l为每条裂缝的长度，单位m；s代表检测区域的面积，单位m2。
[0039]
本发明的有益效果为：
[0040]
本发明利用三维激光扫描仪获取地质露头剖面三维点云数据，对获取的三维点云数据进行滤波去噪、二维转换后进行裂缝轮廓的识别，并对识别出的裂缝进行包括产状、长度、宽度和面密度的关键参数提取；能够准确、高效又安全地对露头剖面中裂缝进行定量识别，具有重要的理论价值和广阔的应用前景，为地下油气目的层的裂缝研究提供可靠依据。
具体实施方式
[0041]
下面对本发明做进一步详细说明：
[0042]
一种基于三维激光扫描技术的地质露头裂缝识别方法，具体步骤如下：
[0043]
一、用三维激光扫描仪获取地质露头剖面三维点云数据；
[0044]
采用三维激光扫描仪对地质露头剖面进行扫描，获取地质露头剖面表面的三维点云数据；以三维激光扫描仪为坐标中心建立三维笛卡尔坐标系，并标定点云数据的笛卡尔坐标；地质露头剖面表面的三维点云数据的坐标标定为p(xi，yi，zi)，i＝1,2,3......n，n为点云数据中点的数量；x、y、z为地质露头剖面表面点的三维坐标，i代表第i个点。
[0045]
二、对步骤一获取的地质露头剖面三维点云数据进行滤波去噪；
[0046]
利用点云处理软件对步骤一获取的三维点云数据中偏离露头主体的孤立点和离散点直接手动删除。
[0047]
三、对步骤二去噪处理后的三维点云数据进行二维转换；
[0048]
将步骤二滤波去噪后的三维点云数据经投影变化后转为二维数据，具体方法如下：
[0049]
第1步：以三维激光扫描仪为坐标中心，建立三维柱坐标系，并建立三维柱坐标系与步骤一中的三维笛卡尔坐标系的转换关系；将笛卡尔空间坐标系中的三维点云对应到柱坐标系中；
[0050]
第2步：将三维柱坐标系的柱面展开，构建二维柱面坐标系，并建立三维柱坐标系与二维柱面坐标系的转换关系；将柱坐标系中的三维点云数据转换为二维数据。
[0051]
四、对步骤三得到的二维点云数据进行裂缝轮廓的识别；
[0052]
利用delaunay三角剖分改进的alpha shapes模型对步骤三得到的二维点云数据进行裂缝轮廓识别，具体实现步骤如下：
[0053]
第1步：建立步骤三中点集的delaunay三角网；
[0054]
第2步：对delaunay三角网的各边依次进行是否符合alpha shape模型判断，当delaunay三角网中某条边的长度大于da，则该边不符合alpha shapes要求，删除该边；
[0055]
第3步：对其余边依次进行是否符合alpha shape模型判断，求出每条边以da为直径的两个外接圆圆心位置；建立该边对应的检索格网，对检索格网内的任意点到两个外接圆圆心的距离进行判断，识别出符合alpha shapes要求的边，边长≤da即为符合，直至所有
边均判断完毕；
[0056]
第4步：所有符合alpha shapes要求的边，构成的轮廓线即为裂缝的边缘线。
[0057]
所述步骤四第3步中da值的选取计算公式为：
[0058][0059]
其中，w是所选取点云的宽度，h是所选取点云的长度，np是点云中点的数量，3是通过试验获取的经验值。
[0060]
五、对步骤四识别出的裂缝进行包括产状、长度、宽度和面密度的关键参数提取；
[0061]
第1步：提取裂缝产状
[0062]
用最小二乘法对平面进行拟合，空间平面的方程用ax by cz 1＝0来表示，由平面方程可以确定裂缝面的法向量为n(a，b，c)，当法向量的z坐标小于0，则给法向量乘以
‑
1将其z坐标变为正值，即保证法向量是垂直于裂缝面向上；根据法向量的三个坐标判断裂缝面的方位，排除多点共线、平面水平以及平面垂直这三种特殊情况，确定裂缝面为倾斜平面，即可以计算倾向和倾角之后计算裂缝面的空间方位，倾向为法向量在水平面上的投影向量与y轴正方向之间的夹角，倾角为法向量与法向量在水平面上的投影向量之间的夹角；
[0063]
第2步：分段测量提取裂缝的长度和宽度
[0064]
先根据裂缝骨架即裂缝内全部最大内切圆的圆心集合与裂缝上下边界点的距离，将整个裂缝分若干段，再针对每段裂缝，确定该段裂缝的最小外接矩形，再根据矩形边界，完成裂缝长度和宽度的计算；
[0065]
第3步:裂缝面密度的计算。
[0066]
面密度指单位面积内裂缝的累积长度，在计算出裂缝的长度后，可由下式阶段裂缝面密度：
[0067][0068]
其中f为面密度，单位l/m，l为每条裂缝的长度，单位m；s代表检测区域的面积，单位m2。
[0069]
所述步骤二中点云处理软件包括cyclone，realworks，pointtools。