Simulation method of stratigraphic uncertainty using a boundary model and generalized coupled Markov chain model
-
摘要: 地层变异性对岩土结构物的性能评价影响显著, 地层变异性的准确表征对工程实际具有重要意义。为此, 提出了一种有效的地层变异性模拟方法, 在概率框架内, 将边界模型和广义耦合马尔可夫链模型相结合形成一种组合模型, 以综合利用两者的优势。首先, 通过贝叶斯方法识别边界模型参数, 进而采用条件随机场对地层边界进行模拟。然后, 将边界模型模拟结果作为背景信息用于广义耦合马尔可夫链模型中, 实现两种模型的结合。最后, 以香港某建筑场地地层为例, 对3种不同模型的地层变异性模拟结果进行了比较, 以阐明该组合模型的优势, 并探讨钻孔数量及其位置对所提组合模型地层变异性模拟的影响。结果表明: 相比于边界模型与广义耦合马尔可夫链模型, 该组合模型不仅能够模拟复杂的地层分布情况, 而且能够考虑边界的空间分布趋势, 能够有效地避免地质异常现象; 钻孔布置方案对地层模拟的不确定性及其实现均具有较大影响。Abstract: Stratigraphic uncertainty has a significant impact on the performance evaluation of geotechnical structures, and it is important for engineering practice to accurately characterize the uncertainty. Hence, an effective simulation method of geological uncertainty is proposed. In the framework of probability, the boundary model and the generalized coupled Markov chain model are coupled to make full use of their advantages. First, the parameters of the boundary model are identified by the Bayes method, and then the boundary of the rock-soil material is simulated by a conditional random field. Then, the boundary model simulation results are used as background information in the generalized coupled Markov chain model to realize the coupling of the two models. Finally, taking a construction site in Hong Kong as an example, the simulation results of stratigraphic uncertainties for three different models are compared, and the advantages of the combination model proposed in this paper are clarified. The effects of the number and location of boreholes on the stratigraphic uncertainty simulation are discussed. The results show that compared with the boundary model and the generalized coupled Markov chain model, the coupling model can not only simulate the complicated stratum distribution but also consider the spatial distribution trend of the boundary and effectively avoid the geological anomaly phenomenon. The borehole layout scheme has a great impact on the stratigraphic uncertainty and its realization.
-
图 1 已知样本点与未知样本点之间的两种可能子方向转移
Figure 1. Two possible transfer modes between a known sample point and an unknown sample point
图 2 基于边界模型和广义耦合马尔可夫链模型的地层模拟过程
a.钻孔数据；b.背景信息(边界模型)；c.地层实现；d.岩土体类型转移方式。Sb, So, Sr, Si, Sp, Sq为相应单元的土体类型
Figure 2. Stratigraphic uncertainty simulation process using the boundary model and generalized coupled Markov charn model
图 3 香港某建筑场地钻孔位置及其揭示的地层
Figure 3. Location and soil strata of the boreholes at a construction site in Hongkong
图 4 基于不同模型的两次地层实现
a1.基于边界模型的地层实现1；b1.基于广义耦合马尔可夫链的地层实现1；c1.基于组合模型的地层实现1；a2.基于边界模型的地层实现2；b2.基于广义耦合马尔可夫链的地层实现2；c2.基于组合模型的地层实现2
Figure 4. Two stratigraphic realizations based on different models
图 5 基于不同模型的信息熵图
a.基于广义耦合马尔可夫链模型的信息熵图；b.基于组合模型的信息熵图
Figure 5. Information entropy maps based on different models
图 6 不同钻孔布置方案对应的信息熵图
Figure 6. Information entropy maps corresponding to different borehole layout schemes
图 7 不同钻孔布置方案对应的一次地层实现
Figure 7. One stratigraphic realization corresponding to different borehole layout schemes
表 1 竖直方向转移概率矩阵
Table 1. Probability matrix of vertical transfer
(a1)竖直向下转移概率矩阵 (a2)竖直向上转移概率矩阵 材料1 材料2 材料3 材料4 材料5 材料1 材料2 材料3 材料4 材料5 材料1 0.773 0.045 0.182 0 0 1.000 0 0 0 0 材料2 0 0.855 0.145 0 0 0.018 0.855 0.127 0 0 材料3 0 0.048 0.910 0.042 0 0.028 0.056 0.910 0.007 0 材料4 0 0 0.006 0.955 0.039 0 0 0.033 0.950 0.017 材料5 0 0 0 0.031 0.969 0 0 0 0.069 0.931 
下载: 导出CSV
表 2 水平方向转移概率矩阵
Table 2. Probability matrix of horizontal transfer
(a1)水平向右转移概率矩阵 (a2)水平向左转移概率矩阵 材料1 材料2 材料3 材料4 材料5 材料1 材料2 材料3 材料4 材料5 材料1 0.946 0.011 0.043 0 0 0.949 0.010 0.041 0 0 材料2 0 0.968 0.032 0 0 0 0.970 0.030 0 0 材料3 0 0.010 0.981 0.009 0 0 0.009 0.982 0.008 0 材料4 0 0 0.001 0.991 0.008 0 0 0.001 0.992 0.007 材料5 0 0 0 0.006 0.994 0 0 0 0.006 0.994
下载: 导出CSV
表 3 已知样本位置及边界1和2对应的深度
Table 3. Known sample locations and corresponding depths of boundaries 1 and 2
钻孔位置x/m 5.32 14.32 26.54 34.25 41.32 边界1深度z/m -19.48 -21.56 -17.10 -14.47 -17.55 边界2深度z/m -42.49 -39.51 -38.12 -34.48 -37.51
下载: 导出CSV
表 4 边界模型中模型参数贝叶斯方法识别结果
Table 4. Bayesian identification results of model parameters in the boundary model
边界 线性趋势函数 标准差σ/m 波动范围δ/m 边界1 Z=0.15x-21.71 1.03 6.32 边界2 Z=0.21x-43.70 0.95 6.78
下载: 导出CSV
表 5 不同的钻孔布置方案
Table 5. Different borehole layout schemes considered in this study
钻孔布置方案 钻孔1 钻孔2 钻孔3 钻孔4 钻孔5 
3 √ √ √ 4A √ √ √ √ 4B √ √ √ √ 5 √ √ √ √ √ 注：“√”代表该钻孔布置方案包含对应的钻孔
下载: 导出CSV
-
[1] Fenton G A, Griffiths D V. Reliability-based deep foundation design[M]. [S. l. ]: Probabilistic Applications in Geotechnical Engineering, 2007: 1-12. [2] 邓志平, 李典庆, 曹子君, 等. 考虑地层变异性和土体参数变异性的边坡可靠度分析[J]. 岩土工程学报, 2017, 39(6): 986-995. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201706005.htmDeng Z P, Li D Q, Cao Z J, et al. Slope reliability analysis considering geological uncertainty and spatial variability of soil parameters[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(6): 986-995(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201706005.htm [3] Li D Q, Qi X H, Cao Z J, et al. Evaluating slope stability uncertainty using coupled Markov chain[J]. Computers and Geotechnics, 2016, 73: 72-82. doi: 10.1016/j.compgeo.2015.11.021 [4] 于正, 杨龙才, 张勇, 等. 考虑地层变异特征一致性的围岩变形不确定性分析[J]. 岩土力学, 2019, 40(5): 1947-1956. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201905038.htmYu Z, Yang L C, Zhang Y, et al. Uncertainty analysis of tunnel surrounding rock deformation considering consistency of geological heterogeneity features[J]. Rock and Soil Mechanics, 2019, 40(5): 1947-1956(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201905038.htm [5] 连志鹏, 徐勇, 付圣, 等. 采用多模型融合方法评价滑坡灾害易发性: 以湖北省五峰县为例[J]. 地质科技通报, 2020, 39(3): 178-186. doi: 10.19509/j.cnki.dzkq.2020.0319Lian Z P, Xu Y, Fu S, et al. Landslide susceptibility assessment based on multi-model fusion method: A case study in Wufeng County, Hubei Province[J]. Bulletin of Geological Science and Technology, 2020, 39(3): 178-186(in Chinese with English abstract). doi: 10.19509/j.cnki.dzkq.2020.0319 [6] Phoon K K, Kulhawy F H. Characterization of geotechnical variability[J]. Canadian Geotechnical Journal, 1999, 36(4): 612-624. doi: 10.1139/t99-038 [7] 官东林, 文国军, 王玉丹, 等. 基于线激光扫描的岩石激光钻孔的三维重建和可视化[J]. 地质科技通报, 2021, 40(3): 173-183. doi: 10.19509/j.cnki.dzkq.2021.0310Guan D L, Wen G J, Wang Y D, et al. 3D reconstruction and visualization for laser drilling hole on rock based on linelaser scanning[J]. Bulletin of Geological Science and Technology, 2021, 40(3): 173-183(in Chinese with English abstract). doi: 10.19509/j.cnki.dzkq.2021.0310 [8] Gong W P, Tang H M, Wang H, et al. Probabilistic analysis and design of stabilizing piles in slope considering stratigraphic uncertainty[J]. Engineering Geology, 2019, 259: 105-162. [9] Zhu H, Zhang L M. Characterizing geotechnical anisotropic spatial variations using random field theory[J]. Canadian Geotechnical Journal, 2013, 50(7): 723-734. doi: 10.1139/cgj-2012-0345 [10] 谭晓慧, 董小乐, 费锁柱, 等. 基于KL展开的可靠度分析方法及其应用[J]. 岩土工程学报, 2020, 42(5): 808-816. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202005005.htmTan X H, Dong X L, Fei S Z, et al. Reliability analysis method based on KL expansion and its application[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(5): 808-816(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202005005.htm [11] Jiang S H, Papaioannou I, Straub D. Bayesian updating of slope reliability in spatially variable soils with in-situ measurements[J]. Engineering Geology, 2018, 239: 310-320. doi: 10.1016/j.enggeo.2018.03.021 [12] 李典庆, 蒋水华, 周创兵, 等. 考虑参数空间变异性的边坡可靠度分析非侵入式随机有限元法[J]. 岩土工程学报, 2013, 35(8): 1413-1422. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201308008.htmLi D Q, Jiang S H, Zhou C B, et al. Reliability analysis of slope using non-intrusive stochastic finite element method considering spatial variability of soil parameters[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(8): 1413-1422(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201308008.htm [13] Wu Y, Zhou X, Gao Y, et al. Effect of soil variability on bearing capacity accounting for non-stationary characteristics of undrained shear strength[J]. Computers and Geotechnics, 2019, 110: 199-210. doi: 10.1016/j.compgeo.2019.02.003 [14] Phoon K K, Retief J V, Ching J, et al. Some observations on ISO2394: 2015 Annex D(reliability of geotechnical structures)[J]. Structural Safety, 2016, 62: 24-33. doi: 10.1016/j.strusafe.2016.05.003 [15] Li Z, Wang X, Wang H, et al. Quantifying stratigraphic uncertainties by stochastic simulation techniques based on Markov random field[J]. Engineering Geology, 2016, 201: 106-122. doi: 10.1016/j.enggeo.2015.12.017 [16] 邓志平, 牛景太, 潘敏, 等. 考虑地层变异性和土体参数空间变异性的边坡可靠度全概率设计方法[J]. 岩土工程学报, 2019, 41(6): 1083-1090. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201906014.htmDeng Z P, Niu J T, Pan M, et al. Full probabilistic design method of slopes considering geological uncertainty and spatial variability of soil parameters[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(6): 1083-1090(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201906014.htm [17] Elkateb T, Chalaturnyk R, Robertson P K. An overview of soil heterogeneity: Quantification and implications on geotechnical field problems[J]. Canadian Geotechnical Journal, 2003, 40(1): 1-15. doi: 10.1139/t02-090 [18] 宋友桂, 兰敏文, 刘慧芳, 等. 关中盆地新生界地层划分对比与第四纪下限[J]. 地质科技通报, 2021, 40(2): 24-35. doi: 10.19509/j.cnki.dzkq.2021.0204Song Y G, Lan M W, Liu H F, et al. Cenozoic stratigraphic correlation and the lower limit of Quaternary in Guanzhong Basin[J]. Bulletin of Geological Science and Technology, 2021, 40(2): 24-35(in Chinese with English abstract). doi: 10.19509/j.cnki.dzkq.2021.0204 [19] Liao T, Mayne P W. Stratigraphic delineation by three-dimensional clustering of piezocone data[J]. Georisk, 2007, 1(2): 102-119. [20] Cao Z, Wang Y. Bayesian approach for probabilistic site characterization using Cone Penetration Tests[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2013, 139(2): 267-276. doi: 10.1061/(ASCE)GT.1943-5606.0000765 [21] Ching J, Wang J, Juang C H, et al. Cone penetration test(CPT)-based stratigraphic profiling using the wavelet transform modulus maxima method[J]. Canadian Geotechnical Journal, 2015, 52(12): 1993-2007. doi: 10.1139/cgj-2015-0027 [22] Zhao T, Wang Y. Interpolation and stratification of multilayer soil property profile from sparse measurements using machine learning methods[J]. Engineering Geology, 2020, 265: 105430. doi: 10.1016/j.enggeo.2019.105430 [23] Li X Y, Zhang L M, Li J H. Using conditioned random field to characterize the variability of geologic profiles[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2016, 142(4): 04015096. doi: 10.1061/(ASCE)GT.1943-5606.0001428 [24] 邓志平, 牛景太, 潘敏, 等. 基于不同材料分类模型的地层变异性模拟比较研究[J]. 自然灾害学报, 2018, 27(4): 128-136. https://www.cnki.com.cn/Article/CJFDTOTAL-ZRZH201804017.htmDeng Z P, Niu J T, Pan M, et al. Comparative study on geological uncertainty simulation based on different material classification models[J]. Journal of Natural Disasters, 2018, 27(4): 128-136(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-ZRZH201804017.htm [25] Gong W, Zhao C, Juang C H, et al. Stratigraphic uncertainty modelling with random field approach[J]. Computers and Geotechnics, 2020, 125: 103681. doi: 10.1016/j.compgeo.2020.103681 [26] Zhao C, Gong W, Li T, et al. Probabilistic characterization of subsurface stratigraphic configuration with modified random field approach[J]. Engineering Geology, 2021, 288: 106138. doi: 10.1016/j.enggeo.2021.106138 [27] Elfeki A, Dekking M. A Markov Chain Model for subsurface characterization: Theory and applications[J]. Mathematical Geology, 2001, 33(5): 569-589. doi: 10.1023/A:1011044812133 [28] Park E. A multidimensional, generalized coupled Markov chain model for surface and subsurface characterization[J]. Water Resources Research, 2010, 46(11): 6291-6297. [29] Deng Z P, Jiang S H, Niu J T, et al. Stratigraphic uncertainty characterization using generalized coupled Markov chain[J]. Bulletin of Engineering Geology and the Environment, 2020, 79(10): 5061-5078. doi: 10.1007/s10064-020-01883-y [30] 邓志平, 李典庆, 祁小辉, 等. 基于广义耦合马尔可夫链的地层变异性模拟方法[J]. 岩土工程学报, 2018, 40(11): 2041-2050. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201811012.htmDeng Z P, Li D Q, Qi X H, et al. Simulation of geological uncertainty using modified generalized coupled Markov chain[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(11): 2041-2050(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201811012.htm [31] Qi X H, Li D Q, Phoon K K, et al. Simulation of geologic uncertainty using coupled Markov chain[J]. Engineering Geology, 2016, 207: 129-140. doi: 10.1016/j.enggeo.2016.04.017 [32] Li W, Zhang C, Burt J E, et al. Two-dimensional Markov chain simulation of soil type spatial distribution[J]. Soil Science Society of America Journal, 2004, 68(5): 1479-1490. doi: 10.2136/sssaj2004.1479 [33] Li J, Cai Y, Li X, et al. Simulating realistic geological stratigraphy using direction-dependent coupled Markov chain model[J]. Computers and Geotechnics, 2019, 115: 103147. doi: 10.1016/j.compgeo.2019.103147 [34] Wang Y, Cao Z, Li D. Bayesian perspective on geotechnical variability and site characterization[J]. Engineering Geology, 2016, 203: 117-125. doi: 10.1016/j.enggeo.2015.08.017 [35] 吴振君, 王水林, 葛修润. 约束随机场下的边坡可靠度随机有限元分析方法[J]. 岩土力学, 2009, 30(10): 3086-3092. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200910040.htmWu Z J, Wang S L, Ge X R. Slope reliability analysis by random FEM under constraint random field[J]. Rock and Soil Mechanics, 2009, 30(10): 3086-3092(in Chinese with English abstract). https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200910040.htm [36] Forrester A, Sobester A, Keane A. Engineering design via surrogate modelling: A practical guide[M]. [S. l. ]: John Wiley & Sons, 2008. -
下载:
点击查看大图
图(7) / 表(5)
计量
- 文章访问数: 719
- PDF下载量: 70
- 被引次数: 0
