技术摘要：
本发明实施例公开了一种电阻率正演仿真方法、装置及计算机可读存储介质，应用GPU正演仿真平台，包括：获取待仿真处理区域中，与仿真处理相关的地层模型参数；根据所确定的每个仿真节点并行处理区块的份数对待仿真处理区域进行网格剖分，得到多个区块；在每个仿真节点并  全部
背景技术：
目前电磁场数值模拟仿真研究对于测井仪器研制和后期资料处理起到越来越关 键的作用，在测井仪器研制过程中通过模拟仿真各仪器结构参数对仪器测量响应的影响， 帮助优化仪器结构参数，给出仪器的最优探测性能，为仪器研制提供指导性建议，而在资料 解释方面其可模拟仿真各种环境因素对仪器测量响应的影响，帮助测井解释人员掌握仪器 的测井响应规律及建立各种环境影响校正图版等，同时对于实际测井中遇到的特殊储层的 响应规律，也可以通过数值模拟仿真进行验证说明，为测井解释提供一定的理论指导。关于 电阻率测井正演仿真模拟研究，目前常用的方法是有限元素法，该有限元素法的特点为适 合于复杂边界形状形成的几何物体，并适合含有复杂边界条件和复杂媒质的定解问题求 解，但该方法所存在的问题是计算速度比较慢。 因此，针对现有技术中的电阻率测井正演仿真模拟计算中所存在的问题，如何在 保证电阻率测井正演仿真模拟计算精度的情况下，实现一种可提高电阻率正演仿真速度的 仿真方法。
技术实现要素：
为了解决上述技术问题，本发明提供了一种电阻率正演仿真方法、装置及计算机 可读存储介质，可以提高电阻率正演仿真计算的速度。 本发明提供了一种电阻率正演仿真方法，应用GPU正演仿真平台，包括： 获取待仿真处理区域中，与仿真处理相关的地层模型参数； 根据每个仿真节点并行处理区块的份数对待仿真处理区域进行网格剖分，得到多 个区块； 在每个仿真节点并行处理所述多个区块的过程中，对每个区块分别构建区块元 素，并根据所获取的地层模型参数对所述区块元素进行电阻率赋值； 分别对每个区块中电阻率赋值后的区块元素进行安装及消元； 当安装及消元所有区块后，计算得到待仿真处理区域的测量响应； 输出包括测量响应的仿真结果。 一种示例性的实施例中，所述GPU正演仿真平台包括：具有独立显卡的GPU计算机， 该GPU计算机包括多个仿真节点。 一种示例性的实施例中，所述根据每个仿真节点并行处理区块的份数对待仿真处 理区域进行网格剖分之前，还包括： 根据所获取的地层模型参数、待仿真处理区域的大小和预先设置的初始网格的划 分，确定每个仿真节点并行处理区块的份数。 4 CN 111581791 A 说　明　书 2/18 页 一种示例性的实施例中，所述根据每个仿真节点并行处理区块的份数对待仿真处 理区域进行网格剖分之前，还包括： 预先获取每个仿真节点的性能参数； 确定每个仿真节点并行处理区块的份数之后，根据每个仿真节点的性能参数、所 述并行处理区块的份数以及所述待仿真处理区域的数量，确定待仿真处理区域并行处理的 个数； 所述地层模型参数包括：测量点的数量； 所述输出包括测量响应的仿真结果包括： 根据待仿真处理区域并行处理的个数输出多个测量点的仿真结果。 一种示例性的实施例中，所述根据每个仿真节点并行处理区块的份数对待仿真处 理区域进行网格剖分，得到多个区块，包括： 根据每个仿真节点并行处理区块的份数N和所述初始网格的划分对待仿真处理区 域进行剖分； 分别将剖分后N个区块中的每个区块边界的值插入到相应的径向、纵向和周向边 界初始网格中，得到网格剖分后的N个区块。 一种示例性的实施例中，所述在每个仿真节点并行处理所述多个区块的过程中， 对每个区块分别构建区块元素，并根据所获取的地层模型参数对所述区块元素进行电阻率 赋值，包括： 在每个仿真节点并行处理所述多个区块的过程中，针对每个区块分别执行以下操 作： 根据径向、纵向和周向网格剖分的份数，确定每个区块构建区块元素时的循环顺 序； 按照循环顺序依次构建离散的区块元素，并对所构建的区块元素进行编号； 根据每个区块元素的编号和所获取的地层模型参数分别对每个区块元素进行电 阻率赋值。 一种示例性的实施例中，所述分别对每个区块中电阻率赋值后的区块元素进行安 装及消元，包括： 对每个区块中的区块元素进行电导阵的计算； 从区块的边界顶点向区块的边界处进行区块元素安装及消元； 当只剩下该区块边界处的区块元素时，对该边界处的区块元素进行消元。 一种示例性的实施例中，所述当安装及消元所有区块后，计算得到待仿真处理区 域的测量响应包括： 对每个区块完成安装及消元后，得到待仿真处理区域的总电导阵； 根据该待仿真处理区域的总电导阵确定仪器的测量响应。 本发明还提供了一种电阻率正演仿真装置，所述装置包括：存储器和处理器； 所述存储器用于保存用于电阻率正演仿真的程序； 所述处理器用于读取执行所述用于电阻率正演仿真的的程序，执行如下操作： 获取待仿真处理区域中，与仿真处理相关的地层模型参数； 根据每个仿真节点并行处理区块的份数对待仿真处理区域进行网格剖分，得到多 5 CN 111581791 A 说　明　书 3/18 页 个区块； 在每个仿真节点并行处理所述多个区块的过程中，对每个区块分别构建区块元 素，并根据所获取的地层模型参数对所述区块元素进行电阻率赋值； 分别对每个区块中电阻率赋值后的区块元素进行安装及消元； 当安装及消元所有区块后，计算得到待仿真处理区域的测量响应； 输出包括测量响应的仿真结果。 一种示例性的实施例中，所述GPU正演仿真平台包括：具有独立显卡的GPU计算机， 该GPU计算机包括多个仿真节点。 一种示例性的实施例中，所述根据每个仿真节点并行处理区块的份数对待仿真处 理区域进行网格剖分之前，还包括： 根据所获取的地层模型参数、待仿真处理区域的大小和预先设置的初始网格的划 分，确定每个仿真节点并行处理区块的份数。 一种示例性的实施例中，所述根据每个仿真节点并行处理区块的份数对待仿真处 理区域进行网格剖分之前，还包括： 预先获取每个仿真节点的性能参数； 确定每个仿真节点并行处理区块的份数之后，根据每个仿真节点的性能参数、所 述并行处理区块的份数以及所述待仿真处理区域的数量，确定待仿真处理区域并行处理的 个数； 所述地层模型参数包括：测量点的数量； 所述输出包括测量响应的仿真结果包括： 根据待仿真处理区域并行处理的个数输出多个测量点的仿真结果。 一种示例性的实施例中，所述根据每个仿真节点并行处理区块的份数对待仿真处 理区域进行网格剖分，得到多个区块，包括： 根据每个仿真节点并行处理区块的份数N和所述初始网格的划分对待仿真处理区 域进行剖分； 分别将剖分后N个区块中的每个区块边界的值插入到相应的径向、纵向和周向边 界初始网格中，得到网格剖分后的N个区块。 一种示例性的实施例中，所述在每个仿真节点并行处理所述多个区块的过程中， 对每个区块分别构建区块元素，并根据所获取的地层模型参数对所述区块元素进行电阻率 赋值，包括： 在每个仿真节点并行处理所述多个区块的过程中，针对每个区块分别执行以下操 作： 根据径向、纵向和周向网格剖分的份数，确定每个区块构建区块元素时的循环顺 序； 按照循环顺序依次构建离散的区块元素，并对所构建的区块元素进行编号； 根据每个区块元素的编号和所获取的地层模型参数分别对每个区块元素进行电 阻率赋值。 一种示例性的实施例中，所述分别对每个区块中电阻率赋值后的区块元素进行安 装及消元，包括： 6 CN 111581791 A 说　明　书 4/18 页 对每个区块中的区块元素进行电导阵的计算； 从区块的边界顶点向区块的边界处进行区块元素安装及消元； 当只剩下该区块边界处的区块元素时，对该边界处的区块元素进行消元。 一种示例性的实施例中，所述当安装及消元所有区块后，计算得到待仿真处理区 域的测量响应包括： 对每个区块完成安装及消元后，得到待仿真处理区域的总电导阵； 根据该待仿真处理区域的总电导阵确定仪器的测量响应。 本发明还提供了一种计算机可读存储介质，存储有计算机可执行指令，所述计算 机可执行指令被处理器执行如下操作： 获取待仿真处理区域中，与仿真处理相关的地层模型参数； 根据每个仿真节点并行处理区块的份数对待仿真处理区域进行网格剖分，得到多 个区块； 在每个仿真节点并行处理所述多个区块的过程中，对每个区块分别构建区块元 素，并根据所获取的地层模型参数对所述区块元素进行电阻率赋值； 分别对每个区块中电阻率赋值后的区块元素进行安装及消元； 当安装及消元所有区块后，计算得到待仿真处理区域的测量响应； 输出包括测量响应的仿真结果。 与现有技术相比，本发明提供了一种电阻率正演仿真方法、装置及计算机可读存 储介质，应用GPU正演仿真平台，包括：获取待仿真处理区域中，与仿真处理相关的地层模型 参数；根据每个仿真节点并行处理区块的份数对待仿真处理区域进行网格剖分，得到多个 区块；在每个仿真节点并行处理所述多个区块的过程中，对每个区块分别构建区块元素，并 根据所获取的地层模型参数对所述区块元素进行电阻率赋值；分别对每个区块中电阻率赋 值后的区块元素进行安装及消元；当安装及消元所有区块后，计算得到待仿真处理区域的 测量响应；输出包括测量响应的仿真结果。通过本发明的方案，应用GPU正演仿真平台，有效 提高电阻率正演仿真的速度。 本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变 得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利 要求书以及附图中所特别指出的结构来实现和获得。 附图说明 附图用来提供对本发明技术方案的进一步理解，并且构成说明书的一部分，与本 申请的实施例一起用于解释本发明的技术方案，并不构成对本发明技术方案的限制。 图1为本发明实施例的电阻率正演仿真方法的流程图； 图2为一个示例中待仿真地层模型示意图； 图3为一个示例中区块元素形成、安装及消元示意图； 图4为一个示例中基于GPU平台进行电阻率正演仿真计算效率对比示意图； 图5为一个示例中电阻率正演仿真测量响应示意图，其中RA1,RA2,RA3,RA4表示四 种不同电流聚焦模式测量响应曲线。 图6为本发明实施例的电阻率正演仿真装置示意图。 7 CN 111581791 A 说　明　书 5/18 页