后将提取的曲面组成的模型作为层面模型。图12为在正断层的 情况下，离散隐函数场的某一剖面效果示意图。图13为在逆断层的情况下，离散隐函数场的 某一剖面效果示意图。图14为在正断层的情况下，离散隐函数场内第一网格的模型示意图。 图15为在逆断层的情况下，离散隐函数场内第一网格的模型示意图。其中，所述隐函数一般 为泊松方程中的指示函数。
[0088] 在本实施方式中，为了快速求解泊松方程，可W根据有限差分算法，将所述泊松方 程离散为线性方程组，然后采用基于多重网格分裂的共辆梯度法，求解所述线性方程组，得 到指示函数。
[0089] 在本实施方式中，可W依次获取每个指示函数的值为零的第一网格。对于获取的 第一网格，可W基于该第一网格的边的连通状态，判断该第一网格是否为完整网格单元。具 体地，当该第一网格的边的连通状态为连通时，可W判断该第一网格为完整网格单元，当该 第一网格的边的连通状态为不连通时，可W判断该第一网格为非完整网格单元。
[0090] 在本实施方式中，如果第一网格是完整网格单元，那么可W基于经典的移动立方 体(March化be,MC)算法，提取该第一网格单元的曲面片。
[0091] 在本实施方式中，如果第一网格是非完整网格单元，那么，可W基于该第一网格每 个边的连通状态将该第一网格划分为多个部分。每一部分可W记录有该部分的有效节点列 表，那么，对于其中的每一部分，可W按照有效节点列表和该部分网格边的连通状态，从已 处理的完整网格单元或非完整单元中捜索邻近控制点集，并基于Ξ维局部克里金插值算法 外推估算该部分的非有效节点集，然后在估算的基础上生成一个伪有效单元，并基于经典 的移动立方体算法从该伪有效单元中提取指示函数的值为零的Ξ角面，然后利用穿过该第 一网格的断面裁剪提取的Ξ角面。
[0092] 在本实施方式中，可W捜集所有完整网格单元的曲面片，W及非完整网格单元的 裁剪的并提取的Ξ角面，得到层面模型。图16为非完整单元的有效节点和非有效节点划分 示意图。在图16中，浅灰色的面为穿过网格的断面，即与网格相交的断面，深灰色的面为该 断面与网格相交所形成的面。图17为完整单元的示意图。图18为在正断层的情况下，得到的 层面模型示意图。图19为在逆断层的情况下，得到的层面模型示意图。
[0093] 图1所对应的实施例，通过泊松方程，可W将地震数据转换为指示函数标量场，通 过所述指示函数标量场可W快速地得到层面的内外区域。因此，与现有技术相比，图1所对 应的实施例可W在存在大量正断层和/或逆断层的情况下，高效地建立层面模型，从而可W 提高层面模型建立过程的适应能力W及建立效率。
[0094] 进一步地，图1所对应的实施例，在Ξ维空间中使用泊松方程，将地震数据中的层 位数据点转换为指示函数标量场，通过抽取指示函数标量场的等值面数据生成层面模型。 因此，图1所对应的实施例可W提高层面模型建立过程的适应能力W及及建立效率。
[00Μ]进一步地，图1所对应的实施例，与现有技术相比，基于泊松方程建立层面模型，因 此对输入数据具有较强的适应性，可W采用统一的方式处理各种复杂的地质界面。例如可 W统一处理含有大量非通透复杂交切情况下的正逆断裂层面、盐丘层面等地质界面。
[0096] 更进一步地，图1所对应的实施例，具有良好的外推能力，W及在大规模断面存在 的情况下，快速建立层面模型的能力，可W在有数据点的地方忠实于实际数据信息，在没有 数据点的地方具有合理的外推效果，从而可W极大地提高层面模型的建立效率并实现边解 释边建模工作模式。
[0097] 虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和 变化而不脱离本申请的精神，希望所附的权利要求包括运些变形和变化而不脱离本申请的 精神。
【主权项】
1. 一种层面模型建立方法，其特征在于，包括： 根据地震数据建立断面模型； 将预设建模区域进行网格划分，将划分的网格作为第一网格； 基于所述断面模型，设置每个第一网格边的连通状态； 基于所述地震数据，以及预设的分段函数，建立关于指示函数的泊松方程，其中，所述 分段函数用于表示层面的变化趋势，所述指示函数为标量场，用于表示层面的内外区域； 基于所述泊松方程和设置网格边连通状态后的第一网格，建立层面模型。2. 如权利要求1所述的方法，其特征在于，所述泊松方程为：其中，函数Φ为指示函数， 汉为根据所述分段函数生成的梯度向量场， ▽ . P为梯度向量场P的散度， Δ Φ为指示函数Φ梯度的散度。3. 如权利要求1所述的方法，其特征在于，所述根据地震数据建立断面模型，具体包括： 从所述地震数据中获取断层的数据点集； 根据每个断层的数据点集，建立该断层的局部坐标系； 在每个断层的局部坐标系中，根据该断层的数据点集，拟合该断层的断面，并获取该断 面的边界，其中，所述断面采用第二网格表示； 基于断面的第二网格和边界，以及断面之间的接触关系，得到断面模型。4. 如权利要求1所述的方法，其特征在于，当所述第一网格为长方体网格时，所述基于 所述断面模型，设置每个第一网格边的连通状态，具体包括： 将每个长方体网格边的连通状态设置为连通； 基于连接于同一个顶点的并属于同一个长方体网格的网格边，生成第一方向、第二方 向和第三方向的网格线，获取所述断面模型分别与第一方向、第二方向和第三方向的网格 线的交点，基于每个交点，修改每个长方体网格边的连通状态。5. 如权利要求1所述的方法，其特征在于，所述基于所述地震数据，以及预设的分段函 数，建立关于指示函数的泊松方程，具体包括： 从所述地震数据中获取层位数据点，并计算每个层位数据点的法向量； 根据预设的分段函数，以及每个层位数据点的法向量，建立关于指示函数的泊松方程。6. 如权利要求1所述的方法，其特征在于，所述基于所述泊松方程和设置网格边连通状 态后的第一网格，建立层面模型，具体包括： 求解所述泊松方程，得到指示函数； 基于所述指示函数，以及设置网格边连通状态后的第一网格，建立层面模型。7. 如权利要求6所述的方法，其特征在于，所述基于所述指示函数，以及设置网格边连 通状态后的第一网格，建立层面模型，具体包括： 获取指示函数的值为零的第一网格，根据获取的第一网格的边的连通状态，判断获取 的第一网格是否为完整单元； 根据判断结果，基于预设规则从获取的第一网格中提取曲面； 基于提取的曲面，建立层面模型。8. 如权利要求7所述的方法，其特征在于，所述预设规则包括移动立方体算法。9. 如权利要求1所述的方法，其特征在于，所述断面之间的接触关系包括不接触、主主 接触和主辅接触。10. 如权利要求1所述的方法，其特征在于，所述分段函数包括：其中，w为预设的与边界之间的距离； d为分段函数於的自变量，表示层位数据点距离层面边界的距离。
【专利摘要】本申请实施例公开了一种层面模型建立方法。所述方法包括：根据地震数据建立断面模型；将预设建模区域进行网格划分，将划分的网格作为第一网格；基于所述断面模型，设置每个第一网格边的连通状态；基于所述地震数据，以及预设的分段函数，建立关于指示函数的泊松方程，其中，所述分段函数用于表示层面的变化趋势，所述指示函数为标量场，用于表示层面的内外区域；基于所述泊松方程和设置网格边连通状态后的第一网格，建立层面模型。本申请实施例的方法，可以适用于存在正断层和/或逆断层的情况，并且可以提高层面模型的建立效率。
【IPC分类】G01V1/30
【公开号】CN105487116
【申请号】CN201510711499
【发明人】于海生, 白雪莲, 钱宇明, 赵亮, 楚万长, 陈海云, 滕振宇, 李咏梅, 王立松, 鱼福灵
【申请人】中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司
【公开日】2016年4月13日
【申请日】2015年10月28日