(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111440519.8 (22)申请日 2021.11.30 (71)申请人 中海石油 （中国） 有限公司 地址 100010 北京市东城区朝阳门北 大街 25号 申请人 中海石油 （中国） 有限公司北京研究 中心 (72)发明人 李中　何英明　周长所　范白涛　 幸雪松　谢仁军　付兴　张奎　 黄宁曼　孙进　 (74)专利代理 机构 北京纪凯知识产权代理有限 公司 11245 代理人 谢斌 (51)Int.Cl. G06Q 10/04(2012.01)G06Q 10/06(2012.01) (54)发明名称 一种基于非线性梯度结构张量算法的火山 通道检测方法 (57)摘要 本发明公开了一种基于非线性梯度结构张 量算法的火山通道检测方法， 首先对地震数据提 取张量积， 在各向异性扩散滤波函数的约束下构 建系数扩散矩 阵对调整后的张量积进行正则化 后， 再计算正则化后张量积的特征值和特征向 量， 通过特征向量计算得到非线性梯度结构张量 属性， 并利用已钻井的漏失数据对非线性梯度结 构张量属性进行标定， 确定属性门槛值， 实现火 山通道的定量刻画， 最后利用该定量描述后的火 山通道去指导其他设计井位中的火山通道预测。 本发明公开的基于非线性梯度结构张量算法的 火山通道检测方法， 对后续井位 设计和调整具有 非常重要的指导 意义。 权利要求书1页 说明书8页 附图5页 CN 114118589 A 2022.03.01 CN 114118589 A 1.一种基于非线性梯度结构张量 算法的火山 通道检测方法， 其特 征在于， 包括： 先通过线性梯度 结构张量求取地震数据的梯度， 然后利用三维标准高斯核函数对地震 数据进行 滤波构造线性梯度结构张量， 得到地震数据的张量积； 先利用三维标准高斯核函数对地震数据进行滤波， 然后通过非线性梯度 结构张量求取 地震数据的梯度， 得到基于梯度扩散的非线性梯度结构张量； 调整所述基于梯度扩散的非线性梯度结构张量中的滤波函数， 进行正则化计算， 得到 正则化后的张量积； 计算所述正则化后的张量积中的特征值和特征向量， 其中不同的特征值和特征向量分 别对应不同的张量属性。 2.根据权利要求1所述的火山 通道检测方法， 其特 征在于， 在进行所述正则化计算中， 包括以下步骤： 步骤C1： 利用各向异性扩散滤波函数调整非线性梯度结构张量； 步骤C2： 根据非线性梯度结构张量调整所述正则化后的张量积中的特征值和特征向 量。 3.一种基于非线性梯度结构张量 算法的火山 通道检测装置， 其特 征在于， 包括 第一处理单元， 用于先通过线性梯度结构张量求取地震数据的梯度， 然后利用三维标 准高斯核函数对地震数据进行 滤波构造线性梯度结构张量， 得到地震数据的张量积； 第二处理单元， 用于先利用三维标准高斯核函数对地震数据进行滤波， 然后通过非线 性梯度结构张量 求取地震数据的梯度， 得到基于梯度扩散的非线性梯度结构张量； 第三处理单元， 用于调整所述基于梯度扩散的非线性梯度结构张量中的滤波函数， 进 行正则化计算， 得到正则化后的张量积； 第四处理单元， 用于计算所述正则化后的张量积中的特征值和特征向量， 其中不同的 特征值和特 征向量分别对应不同的张量属性。 4.一种计算机可读存储介质， 其上存储有计算机程序， 其特征在于， 所述计算机程序被 处理器执行时实现权利要求1 ‑2任一项所述的方法的步骤。 5.一种计算机设备， 包括存储器、 处理器以及存储在存储器上并可在处理器上运行的 计算机程序， 其特征在于， 所述处理器执行所述计算机程序时实现权利要求 1‑2任一所述的 方法的步骤。权　利　要　求　书 1/1 页 2 CN 114118589 A 2一种基于非线性梯度结构张量算法的火山通道检测方 法 技术领域 [0001]本发明涉及火山通道检测技术领域， 具体涉及一种基于非线性梯度结构张量算法 的火山通道检测方法。 背景技术 [0002]根据火山的结构可以分为火山椎体、 岩浆和 火山通道， 火山通道作为岩浆喷发或 者溢流过程中穿过地下岩层的主要通道， 根据形态主要分为中心式和裂隙式。 火山通道在 发育过程中， 通道内存在大量的未经搬运磨圆的火山角砾岩， 由岩浆 或火山灰胶结而成， 再 加上喷发过程中冲破岩层， 形成一些破碎带， 发育着大量的多尺度裂缝。 因此， 在钻井过程 中钻遇的火山通道由于火山灰胶结物和裂缝 的存在， 使得火山通道成为易漏失型地层。 为 了预防钻井过程中发生大规模漏失， 需要对火山通道发育区进行平面和空间的预测， 以此 为依据在钻井设计过程中尽可能避开火山 通道， 从而预防井漏。 发明内容 [0003]本发明的目的在于提供一种基于非线性梯度结构张量算法的火山通道检测方法， 通过对地震数据计算非线性梯度结构张量属 性， 并结合已钻井的漏失情况， 标定张量属 性 门槛值， 最终通过梯度结构张量属性确定火山 通道的平面和空间发育位置 。 [0004]本发明提供一种基于非线性梯度结构张量算法的火山通道检测方法， 包括以下步 骤： [0005]先通过线性梯度结构张量求取地震数据的梯度， 然后利用三维标准高斯核函数对 地震数据进行 滤波构造线性梯度结构张量， 得到地震数据的张量积； [0006]先利用三维标准高斯核函数对地震数据进行滤波， 然后通过非线性梯度结构张量 求取地震数据的梯度， 得到基于梯度扩散的非线性梯度结构张量； [0007]调整所述基于梯度扩散 的非线性梯度结构张量中的滤波函数， 进行正则化计算， 得到正则化后的张量积； [0008]计算所述正则化后的张量积中的特征值和特征向量， 其中不同的特征值和特征向 量分别对应不同的张量属性。 [0009]进一步地， 在进行 所述正则化计算中， 包括以下步骤： [0010]步骤C1： 利用各向异性扩散滤波函数调整非线性梯度结构张量； [0011]步骤C2： 根据非线性梯度结构张量调整所述正则化后的张量积中的特征值和特征 向量。 [0012]本发明还提供一种基于非线性梯度结构张量 算法的火山 通道检测装置， 包括 [0013]第一处理单元， 用于先通过线性梯度结构张量求取地震数据的梯度， 然后利用三 维标准高斯核函数对地震数据进行 滤波构造线性梯度结构张量， 得到地震数据的张量积； [0014]第二处理单元， 用于先利用三维标准高斯核函数对地震数据进行滤波， 然后通过 非线性梯度结构张量 求取地震数据的梯度， 得到基于梯度扩散的非线性梯度结构张量；说　明　书 1/8 页 3 CN 114118589 A 3