基于元胞自动机原理的土石坝溃决分析方法技术

本发明专利技术公开了一种基于元胞自动机原理的土石坝溃决分析方法，首先在坝基或坝顶的合适位置为原点建立坐标系，将土石坝按二维Von模型规范划分为若干个等长度正方形网格，其中每个网格即代表一个元胞，随后表示出每个元胞的时空状态，处在溃口位置的元胞视为中心元胞；然后建立元胞所处的溃口条件，以及冲蚀条件与流量条件限定下的元胞演化模型，并设置边界条件；最后通过MATLAB软件演算出各个步时内坝体各处元胞的所处状态，判断极限状态元胞数量，随后采用数值分析方法对模拟数据与实测结果进行对比，检验相关性与可靠性。本发明专利技术解决了现有技术中存在的土石坝溃决模型运行时工作量大、精度较低的问题。精度较低的问题。精度较低的问题。

【技术实现步骤摘要】
基于元胞自动机原理的土石坝溃决分析方法


[0001]本专利技术属于水利工程
，具体涉及一种基于元胞自动机原理的土石坝溃决分析方法。

技术介绍

[0002]诱发大坝溃决的因素复杂，后果严重，因此加强对土石坝溃决模型的研究有重要意义。
[0003]堤坝的破坏类型，主要包括漫顶、冲蚀与坝坡失稳。土石坝由于其筑坝材料与结构特点，常见的溃决形式为渐溃——这种溃坝模式的特点是坝体随着洪水演进过程逐渐破坏，稳定性随时间逐渐丧失且破坏持续时间较长。造成土石坝溃坝的原因多种多样，包括但不限于渗流、土料特性、溃口形态、破坏位置、水流特征、溃决程度以及溃坝发展过程。不同土石坝失效的原因也不尽相同。例如相较于均质土坝，面板堆石坝更容易受近源地面运动的影响而发生破坏。由此可见，土石坝的溃决受到多重因素影响，不可控变量较多：一方面，数值模型通常有侧重点，着重于研究解决其中的某一问题，因此不能做到同时考虑所有影响因素；另一方面，变量之间通常相互影响，某个变量的改变会带动一系列的连锁反应。如温度的变化会引起渗透系数与水的粘滞系数发生改变。因此目前有关土石坝溃坝机理与模型的研究还处于发展阶段，计算精度仍有较大提升空间。

技术实现思路

