高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于深度学习的地层沉积正演模拟代理模型构建与应用

刘彦锋 段太忠 龚伟 廉培庆 张文彪 黄渊

刘彦锋, 段太忠, 龚伟, 廉培庆, 张文彪, 黄渊. 基于深度学习的地层沉积正演模拟代理模型构建与应用[J]. 沉积学报, 2023, 41(3): 791-803. doi: 10.14027/j.issn.1000-0550.2021.152
引用本文: 刘彦锋, 段太忠, 龚伟, 廉培庆, 张文彪, 黄渊. 基于深度学习的地层沉积正演模拟代理模型构建与应用[J]. 沉积学报, 2023, 41(3): 791-803. doi: 10.14027/j.issn.1000-0550.2021.152
LIU YanFeng, DUAN TaiZhong, GONG Wei, LIAN PeiQing, ZHANG WenBiao, HUANG Yuan. Construction and Application of a Proxy Model for Stratigraphic Forward Modeling Based on Deep Learning[J]. Acta Sedimentologica Sinica, 2023, 41(3): 791-803. doi: 10.14027/j.issn.1000-0550.2021.152
Citation: LIU YanFeng, DUAN TaiZhong, GONG Wei, LIAN PeiQing, ZHANG WenBiao, HUANG Yuan. Construction and Application of a Proxy Model for Stratigraphic Forward Modeling Based on Deep Learning[J]. Acta Sedimentologica Sinica, 2023, 41(3): 791-803. doi: 10.14027/j.issn.1000-0550.2021.152

基于深度学习的地层沉积正演模拟代理模型构建与应用

doi: 10.14027/j.issn.1000-0550.2021.152
基金项目: 

中国石化科技部项目 P21038-3

中国石化科技部项目 P20077kxjgz

中科院先导A类项目 XDA14010204

详细信息
    作者简介:

    刘彦锋,男,1986年出生,博士,副研究员,开发地质学,E-mail: liuyf.syky@sinopec.com

  • 中图分类号: P628+.3

Construction and Application of a Proxy Model for Stratigraphic Forward Modeling Based on Deep Learning

Funds: 

SINOPEC Science and Technology Department P21038-3

SINOPEC Science and Technology Department P20077kxjgz

