(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210118359.3 (22)申请日 2022.02.08 (71)申请人 南京师范大学 地址 210024 江苏省南京市 鼓楼区宁海路 122号 (72)发明人 毛一鹏　刘先裕　李安波　 (74)专利代理 机构 南京苏高专利商标事务所 (普通合伙) 32204 专利代理师 冯艳芬 (51)Int.Cl. G06V 20/10(2022.01) G06V 10/44(2022.01) G06V 10/774(2022.01) G06K 9/62(2022.01) (54)发明名称 一种穹窿构造的自动识别方法及装置 (57)摘要 本发明公开了一种穹窿构造的自动识别方 法及装置， 方法包括如下步骤： (1)从待识别的地 层数据中提取环数≤2的基岩地层， 存入地层集 合D中； (2)依次获取地层集合D中每个地层的地 层边界环数、 包络矩形及其面积， 并存入集合 Index中； (3)从集合Index中识别出具环状结构 的地层， 并存入地层集合Dtemp； (4)根据地层的 长度面积比特征， 删除集合Dtemp中不符合穹窿 构造特征的地层； (5)根据地层的长短轴比特征 和面积比特征， 删除集合Dtemp中不符合穹窿构 造特征的地层； (6)根据地层的产状特征， 将集合 Dtemp中符合穹窿构造特征的地层存入穹窿集合 Domes。 本发明识别结果较为准确， 自动化程度 高。 权利要求书2页 说明书6页 附图9页 CN 114511782 A 2022.05.17 CN 114511782 A 1.一种穹窿构造的自动识别方法， 包括如下步骤： (1)从待识别的地层数据中提取环数≤2的基岩地层， 存 入地层集 合D中； (2)依次获取地层集合D中每个地层的地层边界环数、 包络矩形及其面积， 并存入集合 Index中， 再根据包络矩形的面积从大到小对I ndex进行排序； (3)从集合Index中识别出 具环状结构的地层， 并存 入地层集 合Dtemp； (4)根据地层的长度面积比特 征， 删除集 合Dtemp中不符合穹窿构造特 征的地层； (5)根据地层的长短轴比特征和面积比特征， 删除集合Dtemp中不符合穹窿构造特征的 地层； (6)根据地层的产状特征， 将集合Dtemp中符合穹窿构造特征的地层存入穹窿集合 Domes， 其 余地层存入水平构造集 合Horizo ntalStructure。 2.根据权利要求1所述的穹窿构造的自动识别方法， 其特 征在于： 步骤(1)包括： (1‑1)将待识别的地层数据读入地层集合D中， 集合中每个地层记录的属性包括ID、 年 代属性、 地层边界环数和地层类别； (1‑2)将地层集 合D中年代属性 为松散层的地层删除； (1‑3)将地层集 合D中地层类别为侵入岩的地层删除； (1‑4)将地层集合D中地层边界环数>2的地层删 除， 最终得到地层集合D＝{di|i＝1, 2,…,SN}， di为D中第i个地层， SN 为地层数目。 3.根据权利要求1所述的穹窿构造的自动识别方法， 其特 征在于： 步骤(2)包括： (2‑1)对于集合D中每一地层di， 获取该地层的序号IDi、 地层边界环数Numi， 并提取每个 地层的包络矩形及其面积 Areai， 存入集合Index＝{si(IDi,Numi,xmaxi,xmini,ymaxi,ymini, Areai)|i＝1,2, …,SN}， 其中， si为集合Index中第i个地层， {xmaxi,xmini,ymaxi,ymini}表 示si的包络矩形的横坐标最大值、 横坐标最小值、 纵坐标最大值、 纵坐标最小值， SN为集合D 中的地层数量； (2‑2)按照包络矩形面积从大到小的顺序， 对集 合Index重新 排序。 4.根据权利要求1所述的穹窿构造的自动识别方法， 其特 征在于： 步骤(3)包括： (3‑1)从集合Index中读取任一未读取 过的地层si； (3‑2)获取si的环数Numi， 若环数为1， 则返回执 行步骤(3 ‑1)； 否则， 执 行步骤(3 ‑3)； (3‑3)从集合Index中读取另一未读取过的地层sj， 且j>i， 根据公式下式， 计算si和sj的 MBR包含关系R(si,sj): (3‑4)若R(si,sj)＝0， 则返回执行步骤(3 ‑3)； 若R(si,sj)＝1， 判断Si的内环和Sj的外环 面积是否相等， 若不相等， 则返回执行步骤(3 ‑3),若相等,则 将si的序号IDi,sj的序号IDj存 入地层集 合Dtemp； (3‑5)返回执行步骤(3 ‑1)， 直至集合Index中所有地层被遍历， 得到地层集合Dtemp＝ {tempk|k＝1,2,…,TN}， tempk为Dtemp中第k个环状地层， TN 为环状地层数量。权　利　要　求　书 1/2 页 2 CN 114511782 A 25.根据权利要求1所述的穹窿构造的自动识别方法， 其特 征在于： 步骤(4)包括： (4‑1)将集合Dtemp中具有邻接关系的地层划分为 一组， 得到若干组地层组合； (4‑2)读取任一组具邻接关系的地层组合； (4‑3)根据下式， 计算该组地层中每一地层的长度面积比k1ID与等面积圆 的长度面积比 k0ID的倍数kID； 式中， LengthID、 SizeID分别为对应地层的外周长和补全面积； (4‑4)根据该组地层中所有地层的kID值， 计算该组地层的kID值的方差值； (4‑5)若该组地层的kID值的方差值大于预设阈值K， 则判定该组地层不属于穹窿构造， 将该组地层从Dtemp中删除； (4‑6)重复执 行步骤(4 ‑2)到(4‑5)， 直至所有地层组合被遍历。 6.根据权利要求1所述的穹窿构造的自动识别方法， 其特 征在于： 步骤(5)包括： (5‑1)将集合Dtemp中具有邻接关系的地层划分为 一组， 得到若干组地层组合； (5‑2)读取任一组具邻接关系的地层组合； (5‑3)计算该组地层中面积最大的地层的长短轴比例xID， 若其值大于预设阈值X， 则判 定该组地层不属于穹窿构造， 将其从集合Dtemp中删除， 并执行步骤(5 ‑5)， 否则执行步骤 (5‑4)； (5‑4)计算该组地层中任意两个相邻地层的补全面积比pID， 若其值大于预设阈值P， 则 判定该组地层不属于穹窿构造， 将其从集 合Dtemp中删除； (5‑5)重复执 行步骤(5 ‑2)‑(5‑4)， 直至所有地层组合被遍历。 7.根据权利要求1所述的穹窿构造的自动识别方法， 其特 征在于： 步骤(6)包括： (6‑1)将集合Dtemp中具有邻接关系的地层划分为 一组， 得到若干组地层组合； (6‑2)读取任一组具邻接关系的地层组合； (6‑3)对于该组地层中的每一个地层， 提取该地层范围内包含的所有产状点的倾角信 息， 并计算该地层的平均倾角avgDipAngle； (6‑4)若该组地层中所有地层的平均倾角avgDipAngle均小于预设阈值a， 则将该组地 层存入Domes中； 否则， 存 入水平构造集 合Horizo ntalStructure中； (6‑5)重复操作(6 ‑2)‑(6‑4)， 直至集合Dtemp中所有地层组合被遍历， 得到穹窿集合 Domes＝{domel|l＝1,2, …,DN}和水平构造集合中HorizontalStructure＝{hsv|v＝1, 2,…,DH}， domel表示第l个穹窿构造， DN为穹窿构造数量， hsv表示第v个水平构造， DH 为水平 构造数量。 8.一种穹窿构造的自动识别装置， 包括处理器及存储在存储器上并可在处理器上运行 的计算机程序， 其特征在于， 所述处理器执行所述程序时实现权利要求 1‑7中任意一项 所述 的方法。权　利　要　求　书 2/2 页 3 CN 114511782 A 3