[0004]本专利技术的目的是提供一种基于元胞自动机原理的土石坝溃决分析方法，解决了现有技术中存在的土石坝溃决模型运行时工作量大、精度较低的问题。
[0005]本专利技术所采用的技术方案是，基于元胞自动机原理的土石坝溃决分析方法，具体按照以下步骤实施：
[0006]步骤1、在坝基或坝顶的合适位置为原点建立坐标系，将土石坝按二维Von模型规范划分为若干个等长度正方形网格，其中每个网格即代表一个元胞，随后表示出每个元胞的时空状态，处在溃口位置的元胞视为中心元胞；
[0007]步骤2、建立元胞所处的溃口条件，以及冲蚀条件与流量条件限定下的元胞演化模型，并设置边界条件；
[0008]步骤3、通过MATLAB软件演算出各个步时内坝体各处元胞的所处状态，判断极限状态元胞数量，随后采用数值分析方法对模拟数据与实测结果进行对比，通过均方根误差RMSE、相对误差Re和方差R2三项参数作为敏感性分析的依据，检验相关性与可靠性。
[0009]本专利技术的特点还在于，
[0010]步骤1具体如下：
[0011]按照二维的Von.Neuman模型将研究对象划分为同等尺寸的正方形元胞，并由唯一确定的中心元胞开始向周围8个邻域元胞迭代扩展。
[0012]步骤2具体如下：
[0013]步骤2.1、建立元胞所处的溃口条件以及冲蚀条件与流量条件限定下的元胞演化模型，并设置边界条件，具体表达式如下：
[0014]溃口条件：
[0015]B
Z
＝0.1803KV
0.32
H
0.19
ꢀꢀ
(1)
[0016][0017]其中，B
Z
为溃口最终宽度，B1为溃口初始宽度，B为任意时间t时间的溃口宽度，V为总库容，H为中心元胞的高程，T为预计的溃坝总历时，估算得到；
[0018]步骤2.2、流量条件：
[0019][0020][0021]式中，Q
X(k)
表示X方向上迭代k次时的单宽流量，Q
Y(k)
表示Y方向上迭代k次时的单宽流量，Q
X(k+1)
表示X方向上迭代k+1次时的单宽流量，Q
Y(k+1)
表示Y方向上迭代k+1次时的单宽流量，μ与v分别表示X与Y两个方向的水流速度，k为迭代次数，Z为水面高程，ΔX与ΔY为前后两次迭代的两个元胞的中心高差，由元胞尺寸决定，g为重力加速度，g取值为9.8，h为前后两次迭代的两个元胞浸没的平均水深，t为时间，n为糙率。
[0022]步骤2.3、按照如下计算公式确定边界流量：
[0023]Q
max
＝CK[k1B(Z
‑
H)+k2(Z
‑
H)]1.5
ꢀꢀ
(5)
[0024]式中，Q
max
为坝体所能承受的极限单宽流量，C为修正系数，K为淹没系数，k1和k2为流量系数；
[0025]步骤2.4、在泥沙冲蚀方面采用Schoklitsch推移质公式，即：
[0026][0027]其中，Q
S
为溃口泥沙冲蚀量，Q为通过溃口的流量，由式(3)与式(4)确定，Q＝MAX{Q
X(k+1)
,Q
Y(k+1)
}，D为边坡比，d
84
、d
16
分别为粒径小于该颗粒总重占总体泥沙重量的84％，16％所占比。
[0028]步骤3根据初始实验数据与筑坝材料参数对模型演化过程进行求解，一溃坝时步内的流量与溃口扩展情况作为敏感性分析的评价指标作为对基于元胞自动机原理的土石坝溃决模型参数敏感性进行分析，具体如下：
[0029]步骤3.1、输入参数：溃口初始宽度B1、总库容V、中心元胞的高程H、X与Y两个方向的水流速度μ与v、迭代次数k、水面高程Z、首次迭代的两个元胞的中心高差ΔX与ΔY、修正系数C、淹没系数K、流量系数k1和k2、边坡比D，粒径小于该颗粒总重占总体泥沙重量的84％，16％所占比d
84
、d
16
；
[0030]步骤3.2、输入所述步骤得到的边界条件；
[0031]步骤3.3、综合研究土石坝破坏过程中坝体各处元胞的所处临界状态，当满足临界
破坏条件的元胞数量占比超过总元胞数量的30％时，即可认为后续元胞均发生破坏。临界条件为：max{Q
X
，Q
Y
}>Q
max
；
[0032]反馈实测数据，与步骤2.1、步骤2.2、步骤2.3和步骤2.4得出的流量与溃口面积模拟结果对比并进行敏感性分析，分别与实测值进行均方根误差、相对误差和方差三项参数的对比。
[0033]步骤3.3中敏感性分析具体如下：
[0034]引入的三个参数分别为：(1)相对误差，用Re表示；(2)均方根误差，用RMSE表示；(3)决定系数；用R2表示，三者计算公式如下：
[0035][0036][0037][0038]在上述三式中：O
i
为i时刻的实测值；S
i
为i时刻的模拟值；n为样本容量；为实测值的平均值，
[0039]相对误差Re的作用是衡量模型的可靠程度，相对误差Re值越接近1，则模型的可靠程度越差；均方根误差RMSE的作用是衡量模拟值与实测值的差值可接受程度，均方根误差RMSE值越大，则说明模型模拟结果与实际测量结果的贴合程度越差，而决定系数R2则用衡量相关性，R2数值越小，则模型预估结果与实测数据的偏差越大，模型的相关性越差。
[0040]本专利技术的有益效果是，基于元胞自动机原理的土石坝溃本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于元胞自动机原理的土石坝溃决分析方法，其特征在于，具体按照以下步骤实施：步骤1、在坝基或坝顶的合适位置为原点建立坐标系，将土石坝按二维Von模型规范划分为若干个等长度正方形网格，其中每个网格即代表一个元胞，随后表示出每个元胞的时空状态，处在溃口位置的元胞视为中心元胞；步骤2、建立元胞所处的溃口条件，以及冲蚀条件与流量条件限定下的元胞演化模型，并设置边界条件；步骤3、通过MATLAB软件演算出各个步时内坝体各处元胞的所处状态，判断极限状态元胞数量，随后采用数值分析方法对模拟数据与实测结果进行对比，通过均方根误差RMSE、相对误差Re和方差R2三项参数作为敏感性分析的依据，检验相关性与可靠性。2.根据权利要求1所述的基于元胞自动机原理的土石坝溃决分析方法，其特征在于，所述步骤1具体如下：按照二维的Von.Neuman模型将研究对象划分为同等尺寸的正方形元胞，并由唯一确定的中心元胞开始向周围8个邻域元胞迭代扩展。3.根据权利要求2所述的基于元胞自动机原理的土石坝溃决分析方法，其特征在于，所述步骤2具体如下：步骤2.1、建立元胞所处的溃口条件以及冲蚀条件与流量条件限定下的元胞演化模型，并设置边界条件，具体表达式如下：溃口条件：B
Z
＝0.1803KV
0.32
H
0.19 (1)其中，B
Z
为溃口最终宽度，B1为溃口初始宽度，B为任意时间t时间的溃口宽度，V为总库容，H为中心元胞的高程，T为预计的溃坝总历时，估算得到；步骤2.2、流量条件：步骤2.2、流量条件：式中，Q
X(k)
表示X方向上迭代k次时的单宽流量，Q
Y(k)
表示Y方向上迭代k次时的单宽流量，Q
X(k+1)
表示X方向上迭代k+1次时的单宽流量，Q
Y(k+1)
表示Y方向上迭代k+1次时的单宽流量，μ与v分别表示X与Y两个方向的水流速度，k为迭代次数，Z为水面高程，ΔX与ΔY为前后两次迭代的两个元胞的中心高差，由元胞尺寸决定，g为重力加速度，g取值为9.8，h为前后两次迭代的两个元胞浸没的平均水深，t为时间，n为糙率；步骤2.3、按照如下计算公式确定边界流量：Q
max
＝CK[k1B(Z
‑
H)+k2(Z
‑
H)]
1.5
ꢀꢀꢀꢀꢀꢀꢀ
(5)式中，Q
max
为坝体所能承受的极限单宽流量，C为修正系数，K为淹没系数，k1和k2为流量
系数；步骤2...

【专利技术属性】
技术研发人员：任杰，睢佳衡，雷亚楠，
申请(专利权)人：西安理工大学，
类型：发明
国别省市：