CAS Pilot Class A XDA14010204

  • 摘要: 地层沉积正演模拟方法能更真实地刻画地下地质体的分布规律,比传统的基于地质统计学的建模更有优势,但是条件化难度大,使其应用在实际油藏地质建模时面临较大挑战。地层沉积反演模拟提高了该方法的实用性,基于模拟结果与观测数据定量比较,地层沉积反演模拟在全局优化算法的驱动下不断修正地层沉积正演模拟输入参数,使模拟结果与观测数据吻合度达到最佳。由于反演系统优化参数多,非线性极强,收敛难度大,需要的迭代次数多,单次正演模拟耗时长,导致该方法效率较低。地层沉积反演模拟和深度学习算法中的生成对抗网络相结合,提出了构建地层沉积正演模拟代理模型的方法。以大量的碳酸盐岩地层沉积模拟的合成模型作为样本,通过神经网络训练,形成地层正演模拟器的代理模型,然后将其代入地层沉积反演模拟系统,避免了直接地层正演模拟的长耗时,加快了反演模拟的收敛速度。通过巴哈马西缘碳酸盐岩地层沉积模拟验证了该方法的可行性,采用学习后的生成网络进行沉积反演模拟,反演效率得到了大幅提升。尽管本文展示的是二维实例,也有望扩展应用在三维模型上。
  • 图  1  基于生成对抗网络的代理模型构建方法

    Figure  1.  Proxy model construction method based on generative adversarial network

    Fig.1

    图  2  训练集随机抽取的9个训练样本

    Figure  2.  Nine samples randomly selected from training set

    Fig.2

    图  3  生成对抗网络的结构(生成器和判别器)

    Figure  3.  Structure of generative adversarial network (generator and discriminator)

    Fig.3

    图  4  不同训练轮次的网络在测试集某一个样本上的模拟效果

    Figure  4.  Simulation results of networks different training rounds on a sample of the test set

    Fig.4

    图  5  测试集中随机抽取的6组神经网络生成器模拟(假)与传统正演模拟(真)比较

    Figure  5.  Comparison between generator simulation (false) and conventional forward simulation (true) on randomly selected 6 samples from test set

    Fig.5

    图  6  合成样本反演结果

    Figure  6.  Inverse result on synthetic sample

    Fig.6

    图  7  真实数据反演模拟结果

    Figure  7.  Inverse result on actual observed data

    Fig.7

    图  8  基于正演模拟器的反演模拟结果

    Figure  8.  Inverse result by conventional forward simulation

    Fig.8

    表  1  沉积反演模拟参数区间

    序号参数名字最小值最大值序号参数名字最小值最大值
    1InitTopoElev1-267.00-67.0818PDY130 000.0099 972.00
    2InitTopoElev2-341.00-141.0819KDX1200.00599.84
    3InitTopoElev3-467.50-267.0820KDY1200.00599.84
    4InitTopoElev4-632.00-450.0721Aprod11.305.00
    5AmpSeaL040.0079.9822Kprod10.040.12
    6PerSeaL0230.00699.8123Wprod115.0044.99
    7PhaseSeaL0-0.501.0024kEng10.501.50
    8AmpSeaL110.0049.9825bEng10.200.60
    9PerSeaL110.0099.9626pPr10.100.90
    10PhaseSeaL1-1.001.0027kV110.0039.99
    11SeaLevelDatum-50.0039.9628kW10.020.20
    12SubsRate1-0.30029kR10.010.06
    13SubsRate2-0.30-0.0530kH10.050.17
    14SubsRate3-0.55-0.0131WaveBase120.0059.98
    15SubsRate4-0.65-0.0532kL100
    16SubsRate5-0.87-0.0533DFlex1E+244.998E+26
    17PDX130 000.0079 980.0034rkM100.10
    下载: 导出CSV

    表  附表2  测试9随机抽取的6个样本

    参数名称样本A样本B样本C样本D样本E样本F
    InitTopoElev1-235.872-211.007-240.939-226.946-216.507-134.497
    InitTopoElev2-294.878-317.912-305.656-283.186-314.672-294.664
    InitTopoElev3-440.537-404.218-405.011-363.598-402.461-429.075
    InitTopoElev4-611.515-604.388-572.029-611.646-603.645-602.346
    AmpSeaL045.907 744.462 854.435 144.457 345.539 245.516 2
    PerSeaL0302.895297.319283.025327.798305.036285.534
    PhaseSeaL0-0.313 27-0.268 81-0.295 00-0.265 74-0.125 55-0.318 12
    AmpSeaL131.055 716.873 314.874 916.619 724.722 816.481 8
    PerSeaL122.732 124.420 021.763 342.429 351.685 638.696 4
    PhaseSeaL1-0.751 93-0.711 99-0.647 79-0.664 00-0.199 81-0.670 46
    SeaLevelDatum-23.082 0-29.752 5-33.811 0-35.430 4-37.047 00.349 351
    SubsRate1-0.236 21-0.122 30-0.224 62-0.231 28-0.232 31-0.252 31
    SubsRate2-0.248 49-0.230 61-0.260 81-0.253 01-0.245 30-0.262 11
    SubsRate3-0.450 76-0.489 44-0.455 16-0.483 23-0.454 18-0.485 23
    SubsRate4-0.572 49-0.528 56-0.519 46-0.463 54-0.536 64-0.538 34
    SubsRate5-0.802 85-0.600 26-0.794 73-0.775 51-0.636 48-0.770 71
    PDX139 367.251 846.150 318.148 834.040 938.037 031.7
    PDY149 483.238 700.841 724.946 108.640 312.340 762.6
    KDX1303.702319.104317.759346.898267.878261.639
    KDY1262.673258.522300.591355.141252.191270.105
    Aprod12.488 891.720 792.116 181.935 731.895 371.905 39
    Kprod10.078 6540.057 4640.058 0980.071 6170.050 4010.054 544
    Wprod118.519 719.082 221.298 227.134 619.764 018.834 8
    kEng10.682 2370.835 3050.670 8500.660 1230.678 8570.826 163
    bEng10.348 7710.268 2140.268 6900.302 3730.271 5500.259 995
    pPr10.345 8560.399 6250.372 7460.233 5660.392 5120.267 308
    kV116.119 317.201 619.910 214.708 316.862 415.730 3
    kW10.063 5140.042 1870.049 6630.048 8180.040 6130.063 422
    kR10.016 3350.018 6690.019 1810.035 4210.023 0360.015 796
    kH10.074 9750.070 5310.065 3820.063 5100.086 6670.063 957
    WaveBase126.208 436.326 034.947 224.898 526.646 733.491 0
    kL10.000 1050.000 1388.20E-054.97E-050.000 1244.43E-05
    DFlex1.34E+261.17E+269.29E+256.42E+258.00E+251.58E+26
    rkM10.039 6950.012 1370.030 7000.020 0770.014 1420.022 487
    下载: 导出CSV

    表  附表1  训练集随机抽取的9个样本实例的参数

    参数名称样本A样本B样本C样本D样本E样本F样本G样本H样本I样本J
    InitTopoElev1-241.740-239.795-240.203-189.967-210.110-237.965-212.366-176.158-213.844-228.109
    InitTopoElev2-304.687-308.444-288.532-302.394-315.569-316.927-311.806-313.322-299.977-290.627
    InitTopoElev3-443.803-420.896-440.692-372.690-426.142-429.549-419.037-445.090-444.184-438.012
    InitTopoElev4-608.549-594.739-606.268-608.710-611.325-600.513-599.063-603.759-597.293-608.271
    AmpSeaL047.993 049.966 849.034 452.207 746.511 847.745 151.338 563.144 546.858 550.374 2
    PerSeaL0415.966333.573350.251319.018288.357427.918302.105338.303314.502308.967
    PhaseSeaL0-0.054 27-0.013 71-0.303 06-0.270 38-0.078 13-0.331 04-0.128 830.177 19-0.295 01-0.309 59
    AmpSeaL119.470 515.821 215.208 014.819 318.243 819.648 729.210 715.980 515.315 417.789 8
    PerSeaL136.966 050.087 929.274 725.920 328.986 729.643 435.675 930.717 837.475 227.578 2
    PhaseSeaL1-0.440 86-0.645 22-0.632 37-0.568 44-0.768 05-0.692 32-0.678 35-0.655 46-0.690 46-0.746 04
    SeaLevelDatum-39.560 5-38.872 9-37.119 8-34.170 6-34.217 7-38.279 3-39.223 3-37.920 3-39.168 8-27.961 8
    SubsRate1-0.243 89-0.247 15-0.253 77-0.258 66-0.258 47-0.228 18-0.243 00-0.236 90-0.242 21-0.195 93
    SubsRate2-0.260 98-0.219 83-0.251 05-0.240 68-0.228 05-0.246 30-0.163 37-0.237 58-0.260 05-0.249 80
    SubsRate3-0.370 92-0.417 32-0.476 49-0.453 03-0.408 79-0.467 76-0.427 09-0.454 19-0.458 59-0.487 05
    SubsRate4-0.563 53-0.496 77-0.527 97-0.487 82-0.444 86-0.581 18-0.539 59-0.519 66-0.311 60-0.544 13
    SubsRate5-0.752 94-0.781 76-0.671 17-0.763 93-0.687 47-0.667 36-0.593 09-0.723 55-0.767 86-0.653 68
    PDX137 175.743 489.039 140.935 975.737 643.636 044.036 179.240 122.340 944.552 077.4
    PDY157 366.054 775.444 088.448 277.344 436.440 488.039 890.340 007.238 728.038 879.8
    KDX1275.378266.399245.489256.212251.925276.658245.710255.417396.307275.302
    KDY1262.749249.093278.053312.542309.070284.470270.730259.157435.309438.423
    Aprod12.617 082.159 542.729 962.279 542.192 181.835 661.766 423.566 482.022 241.734 97
    Kprod10.072 5770.049 2170.049 6030.051 5100.052 6950.053 1920.049 7300.059 0710.049 7130.049 448
    Wprod118.598 819.906 022.985 822.905 626.301 742.079 823.508 324.068 423.353 318.509 9
    kEng10.671 6430.738 1030.727 7840.666 1550.819 8350.799 8260.620 3400.714 9960.923 4370.725 410
    bEng10.271 1130.282 1440.299 4230.336 5460.261 0210.317 6330.246 5010.282 420.257 8690.391 195
    pPr10.206 7960.196 9490.268 6950.204 3560.305 7930.231 9750.261 1850.364 4550.210 0160.258 421
    kV128.296 315.356 414.039 513.641 217.244 115.175 515.189 615.86520.368 916.541 7
    kW10.058 6210.051 8510.061 6310.044 9440.053 1900.125 9660.043 5920.060 9930.050 8900.041 966
    kR10.022 4550.024 6010.025 0260.034 4260.015 8120.016 9740.016 9300.017 060.017 7120.016 261
    kH10.063 8970.065 0080.120 7050.069 2720.063 4450.077 0700.137 2120.083 2970.077 0880.076 075
    WaveBase128.713 327.629 230.370 325.878 631.746 928.013 127.847 924.561930.966 831.201 4
    kL10.000 1074.65E-054.53E-056.79E-055.09E-054.26E-056.15E-054.73E-055.87E-054.86E-05
    DFlex1.19E+267.95E+257.63E+256.58E+251.13E+265.80E+256.33E+255.81E+257.15E+251.94E+26
    rkM10.016 5680.027 5960.031 4030.012 2240.031 8960.017 9190.015 1680.018 8520.014 1580.034 389
    下载: 导出CSV
  • [1] Tetzlaff D M, Harbaugh J W. Simulating clastic sedimentation[M]. New York: Springer, 1989.
    [2] Burgess P M, Allen P A. A forward-modelling analysis of the controls on sequence stratigraphical geometries[J]. Geological Society, London, Special Publications, 1996, 103(1): 9-24.
    [3] Lessenger M A, Cross T A. A stratigraphic inverse simulation model[C]//Proceedings of annual meeting of the American association of petroleum geologists. Dallas: AAPG, 1991.
    [4] Cross T A, Lessenger M A. Construction and application of a stratigraphic inverse model[M]//Harbaugh J W, Watney W L, Rankey E C, et al. Numerical experiments in stratigraphy: Recent advances in stratigraphic and sedimentologic computer simulations.Oklahoma: SEPM, 1999: 69-83.
    [5] Bowman S A, Vail P R. Interpreting the stratigraphy of the Baltimore canyon section, offshore New Jersey with PHil, a stratigraphic simulator[M]//Harbaugh J W, Watney W L, Rankey E C, et al. Numerical experiments in stratigraphy: Recent advances in stratigraphic and sedimentologic computer simulations. Oklahoma: SEPM, 1999: 117-138.
    [6] Tetzlaff D M. Input uncertainty and conditioning in siliciclastic process modelling[J]. Geological Society, London, Special Publications, 2004, 239(1): 95-109.
    [7] Imhof M G, Sharma A K. Quantitative seismostratigraphic inversion of a prograding delta from seismic data[J]. Marine and Petroleum Geology, 2006, 23(7): 735-744.
    [8] Charvin K, Gallagher K, Hampson G L, et al. A Bayesian approach to inverse modelling of stratigraphy, Part 1: Method[J]. Basin Research, 2009, 21(1): 5-25.
    [9] Falivene O, Frascati A, Gesbert S, et al. Automatic calibration of stratigraphic forward models for predicting reservoir presence in exploration[J]. AAPG Bulletin, 2014, 98(9): 1811-1835.
    [10] Sacchi Q, Weltje G J, Verga F. Towards process-based geological reservoir modelling: Obtaining basin-scale constraints from seismic and well data[J]. Marine and Petroleum Geology, 2015, 61: 56-68.
    [11] Sharma A K. Quantitative stratigraphic inversion[D]. Virginia: Virginia Polytechnic Institute and State University, 2006.
    [12] Chan S, Elsheikh A H. Parametrization of stochastic inputs using generative adversarial networks with application in geology[J]. Frontiers in Water, 2020, 2: 5.
    [13] Ren X L, Li X Y, Ren K J, et al. Deep learning-based weather prediction: A survey[J]. Big Data Research, 2021, 23: 100178.
    [14] Gong W, Duan Q Y. An adaptive surrogate modeling-based sampling strategy for parameter optimization and distribution estimation (ASMO-PODE)[J]. Environmental Modelling & Software, 2017, 95: 61-75.
    [15] Pomar L. Ecological control of sedimentary accommodation: Evolution from a carbonate ramp to rimmed shelf, Upper Miocene, Balearic Islands[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2001, 175(1/2/3/4): 249-272.
    [16] Pomar L. Carbonate systems[M]//Scarselli N, Adam J, Chiarella D, et al. Regional geology and tectonics. Amsterdam: Elsevier, 2020: 235-311.
    [17] Pomar L, Haq B U. Decoding depositional sequences in carbonate systems: Concepts vs experience[J]. Global and Planetary Change, 2016, 146: 190-225.
    [18] 段太忠,王光付,廉培庆,等. 油气藏定量地质建模方法与应用[M]. 北京:石油工业出版社,2019.

    Duan Taizhong, Wang Guangfu, Lian Peiqing, et al. Quantitative geological modeling method and application of oil and gas reservoir [M].Beijing: Petroleum Industry Press, 2019.
    [19] Duan Q Y, Sorooshian S, Gupta V K. Optimal use of the SCE-UA global optimization method for calibrating watershed models[J]. Journal of Hydrology, 1994, 158(3/4): 265-284.
    [20] Goodfellow I J, Pouget-Abadie J, Mirza M, et al. Generative adversarial nets[C]//Proceedings of the 27th international conference on neural information processing systems. Montreal, Canada: MIT Press, 2014.
    [21] 刘彦锋,张文彪,段太忠,等. 深度学习油气藏地质建模研究进展[J]. 地质科技通报,2021,40(4):235-241.

    Liu Yanfeng, Zhang Wenbiao, Duan Taizhong, et al. Progress of deep learning in oil and gas reservoir geological modeling[J]. Bulletin of Geological Science and Technology, 2021, 40(4): 235-241.
    [22] Betzler C, Reijmer J J G, Bernet K, et al. Sedimentary patterns and geometries of the Bahamian outer carbonate ramp (Miocene⁃Lower Pliocene, Great Bahama Bank)[J]. Sedimentology, 1999, 46(6): 1127-1143.
    [23] Ginsburg R N. Subsurface geology of a prograding carbonate platform margin, great Bahama bank: Results of the bahamas drilling project[M]. Oklahoma:SEPM, 2001.
    [24] Ginsburg R N. The Bahamas drilling project: Background and acquisition of cores and logs[M]//Ginsburg R N. Subsurface geology of a prograding carbonate platform margin, great Bahama bank: Results of the Bahamas drilling project. Oklahoma:SEPM, 2001.
    [25] Eberli G P, Kendall C G S C, Moore P, et al. Testing a seismic interpretation of Great Bahama Bank with a computer simulation[J]. AAPG Bulletin, 1994, 78(6): 981-1004.
    [26] Eberli G P, Anselmetti F S, Kroon D, et al. The chronostratigraphic significance of seismic reflections along the Bahamas Transect[J]. Marine Geology, 2002, 185(1/2): 1-17.
    [27] Qi G J. Loss-sensitive generative adversarial networks on lipschitz densities[J]. International Journal of Computer Vision, 2020, 128(5): 1118-1140.
  • [1] 齐昆, 龚承林, 金振奎.  块体搬运沉积在大陆边缘深水区的发育——来自晚第四纪珠江陆缘地层正演模拟的启示 . 沉积学报, 2024, 42(4): 1309-1325. doi: 10.14027/j.issn.1000-0550.2023.031
    [2] 邵俊哲, 李阳, 吴楠, 冯文杰, 魏薇.  基于机器学习的海陆过渡环境中陆源有机质分布预测模型——从沉积模拟实验到地质应用【水槽沉积模拟实验专辑】 . 沉积学报, 2024, (): -. doi: 10.14027/j.issn.1000-0550.2024.056
    [3] 胡勇, 高小洋, 何文祥, 李顺利, 朱建斌, 司锦, 陆雨诗.  利用条件生成对抗网络建立曲流河地质模型 . 沉积学报, 2024, 42(1): 201-218. doi: 10.14027/j.issn.1000-0550.2022.046
    [4] 罗歆, 闫建平, 王军, 耿斌, 王敏, 钟广海, 张帆, 李志鹏, 高松洋.  基于FMI图像深度学习的砂砾岩体沉积微相识别方法 . 沉积学报, 2023, 41(4): 1138-1152. doi: 10.14027/j.issn.1000-0550.2021.148
    [5] 马艺萍, 王荣华, 戴霜, 马晓军.  北祁连北大河沉积物碎屑组成及物源正演分析 . 沉积学报, 2022, 40(6): 1525-1541. doi: 10.14027/j.issn.1000-0550.2022.086
    [6] 谈明轩, 朱筱敏, 张自力, 孙浩南.  构造掀斜主导的断陷湖盆缓坡层序“源—汇”正演模拟定量研究 . 沉积学报, 2022, 40(6): 1481-1493. doi: 10.14027/j.issn.1000-0550.2022.061
    [7] 李红斌, 王贵文, 王松, 庞小娇, 刘士琛, 包萌, 彭寿昌, 赖锦.  基于Kohonen神经网络的页岩油岩相测井识别方法 . 沉积学报, 2022, 40(3): 626-640. doi: 10.14027/j.issn.1000-0550.2021.045
    [8] 高志勇, 马建英, 崔京钢, 冯佳睿, 周川闽, 吴昊.  埋藏(机械)压实-侧向挤压地质过程下深层储层孔隙演化与预测模型 . 沉积学报, 2018, 36(1): 176-187. doi: 10.3969/j.issn.1000-0550.2018.019
    [9] 金中国, 邹林, 张力, 郑明泓, 韩英, 谷静, 谢樨, 陈兴龙, 王洪.  贵州务—正—道地区铝土矿成矿模式与找矿模型 . 沉积学报, 2018, 36(5): 914-926. doi: 10.14027/j.issn.1000-0550.2018.081
    [10] 杨伟利, 王毅, 张众琛, 梁晨, 李京昌.  塔里木盆地顺南地区中-下奥陶统正演沉积数值模拟 . 沉积学报, 2017, 35(1): 117-123. doi: 10.14027/j.cnki.cjxb.2017.01.012
    [11] 朱锐.  正弦曲线型地层基准面旋回一维定量理论模型的建立 . 沉积学报, 2009, 27(1): 64-69.
    [12] 赖生华.  断陷盆地沉积体系研究新思路:从古地貌、岩性变化、水体深度到沉积体系 . 沉积学报, 2007, 25(5): 663-670.
    [13] 石广仁.  蒙皂石转化伊利石的数值模拟——溶解沉淀模型与化学动力学模型 . 沉积学报, 2007, 25(5): 693-700.
    [14] 唐为清, 郭荣坤, 王忠东, 王红, 罗安银, 毋学平.  沉积微相测井资料神经网络判别方法研究 . 沉积学报, 2001, 19(4): 581-585.
    [15] 马世忠, 杨清彦.  曲流点坝沉积模式、三维构形及其非均质模型 . 沉积学报, 2000, 18(2): 241-247.
    [16] 覃建雄, 陈洪德, 田景春, 李余生, 杨作升.  川滇黔桂地区二叠纪不同成因盆地层序地层模型 . 沉积学报, 1999, 17(2): 183-187.
    [17] 卫延召, 冯江, 宋海明, 王多云.  陕甘宁盆地鄂托克前旗地区奥陶系风化壳的地震正演模型 . 沉积学报, 1999, 17(S1): 734-741.
    [18] 张福明, 李洪奇, 邵才瑞, 白彦彬, 张巨兴, 曹峰, 王占忠, 焦翠华.  神经网络和分形几何方法在识别测井沉积微相中的应用 . 沉积学报, 1997, 15(3): 62-66.
    [19] 王硕儒, 范德江, 汪丙柱.  岩相识别的神经网络计算 . 沉积学报, 1996, 14(4): 154-160.
    [20] 冉启全, 李士伦, 顾小芸.  神经网络模式识别沉积微相 . 沉积学报, 1996, 14(S1): 186-191.
  • 加载中
图(8) / 表 (3)
计量
  • 文章访问数:  367
  • HTML全文浏览量:  74
  • PDF下载量:  91
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-09-10
  • 修回日期:  2021-11-12
  • 录用日期:  2021-12-22
  • 网络出版日期:  2021-12-22
  • 刊出日期:  2023-06-10

目录

    基于深度学习的地层沉积正演模拟代理模型构建与应用

    doi: 10.14027/j.issn.1000-0550.2021.152
      基金项目:

      中国石化科技部项目 P21038-3

      中国石化科技部项目 P20077kxjgz

      中科院先导A类项目 XDA14010204

      作者简介:

      刘彦锋,男,1986年出生,博士,副研究员,开发地质学,E-mail: liuyf.syky@sinopec.com

    • 中图分类号: P628+.3

    摘要: 地层沉积正演模拟方法能更真实地刻画地下地质体的分布规律,比传统的基于地质统计学的建模更有优势,但是条件化难度大,使其应用在实际油藏地质建模时面临较大挑战。地层沉积反演模拟提高了该方法的实用性,基于模拟结果与观测数据定量比较,地层沉积反演模拟在全局优化算法的驱动下不断修正地层沉积正演模拟输入参数,使模拟结果与观测数据吻合度达到最佳。由于反演系统优化参数多,非线性极强,收敛难度大,需要的迭代次数多,单次正演模拟耗时长,导致该方法效率较低。地层沉积反演模拟和深度学习算法中的生成对抗网络相结合,提出了构建地层沉积正演模拟代理模型的方法。以大量的碳酸盐岩地层沉积模拟的合成模型作为样本,通过神经网络训练,形成地层正演模拟器的代理模型,然后将其代入地层沉积反演模拟系统,避免了直接地层正演模拟的长耗时,加快了反演模拟的收敛速度。通过巴哈马西缘碳酸盐岩地层沉积模拟验证了该方法的可行性,采用学习后的生成网络进行沉积反演模拟,反演效率得到了大幅提升。尽管本文展示的是二维实例,也有望扩展应用在三维模型上。

    English Abstract

    刘彦锋, 段太忠, 龚伟, 廉培庆, 张文彪, 黄渊. 基于深度学习的地层沉积正演模拟代理模型构建与应用[J]. 沉积学报, 2023, 41(3): 791-803. doi: 10.14027/j.issn.1000-0550.2021.152
    引用本文: 刘彦锋, 段太忠, 龚伟, 廉培庆, 张文彪, 黄渊. 基于深度学习的地层沉积正演模拟代理模型构建与应用[J]. 沉积学报, 2023, 41(3): 791-803. doi: 10.14027/j.issn.1000-0550.2021.152
    LIU YanFeng, DUAN TaiZhong, GONG Wei, LIAN PeiQing, ZHANG WenBiao, HUANG Yuan. Construction and Application of a Proxy Model for Stratigraphic Forward Modeling Based on Deep Learning[J]. Acta Sedimentologica Sinica, 2023, 41(3): 791-803. doi: 10.14027/j.issn.1000-0550.2021.152
    Citation: LIU YanFeng, DUAN TaiZhong, GONG Wei, LIAN PeiQing, ZHANG WenBiao, HUANG Yuan. Construction and Application of a Proxy Model for Stratigraphic Forward Modeling Based on Deep Learning[J]. Acta Sedimentologica Sinica, 2023, 41(3): 791-803. doi: 10.14027/j.issn.1000-0550.2021.152
      • 地层沉积过程正演模拟比传统的基于统计学的建模能更真实地刻画地下地质体的分布规律,但条件化方面仍需改善。地层沉积过程正演模拟以质量守恒为基础,考虑了沉积物的供给、产生、搬运、剥蚀和再沉积等地质作用,模拟结果符合地质规律[1]。正演模拟的输入参数有初始地形、海平面曲线、沉积物搬运系数、沉积物供给曲线、基底沉降曲线、碳酸盐岩生产的相关系数,以及这些系数随时间和空间的变化等[2]。由于这些输入参数很难通过现有的观测手段获得,使得模型存在较大不确定性,与井数据和地震解释数据吻合难度大,导致该建模技术在实际油藏地质建模中应用推广缓慢。

        地层沉积反演模拟技术提高了地层沉积过程模拟的实用性。Lessenger et al.[3],Cross et al.[4]提出了地层沉积反演模拟的基本框架,认为地层沉积反演模拟包括地层沉积正演、模拟结果和观测数据之间的比较,以及自动调整正演参数的最优化算法。基本思路是:以有限的、能反映地层分布特征的观测数据,比如井上钻遇的沉积微相、地震层位解释等,反推地层沉积正演模拟的输入参数,实现模拟结果与观测数据吻合,然后进行地层空间展布的预测。地层沉积模拟的基本过程是从一组初始的过程参数开始执行正演模拟,将模拟结果与观测数据比较,优化算法按照一定的策略调整过程参数,再次执行正演模拟,直到模拟结果与观测数据的吻合度达到门槛值。前人在这些方面做了很多研究,涵盖正演模拟技术、比较技术和优化算法等[511]

        在反演系统中,通常正演模拟的计算量很大,是最耗时的计算部分。代理模型是工程问题中常用的一个优化方法。当实际问题计算量很大、不容易求解时,可以使用计算量较小、求解速度快的代理模型替代原模型,加速优化过程,比如多项式法、二次响应曲面和克里金方法等。尽管代理模型是复杂系统反演的常见处理方式[1213],但针对沉积模拟代理模型的相关研究未见报道。

        与油气藏渗流、地下水流动、数值天气预报等模型类似,地层沉积反演模拟也是一个复杂的非线性系统,本质上求解这些复杂系统都是大规模高度非线性问题的优化问题,需要借助高效的优化算法反复进行迭代求解,其中单次正演模拟求解耗时长、非线性强是影响问题快速求解的主要瓶颈。对于这类问题,代理模型是用另外一种更快的数学方法产生与原模型误差最小的结果,以达到快速迭代。目前,已经有一些替代正演模拟的方法,这些模型替代非线性较弱的模型时效果较好,对复杂性强的模型仍无法替代[14],深度神经网络具有很强的模式表达能力,为代理模型构建提供了新的途径,也为快速的地层沉积反演模拟提供了可能。

        最近几年发展迅速的生成对抗人工神经网络方法具有很强的模型生成能力,已经用于模拟复杂地质模型的生成和油气藏数值模拟中替代模型的建立。神经网络主要分为判别网络和生成网络,判别网络建立高维参数向低维参数的映射,常用于分类、聚类和判别问题等;生成网络建立低维参数向高维参数的映射,常用于回归和模型生成问题。复杂三维地质建模可以通过生成网络实现。常见的生成网络有贝叶斯网络、变分自编码器、玻尔兹曼机、生成对抗网络等,其中生成对抗网络的应用最广泛。

        相比于其他方法,生成对抗网络具有更强的模式生成能力,可以生成更逼真的模型,但存在样本需求量的问题,不适合样本少的情况。对于地下地质类问题,难以获取真实的三维模型作为样本,通过正演模拟方法合成样本是深度学习技术在该领域落地应用的主要方式。在沉积反演模拟系统中,样本通过正演模拟器生成,不存在样本不足的问题。使用生成对抗网络,有望提高反演的收敛速度,解决目前的效率问题。

        本文在地层沉积反演模拟基本框架、生成对抗网络基本原理介绍的基础上,提出了基于生成对抗网络的地层沉积代理模型构建方法,并通过实例研究说明该方法的可行性,最后讨论了该方法存在的局限性和未来发展的方向。

      • 地层沉积反演模拟包括生成三维模型的正演模拟器、模拟结果与观测数据比较的误差计算器,调整正演模拟输入参数的全局优化器。优化器是整个系统的核心驱动力,它根据误差计算器的响应不断调整正演模拟输入的参数,使模拟结果与观测数据吻合度最高。

        地层正演模型的输入参数考虑了边界条件和初始条件,它们描述了可容空间、沉积物剥蚀、供给、生产、搬运、堆积以及压实作用等。输出结果是模拟的地层和一系列的古环境条件,比如整个演化历程的古地貌、沉积间断等。区别于碎屑岩地层,碳酸盐岩沉积体系的形成、演化和消亡,除了受构造活动、海平面变化、气候条件、海洋环境和水动力条件等多种作用于碎屑岩沉积体系的共同因素的直接或间接的控制外,还受体系中生物与生态因素的重要控制[1517]。此次采用了自主研发的碳酸盐岩地层沉积正演模拟方法,它考虑了最新的碳酸盐岩工厂、生态可容空间和层序地层学等基本原理,建立了基于环境能量和生物碳酸盐产能的地层沉积正演模拟方法和软件系统,可以模拟台缘带、缓坡、孤立台地以及复杂地形下的碳酸盐岩沉积体系等[18],限于篇幅问题,具体的正演模拟方法不再展开。

        模拟结果与观测数据的比较方法依赖于观测数据的多少及类型,往往是较少的观测资料,比如少量的几口钻井资料、分辨率较低的地震解释资料等。为了降低地震资料引起的不确定性,本文选择井上的解释的沉积相与实际模拟的过井处沉积相的均方根误差作为反演的目标函数。对于复杂非线性且不可求导的非线性系统,需要采用启发式优化算法,比如模拟退火、差异演化、粒子群算法等,通过多参数下的优化算法收敛性能对比,基于数值梯度的复合多重形优化算法(SCE-UA)在多参数反演方面具备更多优势[19]

        由于整个地层沉积反演模拟系统非线性极强,收敛难度大,需要的迭代次数很多,导致效率较低。提出了基于生成对抗网络沉积模拟代理模型构建方法,替代其中的正演模拟过程,充分利用神经网络全程可微分的性质,提高整个反演系统的收敛性能。

      • 生成对抗神经网络为地质建模提供了新方法。随着深度学习等人工智能技术的发展,地质建模技术也朝着智能化方向快速发展。由于生成神经网络与地质建模之间的高度相似性,这项技术在地质建模领域快速落地应用。在多层神经网络框架下,地质建模是一个生成问题,油气藏地质建模通常是输入稀疏的、高精度的井点数据,在稠密的、相对低精度的地震数据约束条件下建立高精度三维网格模型,可以视为生成神经网络问题,其中由Goodfellow et al.[20]提出的生成对抗神经网络(Generative Adversarial Networks,GAN)是解决这类问题的主要方式。

        生成对抗网络是一种适用于任意形式复杂分布的无监督对抗式深度学习方法。它一般由两个模块构成:生成模型(Generative Model,简写为G)和判别模型(Discriminative Model,简写为D)组成,生成模型也称为生成器,判别模型也称为判别器。生成模型G根据一系列特征,可以是数值或者字符串等,来生成一组数据,生成的数据可以是图片、影像、声音等,也可以是地质模型。生成的数据和真实数据相比要尽可能相似。判别模型D判断生成模型输出的结果是否和真实数据相似,其目标是将生成数据和真实数据尽可能分开。G和D两个模型在对抗性训练中共同进步,G生成的数据越来越逼真,D的辨别能力越来越强。GAN训练终止的条件和一般的网络有所区别,并非以损失函数达到某一阈值作为终止条件,而是人工判断生成模型的结果逼真程度是否还在持续改善。如果持续多轮迭代不再改善则终止训练。

        生成对抗网络是目前深度学习地质建模中主要采用的神经网络[21],为替代模型建立提供了新途径。对于标准的生成对抗网络,需要输入一个低维的随机数序列得到高维的网络层,与非条件地质建模过程较为相似,即通过变差函数、训练图像或目标体长宽高等输入参数得到一系列三维模型实现,这些实现不考虑与井点数据或地震约束数据的吻合情况,但能够体现较为合理的地质模式。生成对抗网络深度学习体现了类似过程,通过学习地质模式,利用训练后的生成器可以得到满足地质认识的模拟结果。由于神经网络固有的线性化性质,按照这种方式建立的代理模型,具备更好的反向求导和快速收敛特征。

      • 深度学习模型的训练需要大量的样本,不断扰动沉积正演模拟器的输入参数,可以获得大规模沉积模拟模型作为深度学习的样本。沉积模拟的输入包括初始地形、碳酸盐岩产率、沉积物搬运、海平面曲线、构造沉降等相关参数,改变任意的输入参数都可能产生不同的结果。采用生成对抗神经网络框架进行训练,其中生成器接收与沉积模拟一样的输入参数,产生二维或三维模拟结果,判别器输入二维或三维模型,输出真或假,其中训练后的生成器即为最终需要的沉积模拟代理模型(图1)。具体构建步骤如下。

        图  1  基于生成对抗网络的代理模型构建方法

        Figure 1.  Proxy model construction method based on generative adversarial network

        步骤一:分析具体的地质问题,搭建反演模拟基本框架。根据钻井和地震解释的地层特征,结合区域地质特征,确定目标地区可能的地层特征和沉积相展布规律。通过反复的正演模拟测试,确定初始地形和构造沉降的基本特征,以及沉积物搬运、碳酸盐岩产率、水体能量相关参数的数值区间,为样本生成提供基础。

        步骤二:基于沉积正演模拟器生产大规模样本。选择对模拟结果敏感的参数,且尽量涵盖更大的解空间,通过超立方采样的方式,产生不同的参数组合,开展正演模拟,建立样本库。一个样本包括输入的参数和正演模拟的结果,其中模拟结果是二维或三维的相模型,通过离散数值(比如a1、b1、b2、c1、c2、c3,c4共7种相类型)的矩阵表达。随机选择其中的80%作为训练集,20%作为测试集。

        步骤三:基于生成对抗网络的代理模型构建。根据样本数据特征及生成对抗网络的要求设计生成网络和判别网络的结构,为了更好地捕捉地质模式并增加效率,采用全卷积式网络,即网络的核心处理单元均采用卷积和转置卷积层,每个卷积层和转置卷积层配套标准化层和非线性层,不采用全连接层。在训练集上对网络进行训练,在测试集上对训练后的生成器和判别器进行测试,如果生成器或判别器的效果不再改进,则停止训练。如果训练后网络无法生成合理的地质模型则调整网络结构,甚至回到步骤二生成更多的样本,直到生成器满足要求。

        步骤四:基于深度学习代理模型的地层沉积反演模拟。将步骤三训练得到的代理模型(即对抗神经网络中的生成器)带入地层沉积反演模拟系统,以计算结果与观测数据的误差最小为目标,采用全局优化器不断调整代理模型的输入参数,直到得到最优结果。将反演得到的最优参数带入代理模型,得到最终的地质模型。

        地层沉积正演模拟器有很强的灵活性。理论上,可以构建代表任意参数个数和参数区间的代理模型,但它需要样本太多、网络的深度也太大,训练难度大,甚至不可能。但是,针对一定的情况,参数空间会大大减小,需要的样本也会减少。针对常见的几种碳酸盐岩沉积环境类型,比如斜坡、镶边台地、孤立台地等类型,可以训练出有针对性、同时又不失一般性的神经网络代理模型。下面将以巴哈马滩碳酸盐岩沉积为例,介绍地层沉积模拟代理模型构建方法及其在反演模拟系统中的应用。

      • 巴哈马滩是国内外研究碳酸盐岩地层沉积的热点地区,在其西缘斜坡有几口科学钻探井和实测的地震剖面,前人在该地区开展了大量研究[2224]。台地整体形态、初始地形、构造沉降速度、沉积演化过程等都有比较清楚的认识,也有沉积过程定量模拟的相关研究,但都侧重于沉积过程的控制因素分析,缺少精细的沉积过程建模研究。结合地质分析和前人研究成果[2526],用地层沉积正演的方法检测地区地震解释剖面,初步获取该地区碳酸盐岩海平面变化、构造沉降、碳酸盐岩产率等参数,搭建了该地区的地层沉积过程模拟基本框架。沉积模拟范围为32 km×32 km,时间跨度为5.3 Ma,平面网格数为129×129,总时间步数为100。

        通过正演模拟的参数尝试,获取了主要参数的分布区间,为生成大量样本和代理模型构建奠定了基础。地层沉积正演模拟要输入的参数有初始地形、构造沉降、海平面变化曲线、碳酸盐岩产率相关参数、沉积物搬运参数等34个。根据对该地区前期的资料分析认为初始地形和构造沉降的不确定性较低,而海平面曲线、碳酸盐岩产率和沉积物搬运的相关参数不确定性较大,建立反演模拟系统时确定性较强的参数给了较小的变化区间,不确定性较大的参数给了较大的变化区间(表1)。这34个参数中,第1~4个参数与初始地形有关(沿着斜坡方向均匀分布的个4点,用来插值生成初始地形),第5~11个参数与海平面曲线有关(分别是表示三级和四级海平面旋回的正弦函数的振幅、周期和相位,以及海平面基准值),第12~16个参数与构造沉降有关(沿着斜坡方向均匀分布的4个点,用来插值生成构造沉降面),第17~20个参数与沉积物搬运有关(分别是势能和动能在X和Y方向的搬运系数),第21~26个参数与产率有关(分别是势能产率最大幅度、势能产率递减系数、透光带厚度,动能产率系数,动能产率基准值,势能产率权重),第27~32个参数与水体动能有关(动能幅度、波浪能下降系数、地形消浪能系数、浪基面、风能系数),第33个参数为地层挠曲系数,第34个参数为生物群落系数,各个参数更详细的介绍可以参考文献[18]。

        表 1  沉积反演模拟参数区间

        序号参数名字最小值最大值序号参数名字最小值最大值
        1InitTopoElev1-267.00-67.0818PDY130 000.0099 972.00
        2InitTopoElev2-341.00-141.0819KDX1200.00599.84
        3InitTopoElev3-467.50-267.0820KDY1200.00599.84
        4InitTopoElev4-632.00-450.0721Aprod11.305.00
        5AmpSeaL040.0079.9822Kprod10.040.12
        6PerSeaL0230.00699.8123Wprod115.0044.99
        7PhaseSeaL0-0.501.0024kEng10.501.50
        8AmpSeaL110.0049.9825bEng10.200.60
        9PerSeaL110.0099.9626pPr10.100.90
        10PhaseSeaL1-1.001.0027kV110.0039.99
        11SeaLevelDatum-50.0039.9628kW10.020.20
        12SubsRate1-0.30029kR10.010.06
        13SubsRate2-0.30-0.0530kH10.050.17
        14SubsRate3-0.55-0.0131WaveBase120.0059.98
        15SubsRate4-0.65-0.0532kL100
        16SubsRate5-0.87-0.0533DFlex1E+244.998E+26
        17PDX130 000.0079 980.0034rkM100.10
      • 纳入全部34个可调整的物理参数,并通过随机抽样生成了大量的样本,构建了深度学习所需的数据库。样本为沉积相剖面相模型及相对应的34个输入参数,沉积相取值0~7,其中0表示背景值,数字1~7对应图2中图例中的a1、b1、b2、c1、c2、c3、c4,表示不同的相类型,该相类型的划分方案同时考虑了沉积物的水体深度和水体能量。模型以二维矩阵形式表示,矩阵尺寸为100×129,建立了地层沉积过程正演模型样本库。通过人工和自动筛选,去除模拟结果明显不真实不合理的样本后,该数据库共计纳入训练样本416 079个,约占总样本的83%;测试样本83 167个,约占总样本的17%。训练样本和测试样本相互独立,保证深度学习模型不会陷入过拟合。从训练样本中随机抽取几个样本(图2),可以看出它们具备不同的沉积相结构样式,体现了样本包含的沉积模式的多样性,这9个样本采用的正演模拟参数具体见附表1

        图  2  训练集随机抽取的9个训练样本

        Figure 2.  Nine samples randomly selected from training set

      • 通过测试,选择了改进后的卷积对抗神经网络GLS-GAN[27]。原始GAN要最小化生成分布与真实分布的KL散度,同时要最大化两者的JS散度,在数值上会导致梯度不稳定(通常存在梯度消失和梯度爆炸的可能),并且KL散度具有不对称性,降低了生成器的多样性,会造成模式坍缩问题。GLS-GAN具备损失敏感的自适应损失函数,在训练过程中自动调整损失函数,如果某些区域的真实样本和生成样本已经很接近,则生成器的优化重点就转移到真实样本和生成样本依然差异很大的区域,给GAN提供“按需分配”的建模能力,解决了原版GAN及其衍生算法的梯度消失和模式坍缩问题。

        为了保持物理参数和正演结果的一一对应关系,我们在训练生成器时还需要考虑逐个像素的均方误差(Mean Square Error, MSE)。整个网络均由卷积Conv2D单元和转置卷积ConvTranspose2d单元组成,不包含全连接层,每个单元配套一个批量标准化层BatchNorm2d和非线性化层ReLU,通过大量测试推荐采用5~8个卷积层,生成网络和判别网络的具体结构如图3

        图  3  生成对抗网络的结构(生成器和判别器)

        Figure 3.  Structure of generative adversarial network (generator and discriminator)

        利用上述416 079个训练样本,进行56轮次的迭代训练,取得了比较理想的效果。将不同训练阶段的网络应用在测试集的同一套参数上,可以看出不断优化的模拟效果,相序关系变得更加合理(图4)。把最终训练后的网络应用在测试集,随机抽取几个模型,都具备合理的沉积相分布模式(图5),这6个样本的输入参数见附表2。可见相同输入参数下LS-GAN模型的输出与地层沉积过程模型的输出一致性很强,虽然没有达到完全相同,但是基本实现了合理的相带组合特征。通过34个参数生成符合条件的相模型本身是个难度很大问题,GLS-GAN很大程度上改善了模式坍缩问题,但训练难度大、收敛过程不稳定是GAN类方法自身的属性,无法完全避免。另外,网络结构、样本数量,样本质量等也会影响最终的效果。目前认为形成的替代模型可以用于下一步的反演模拟。

        图  4  不同训练轮次的网络在测试集某一个样本上的模拟效果

        Figure 4.  Simulation results of networks different training rounds on a sample of the test set

        图  5  测试集中随机抽取的6组神经网络生成器模拟(假)与传统正演模拟(真)比较

        Figure 5.  Comparison between generator simulation (false) and conventional forward simulation (true) on randomly selected 6 samples from test set

        表 附表2  测试9随机抽取的6个样本

        参数名称样本A样本B样本C样本D样本E样本F
        InitTopoElev1-235.872-211.007-240.939-226.946-216.507-134.497
        InitTopoElev2-294.878-317.912-305.656-283.186-314.672-294.664
        InitTopoElev3-440.537-404.218-405.011-363.598-402.461-429.075
        InitTopoElev4-611.515-604.388-572.029-611.646-603.645-602.346
        AmpSeaL045.907 744.462 854.435 144.457 345.539 245.516 2
        PerSeaL0302.895297.319283.025327.798305.036285.534
        PhaseSeaL0-0.313 27-0.268 81-0.295 00-0.265 74-0.125 55-0.318 12
        AmpSeaL131.055 716.873 314.874 916.619 724.722 816.481 8
        PerSeaL122.732 124.420 021.763 342.429 351.685 638.696 4
        PhaseSeaL1-0.751 93-0.711 99-0.647 79-0.664 00-0.199 81-0.670 46
        SeaLevelDatum-23.082 0-29.752 5-33.811 0-35.430 4-37.047 00.349 351
        SubsRate1-0.236 21-0.122 30-0.224 62-0.231 28-0.232 31-0.252 31
        SubsRate2-0.248 49-0.230 61-0.260 81-0.253 01-0.245 30-0.262 11
        SubsRate3-0.450 76-0.489 44-0.455 16-0.483 23-0.454 18-0.485 23
        SubsRate4-0.572 49-0.528 56-0.519 46-0.463 54-0.536 64-0.538 34
        SubsRate5-0.802 85-0.600 26-0.794 73-0.775 51-0.636 48-0.770 71
        PDX139 367.251 846.150 318.148 834.040 938.037 031.7
        PDY149 483.238 700.841 724.946 108.640 312.340 762.6
        KDX1303.702319.104317.759346.898267.878261.639
        KDY1262.673258.522300.591355.141252.191270.105
        Aprod12.488 891.720 792.116 181.935 731.895 371.905 39
        Kprod10.078 6540.057 4640.058 0980.071 6170.050 4010.054 544
        Wprod118.519 719.082 221.298 227.134 619.764 018.834 8
        kEng10.682 2370.835 3050.670 8500.660 1230.678 8570.826 163
        bEng10.348 7710.268 2140.268 6900.302 3730.271 5500.259 995
        pPr10.345 8560.399 6250.372 7460.233 5660.392 5120.267 308
        kV116.119 317.201 619.910 214.708 316.862 415.730 3
        kW10.063 5140.042 1870.049 6630.048 8180.040 6130.063 422
        kR10.016 3350.018 6690.019 1810.035 4210.023 0360.015 796
        kH10.074 9750.070 5310.065 3820.063 5100.086 6670.063 957
        WaveBase126.208 436.326 034.947 224.898 526.646 733.491 0
        kL10.000 1050.000 1388.20E-054.97E-050.000 1244.43E-05
        DFlex1.34E+261.17E+269.29E+256.42E+258.00E+251.58E+26
        rkM10.039 6950.012 1370.030 7000.020 0770.014 1420.022 487
      • 在开始真实观测数据反演模拟之前,先从测试集中随机选择样本进行参数反演。该样本的参数及其沉积过程模型模拟结果都是已知的,可以用于验证参数反演是否有效。在水平方向等间距设置了6个虚拟测井,作为反演系统的观测数据。用观测数据与模拟结果的均方根误差作为目标函数,优化器采用上述提到的SCE-UA算法[19]

        用SCU-UA进行了500 000次迭代优化,得到了预期的效果(图6)。图6a展示了参数的真实值(灰色虚线)和SCE-UA算法得到的最优参数值(蓝色实线),图6b展示了SCE-UA优化过程中目标函数的变化,图6c展示了替代模型LS-GAN给出的地层沉积相,其中虚拟测井位置上的沉积相用真实值(即沉积过程模型模拟值)代替。模拟结果表明,随着优化过程的进行,目标函数稳定下降并不再降低,表明SCE-UA算法已经收敛。然而反演所得到的参数值与参数的真实值有较大差异,且互不相同。这一结果表明地层沉积过程的参数反演是一个高度复杂的问题,存在多个局部极小值。

        图  6  合成样本反演结果

        Figure 6.  Inverse result on synthetic sample

        采用类似的方法进行真实观测数据的反演。将代理模型加入地层沉积反演模拟系统,替换传统的地层正演模拟器,实现参数反演。选择4口科考钻井的人工解释的沉积数据为观测数据,将钻井穿过网格的像素均方根误差为目标函数,SCE-UA为反演优化算法,进行参数反演。对于真实数据,我们并不知道测井以外区域的沉积相,也不知道真实的物理参数的值,因此只能靠专家经验判断结果是否合理。从反演结果可以看出,使用替代模型反演得到的地层沉积相的层序结构明显,分布较为合理(图7)。与基于正演模拟器的反演方法相比,结果一致性较强,效率有大幅提升(图8)。

        图  7  真实数据反演模拟结果

        Figure 7.  Inverse result on actual observed data

        图  8  基于正演模拟器的反演模拟结果

        Figure 8.  Inverse result by conventional forward simulation

      • 相比传统的SCE-UA算法驱动下调直接采用正演模拟器反演,基于深度学习代理模型的反演在效率方面有大幅提升。选取全部34个参数参与优化,为了使SCE-UA算法充分收敛,使用SCE-UA直接优化地层沉积过程模型往往需要不间断执行一个月。相比之下,基于神经网络代理模型的优化过程在普通计算机上仅需数分钟即可完成。基于深度学习代理模型的地层沉积反演模拟在效率方面有了大大提高,但也存在一定的局限性。

        代理模型依赖于正演模拟器,而且需要一定的限制条件。代理模型缺乏正演模拟器所拥有的模拟能力完备性,尽管生成样本时已经考虑了模型的泛化性能,但是由于沉积模拟的参数根本无法穷尽,特别是初始地形和构造沉降参数,因此只能构建限定条件下的模型代理。这也是代理模型方法逊于传统反演系统的地方。

        基于深度学习代理模型构建可以进一步推广应用。本文实例研究展示的是斜坡和镶边类型的碳酸盐岩沉积在二维剖面上的代理模型,其实也可以进行孤立台地、潟湖等更复杂地形条件的下的代理模型构建。也有望推广到三维沉积模拟,不过神经网络的结构、需要的样本数量等具体参数也需要进一步探索。本文展示的是针对碳酸盐岩地层沉积正演模拟的深度学习代理模型构建方法,理论上该方法可以用于其他领域的正演模拟替代模型构建,比如油气藏渗流模型、地下水流动模型等。

      • (1) 地层沉积过程反演模拟系统是在优化器的驱动下,不断调整正演模拟输入参数,提高模拟结果与观测数据的吻合程度,大大增加了沉积过程模拟建模方法的实用性,但是存在反演系统非线性强、耗时长、效率较低的问题。

        (2) 生成对抗网络具有较强的模式生成能力,在一定条件下,利用该方法建立的碳酸盐岩地层沉积过程模拟系统的代理模型,具有较强的泛化能力。

        (3) 将代理模型带入沉积模拟反演系统,替换正演模拟器,可以大大提高反演的效率。通过巴哈马滩西缘斜坡现代沉积的反演模型实例,验证了该方法的可行性。

        表 附表1  训练集随机抽取的9个样本实例的参数

        参数名称样本A样本B样本C样本D样本E样本F样本G样本H样本I样本J
        InitTopoElev1-241.740-239.795-240.203-189.967-210.110-237.965-212.366-176.158-213.844-228.109
        InitTopoElev2-304.687-308.444-288.532-302.394-315.569-316.927-311.806-313.322-299.977-290.627
        InitTopoElev3-443.803-420.896-440.692-372.690-426.142-429.549-419.037-445.090-444.184-438.012
        InitTopoElev4-608.549-594.739-606.268-608.710-611.325-600.513-599.063-603.759-597.293-608.271
        AmpSeaL047.993 049.966 849.034 452.207 746.511 847.745 151.338 563.144 546.858 550.374 2
        PerSeaL0415.966333.573350.251319.018288.357427.918302.105338.303314.502308.967
        PhaseSeaL0-0.054 27-0.013 71-0.303 06-0.270 38-0.078 13-0.331 04-0.128 830.177 19-0.295 01-0.309 59
        AmpSeaL119.470 515.821 215.208 014.819 318.243 819.648 729.210 715.980 515.315 417.789 8
        PerSeaL136.966 050.087 929.274 725.920 328.986 729.643 435.675 930.717 837.475 227.578 2
        PhaseSeaL1-0.440 86-0.645 22-0.632 37-0.568 44-0.768 05-0.692 32-0.678 35-0.655 46-0.690 46-0.746 04
        SeaLevelDatum-39.560 5-38.872 9-37.119 8-34.170 6-34.217 7-38.279 3-39.223 3-37.920 3-39.168 8-27.961 8
        SubsRate1-0.243 89-0.247 15-0.253 77-0.258 66-0.258 47-0.228 18-0.243 00-0.236 90-0.242 21-0.195 93
        SubsRate2-0.260 98-0.219 83-0.251 05-0.240 68-0.228 05-0.246 30-0.163 37-0.237 58-0.260 05-0.249 80
        SubsRate3-0.370 92-0.417 32-0.476 49-0.453 03-0.408 79-0.467 76-0.427 09-0.454 19-0.458 59-0.487 05
        SubsRate4-0.563 53-0.496 77-0.527 97-0.487 82-0.444 86-0.581 18-0.539 59-0.519 66-0.311 60-0.544 13
        SubsRate5-0.752 94-0.781 76-0.671 17-0.763 93-0.687 47-0.667 36-0.593 09-0.723 55-0.767 86-0.653 68
        PDX137 175.743 489.039 140.935 975.737 643.636 044.036 179.240 122.340 944.552 077.4
        PDY157 366.054 775.444 088.448 277.344 436.440 488.039 890.340 007.238 728.038 879.8
        KDX1275.378266.399245.489256.212251.925276.658245.710255.417396.307275.302
        KDY1262.749249.093278.053312.542309.070284.470270.730259.157435.309438.423
        Aprod12.617 082.159 542.729 962.279 542.192 181.835 661.766 423.566 482.022 241.734 97
        Kprod10.072 5770.049 2170.049 6030.051 5100.052 6950.053 1920.049 7300.059 0710.049 7130.049 448
        Wprod118.598 819.906 022.985 822.905 626.301 742.079 823.508 324.068 423.353 318.509 9
        kEng10.671 6430.738 1030.727 7840.666 1550.819 8350.799 8260.620 3400.714 9960.923 4370.725 410
        bEng10.271 1130.282 1440.299 4230.336 5460.261 0210.317 6330.246 5010.282 420.257 8690.391 195
        pPr10.206 7960.196 9490.268 6950.204 3560.305 7930.231 9750.261 1850.364 4550.210 0160.258 421
        kV128.296 315.356 414.039 513.641 217.244 115.175 515.189 615.86520.368 916.541 7
        kW10.058 6210.051 8510.061 6310.044 9440.053 1900.125 9660.043 5920.060 9930.050 8900.041 966
        kR10.022 4550.024 6010.025 0260.034 4260.015 8120.016 9740.016 9300.017 060.017 7120.016 261
        kH10.063 8970.065 0080.120 7050.069 2720.063 4450.077 0700.137 2120.083 2970.077 0880.076 075
        WaveBase128.713 327.629 230.370 325.878 631.746 928.013 127.847 924.561930.966 831.201 4
        kL10.000 1074.65E-054.53E-056.79E-055.09E-054.26E-056.15E-054.73E-055.87E-054.86E-05
        DFlex1.19E+267.95E+257.63E+256.58E+251.13E+265.80E+256.33E+255.81E+257.15E+251.94E+26
        rkM10.016 5680.027 5960.031 4030.012 2240.031 8960.017 9190.015 1680.018 8520.014 1580.034 389
    参考文献 (27)

    目录

      /

      返回文章
      返回