地层分布模式约束的钻孔对比剖面模型构建方法及装置制造方法及图纸

本发明专利技术公开了一种地层分布模式约束的钻孔对比剖面模型构建方法及装置，方法包括：(1)读取剖面线上实测钻孔的钻孔集合D和钻孔地层厚度集合T；(2)从钻孔集合D中读取任一相邻实测钻孔组合，解析地层分布模式，形成地层分布模式标志列表；(3)插值生成所述相邻实测钻孔组合之间的虚拟钻孔，与所述相邻实测钻孔组合形成钻孔集合PD；(4)基于地层尖灭规律，根据地层分布模式标志列表M初步推断钻孔集合PD中各钻孔的各地层厚度，存入地层厚度矩阵PT；(5)优化地层厚度，构建所述相邻实测钻孔组合间所有地层面；(6)根据构建的地层面和钻孔集合PD生成对比剖面图。本发明专利技术建模质量更高，效率更高，自动化程度更高。自动化程度更高。自动化程度更高。

【技术实现步骤摘要】
地层分布模式约束的钻孔对比剖面模型构建方法及装置


[0001]本专利技术涉及地理信息技术，尤其涉及一种地层分布模式约束的钻孔对比剖面模型构建方法及装置。

技术介绍

[0002]基于钻孔构建地层对比剖面是地层学家对比地层的主要手段。其主要构建方法为将钻孔和地层的属性数据转换为用于绘制的几何要素，进而实现剖面数据的可视化。地层在空间上分布不连续或者厚度的不均匀，往往增加了构建地层剖面的复杂性。因此，地层分布模式的科学认知与推理规则的正确运用是在复杂情况下，合理构建地层对比剖面的基础和前提。
[0003]现有方法多是采用了简单的地层尖灭处理规则，导致多地层连续尖灭、多地层间断缺失等复杂情况下的剖面建模质量不够理想，与地层发育规律有所冲突。此外，所生成的剖面虽可进一步依靠专家手工编辑和优化，但效率低下，且剖面建模质量不够稳定。

技术实现思路

