基于petrel插件的SIMPAT多点地质统计学建模算法实现

Ocean为我们提供了对petrel的二次开发平台，开放了petrel中几乎所有数据接口，可以根据需要编写新的功能。最近把SIMPAT建模算法以plugin的形式实现了。

Petrel是各大建模软件中最早引入的多点地质统计学建模算法的，是SNESIM算法，而且到目前似乎其它软件也都没引入。但从2009年加入SNESIM之后这个模块就没什么更新，最近几年文献中新提出的其它多点地质统计学建模算法都没有加进来，比如基于模式的SIMPAT算法。

自从不久前加入中石化石勘院Duan团队，他大力支持多点地质统计学新算法在petrel中实现，经过一段时间的探索终于把测试版的SIMPAT做出来了。在这里不在赘述SIMPAT算法的具体内容，感兴趣的可以参考斯坦福大学Arpat的论文。

先看界面。安装插件之后在FaciesModeling模块下有个simpatAlgorith：

 

插件使用包括两个步骤：

（1）输入参数，然后点makePattern按钮，生成模式库。

（2）点击Apply或OK进行计算，有进度条提示计算。测试版，计算过程比较慢，尽管已经用了并行计算。

参数输入很简单：

（1）训练图像

（2）模式库名称（任意）

（3）模板在I、J、K方向的网格数

（4）自定义模式个数（可选）

尽管输入的参数很简单，但具体应该输多少却不容易。比如什么样的训练图像合适，多大的模板大小合适等，需要自己摸索，根据实际情况调整。

相比于文献中的SIMPAT算法，该插件提供了自定义模式个数的可选选项，给的个数越多模拟结果更可靠，但计算也越慢。给的太多也许训练图像中根本没有那么多模式。而不选的话它会自动进行全部扫描，尽管内置了个不大不小的步长。

下图是一个曲流河构型模型的简单测试。为了提高计算速度设置的模式个数比较小，但已经可以看出SIMPAT比SNESIM的优势，SIMPAT计算的侧积层更加符合规律。

 

插件是个叫simpat_installer.pip的15kb的文件，已上传在“油藏地质和开发交流63231398”群共享，欢迎加群下载使用，也可以发邮件给yanfeng1022@126.com索取。具体安装请参考petrel帮助文件。