基于地面激光扫描的碳酸盐岩白云石含量模拟方法

技术领域

本发明涉及一种基于地面激光扫描的碳酸盐岩白云石含量模拟方法，属于遥感技术在碳酸盐岩精细地质中的应用技术领域。

背景技术

碳酸盐岩是重要的储油岩，全世界50％的石油和天然气储存于碳酸盐岩中。碳酸盐岩主要是由方解石和白云石等组成的沉积岩。以方解石为主的岩石称为石灰石，以白云石为主的岩石称为白云岩。方解石的主要成分是碳酸钙，这种沉积岩比较坚实，孔隙度差，渗透率低，地下的流体往往不能很好的渗透其中。碳酸钙在沉积的地质构造变化过程中，由于古代海水的侵蚀，海水中大量的镁离子和碳酸钙中的钙离子进行交换，这个过程即是白云岩化，即形成白云石。白云石的裂缝和孔洞比较发育，渗透性好，对油气成藏提供了良好的环境。因此，白云石的含量对于研究碳酸盐岩储层至关重要。

方解石和白云石用肉眼是难以区分的，一般要通过薄片镜下观察或成分分析来区分。只能对野外地质剖面的局部样品进行测定。

因此，提供一种新型的基于地面激光扫描的碳酸盐岩白云石含量模拟方法已经成为本领域亟需解决的技术问题。

发明内容

为了解决上述的缺点和不足，本发明的目的在于提供一种基于地面激光扫描的碳酸盐岩白云石含量模拟方法。

为了实现以上目的，本发明提供了一种基于地面激光扫描的碳酸盐岩白云石含量模拟方法，其中，所述基于地面激光扫描的碳酸盐岩白云石含量模拟方法包括：

步骤1，获取野外碳酸盐岩储层地质剖面数据；

步骤2，建立碳酸盐岩岩石样品激光强度值和白云石含量之间的关系模型，并将其作为碳酸盐岩岩石样品白云石含量的预测模型；

步骤3，根据所述野外碳酸盐岩储层地质剖面数据，建立碳酸盐岩储层数字地质剖面表层模型；

步骤4，利用所述碳酸盐岩岩石样品白云石含量的预测模型对碳酸盐岩储层数字地质剖面表层模型所有节点进行白云石含量预测表征，得到碳酸盐岩储层地质剖面白云石含量的平面分布。

根据本发明具体实施方案，在所述的方法中，优选地，步骤1中所述获取野外碳酸盐岩储层地质剖面数据，包括：

利用地面激光雷达扫描仪对野外碳酸盐岩储层进行扫描，以获得野外碳酸盐岩储层地质剖面表层三维点云，从而精确描述地质剖面表层的相对空间几何信息；

利用高分辨率数码相机获取野外碳酸盐岩储层地质剖面的高精度纹理影像。

根据本发明具体实施方案，在所述的方法中，优选地，步骤2中所述建立碳酸盐岩岩石样品激光强度值和白云石含量之间的关系模型，并将其作为碳酸盐岩岩石样品白云石含量的预测模型，包括：

利用地面激光雷达扫描仪对碳酸盐岩岩石样品进行扫描，以获得碳酸盐岩岩石样品的点云数据，所述点云数据包括三维空间坐标(X、Y、Z)和激光强度值I；

计算碳酸盐岩岩石样品的激光强度值和白云石含量数据之间的相关系数；

在碳酸盐岩岩石样品的激光强度值和白云石含量具有显著相关性的基础上，建立碳酸盐岩岩石样品白云石含量的预测模型。

根据本发明具体实施方案，在所述的方法中，优选地，当碳酸盐岩岩石样品的激光强度值和白云石含量的相关系数为0.8以上时，确定二者具有显著相关性。

所述步骤2中，利用地面激光雷达扫描仪对多个岩石样品进行扫描，获得岩石样品的点云数据，包括三维空间坐标(X、Y、Z)和激光强度值I；同时进行岩心分析，测定岩石样品的白云石含量。

在接收机内部，I通常被表示为接收信号的峰值振幅。假设强度值与接收信号功率之间存在简单的线性关系，如式(1)所示：

式(1)中，C为常数。

从式(1)中可以看出，强度值的影响因素有激光发射功率P

因此，式(1)可化简为式(2)：

式(2)中，

由式(2)可以看出，地面激光雷达强度值与目标反射率，激光测距值和激光入射角有关。在地面激光扫描时，设置合理的扫描距离和角度，可减弱扫描距离和角度对激光强度值的影响，使激光强度值更大程度上反映目标反射特性。

不同白云石含量的碳酸盐岩岩石样品具有不同的成分、颜色等特征，因而具有不同的反射特性，使得利用激光强度值预测碳酸盐岩白云石含量成为可能。利用碳酸盐岩岩石样品的激光强度值和白云石含量数据，并计算两者之间的相关系数，在两者具有显著相关性的基础上，建立碳酸盐岩岩石样品激光强度值和白云石含量的关系模型，作为白云石含量的预测模型。预测模型可包括线性模型和非线性模型。

