层面模型建立方法
层面模型建立方法
【技术领域】
[0001] 本申请设及地球物理勘探地震数据处理技术领域，特别设及一种层面模型建立方 法。
【背景技术】
[0002] 层面，又称地层面，一般为地层与地层的分界面。层面模型一般是对地层层面的数 学描述，是研究地质参数空间分布的重要基础，也是研究Ξ维地层模型的重要方法。层面模 型的建立为地质构造图和储层参数分布的绘制和显示提供了方法基础。
[0003] 现有技术中，建立层面模型的方法主要包括基于断面分块的方法和基于连续曲面 的方法。
[0004][0005]
【发明内容】
[0006] 本申请实施例的目的是提供一种层面模型建立方法。所述方法可W适用于存在正 断层和/或逆断层的情况，并且可W提高层面模型的建立效率。
[0007] 为解决上述技术问题，本申请实施例提供一种层面模型建立方法是运样实现的： [000引一种层面模型建立方法，包括：
[0009] 根据地震数据建立断面模型；
[0010] 将预设建模区域进行网格划分，将划分的网格作为第一网格；
[0011] 基于所述断面模型，设置每个第一网格边的连通状态；
[0012] 基于所述地震数据，W及预设的分段函数，建立关于指示函数的泊松方程，其中， 所述分段函数用于表示层面的变化趋势，所述指示函数为标量场，用于表示层面的内外区 域；
[0013] 基于所述泊松方程和设置网格边连通状态后的第一网格，建立层面模型。
[0014] 由w上本申请实施例提供的技术方案可见，本申请实施例通过泊松方程，可w将 地震数据转换为指示函数标量场，通过所述指示函数标量场可W快速地得到层面的内外区 域。因此，与现有技术相比，本申请实施例的方法可W在存在大量正断层和/或逆断层的情 况下，高效地建立层面模型，从而可W提高层面模型建立过程的适应能力W及建立效率。
【附图说明】
[0015] 为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提 下，还可W根据运些附图获得其他的附图。
[0016] 图1为本申请实施例层面模型建立方法的流程图；
[0017] 图2为本申请实施例断层局部坐标系示意图；
[0018] 图3为本申请实施例断面接触关系为主辅关系的示意图；
[0019] 图4为本申请实施例断面接触关系为主主关系的示意图；
[0020] 图5为本申请实施例在正断层的情况下，断面模型的示意图；
[0021] 图6为本申请实施例在逆断层的情况下，断面模型的示意图；
[0022] 图7为本申请实施例在对建模区域进行网格划分后，标记穿过图5所示的断面模型 的网格边的示意图；
[0023] 图8为本申请实施例在对建模区域进行网格划分后，标记穿过图6所示的断面模型 的网格边的示意图；
[0024] 图9为本申请实施例图5所示的断面模型和层位数据点的分布示意图；
[0025] 图10为本申请实施例图6所示的断面模型和层位数据点的分布示意图；
[0026] 图11为本申请实施例层位数据点的最优邻域点集和法向量估计示意图；
[0027] 图12为本申请实施例在正断层的情况下，离散隐函数场的某一剖面效果示意图；
[0028] 图13为本申请实施例在逆断层的情况下，离散隐函数场的某一剖面效果示意图；
[0029] 图14为本申请实施例在正断层的情况下，离散隐函数场中第一网格的模型示意 图；
[0030] 图15为本申请实施例在逆断层的情况下，离散隐函数场中第一网格的模型示意 图；
[0031] 图16为本申请实施例非完整单元的有效节点和非有效节点划分示意图；
[0032] 图17为本申请实施例完整单元的示意图；
[0033] 图18为本申请实施例在正断层的情况下，建立的层面模型示意图；
[0034] 图19为本申请实施例在逆断层的情况下，建立的层面模型示意图。
【具体实施方式】
[0035] 为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实 施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施 例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通 技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护 的范围。
[0036] 在石油勘探开发过程中，常常需要依据地震数据或测井数据建立地下地质模型， W模拟和分析油气分布、计算剩余开采储量、W及对开采设计方案进行数值模拟计算和分 析，从而确定最佳的开发设计方案。因此，地质模型的建立是勘探开发过程中的一项重要工 作。
[0037] 在建立地质模型的过程中，通常需要建立描述地层层面的层面模型。快速准确的 建立层面模型是建立地质模型的基础和关键，层面模型的不确定性对油气剩余资源的分析 影响很大。
[0038] 建立层面模型的方法主要包括基于断面分块的方法和基于连续曲面的方法。上述 两种方法目前均广泛用于地质建模系统中。但面对越来越复杂的构造地质情况，比如存在 大量正逆断层的情况或上千个断面的情况，上述层面模型建立方法常常需要花费较长的时 间，甚至无法完成建立，严重制约着层面模型的建立效率。
[0039] 本申请实施例提供一种基于泊松方程的自动地质层面建模方法，W解决现有技术 中的建模方法无法适应复杂构造情况下的正断层和逆断层问题，W及层面模型建立效率低 下的问题，从而为后续的地层网格模型和属性模型的分析和模拟提供了保证。所述方法可 W包括如下的步骤：
[0040] S101:根据地震数据建立断面模型。
[0041] 具体地，所述地震数据可W是根据地震观测系统采集得到的地震数据。所述断面 又称断层面，是因受挤压或拉伸等作用过度，使地层发生断裂而形成的面。每个断层一般均 具有断面。通常地，断面与层面之间可W相交，断面与断面之间也可W相交。
[0042] -般地，可W通过如下的方式，根据地震数据建立断面模型：
[0043] a:从所述地震数据中获取每个断层的数据点集。
[0044] -般地，地震数据中可W包括多个断层的数据。地震数据一般是通过采集得到的， 因此，每个断层的数据可W包括多个数据点。那么，可W将断层的多个数据点作为该断层的 数据点集。
[0045] b、根据每个断层的数据点集，建立该断层的局部坐标系。
[0046] 具体地，可W根据每个断层的数据点集，基于主分量分析法，估计该断层的最大投 影局部坐标系化ocal Coordinate)，并将该断层的数据点集转换至该局部坐标系。其中，所 述局部坐标系，一般为W物体的中屯、为坐标原点建立的坐标系。物体的旋转或平移等操作 都是围绕该物体的局部坐标系进行的。当物体模型进行旋转或平移等操作时，局部坐标系 也执行相应的旋转或平移操作。所述局部坐标系可W包括局部笛卡尔直角坐标系。图2为断 层的局部坐标系示意图。
[0047] C:在每个断层的局部坐标系中，根据该断层的数据点集，拟合该断层的断面，并获 取该断面的边界，其中，所述断面采用第二网格表示。
[0048] 所述第二网格可W包括矩形网格。所述断面的边界的形状可W为多边形。
[0049] 具体地，在每个断层的局部坐标系中，可W根据该断层的数据点集，基于最小曲率 插值算法，拟合该断层的断面，然后基于A1地a化apes算法，计算该断面的边界。其中，断层 的断面可W由多个第二网格组 成。
[0050] d:基于断面的第二网格和边界，W及断面之间的接触关系，得到断面模型。
[0051] 断面之间的接触关系可W包括无关系、主辅关系和主主关系。其中，无关系表示断 面之间不接触;主辅关系表示主断面裁剪辅断面;主主关系表示断面之间存在交切，但求交 后断面之间互相不裁剪。图3为断面接触关系为主辅关系的示意图，图4为断面接触关系为 主主关系的示意图。
[0052] 具体地，对于每个断层的断面，可W在该断层的局部坐标系中，利用该断面的边界 裁剪该断面的第二网格，得到该断层的初始断面模型。例如，当第二网格为矩形网格时，可 W利用断面的边界裁剪该断面的矩形网格，保留边界内部的矩形网格，并将边界内部的矩 形网格转换为预设的Ξ角形网格，最后将所述Ξ角形网格变换至世界坐标系，从而得到与 该断面对应的断层的初始断面模型。然后根据各个断面之间的接触关系，计算断面之间的 交线，并沿所述交线对所述初始断面模型进行裁剪和网格划分，得到最终的断面模型。图5 为在正断层的情况下，断面模型的示意图。图6为本申请实施例在逆断层的情况下，断面模 型的示意图。
[0053] S102:将预设建模区域进行网格划分，将划分的网格作为第一网格。
[0054] 所述第一网格可W为Ξ维计算网格，具体可W包括四面体网格和长方体网格。
[0055] 具体地，可W根据预设层面网格的大小，将预设建模区域划分为多个第一网格。所 述预设建模区域可W为预设的建立层面模型的区域。
[0056] S103:基于所述断面模型，设置每个第一网格边的连通状态。
[0057] 在划分第一网格后，为了处理断面引起的层面不连续性，可W基于所述断面模型， 设置每个第一网格的边的连通状态。下面W长方体网格为例，详细说明设置每个第一网格 边连通状态的过程。
[0058] f:将每个长方体网格边的连通状态设置为连通。
[0059] 具体地，对于每个长方体网格顶点，可W将与该顶点相连的网格边的连通状态设 置为连通。一般地，在所述建模区域内，除边界顶点外，每个长方体网格顶点可W有六个网 格边与该顶点相连接，那么，可W将与该顶点相连的六个网格边的连通状态设置为连通。例 如可W将与该顶点相连的六个网格边的连通状态标记为0。其中，0表示网格边连通。
[0060] g:基于连接于同一个顶点的并属于同一个长方体网格的网格边，生成第一方向、 第二方向和第Ξ方向的网格线，获取所述断面模型分别与第一方向、第二方向和第Ξ方向 的网格线的交点，基于每个交点，修改每个长方体网格边的连通状态。
[0061] 具体地，在所述建模区域内，除边界顶点外，与同一个长方体网格顶点相连接的并 且属于同一个长方体网格的网格边有Ξ个，该Ξ个网格边相互垂直，那么，可W将该Ξ个网 格边的方向分别作为第一方向、第二方向和第Ξ方向。然后根据所述第一方向、第二方向和 第Ξ方向，分别生成第一方向、第二方向和第Ξ方向的网格线。然后可W将所述断面模型分 别与第一方向、第二方向和第Ξ方向的网格线求交，得到断面模型与每个网格线的交点，然 后利用交点判断每个长方体网格边是否与断面模型的断面相交，如果相交则修改该网格边 的连通状态为不连通。例如，可W将该网格边的连通状态标记修改为1。例如，在图2中，X方 向为第一方向，y方向为第二方向，Z方向为第Ξ方向。
[0062] 图7为在对建模区域进行网格划分后，标记穿过图5所示的断面模型的网格边的示 意图。图8为在对建模区域进行网格划分后，标记穿过图6所示的断面模型的网格边的示意 图。在图7和图8中，若网格边穿过断面模型，则修改该网格边的连通状态为不连通。图7和图 8中的标记过程可W在求解泊松方程时保证断面两侧的网格计算互不影响，从而巧妙的解 决了隐函数场的不连续计算问题。
[0063] S104:基于所述地震数据，W及预设的分段函数，建立关于指示函数的泊松方程， 其中，所述分段函数用于表示层面的变化趋势，所述指示函数为标量场，用于表示层面的内 外区域。
[0064] 泊松方程是一种偏微分方程，主要用于连续闭合曲面的重构。泊松方程的基本原 理为：已知计算定义域A，W及A上的向量场京，求出函数Φ使其梯度最佳逼近向量场K，即 满足算式
。依据变分原理等相关定理，上述求解过程等价于求解泊松方程 这分=矿护。其中，
[0065] V表示梯度向量场r的散度，具体为化尿；
[0066] Δ Φ =div grad Φ，其中，grad Φ表示标量Φ的梯度，div grad Φ表示标量Φ梯 度的散度。
[0067] 泊松方程可W用于工业产品设计中逆向工程的闭合曲面重构，并能够取得很好的 效果。该闭合曲面重构通常情况下是一种隐式闭合曲面重构方法。闭合曲面重构的基本思 路为:依据给定测量散点的空间位置和法矢信息，利用分段函数定义闭合曲面的内外部分， 并基于分段函数建立泊松方程，然后通过求解泊松方程提取指示函数值为零的等值面。
[0068] 一般地，表示层面变化趋势的函数矿在层面边界处是不连续的，因此，为了数值计 算的稳定性，可^在距边界的距离*内设置为线性过渡形式。分段函数^^的形式可^如下：
[0069]
[0070] 其中，W为预设的与边界之间的距离；
[0071] d为分段函数的自变量，表示层位数据点距离层面边界的距离。
[0072] 具体地，可W从所述地震数据中获取层位数据点，并计算每个层位数据点的法向 量，然后根据所述分段函数W及每个层位数据点的法向量，建立关于指示函数的泊松方程。 图9为图5所示的断面模型和层位数据点的分布示意图。图10为图6所示的断面模型和层位 数据点的分布示意图。
[0073] 在本实施方式中，可W基于自然邻域捜索算法建立每个层位数据点的最优邻域点 集，然后基于最优邻域点集进行最小二乘法平面拟合，得到拟合平面，然后获取该拟合平面 的法向量，并将该拟合平面的法向量的方向调整为竖直向上的方向，将调整方向后的法向 量作为该层位数据点的法向量。图11为层位数据点的最优邻域点集和法向量估计示意图。 其中，Z方向为拟合平面的法向量方向。
[0074] 在本实施方式中，可W基于所述分段函数，W及每个层位数据点的法向量，建立梯 度向量场戸。具体地，可W对所述分段函数进行平滑滤波卷积，然后根据平滑滤波卷积后的 分段函数，W及每个层位数据点的法向量，通过如下的公式建立梯度向量场r:
[0075]
[0076] 其中，Μ为层面包围的下半部分区域；
[0077] 嫌?为层面；
[007引擇Μ为表示层面上下区域的指示函数；
[0079] F表示平滑滤波算子，包括高斯滤波算子；
[0080] 而,^/,：)表示层面数据点口的法向量。
[0081] 在本实施方式中，根据梯度向量场f，可W建立关于指示函数Φ的泊松方程。所述 泊松方程具体为：A# = V节。其中，
[0082] 函数Φ为指示函数，
[0083] F为根据所述分段函数生成的梯度向量场，
[0084] V .F为梯度向量场F的散度，
[00化]Δ Φ为指示函数Φ梯度的散度。
[0086] S105:基于所述泊松方程和设置网格边连通状态后的第一网格，建立层面模型。
[0087] 具体地，可W求解所述泊松方程，得到指示函数。然后基于所述指示函数，W及设 置网格边连通状态后的第一网格，建立层面模型。所述指示函数一般为标量场，其内部可W 包括多个第一网格。一般地，网格顶点的指示函数的值可W作为该网格的指示函数的值。那 么，可W基于指示函数，获取指示函数的值为零的第一网格，并根据获取的第一网格的边的 连通状态，判断获取的第一网格是否为完整单元，然后根据判断结果，基于预设规则从获取 的第一网格中提取曲面，最 后将提取的曲面组成的模型作为层面模型。图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日
【技术领域】
[0001] 本申请设及地球物理勘探地震数据处理技术领域，特别设及一种层面模型建立方 法。
【背景技术】
[0002] 层面，又称地层面，一般为地层与地层的分界面。层面模型一般是对地层层面的数 学描述，是研究地质参数空间分布的重要基础，也是研究Ξ维地层模型的重要方法。层面模 型的建立为地质构造图和储层参数分布的绘制和显示提供了方法基础。
[0003] 现有技术中，建立层面模型的方法主要包括基于断面分块的方法和基于连续曲面 的方法。
[0004][0005]
【发明内容】
[0006] 本申请实施例的目的是提供一种层面模型建立方法。所述方法可W适用于存在正 断层和/或逆断层的情况，并且可W提高层面模型的建立效率。
[0007] 为解决上述技术问题，本申请实施例提供一种层面模型建立方法是运样实现的： [000引一种层面模型建立方法，包括：
[0009] 根据地震数据建立断面模型；
[0010] 将预设建模区域进行网格划分，将划分的网格作为第一网格；
[0011] 基于所述断面模型，设置每个第一网格边的连通状态；
[0012] 基于所述地震数据，W及预设的分段函数，建立关于指示函数的泊松方程，其中， 所述分段函数用于表示层面的变化趋势，所述指示函数为标量场，用于表示层面的内外区 域；
[0013] 基于所述泊松方程和设置网格边连通状态后的第一网格，建立层面模型。
[0014] 由w上本申请实施例提供的技术方案可见，本申请实施例通过泊松方程，可w将 地震数据转换为指示函数标量场，通过所述指示函数标量场可W快速地得到层面的内外区 域。因此，与现有技术相比，本申请实施例的方法可W在存在大量正断层和/或逆断层的情 况下，高效地建立层面模型，从而可W提高层面模型建立过程的适应能力W及建立效率。
【附图说明】
[0015] 为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提 下，还可W根据运些附图获得其他的附图。
[0016] 图1为本申请实施例层面模型建立方法的流程图；
[0017] 图2为本申请实施例断层局部坐标系示意图；
[0018] 图3为本申请实施例断面接触关系为主辅关系的示意图；
[0019] 图4为本申请实施例断面接触关系为主主关系的示意图；
[0020] 图5为本申请实施例在正断层的情况下，断面模型的示意图；
[0021] 图6为本申请实施例在逆断层的情况下，断面模型的示意图；
[0022] 图7为本申请实施例在对建模区域进行网格划分后，标记穿过图5所示的断面模型 的网格边的示意图；
[0023] 图8为本申请实施例在对建模区域进行网格划分后，标记穿过图6所示的断面模型 的网格边的示意图；
[0024] 图9为本申请实施例图5所示的断面模型和层位数据点的分布示意图；
[0025] 图10为本申请实施例图6所示的断面模型和层位数据点的分布示意图；
[0026] 图11为本申请实施例层位数据点的最优邻域点集和法向量估计示意图；
[0027] 图12为本申请实施例在正断层的情况下，离散隐函数场的某一剖面效果示意图；
[0028] 图13为本申请实施例在逆断层的情况下，离散隐函数场的某一剖面效果示意图；
[0029] 图14为本申请实施例在正断层的情况下，离散隐函数场中第一网格的模型示意 图；
[0030] 图15为本申请实施例在逆断层的情况下，离散隐函数场中第一网格的模型示意 图；
[0031] 图16为本申请实施例非完整单元的有效节点和非有效节点划分示意图；
[0032] 图17为本申请实施例完整单元的示意图；
[0033] 图18为本申请实施例在正断层的情况下，建立的层面模型示意图；
[0034] 图19为本申请实施例在逆断层的情况下，建立的层面模型示意图。
【具体实施方式】
[0035] 为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实 施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施 例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通 技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护 的范围。
[0036] 在石油勘探开发过程中，常常需要依据地震数据或测井数据建立地下地质模型， W模拟和分析油气分布、计算剩余开采储量、W及对开采设计方案进行数值模拟计算和分 析，从而确定最佳的开发设计方案。因此，地质模型的建立是勘探开发过程中的一项重要工 作。
[0037] 在建立地质模型的过程中，通常需要建立描述地层层面的层面模型。快速准确的 建立层面模型是建立地质模型的基础和关键，层面模型的不确定性对油气剩余资源的分析 影响很大。
[0038] 建立层面模型的方法主要包括基于断面分块的方法和基于连续曲面的方法。上述 两种方法目前均广泛用于地质建模系统中。但面对越来越复杂的构造地质情况，比如存在 大量正逆断层的情况或上千个断面的情况，上述层面模型建立方法常常需要花费较长的时 间，甚至无法完成建立，严重制约着层面模型的建立效率。
[0039] 本申请实施例提供一种基于泊松方程的自动地质层面建模方法，W解决现有技术 中的建模方法无法适应复杂构造情况下的正断层和逆断层问题，W及层面模型建立效率低 下的问题，从而为后续的地层网格模型和属性模型的分析和模拟提供了保证。所述方法可 W包括如下的步骤：
[0040] S101:根据地震数据建立断面模型。
[0041] 具体地，所述地震数据可W是根据地震观测系统采集得到的地震数据。所述断面 又称断层面，是因受挤压或拉伸等作用过度，使地层发生断裂而形成的面。每个断层一般均 具有断面。通常地，断面与层面之间可W相交，断面与断面之间也可W相交。
[0042] -般地，可W通过如下的方式，根据地震数据建立断面模型：
[0043] a:从所述地震数据中获取每个断层的数据点集。
[0044] -般地，地震数据中可W包括多个断层的数据。地震数据一般是通过采集得到的， 因此，每个断层的数据可W包括多个数据点。那么，可W将断层的多个数据点作为该断层的 数据点集。
[0045] b、根据每个断层的数据点集，建立该断层的局部坐标系。
[0046] 具体地，可W根据每个断层的数据点集，基于主分量分析法，估计该断层的最大投 影局部坐标系化ocal Coordinate)，并将该断层的数据点集转换至该局部坐标系。其中，所 述局部坐标系，一般为W物体的中屯、为坐标原点建立的坐标系。物体的旋转或平移等操作 都是围绕该物体的局部坐标系进行的。当物体模型进行旋转或平移等操作时，局部坐标系 也执行相应的旋转或平移操作。所述局部坐标系可W包括局部笛卡尔直角坐标系。图2为断 层的局部坐标系示意图。
[0047] C:在每个断层的局部坐标系中，根据该断层的数据点集，拟合该断层的断面，并获 取该断面的边界，其中，所述断面采用第二网格表示。
[0048] 所述第二网格可W包括矩形网格。所述断面的边界的形状可W为多边形。
[0049] 具体地，在每个断层的局部坐标系中，可W根据该断层的数据点集，基于最小曲率 插值算法，拟合该断层的断面，然后基于A1地a化apes算法，计算该断面的边界。其中，断层 的断面可W由多个第二网格组 成。
[0050] d:基于断面的第二网格和边界，W及断面之间的接触关系，得到断面模型。
[0051] 断面之间的接触关系可W包括无关系、主辅关系和主主关系。其中，无关系表示断 面之间不接触;主辅关系表示主断面裁剪辅断面;主主关系表示断面之间存在交切，但求交 后断面之间互相不裁剪。图3为断面接触关系为主辅关系的示意图，图4为断面接触关系为 主主关系的示意图。
[0052] 具体地，对于每个断层的断面，可W在该断层的局部坐标系中，利用该断面的边界 裁剪该断面的第二网格，得到该断层的初始断面模型。例如，当第二网格为矩形网格时，可 W利用断面的边界裁剪该断面的矩形网格，保留边界内部的矩形网格，并将边界内部的矩 形网格转换为预设的Ξ角形网格，最后将所述Ξ角形网格变换至世界坐标系，从而得到与 该断面对应的断层的初始断面模型。然后根据各个断面之间的接触关系，计算断面之间的 交线，并沿所述交线对所述初始断面模型进行裁剪和网格划分，得到最终的断面模型。图5 为在正断层的情况下，断面模型的示意图。图6为本申请实施例在逆断层的情况下，断面模 型的示意图。
[0053] S102:将预设建模区域进行网格划分，将划分的网格作为第一网格。
[0054] 所述第一网格可W为Ξ维计算网格，具体可W包括四面体网格和长方体网格。
[0055] 具体地，可W根据预设层面网格的大小，将预设建模区域划分为多个第一网格。所 述预设建模区域可W为预设的建立层面模型的区域。
[0056] S103:基于所述断面模型，设置每个第一网格边的连通状态。
[0057] 在划分第一网格后，为了处理断面引起的层面不连续性，可W基于所述断面模型， 设置每个第一网格的边的连通状态。下面W长方体网格为例，详细说明设置每个第一网格 边连通状态的过程。
[0058] f:将每个长方体网格边的连通状态设置为连通。
[0059] 具体地，对于每个长方体网格顶点，可W将与该顶点相连的网格边的连通状态设 置为连通。一般地，在所述建模区域内，除边界顶点外，每个长方体网格顶点可W有六个网 格边与该顶点相连接，那么，可W将与该顶点相连的六个网格边的连通状态设置为连通。例 如可W将与该顶点相连的六个网格边的连通状态标记为0。其中，0表示网格边连通。
[0060] g:基于连接于同一个顶点的并属于同一个长方体网格的网格边，生成第一方向、 第二方向和第Ξ方向的网格线，获取所述断面模型分别与第一方向、第二方向和第Ξ方向 的网格线的交点，基于每个交点，修改每个长方体网格边的连通状态。
[0061] 具体地，在所述建模区域内，除边界顶点外，与同一个长方体网格顶点相连接的并 且属于同一个长方体网格的网格边有Ξ个，该Ξ个网格边相互垂直，那么，可W将该Ξ个网 格边的方向分别作为第一方向、第二方向和第Ξ方向。然后根据所述第一方向、第二方向和 第Ξ方向，分别生成第一方向、第二方向和第Ξ方向的网格线。然后可W将所述断面模型分 别与第一方向、第二方向和第Ξ方向的网格线求交，得到断面模型与每个网格线的交点，然 后利用交点判断每个长方体网格边是否与断面模型的断面相交，如果相交则修改该网格边 的连通状态为不连通。例如，可W将该网格边的连通状态标记修改为1。例如，在图2中，X方 向为第一方向，y方向为第二方向，Z方向为第Ξ方向。
[0062] 图7为在对建模区域进行网格划分后，标记穿过图5所示的断面模型的网格边的示 意图。图8为在对建模区域进行网格划分后，标记穿过图6所示的断面模型的网格边的示意 图。在图7和图8中，若网格边穿过断面模型，则修改该网格边的连通状态为不连通。图7和图 8中的标记过程可W在求解泊松方程时保证断面两侧的网格计算互不影响，从而巧妙的解 决了隐函数场的不连续计算问题。
[0063] S104:基于所述地震数据，W及预设的分段函数，建立关于指示函数的泊松方程， 其中，所述分段函数用于表示层面的变化趋势，所述指示函数为标量场，用于表示层面的内 外区域。
[0064] 泊松方程是一种偏微分方程，主要用于连续闭合曲面的重构。泊松方程的基本原 理为：已知计算定义域A，W及A上的向量场京，求出函数Φ使其梯度最佳逼近向量场K，即 满足算式
。依据变分原理等相关定理，上述求解过程等价于求解泊松方程 这分=矿护。其中，
[0065] V表示梯度向量场r的散度，具体为化尿；
[0066] Δ Φ =div grad Φ，其中，grad Φ表示标量Φ的梯度，div grad Φ表示标量Φ梯 度的散度。
[0067] 泊松方程可W用于工业产品设计中逆向工程的闭合曲面重构，并能够取得很好的 效果。该闭合曲面重构通常情况下是一种隐式闭合曲面重构方法。闭合曲面重构的基本思 路为:依据给定测量散点的空间位置和法矢信息，利用分段函数定义闭合曲面的内外部分， 并基于分段函数建立泊松方程，然后通过求解泊松方程提取指示函数值为零的等值面。
[0068] 一般地，表示层面变化趋势的函数矿在层面边界处是不连续的，因此，为了数值计 算的稳定性，可^在距边界的距离*内设置为线性过渡形式。分段函数^^的形式可^如下：
[0069]
[0070] 其中，W为预设的与边界之间的距离；
[0071] d为分段函数的自变量，表示层位数据点距离层面边界的距离。
[0072] 具体地，可W从所述地震数据中获取层位数据点，并计算每个层位数据点的法向 量，然后根据所述分段函数W及每个层位数据点的法向量，建立关于指示函数的泊松方程。 图9为图5所示的断面模型和层位数据点的分布示意图。图10为图6所示的断面模型和层位 数据点的分布示意图。
[0073] 在本实施方式中，可W基于自然邻域捜索算法建立每个层位数据点的最优邻域点 集，然后基于最优邻域点集进行最小二乘法平面拟合，得到拟合平面，然后获取该拟合平面 的法向量，并将该拟合平面的法向量的方向调整为竖直向上的方向，将调整方向后的法向 量作为该层位数据点的法向量。图11为层位数据点的最优邻域点集和法向量估计示意图。 其中，Z方向为拟合平面的法向量方向。
[0074] 在本实施方式中，可W基于所述分段函数，W及每个层位数据点的法向量，建立梯 度向量场戸。具体地，可W对所述分段函数进行平滑滤波卷积，然后根据平滑滤波卷积后的 分段函数，W及每个层位数据点的法向量，通过如下的公式建立梯度向量场r:
[0075]
[0076] 其中，Μ为层面包围的下半部分区域；
[0077] 嫌?为层面；
[007引擇Μ为表示层面上下区域的指示函数；
[0079] F表示平滑滤波算子，包括高斯滤波算子；
[0080] 而,^/,：)表示层面数据点口的法向量。
[0081] 在本实施方式中，根据梯度向量场f，可W建立关于指示函数Φ的泊松方程。所述 泊松方程具体为：A# = V节。其中，
[0082] 函数Φ为指示函数，
[0083] F为根据所述分段函数生成的梯度向量场，
[0084] V .F为梯度向量场F的散度，
[00化]Δ Φ为指示函数Φ梯度的散度。
[0086] S105:基于所述泊松方程和设置网格边连通状态后的第一网格，建立层面模型。
[0087] 具体地，可W求解所述泊松方程，得到指示函数。然后基于所述指示函数，W及设 置网格边连通状态后的第一网格，建立层面模型。所述指示函数一般为标量场，其内部可W 包括多个第一网格。一般地，网格顶点的指示函数的值可W作为该网格的指示函数的值。那 么，可W基于指示函数，获取指示函数的值为零的第一网格，并根据获取的第一网格的边的 连通状态，判断获取的第一网格是否为完整单元，然后根据判断结果，基于预设规则从获取 的第一网格中提取曲面，最 后将提取的曲面组成的模型作为层面模型。图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日
最新回复(0)