[0004]专利技术目的：本专利技术针对现有技术存在的问题，提供一种效率更高、质量更好的地层分布模式约束的钻孔对比剖面模型构建方法及装置。
[0005]技术方案：本专利技术所述的地层分布模式约束的钻孔对比剖面模型构建方法包括：
[0006](1)读取剖面线上实测钻孔的钻孔集合D和钻孔地层厚度集合T；
[0007](2)从钻孔集合D中读取任一相邻实测钻孔组合，基于钻孔地层厚度集合T建立相邻实测钻孔组合间的缺失提取地层的缺失标志列表SL，并根据SL构建地层缺失区间集合DS，再基于DS解析地层分布模式，形成地层分布模式标志列表M；
[0008](3)插值生成所述相邻实测钻孔组合之间的虚拟钻孔，并推断各虚拟钻孔的底面标高和地面标高，将所述相邻实测钻孔组合和之间的虚拟钻孔存入钻孔集合PD；
[0009](4)基于地层尖灭规律，根据地层分布模式标志列表M初步推断钻孔集合PD中各钻孔的各地层厚度，存入地层厚度矩阵PT；
[0010](5)根据各虚拟钻孔的底面标高和地面标高优化地层厚度，并构建所述相邻实测钻孔组合间所有地层面；
[0011](6)根据构建的地层面和钻孔集合PD生成所述相邻实测钻孔组合间的对比剖面图；
[0012](7)循环执行步骤(2)
‑
(6)，直至钻孔集合D中所有相邻实测钻孔组合被遍历，得到所有对比剖面图。
[0013]进一步的，步骤(1)包括；
[0014](1
‑
1)读取剖面线上实测钻孔的钻孔信息，得到钻孔集合D＝{d
i
|i＝1，2，...，dn}；其中，d
i
表示第i个钻孔，dn表示钻孔数量；
[0015](1
‑
2)读取剖面线上实测钻孔的钻孔地层信息，得到地层厚度集合T＝{t
i，j
|i＝1，
2，...，dn；j＝1，2，...，sni}；其中，t
i，j
表示钻孔d
i
的第j个地层的厚度，sni表示钻孔d
i
钻遇的地层总数。
[0016]进一步的，步骤(2)包括；
[0017](2
‑
1)从钻孔集合D中读取任意一组相邻实测钻孔d
i
和d
i+1
；
[0018](2
‑
2)遍历钻孔地层厚度集合T，根据下式匹配钻孔d
i
和钻孔d
i+1
间的缺失提取地层的缺失标志列表SL＝{sl
j
|j＝1，2，...，sni}；
[0019][0020]其中，sl
j
表示地层j的地层缺失标志，t
i，j
、t
i+1，j
表示钻孔d
i
、d
i+1
中地层j的厚度，sni表示钻孔d
i
钻遇的地层总数；
[0021](2
‑
3)基于地层缺失标志列表SL
i
，构建钻孔d
i
和钻孔d
i+1
间地层缺失区间集合DS＝{ds
u
|u＝1，2，...，dsn}；其中，ds
u
＝(start，end，count)为一个三元组，表示第u个地层缺失区间，dsn为缺失区间数量，start为区间ds
u
的起始地层编号，end为区间ds
u
的终止地层编号，count为区间ds
u
中的地层个数；
[0022](2
‑
4)对于地层缺失区间集合DS中每一地层缺失区间ds
u
，遍历区间ds
u
内的所有地层的缺失情况，并按照下式计算区间ds
u
内的地层分布模式标志m
u
，存入地层分布模式标志列表M＝{m
u
|u＝1，2，...，dsn}；
[0023][0024]进一步的，步骤(2
‑
3)包括；
[0025](2
‑3‑
1)顺序遍历缺失标志列表SL，直至sl
j
不等于0，将j的值赋予start，并将count值赋为1；
[0026](2
‑3‑
2)继续遍历缺失标志列表SL，当j小于等于sni且|sl
j
|等于|sl
start
|时，令count值递增1；否则，将end值赋为j
‑
1，新建区间元素ds
u
＝(start，end，count)，并添加至地层缺失区间集合DS，之后返回至步骤(2
‑3‑
2)；其中，||为绝对值符号；
[0027](2
‑3‑
3)循环执行步骤(2
‑3‑
1)至(2
‑3‑
2)，直至缺失标志列表SL被遍历完，得到地层缺失区间集合DS。
[0028]进一步的，步骤(3)包括；
[0029](3
‑
1)遍历集合DS，获取DS中所有ds
u
(count)的最大值，记为最大区间缺失地层数max，作为待插入的虚拟钻孔数量，其中，ds
u
(count)为区间ds
u
中的地层个数；
[0030](3
‑
2)新建空的钻孔集合PD＝{pd
k
|k＝1，2，...，max+2}，存储虚拟钻孔及读取的相邻实测钻孔d
i
和d
i+1
；其中，pd1值赋为实测钻孔d
i
，pd
max+2
值赋为实测钻孔d
i+1
，pd2至pd
max+1
依次为第1至max个虚拟钻孔；
[0031](3
‑
3)顺序遍历钻孔集合PD，根据下式通过线性插值，计算待插入虚拟钻孔pd...

【技术保护点】