根据本发明具体实施方案，在所述的方法中，优选地，步骤3中所述根据所述野外碳酸盐岩储层地质剖面数据，建立碳酸盐岩储层数字地质剖面表层模型，包括：

将多站扫描所得的点云数据进行拼接，形成完整的地质剖面点云数据；

人工手动选择高精度纹理影像与点云的匹配点，计算配准参数并纠正高精度纹理影像，使得每个三维点云都有RGB颜色信息；

采用基于最优趋势面构建三角网的方法对不规则海量点云进行建模，得到碳酸盐岩储层数字地质剖面表层模型。

根据本发明具体实施方案，在所述的方法中，优选地，所述采用基于最优趋势面构建三角网的方法对不规则海量点云进行建模，得到碳酸盐岩储层数字地质剖面表层模型，包括：

将所有点云向各个方向进行投影，选择投影面积最大的方向的平面作为最优趋势面，将所有点云投影到此最优趋势面上，在平面上建立不规则三角网模型，再把平面三角网通过高程值还原到三维空间中，从而形成具有颜色纹理信息的数字地质剖面模型。

根据本发明具体实施方案，优选地，该方法具体包括以下步骤：

步骤1，利用地面激光雷达扫描仪对野外碳酸盐岩储层进行扫描，以获得野外碳酸盐岩储层地质剖面表层三维点云；

利用高分辨率数码相机获取野外碳酸盐岩储层地质剖面的高精度纹理影像；

步骤2，利用地面激光雷达扫描仪对碳酸盐岩岩石样品进行扫描，以获得碳酸盐岩岩石样品的点云数据，所述点云数据包括三维空间坐标(X、Y、Z)和激光强度值I；

计算碳酸盐岩岩石样品的激光强度值和白云石含量数据之间的相关系数；

在碳酸盐岩岩石样品的激光强度值和白云石含量具有显著相关性的基础上，建立碳酸盐岩岩石样品白云石含量的预测模型；

步骤3，将多站扫描所得的点云数据进行拼接，形成完整的地质剖面点云数据；

人工手动选择高精度纹理影像与点云的匹配点，计算配准参数并纠正高精度纹理影像，使得每个三维点云都有RGB颜色信息；

采用基于最优趋势面构建三角网的方法对不规则海量点云进行建模，得到碳酸盐岩储层数字地质剖面表层模型；

步骤4，利用所述碳酸盐岩岩石样品白云石含量的预测模型对碳酸盐岩储层数字地质剖面表层模型所有节点进行白云石含量预测表征，得到碳酸盐岩储层地质剖面白云石含量的平面分布。

根据本发明具体实施方案，在所述的方法中，优选地，所述采用基于最优趋势面构建三角网的方法对不规则海量点云进行建模，得到碳酸盐岩储层数字地质剖面表层模型，包括：

将所有点云向各个方向进行投影，选择投影面积最大的方向的平面作为最优趋势面，将所有点云投影到此最优趋势面上，在平面上建立不规则三角网模型，再把平面三角网通过高程值还原到三维空间中，从而形成具有颜色纹理信息的数字地质剖面模型。

根据本发明具体实施方案，在所述的方法中，优选地，当碳酸盐岩岩石样品的激光强度值和白云石含量的相关系数为0.8以上时，确定二者具有显著相关性。

目前，本领域常用的传统三角网建模都是将点投影到水平面建网，但是该种方法不适用于地质剖面垂向信息丰富的建网需求，因此，在本申请的该方法中采用基于最优趋势面构建三角网的方法对不规则海量点云进行建模，以得到碳酸盐岩储层数字地质剖面表层模型。

本发明所提供的该方法包括：对野外碳酸盐岩储层地质剖面进行地面激光扫描、高精度数码照片采集以及碳酸盐岩样品采样；室内测量岩样的地面激光强度值和白云石含量，建立激光强度值与白云石含量的关系模型，作为白云石含量的预测模型；利用地面激光点云数据和高精度数码照片构建高精度储层数字地质剖面表层模型；应用预测模型，预测碳酸盐岩储层地质剖面白云岩含量分布。本方法具有预测快速、准确等优点，为碳酸盐岩储层剖面的白云岩含量精细与定量表征提供一种新思路。

本发明所提供的该基于地面激光扫描的碳酸盐岩白云石含量模拟方法可以快速模拟地质剖面的白云石含量的展布，该方法中通过地面激光雷达扫描获得的激光强度值数据反应了扫描目标对激光的反射能力，利用不同白云石含量具有不同的反射能力的特点，可以建立激光强度值与白云石含量的关系模型，从而实现数字地质剖面上白云石含量的平面预测，为碳酸盐岩储层地质剖面白云石含量精细与定量表征提供一种新方法，并且该方法在油气勘探开发及地质应用中具有较好的应用前景。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例中所提供的基于地面激光扫描的碳酸盐岩白云石含量模拟方法的具体工艺流程图。

