-
储层建模包括多个步骤，不同学者有不同的划分方法，如二步建模[1]、三步建模以及多步建模[2]，沉积相建模是其中的一个重要环节。通常会先建立沉积相模型，然后在相模型基础上建立储层物性参数模型。为了评估建模参数不确定性的影响，通常还需要建立概率储量[3-5]。一般来说，不同相带物性参数的分布特征是不一样的，在相模型基础之上建立物性参数模型能够更加合理地刻画物性参数的空间分布，中国学者称之为相控建模[6]，已得到了广泛的应用[7-8]，并且这个概念及应用领域也被进一步扩大，如多级相控建模、碳酸盐岩相控建模、地质模式约束建模等[9-11]。有多种建立相模型的技术，包括基于两点地质统计学的方法（序贯指示模拟、截断高斯模拟），基于目标的建模方法、多点地质统计学方法（MPS，Multiple-points Statistics Methods） 和基于沉积事件（过程）的建模方法[12]。自 1993 年提出MPS方法，在相建模方面展示出良好的应用前景，该方法综合了两点地质统计学建模方法和基于目标建模方法的优点，既能较好地刻画复杂地质体的几何形态，又容易满足已有的条件数据。MPS方法逐渐成为研究的热点，中外学者对经典的MPS方法提出了很多改进[13-17]。MPS方法基于训练图像进行模拟，例如经典的基于概率的单一正规方程模拟 （Snesim）方法，其核心技术包括：一是对训练图像建立搜索树，二是从搜索树中获取条件概率[18]；基于样式的 Simpat方法，通过计算数据事件与训练图像中样式的差异进行模拟[19]，能够更好地再现地质模式。针对一个工区通常基于一个训练图像进行模拟，除非是对该工区地质模式认识上存在较大不确定性，可以采取几张训练图像代表不同的地质模式[20]。基于训练图像 MPS 方法的重点是提取训练图像的模式，但是这些算法不能生成真实的非线性模式，也不能展现地质推断的可变性和不确定性， DUPONT认为基于GAN模拟可以解决这些问题[21]。
-
1968年，IVAKHNENKO公开了感知器（Perception）的概念，被认为是深度学习的起源[22]。2014 年，IAN 等提出了 GAN 方法[23]，主要包含生成器和判别器，两者相互博弈，该方法受到广大学者的关注，并在其基础上做了大量的改进研究。近年来，人工智能方法在油气领域的应用越来越广泛，例如以支持向量机回归和XGBoost回归算法为基础的储层裂缝开度预测，结合卡尔曼滤波和人工神经网络进行油藏井间连通性研究，基于测井的岩相自动识别，应用多次数据吸收集合平滑算法进行自动油藏历史拟合[24-28]。ERIC 等提出基于训练图像的空间生成对抗神经网络地质统计学反演方法，该方法能够快速生成二维和三维的非条件模拟结果[29]。 COIFFIER 等使用生成对抗神经网络将二维样本合成为三维地质图像，该方法实现简单，在标准 GAN 的生成器和判别器之间引入随机截取采样步骤[30]。ZHANG 等使用基于沉积结构的生成对抗神经网络方法生成逼真的三维相模型，该方法可以再现各种概念地质模型，同时具有满足约束条件（如井资料） 所必需的灵活性[31]。目前利用 GAN 方法建模需要大量的训练图像，基于 15 000张训练图像学习开展模拟预测，尽管可以得到与训练图像非常逼真的模拟结果，但是大量训练图像的生成与学习训练增加了该方法应用的难度。SinGAN 方法能够从单个自然景观图像生成新的模型，克服了传统的生成对抗神经网络需要大量的图像进行训练的限制[32]。为此，笔者尝试将 SinGAN 方法应用到相建模中，以辫状河储层为例，并与经典的基于样式的 Simpat方法进行对比，SinGAN方法展示出良好的应用前景。
-
1 SinGAN方法基本原理
-
GAN 方法可以生成与期望地质模式高度匹配并且满足条件数据的模型，但是训练数据必须基于大量的特定训练集[24]。SinGAN 方法能够从单个自然图像中学习并生成模型[26]。
-
SinGAN 方法是一种金字塔式的全卷积轻量权重的 GAN 方法，与常见的 GAN 方法类似，其重要改进之处是训练样本采用的是单个训练图像，核心在于通过不同尺度来捕获复杂图像结构的统计特征。模型包含金字塔式的生成器和训练对应的金字塔式图像。降采样是一种图像处理技术，降采样的结果是原始图像的缩略图。训练过程是从最粗尺度开始，多次迭代到最细尺度。所有的生成器和判别器具有相同的感受野，因此随着尺度的逐步增加能够捕获到更小尺寸的图像结构。在最粗尺度上由生成器根据高斯分布的噪声生成的图像，其表达式为：
-
在最粗尺度下，有效感受野大约是训练图像尺寸的一半，因此生成器生成的是图像的布局和全局结构。每一个更细尺度的生成器会增添上一个生成器没有生成的细节部分。因此，除了空间噪声，每一个生成器还要接受上一个尺度的上采样图像数据，其表达式为：
-
所有的生成器结构类似，特别之处在于在噪声进入卷积层之前，噪声会被添加到上一个尺度生成的图像数据中，这既保证 GAN 方法未忽视噪声，又避免完全随机。卷积层目的是生成上一个尺度所缺乏的细节，其表达式为：
-
卷积核大小为 3×3，在最粗尺度上每个卷积块使用 32 个卷积核，每增加一个尺度，卷积核数量增加16个。因为生成器是全连接的卷积网络，通过修改噪声的维度，可以生成任意大小和长宽比例的图像。
-
训练过程采用从最粗尺度到最细尺度的顺序模型。一个 GAN 方法训练完毕，就保持固定状态。训练损失包括对抗损失和重建损失，其表达式为：
-
对抗损失用来减少训练图像采样后的图像中块的分布和生成样本中块的分布之间的差异，重建损失确保存在可以产生训练图像降采样后的图像的特定噪声图集。每一个生成器与一个马尔科夫判别器耦合，判别器将输入的每个重叠块分为真和假两类。SinGAN 方法定义的损失是整个图像的损失，这样网络可以学习边界条件。马尔科夫判别器的结构和生成器中的由 BatchNorm-LeakyReLU 的 5 个卷积块组成的全连接卷积网络相同，因此块的大小也相同。要生成原始的图像 x，必须保证存在一组特定的输入噪声集。重建损失的表达式为：
-
当尺度数量为 N 时，重建的图像在训练的过程中还有一个重要的作用，重建的图像决定了每一个尺度上噪声的标准差。
-
SinGAN方法的结构尺度数量会影响模拟结果。当尺度数量较少时，粗尺度的有效感受野更小，只能捕获细微的特征。随着尺度数量增加，更高质量的结构和整体目标布局会更好地保存下来，图像结构随着尺度数量增加而变化。
-
2 实例应用与分析
-
N气田位于中国东海陆架盆地某凹陷的中央反转构造带中北部，为中国东海大型低渗透气田，已经勘探的3套气藏储层的厚度均大于100 m，总气藏的厚度大于 400 m。其主力层段花港组上段主要发育灰色泥岩、粉砂质泥岩和泥质粉砂岩等，且夹有薄层细砂岩；花港组下段主要发育厚层细砂岩，局部发育中粗砂岩，并夹杂少量粉砂岩及泥岩。
-
为了检验提出方法的合理性，选取了 N 气田 2 个有代表性的小层沉积微相图作为训练图像，包括河道、心滩和泥岩，训练图像的大小为256×159个网格单元。其中图1a平稳性较好，图1b较差，在训练图像上部与下部具有较大差异，上部的河道间泥岩更为发育。为对比该方法的适用性，采用 2 组不同的训练图像。测试硬件环境包括：CPU 采用第八代英特尔酷睿 i7处理器，频率为 3.7 GHz，运行内存为 16 G，显卡采用英伟达GeForce GTX1 080 Ti。
-
图1 2个小层的沉积微相图
-
Fig.1 Microfacies maps of two layers
-
2.1 利用MPS方法建立沉积微相模型
-
Simpat方法由ARPAT提出[19]，与基于概率统计的 Snesim方法不同，该方法将图像重建的思想引入到储层建模，主要模拟步骤包括：①以数据样板扫描训练图像，构建样式库。②在模拟过程中，计算数据事件与样式库中所有样式的相似性，Simpat 方法采用 Manhattan 距离计算数据事件与样式的相似度。③将最相似的样式整体粘贴到模拟网格中。 ④继续模拟，直至所有网格都被模拟。
-
利用 Simpat 方法分别采用 3×3，5×5，7×7，9×9， 11×11，13×13，15×15，17×17共8组数据样板进行测试。由基于相对平稳训练图像 Simpat 方法模拟结果（图2）和基于非平稳训练图像的模拟结果（图3） 可以直观地看出，数据样板大小对于模拟结果的影响很大，数据样板过小模拟结果很离散，过大容易失真，且计算量增加。对于 MPS 方法来说，合理确定数据样板大小对于建模质量的好坏具有重要的影响。另外，相对平稳的训练图像模拟效果比非平稳的训练图像好。
-
2.2 利用SinGAN方法建立沉积微相模型
-
基于 SinGAN 方法的相关理论，搭建了相应的 PyTorch 环境。初始化参数设置包括噪声、尺度因子、输出路径、种子数等，本次设置尺度因子为0.75，种子数使用程序随机产生。训练的过程会根据尺度数量以金字塔式的顺序依次从最粗尺度模拟至最细尺度，每个尺度迭代次数可以人为设置，本次设置为 2 000 次，目标函数重建损失前的系数设置为10。SinGAN方法除了具有GAN方法共有的核心生成器和判别器，还需要进行尺度数量的选择。尺度数量的选取一种是自适应的方式，主要是根据图像的像素大小计算选取，另一种方式是人工指定尺度数量。
-
图2 基于相对平稳训练图像的Simpat方法模拟结果
-
Fig.2 Simulation results of Simpat based on stationary training image
-
图3 基于非平稳训练图像的Simpat方法模拟结果
-
Fig.3 Simulation results of Simpat based on non-stationary training image
-
首先利用 SinGAN 方法对相对平稳的训练图像进行模拟，选择图1a相对平稳的辫状河训练图像进行测试。分别采用尺度数量为 3，5和 7进行模拟计算，设置成图的数量为 50 幅，部分模拟结果如图4 所示。
-
从图4 中可以看出，N 值对模拟结果的影响很大，N 值较小的时候，会出现大量的过渡色，因为 SinGAN方法是针对自然景观图像，是连续变量。当 N为 7时，效果良好，能够刻画不同沉积微相的结构和空间展布规律，直观对比发现，优于经典 Simpat 方法模拟结果。
-
为进一步验证 SinGAN 方法对非平稳训练图像的模拟效果，选择图1b非平稳的辫状河训练图像进行测试。同样采取N为3，5和7进行模拟计算，模拟结果如图5所示。
-
图4 基于相对平稳训练图像的SinGAN方法模拟结果
-
Fig.4 Simulation results of SinGAN based on stationary training image
-
由图5 可以看出，N 值对模拟结果的影响更加明显，N为3时，明显失真，N为5变形也比较严重，N 为 7 效果较好。进一步增加 N 值容易导致系统崩溃，一般不要超过 7 [26]。对比 Simpat 方法模拟结果可以看出，SinGAN 方法优于经典的 MPS方法，与训练图像刻画的沉积微相空间结构更相似。因为Sin⁃ GAN方法是基于自然景观图像，利用该方法建立沉积相模型会出现一些过渡色，即使当采用 N 为 7 时也会在个别地方出现，且非平稳的训练图像这种现象更明显一些。这个问题可以采用聚类的方法进行后处理，也就是对带过渡色的图像根据RGB值进行聚类，判断每个像素的颜色与 3 种颜色的差异，将最接近的那种颜色赋值给过渡色的像素，从而得到只有3种沉积相类型的模拟结果。非平稳训练图像基于 N 为 6时得到模拟结果 3的聚类后处理结果 （图6），在不同沉积相类型的接触边界，原始模拟结果（图6a）呈现过渡色，通过聚类后处理的结果（图6b）有效地消除这种过渡色，模型只保留了3种沉积相类型。
-
图5 基于非平稳训练图像的SinGAN方法模拟结果
-
Fig.5 Simulation results of SinGAN based on non-stationary training image
-
图6 聚类前后结果对比
-
Fig.6 Comparison of results before and after clustering
-
2.3 模拟结果定量对比
-
为了进一步比较 SinGAN 方法与经典 MPS方法的样式再现能力，选取了基于 Simpat方法模拟效果较好的 2 个模拟结果（分别对应数据样板为 13×13 和15×15），N为7的SinGAN方法模拟结果和原始训练图像（相对平稳）进行样式相似性对比分析。将模拟结果中的样式全部提取出来，与训练图像所包含的样式进行比对，采用逐网格属性比对的方式刻画样式的相似性。其表达式为：
-
通过计算模拟结果的每个样式的相似性，即可统计比较模拟结果与训练图像的样式相似性差异。
-
由对比结果（图7）可知，SinGAN 方法模拟结果与原始训练图像的样式相似性最好，其次是基于 13×13 样板模拟的结果。通过对比分析表明，SinGAN 方法在再现训练图像内部结构方面具有强大的能力，优于经典的 MPS 方法，展示了良好的应用前景。本次仅开展了二维非条件的沉积相模拟，进一步实现条件化的三维模拟将极大推进该方法在储层建模中的实际应用。
-
图7 数据样式相似性对比
-
Fig.7 Comparison of similarity of data patterns
-
3 结论
-
基于单幅自然景观图像的 SinGAN 方法应用于沉积相建模，取得了良好的效果，通过直观分析及数据样式再现的定量化对比，SinGAN方法模拟的结果优于 Simpat 方法。SinGAN 方法模拟结果不论是相对平稳的沉积微相图还是非平稳的沉积微相图受模拟尺度数量的影响明显。相对而言，对于非平稳的沉积微相图影响更大，需要选取更大的尺度数量才能够达到较好的模拟效果。对于沉积相建模一般选取最粗尺度数量为 7，过大可能会导致系统崩溃。当 SinGAN 方法用于建立离散的沉积相模型时，即使采用最粗尺度数量为 7 也会在个别地方出现一些过渡色，非平稳情况下更为明显。针对这种情况可以采取聚类的方式处理，消除过渡类型，达到模拟离散模型的目的。SinGAN 方法是针对连续变量的，可以直接应用于连续变量的模拟。目前主要进行了二维非条件的沉积相模拟测试，三维空间的模拟可以借鉴已有做法，需要考虑不同尺度相互匹配的问题。如何将条件数据、地震属性等更多的约束信息加入到 SinGAN 方法仍是一个难点，但该方法在储层建模领域已经展示出良好的应用前景。
-
符号解释
-
Dist ——计算模拟结果的第i个样式与PatDBT中的每个样式的距离值，无量纲；
-
Dn ——马尔科夫判别器；
-
Gn ——生成器；
-
GN ——金字塔式的生成器；
-
i ——编号；
-
L adv ——对抗损失，无量纲；
-
L rec——重建损失，无量纲；
-
n ——当前尺度数量；
-
N ——最粗尺度数量；
-
Pat（ i）——模拟结果的第i个样式；
-
PatDBT ——训练图像里所有的样式；
-
S（ i）——训练图像中与模拟结果的第i个样式最相似样式的相似度；
-
xn ——当前尺度数量为n时的图像；
-
——对应xn生成器生成的图像；
-
——当前尺度数量为n时噪声集生成的图像；
-
——最粗尺度数量为N时生成的图像；
-
zn ——当前尺度数量为n时的空间噪声图像；
-
zN ——最粗尺度数量为N时的空间噪声图像；
-
α——目标函数重建损失前的系数，无量纲；
-
↑ r ——重采样；
-
ψn——BatchNorm-LeakyReLU的5个卷积块组成的全连接卷积网络。
-
参考文献
-
[1] HALDORSEN H H，DAMSLETH E.Stochastic modeling［J］.Jour⁃ nal of Petroleum Technology，1990，42（48）：404-412.
-
[2] CANNON Stephen.Reservoir modelling：A practical guide［M］.chichester：John Wiley & Sons，2018.
-
[3] 霍春亮，刘松，古莉，等.一种定量评价储集层地质模型不确定性的方法［J］.石油勘探与开发，2007，34（5）：574-579.HUO Chunliang，LIU Song，GU Li，et al.A quantitative method for appraising reservoir geological model uncertainty［J］.Petroleum Exploration and Development，2007，34（5）：574-579.
-
[4] 叶小明，王鹏飞，霍春亮，等.海上复杂碎屑岩储层油气藏地质建模关键技术［J］.中国海上油气，2018，30（3）：110-115． YE Xiaoming，WANG Pengfei，HUO Chunliang，et al.Key tech⁃ niques for geological modelingof offshore complex clastic rock res⁃ ervoirs［J］.China Offshore Oil and Gas，2018，30（3）：110-115.
-
[5] 王鸣川，商晓飞，段太忠，等.海外油气田开发新项目储量评价方法［J］.石油实验地质，2020，42（2）：296-301.WANG Mingchuan，SHANG Xiaofei，DUAN Taizhong，et al.Re⁃ serves evaluation for new investment projects in overseas oil and gas field development［J］.Petroleum Geology & Experiment，2020，42（2）：296-301.
-
[6] 蒋建平，康贤，邓礼正.储层物性参数展布的相控模型［J］.成都理工学院学报，1995，22（1）：12-17.JIANG Jianping，KANG Xian，DENG Lizheng.Phase control mod⁃ el for the distribution of reservoir parameters［J］.Journal of Cheng⁃ du Institute of Technology，1995，22（1）：12-17.
-
[7] 桂峰，黄智辉，马正.利用相控模型进行井间参数预测［J］.地球科学，2001，26（1）：49-53.GUI Feng，HUANG Zhihui，MA Zheng.Application of lithofaciescontrolled model to prediction of physical parameters between wells［J］.Earth Science，2001，26（1）：49-53.
-
[8] 李少华，张昌民，张尚锋，等.沉积微相控制下的储层物性参数建模［J］.江汉石油学院学报，2003，25（1）：24-26.LI Shaohua，ZHANG Changmin，ZHANG Shangfeng，et a1.Reser⁃ voir petrophysical modeling under control of sedimentary facies ［J］.Journal of Jiang han Petroleum Institute，2003，25（1）：24-26.
-
[9] 廖保方，龙国清，刘卓，等.非均质河道砂体多级相控建模方法探讨［J］.石油天然气学报，2010，32（3）：32-37.LIAO Baofang，LONG Guoqing，LIU Zhuo，et al.Multiple faciescontrolled modeling for anisotropic channel sands［J］.Journal of Oil and Gas Technology，2010，32（3）：32-37.
-
[10] 李辉，林承焰，马存飞，等.基于地质模式约束的致密砂岩储层断缝体地质建模［J］.大庆石油地质与开发，2021，40（1）：38-46.LI Hui，LIN Chengyan，MA Cunfei，et al.Geological modeling of the fault-fracture body in the tight sandstone reservoir based on the geological model constraint［J］.Petroleum Geology & Oilfield Development in Daqing，2021，40（1）：38-46.
-
[11] 兰义飞，王东旭，樊友宏.相控建模在剥蚀残留碳酸盐岩储层中的应用［J］.天然气工业，2007，27（8）：52-53.LAN Yifei，WANG Dongxu，FAN Youhong.Application of faciescontrolling modeling in the erosion carbonate reservoir［J］.Natu⁃ ral Gas Industry，2007，27（8）：52-53.
-
[12] 李少华，张昌民，尹艳树.储层建模算法剖析［M］.北京：石油工业出版社，2012.LI Shaohua，ZHANG Changmin，YIN Yanshu.Analysis of reser⁃ voir modeling algorithm［M］.Beijing：Petroleum Industry Press，2012.
-
[13] GREGOIRE Mariethoz，JEF Caers.Multiple-point geostatistics：Stochastic modeling with training images［M］.Hoboken：WileyBlackwell，2014.
-
[14] 喻思羽，李少华，何幼斌，等.基于样式降维聚类的多点地质统计建模算法［J］.石油学报，2016，37（11）：1 403-1 409.YU Siyu，LI Shaohua，HE Youbin，et al.Multiple-point geostatis⁃ tics algorithm based on pattern scale-down cluster［J］.Acta Petro⁃ lei Sinica，2016，37（11）：1 403-1 409.
-
[15] 喻思羽，李少华，段太忠，等.基于局部各向异性的非平稳多点地质统计学算法［J］.物探与化探，2017，41（2）：262-269.YU Siyu，LI Shaohua，Duan Taizhong，et al.Non-stationary multi⁃ ple-point geostatistics algorithm base on local anisotropy［J］.Geo⁃ physical & Geochemical Exploration，2017，41（2）：262-269.
-
[16] 喻思羽，李少华，王端平，等.基于 p-stable LSH 的多点地质统计建模算法.石油学报，2017，38（12）：1 425-1 433.YU Siyu，LI Shaohua，WANG Duanping，et al.Multipoint geosta⁃ tistical modeling algorithm based on p-stable LSH［J］.Acta Petro⁃ lei Sinica，2017，38（12）：1 425-1 433.
-
[17] 麻平山，李少华，卢昌盛，等.基于众数法聚类的多点地质统计学方法［J］.天然气地球科学，2020，31（6）：809-817.MA Pingshan，LI Shaohua，LU Changsheng，et al.Multi-point geo⁃ statistical method based on mode method clustering［J］.Natural Gas Geoscience，2020，31（6）：809-817.
-
[18] SEBASTIEN Strebelle.Conditional simulation of complex geologi⁃ cal structures using multiple-point statistics［J］.Mathematical Ge⁃ ology，2002，34（1）：1-21.
-
[19] ARPAT G B，CAERS J.Conditional simulation with patterns［J］.Mathematical Geoscience，2007，39（2）：177-203.
-
[20] 王立鑫，尹艳树，文杰，等.多点地质统计学中训练图像优选方法及其在地质建模中的应用［J］.石油勘探与开发，2019，46（4）：703-709.WANG Lixin，YIN Yanshu，WEN Jie，et al.A training image opti⁃ mization method in multiple-point geostatistics and its applica⁃ tion in geological modeling［J］.Petroleum Exploration and Devel⁃ opment，2019，46（4）：703-709.
-
[21] DUPONT E，ZHANG T F，TILKE P，et al.Generating realistic ge⁃ ology conditioned on physical measurements with generative ad⁃ versarial networks［J/OL］.2018.http：//arxiv.org/abs/1802.03065v3.
-
[22] IVAKHNENKO Alexey Grigorevich.The group method of data of handling-a rival of the method of stochastic approximation［J］.So⁃ viet Automatic Control，1968，13（3）：43-55.
-
[23] IAN Goodfellow，JEAN Pouget Abadie，MEHDI Mirza，et al.In ad⁃ vances in neural information processing systems［C］.Montreal：An⁃ nual Conference on Neural Information Processing Systems，2014：2 672-2 680.
-
[24] 孙致学，姜宝胜，肖康，等.基于新型集成学习算法的基岩潜山油藏储层裂缝开度预测算法［J］.油气地质与采收率，2020，27（3）：32-38.SUN Zhixue，JIANG Baosheng，XIAO Kang，et al.Prediction of fracture aperture in bedrock buried hill oil reservoir based on nov⁃ el ensemble learning algorithm［J］.Petroleum Geology and Recov⁃ ery Efficiency，2020，27（3）：32-38.
-
[25] 刘巍，刘威，谷建伟，等.利用卡尔曼滤波和人工神经网络相结合的油藏井间连通性研究［J］.油气地质与采收率，2020，27（2）：118-124.LIU Wei，LIU Wei，GU Jianwei，et al.Research on interwell con⁃ nectivity of oil reservoirs based on Kalman filter and artificial neu⁃ ral network［J］.Petroleum Geology and Recovery Efficiency，2020，27（2）：118-124.
-
[26] 李昌，沈安江，常少英，等.机器学习法在碳酸盐岩岩相测井识别中应用及对比——以四川盆地MX地区龙王庙组地层为例 ［J］.油气藏评价与开发，2021，11（4）：586-596.LI Chang，SHEN Anjiang，CHANG Shaoying，et al.Application and contrast of machine learning in carbonate lithofacies log iden⁃ tification：A case study of Longwangmiao Formation of MX area in Sichuan Basin［J］.Reservoir Evaluation and Development，2021，11（4）：586-596.
-
[27] 石玉江，刘国强，钟吉彬，等.基于大数据的测井智能解释系统开发与应用［J］.中国石油勘探，2021，26（2）：113-126.SHI Yujiang，LIU Guoqiang，ZHONG Jibin，et al.Development and application of intelligent logging interpretation system based on big data［J］.China Petroleum Exploration，2021，26（2）：113-126.
-
[28] 王泽龙，刘先贵，唐海发，等.基于多次数据吸收集合平滑算法的自动油藏历史拟合研究［J］.特种油气藏，2021，28（3）：99-105.WANG Zelong，LIU Xiangui，TANG Haifa，et al.Study on auto⁃ matic reservoir history matching based on ES-MDA algorithm［J］.Special Oil & Gas Reservoirs，2021，28（3）：99-105.
-
[29] ERIC Laloy，ROMAIN Herault，DIEDERIK Jacques，et al.Train⁃ ing-image based geostatistical inversion using a spatial generative adversarial neural network［J］.Water Resources Research，2018，54（1）：381-406.
-
[30] COIFFIER G，RENARD P.3D Geological image synthesis from 2D examples using generative adversarial networks［C］.Florence：Fourth EAGE Conference on Petroleum Geostatistics，2019：2-6.
-
[31] ZHANG Tuanfeng，PETER Tilke，EMILIEN Dupont，et al.Gener⁃ ating geologically realistic 3D reservoir facies models using deep learning of sedimentary architecture with generative adversarial networks［J］.Petroleum Science，2019，16（3）：541-549.
-
[32] SHAHAM T R，DEKEL T，MICHAELI T.SinGAN：learning a gen⁃ erative model from a single natural image［C］.Seoul：IEEE/CVF International Conference on Computer Vision（ICCV），2019：4 569-4 579.
-
摘要
沉积相建模是储层建模中的一个重要环节，有多种方法可以用来建立沉积相模型。传统的建模方法需要利用各种参数对变量的空间结构信息进行刻画，如变差函数、数据样式等，在模拟中再现这种空间结构。利用生成对抗神经网络方法（GAN，Generative Adversarial Nets）建模采用了不同的策略，通过对大量图像（模型）的学习，生成与学习样本具有高度相似特征的模型。基于单一图像生成对抗神经网络方法（SinGAN，Generative Adversarial Nets based on single image）对传统的GAN方法进行改进，仅需一张图像进行训练就能够生成高度相似的图像。以N气田 2个小层的沉积微相图为例，建立了相应的沉积相模型，并与经典的基于样式的多点地质统计学建模方法（Simpat） 对比可以看出，SinGAN方法与训练图像刻画的沉积微相空间结构更相似，具有良好的应用前景。
Abstract
Sedimentary facies modeling is an important part of reservoir modeling and there are many methods to build sedi- mentary facies models. Traditional modeling methods need to describe the spatial structure information of variables by vari- ous parameters such as variogram and data patterns，and then reproduce the spatial structure in the realizations. With differ- ent strategies，the reservoir modeling based on Generative Adversarial Nets（GANs）learns a large number of images（mod- els）to generate the model possessing highly similar characteristics with the learning samples. Generative Adversarial Nets based on the single image（SinGAN）only need one image for training to generate highly similar images，improving the tradi- tional GANs. With the sedimentary microfacies diagram of two layers in N gas field as an example，the corresponding sedi- mentary facies model is built. Compared with the classical multiple-point geostatistics method Simpat，the SinGAN method obtains more similar spatial structure of sedimentary microfacies with that depicted by training images and has a good appli- cation prospect.