【技术特征摘要】
1.一种地层分布模式约束的钻孔对比剖面模型构建方法，其特征在于该方法具体包括：(1)读取剖面线上实测钻孔的钻孔集合D和钻孔地层厚度集合T；(2)从钻孔集合D中读取任一相邻实测钻孔组合，基于钻孔地层厚度集合T建立相邻实测钻孔组合间的缺失提取地层的缺失标志列表SL，并根据SL构建地层缺失区间集合DS，再基于DS解析地层分布模式，形成地层分布模式标志列表M；(3)插值生成所述相邻实测钻孔组合之间的虚拟钻孔，并推断各虚拟钻孔的底面标高和地面标高，将所述相邻实测钻孔组合和之间的虚拟钻孔存入钻孔集合PD；(4)基于地层尖灭规律，根据地层分布模式标志列表M初步推断钻孔集合PD中各钻孔的各地层厚度，存入地层厚度矩阵PT；(5)根据各虚拟钻孔的底面标高和地面标高优化地层厚度，并构建所述相邻实测钻孔组合间所有地层面；(6)根据构建的地层面和钻孔集合PD生成所述相邻实测钻孔组合间的对比剖面图；(7)循环执行步骤(2)
‑
(6)，直至钻孔集合D中所有相邻实测钻孔组合被遍历，得到所有对比剖面图。2.根据权利要求1所述的地层分布模式约束的钻孔对比剖面模型构建方法，其特征在于：步骤(1)包括；(1
‑
1)读取剖面线上实测钻孔的钻孔信息，得到钻孔集合D＝{d
i
|i＝1，2，...，dn}；其中，d
i
表示第i个钻孔，dn表示钻孔数量；(1
‑
2)读取剖面线上实测钻孔的钻孔地层信息，得到地层厚度集合T＝{t
i，j
|i＝1，2，...，dn；j＝1，2，...，sni}；其中，t
i，j
表示钻孔d
i
的第j个地层的厚度，sni表示钻孔d
i
钻遇的地层总数。3.根据权利要求1所述的地层分布模式约束的钻孔对比剖面模型构建方法，其特征在于：步骤(2)包括；(2
‑
1)从钻孔集合D中读取任意一组相邻实测钻孔d
i
和d
i+1
；(2
‑
2)遍历钻孔地层厚度集合T，根据下式匹配钻孔d
i
和钻孔d
i+1
间的缺失提取地层的缺失标志列表SL＝{sl
j
|j＝1，2，...，sni}；其中，sl
j
表示地层j的地层缺失标志，t
i，j
、t
i+1，j
表示钻孔d
i
、d
i+1
中地层j的厚度，sni表示钻孔d
i
钻遇的地层总数；(2
‑
3)基于地层缺失标志列表SL
i
，构建钻孔d
i
和钻孔d
i+1
间地层缺失区间集合DS＝{ds
u
|u＝1，2，...，dsn}；其中，ds
u
＝(start，end，count)为一个三元组，表示第u个地层缺失区间，dsn为缺失区间数量，start为区间ds
u
的起始地层编号，end为区间ds
u
的终止地层编号，count为区间ds
u
中的地层个数；(2
‑
4)对于地层缺失区间集合DS中每一地层缺失区间ds
u
，遍历区间ds
u
内的所有地层的
缺失情况，并按照下式计算区间ds
u
内的地层分布模式标志m
u
，存入地层分布模式标志列表M＝{m
u
|u＝1，2，...，dsn}；4.根据权利要求3所述的地层分布模式约束的钻孔对比剖面模型构建方法，其特征在于：步骤(2
‑
3)包括；(2
‑3‑
1)顺序遍历缺失标志列表SL，直至sl
j
不等于0，将j的值赋予start，并将count值赋为1；(2
‑3‑
2)继续遍历缺失标志列表SL，当j小于等于sni且|sl
j
|等于|sl
start
|时，令count值递增1；否则，将end值赋为j
‑
1，新建区间元素ds
u
＝(start，end，count)，并添加至地层缺失区间集合DS，之后返回至步骤(2
‑3‑
2)；其中，||为绝对值符号；(2
‑3‑
3)循环执行步骤(2
‑3‑
1)至(2
‑3‑
2)，直至缺失标志列表SL被遍历完，得到地层缺失区间集合DS。5.根据权利要求1所述的地层分布模式约束的钻孔对比剖面模型构建方法，其特征在于：步骤(3)包括；(3
‑
1)遍历集合DS，获取DS中所有ds
u
(count)的最大值，记为最大区间缺失地层数max，作为待插入的虚拟钻孔数量，其中，ds
u
(count)为区间ds
u
中的地层个数；(3
‑
2)新建空的钻孔集合PD＝{pd
k
|k＝1，2，...，max+2}，存储虚拟钻孔及读取的相邻实测钻孔d
i
和d
i+1
；其中，pd1值赋为实测钻孔d
i
，pd
max+2
值赋为实测钻孔d
i+1
，pd2至pd
max+1
依次为第1至max个虚拟钻孔；(3
‑
3)顺序遍历钻孔集合PD，根据下式通过线性插值，计算待插入虚拟钻孔pd
k
的坐标和底面标高pd
k
(bottom)；其中，pd
k
(x)和pd
k
(y)为虚拟钻孔pd
k
的横坐标和纵坐标，d
i
(x)和d
i
(y)为钻孔d
i
的横坐标和纵坐标，d
i+1
(x)和d
i+1
(y)为钻孔d
i+1
的横坐标和纵坐标，d
i
(bottom)、d
i+1
(bottom)和pd
k
(bottom)分别表示钻孔d
i
、d
i+1
和pd
k
的孔底标高；(3
‑
4)根据虚拟钻孔pd
k
坐标从DEM数据中提取虚拟钻孔pd
k
的地面标高pd
k
(h)，k＝2，...，max+1，完成钻孔集合PD中的各钻孔赋值...

【专利技术属性】
技术研发人员：李安波，徐诗宇，黄健初，陈浩，沈言根，
申请(专利权)人：南京师范大学，
类型：发明
国别省市：