图2为本发明实施例中利用地面激光雷达扫描仪采集到的四川盆地盘龙洞碳酸盐岩剖面表层三维点云带激光强度信息图。

图3为本发明实施例中四川盆地盘龙洞碳酸盐岩剖面的高精度纹理影像。

图4为本发明实施例中所建立得不规则三角网模型的示意图。

图5为本发明实施例中最优趋势面建模示意图。

图6为本发明实施例中线性预测模型得到的预测白云石含量与实测白云石含量的对比图。

具体实施方式

为了对本发明的技术特征、目的和有益效果有更加清楚的理解，现结合以下具体实施例对本发明的技术方案进行以下详细说明，但不能理解为对本发明的可实施范围的限定。

图1显示了本实施例所提供的基于地面激光扫描的碳酸盐岩白云石含量模拟方法的具体工艺流程图，如图1所示，该方法包括以下步骤：

1.野外地质剖面数据采集及岩石样品采样

野外地质剖面数据采集系统包括地面激光雷达扫描仪和高分辨率数码相机。本实施例中，选用的地面激光雷达设备是奥地利rigel-vz400型地面激光雷达扫描仪，本实施例中利用该地面激光雷达扫描仪采集地质剖面表层三维点云(如图2所示)，设定扫描距离和扫描点间距；其中，扫描距离为10-15m，点间距为1mm；

本实施例中选用的高分辨率数码相机为宾得645D，其像素可达到4000万，用于获取地质剖面(本实施例中为四川盆地盘龙洞碳酸盐岩剖面)的高精度纹理影像(如图3所示)；

本实施例中，岩石采样是对不同白云石含量的碳酸盐岩进行采样，共采集10块样品，并通过实验分析获得岩样中白云石和方解石的含量。

2.建立碳酸盐岩白云石含量的预测模型

在室内，对步骤1中所采集的10块碳酸盐岩样品利用地面激光雷达扫描仪进行扫描，扫描距离设定约为15米，点云间距为1mm。扫描后获得岩石样品的点云数据，该点云数据包括三维空间坐标(X、Y、Z)和激光强度值I；

计算10块岩样的激光强度值与白云石含量的相关系数，计算得到R

在此基础上，以传统线性模型为例，建立碳酸盐岩岩石样品白云石含量的预测模型，该模型为：C＝-35.595×I-34.262，其中，C为白云石含量，相关关系图如图6所示。

3.建立碳酸盐岩储层数字地质剖面表层模型

(1)点云处理：对地面激光雷达扫描仪和高分辨率数码相机采集的数据进行处理，具体地，利用该地面激光雷达扫描仪自带的处理软件将多站扫描的点云数据进行拼接，形成完整的地质剖面点云数据；然后利用该地面激光雷达扫描仪自带的处理软件，人工手动选择照片(地质剖面的高精度纹理影像)与点云的匹配点，计算配准参数并纠正照片，以实现每个三维点云都有RGB颜色信息；

(2)模型建立：采用基于最优趋势面构建三角网的方法对不规则海量点云进行建模，具体地，将所有点云向各个方向进行投影，选择投影面积最大的方向的平面作为最优趋势面，将所有点云投影到此最优趋势面(如图5所示)上，在平面上建立不规则三角网模型(如图4所示)，再把平面三角网通过高程值还原到三维空间中，从而形成具有颜色纹理信息的数字地质剖面模型。

4.数字地质剖面模型上白云石含量表征

利用所述碳酸盐岩岩石样品白云石含量的预测模型对碳酸盐岩储层数字地质剖面表层模型所有节点进行白云石含量预测表征，得到碳酸盐岩储层地质剖面白云石含量的平面分布，具体地，碳酸盐岩储层数字地质剖面表层模型上的所有节点均具有一个激光强度值，应用所述碳酸盐岩岩石样品白云石含量的预测模型，可得到所有节点处的白云石含量，进而实现对碳酸盐岩储层地质剖面上白云石含量的平面预测表征，实验数据如下表1所示。

表1

表1中，白云石含量、方解石含量及石英含量为在显微镜下观察岩石样品薄片后所计算得到的面积百分含量。

以上所述，仅为本发明的具体实施例，不能以其限定发明实施的范围，所以其等同组件的置换，或依本发明专利保护范围所作的等同变化与修饰，都应仍属于本专利涵盖的范畴。另外，本发明中的技术特征与技术特征之间、技术特征与技术发明之间、技术发明与技术发明之间均可以自由组合使用。