DOI:10.19476/j.ysxb.1004.0609.2002.02.037
基于胞映射理论的岩体动力系统可预测尺度模型
周辉 冯夏庭 谭云亮 王泳嘉
中国科学院武汉岩土力学研究所岩土力学重点实验室
山东科技大学矿压研究所
东北大学资源与土木工程学院 武汉430071
武汉430071.东北大学资源与土木工程学院
沈阳110004
泰安271019
摘 要:
岩体混沌动力系统演化过程的可预测尺度是非线性预测理论研究的重要内容。在利用传统的理论和方法分析实际问题时 , 常常由于观测资料不足而无法计算岩体动力系统的可预测尺度。为了解决这一问题 , 基于胞映射理论提出了一种新的岩体动力系统可预测尺度的计算模型。该模型基于系统的非线性动力学方程和随机过程理论 , 将点到点的映射转化成胞到胞的映射 , 以系统演化达到极限状态概率分布前的时段作为系统的可预测尺度 , 消除了初值观测误差和系统离散化对系统造成的影响 , 从而可以更精确地描述系统的混沌性本质。实例的应用也表明了这一点
关键词:
胞映射理论 ;岩体动力系统 ;可预测尺度 ;
中图分类号: TD315
收稿日期: 2001-07-09
基金: 国家自然科学基金资助项目 (5 0 1790 3 4;5 980 40 0 5 );
Study on predictable time scale of rock mass dynamic system using cell mapping theory
Abstract:
Predictable time scale of the evolution of chaotic rock mass dynamic systems is one of the most important parameters in nonlinear forecasting theory. However, the predictable time scale of most rock mass dynamic systems can't be calculated because of lack of the observing data. In order to solve this problem, a new model to calculate the predictive time scale of rock mass dynamic system was proposed using the cell mapping theory. The model transforms the point to point mapping to cell to cell mapping using the nonlinear dynamic equation of system and stochastic process theory, and considers period of time before evolution of system reaches to the limit probability distribution as the predictable time scale of system. The model eliminates the influences of error of initial conditions and dispersal on the evolution of system, so chaotic essence of system can be characterized more accurately which is also proved by the example. Case study indicates these advantages.
Keyword:
cell mapping theory; rock mass dynamic system; predictable time scale;
Received: 2001-07-09
大量研究表明
[1 ,2 ,3 ,4 ,5 ,6 ,7 ,8 ,9 ]
, 大多数岩体动力系统的演化过程都会表现出明显的混沌性, 从而导致系统演化的长期不可预测性
[8 ]
。 因此, 对于岩体混沌动力系统而言, 其预测理论应包括两个方面的内容, 可预测尺度理论及预测理论和方法。 其中, 可预测尺度问题是岩体混沌动力系统非线性预测理论研究的前提和重要内容。 在传统理论中, 混沌系统的可预测尺度通常采用最大Lyapunov指数LE 1 和所有正的Lyapunov指数之和LE + 的倒数来衡量, 并分别将其定义为混沌系统的最大可预测尺度T max 和平均可预测尺度
Τ
?
。 然而, 在大多数情况下, 要由实际岩体动力系统的观测资料来计算其Lyapunov 指数并不容易。 因为要得到一个较好的最大Lyapunov 指数的估计值, 要求观测资料的个数最少不小于10D0 (D0 为系统混沌吸引子的分维)
[10 ]
; 且由于至少要在大于2维的空间中才会产生混沌, 所以D0 应大于2。 因此, 计算最大Lyapunov 指数所需的观测资料的个数至少应大于102 。 当实际观测资料的个数小于102 时, 所得到的Lyapunov 指数就很可能存在较大误差, 这时, 我们就无法准确计算系统的可预测尺度。 为了解决这一问题, 作者基于胞映射理论的基本原理, 提出一种新的岩体混沌动力系统可预测尺度的计算模型, 该模型可以更深刻地描述系统的混沌性和可预测尺度。
1胞映射理论的基本原理
胞映射理论由Hsu 提出
[11 ,13 ,14 ,15 ,16 ]
, 是非线性动力系统特别是强混沌性系统全局分析的一种有效工具。 胞映射将点到点的映射转化为胞到胞的映射, 可以“平滑”并削弱观测初值的随机性和系统的内在随机性对系统动力学模式的影响。 其基本原理如下。
设N维动力系统的非线性动力学方程可以表示为
d xi /d t=fi (x1 , x2 , …, xN ; αi ) (1)
式中 xi 为系统的状态变量, αi 为控制变量。
把N维状态空间的每个坐标都分成许多小段, 每段用一个整数Zi 表示, 离散化状态空间的N维有限单元称为胞, 以N个整数Zi (i=1, 2, …, N) 构成向量Z (称为胞向量) 来表征。
将系统的动力学行为 (从n 时刻到n +1时刻) 用胞对胞的映射c 表示为离散化的整数方程
Z (n +1) =c Z (n ) (2)
对于胞映射, 需要计算出每个胞在下一时刻映射到多个胞i 中的概率, 记为ξ i 。 把这些概率汇总成胞概率向量, 记为ξ (t ) (t 为时间) 。 在实际计算中, ξ (t ) 可按下述的样本法进行计算:
将第j 个胞离散细化为若干个子胞, 设第k 个子胞的中心点状态向量为X d (n ) , 则其一次映射后的状态向量可由式 (1) 经差分运算求得
X d (n +1) =F (X d (n ) ) (3)
记胞j 的所有子胞经上述映射后落入胞i 的个数为M ij , 则从胞j 到胞i 的一次转移概率为
P ij =M ij /N s (4)
式中 N s 为一个胞内子胞的总数。
由P ij 构成Markov一次转移概率矩阵P , 则系统的动力学行为可用概率方程表示为
ξ (n +1) =P ξ (n ) =P n ξ (0) (5)
式中 ξ (0) 为系统初始时刻的状态概率向量; ξ (n +1) 为系统第n +1时刻的状态概率向量。 利用上式, 可以分析非线性动力系统的各种演化特性。
2基于胞映射理论的岩体动力系统可预测尺度模型
胞映射和随机过程理论已经证明
[11 ,13 ,14 ,15 ,16 ]
, 当n →∞时, 由式 (5) 求得的系统状态概率的极限分布ξ * 与初始状态概率分布ξ (0) (即初始信息) 无关, 初始信息已不再对系统的继续演化起作用。 也就是说, 对于达到极限状态概率分布ξ * 前的系统状态, 可利用初始信息进行预测, 超过了这一范围, 就不能再利用初始信息进行预测。 因此, 达到极限状态概率分布ξ * 前的时段即是系统的可预测尺度。
由随机过程理论可知, 极限状态概率分布ξ * 不依赖于初始状态的信息ξ (0) , 系统的状态概率分布将随演化进程而逐渐趋近于ξ * 。 本文作者采用以下两种指标来判别系统趋近于极限状态概率分布的程度。
1) Shannon信息量:
若事件每个可能结果发生的概率为q i , 则定义与该结果相对应的信息量为
Ι
1
=
-
∑
i
=
1
n
q
i
lg
qi (6)
式中 n 为事件所有可能结果的个数。
在胞映射中, 当迭代到第t 步时, 状态概率不为0的胞的集合对应于系统状态所有的可能结果, 根据式 (6) 可以计算出系统状态在不同迭代时步的信息量。 根据该信息量的演化即可识别系统的状态概率分布随时间而逐渐趋近于极限分布ξ * (某一稳定值) 的程度。
2) 均方误差Δ:
t 时刻系统状态概率分布ξ (t ) 与t -1时刻状态概率分布ξ (t -1) 之间的均方误差为
Δ
=
∑
i
=
1
n
(
ξ
i
(
t
)
-
ξ
i
(
t
-
1
)
)
2
?
?
?
(
7
)
式中 i 为状态胞编号。
若在时刻T 时, Δ≤Δ* (其中, Δ* 为给定的判别阈值) , 就可以近似判定系统的状态在T 时刻已处于极限概率分布状态, 从而系统的可预测尺度即为T 。
采用上述两种判别指标可以尽可能地消除由数值相对量的大小而引起的主观判别误差。
3计算实例
表1为门头沟煤矿某工作面的试验观测数据, 由于只有少量的观测数据, 因此采用传统的方法无法计算系统的可预测尺度。 利用本文作者提出的方法, 则可以计算该系统的可预测尺度, 其方法和步骤如下。
表1 门头沟煤矿某工作面系统变量的试验观测数据
Table 1 Variable data of a working face in Mentougou Coal Mine
Date
y 1
y 2
y 3
1
13
856
9
2
119
6 294
17
3
18
1 283
0
4
10
454
36
5
24
2 068
5
6
30
2 220
5
7
28
1 742
55
8
3
45
37
9
32.3
176
-19
10
35
2 900
-18
11
391
4 129
23
12
349
1 265
41
13
3.6
229
29
14
21
404
-19
15
224
2 313
-10
y 1 is the frequency of acoustic emission, times/minute; y 2 is the energy rate of acoustic emission; y 3 is strain difference of photo-elasticity.
3.1反演系统的非线性动力学方程
为了削弱观测数据的随机性, 消除变量之间的非独立性, 提高反演非线性动力学方程对原始系统的描述精度, 对传统的非线性动力学方程反演方法作如下改进处理。
1) 为了消弱观测数据中的“噪音”成分, 首先对观测数据进行五点三次平滑滤波; 为了消除滤波后变量之间的非独立性, 将滤波值序列生成独立主变量 (具体方法见文献
[
8 ]
) ; 为了加强独立主变量序列中的确定性信息, 削弱其随机成分, 将独立主变量序列进行AGO生成:
y
i
(
1
)
(
k
)
=
∑
j
=
1
k
Y
?
i
(
j
)
?
(
i
=
1
,
2
,
3
)
?
?
?
(
8
)
式中
Y
?
i (j ) 为由滤波序列生成的独立主变量。 利用新的生成序列y
i
(
1
)
(i =1, 2, 3) 作为系统的状态变量来反演其动力学方程。
2) 将系统的非线性动力学方程设为四次多项式 (见文献
[
8 ]
) 。
对表1中3日~9日的观测数据进行上述处理, 可以反演得到系统在3日~9日的非线性动力学方程 (此过程由计算机完成) 。 经过上述改进, 可以显著提高反演非线性动力学方程对系统的描述精度
[8 ]
。
3.2生成胞向量空间和初始胞向量
首先确定系统的可能状态空间。 将原始观测数据进行3次滤波, 进而生成独立主变量序列, 再经累加处理后, 得到系统的可能状态空间为
-1.842≤y
1
(
1
)
≤0.313×10-7 ,
-2.203≤y
2
(
1
)
≤3.58×10-7 ,
-1.402≤y
3
(
1
)
≤1.656 (9)
将上述区间均匀地划分成10×10×10=1 000个胞 (胞向量编号为1至1 000) , 连同沉没胞 (即由系统的所有不可能状态空间组成的胞, 编号为0) 共1 001个胞; 将累加生成序列中第“3日”的数据所对应的胞向量
(
0
?
?
?
0
8
1
3
个
0
?
1
?
0
?
?
?
0
1
8
7
个
0
)
作为初始胞向量。
3.3确定系统的可预测尺度
利用上文所建立的系统的非线性动力学方程和划分的胞向量空间建立一次转移概率矩阵 (用样本法计算一次转移概率矩阵时, 在每个胞内各取1000个样本点进行计算) ; 将初始状态概率向量与一次转移概率矩阵进行迭代, 分别得到系统演化的Shannon 信息量和第i步与第i+1步的系统状态概率向量间的均方误差随迭代步数的变化曲线 (如图1所示) 。 由图1 (a ) 可见, 当迭代至第6步时, Shannon 信息量基本趋于稳定 (第6步与第7步的Shannon 信息量差值为0.003) , 表明在第6步以后, 系统的Shannon 信息量已不再随迭代步数的增加而发生有意义的变化; 同时, 由图1 (b ) 可见, 第6步与第7步的系统状态概率分布的均方误差为8.873×10-5 , 表明第6步的系统状态概率分布已充分接近于系统的极限状态概率分布。 综合上述两种不同指标的识别结果可以判定: 由3日~9日的观测数据所确定的动力系统在演化至第6步时已经进入混沌态, 初始信息不再对系统的继续演化起作用, 因此在该时段内, 系统的可预测尺度为6 d 。
为了比较不同初始状态概率分布 (初始信息) 对可预测尺度的影响, 分别用另外两种初始条件, 条件Ⅰ: 初始状态概率分布—
(
0
?
?
?
0
5
0
0
个
0
?
1
?
0
?
?
?
0
5
0
0
个
0
)
; 条件Ⅱ: 初始状态概率分布—
(
1
1
0
0
1
?
1
1
0
0
1
?
?
?
1
1
0
0
1
1
0
0
1
个
)
, 与上述动力系统的一次状态转移概率矩阵进行迭代, 相应得到这两种初始条件下的Shannon信息量和均方误差随迭代步数的变化曲线分别如图2和图3所示。 可见, 在条件Ⅰ和条件Ⅱ两种不同的初始条件下, 系统的可预测尺度均为6 d, 与图1相同。 因此可以得到如下结论:
图1 Shannon信息量 (a) 和状态概率向量的均方误差 (b) 随迭代步数的变化曲线
Fig.1 Curves of Shannon information (a) and mean squared error of state probability vector (b) according to iteration steps
1) 用胞映射理论可将动力系统的演化描述为: 相对于初始信息而言, 动力系统的演化过程即是状态概率向量的渐进稳定过程;
2) 动力系统的极限状态概率分布仅由其非线性动力学方程完全确定, 而与初始信息无关。
利用本文作者提出的方法计算表1所示的动力系统在不同时段的可预测尺度 (如表2所列) , 可得到该系统的可预测尺度随时间的变化曲线 (如图4所示) 。
需要进一步说明的是: 1) 由胞映射理论得到的动力系统的可预测尺度不是对个别变量的反映, 而是构成系统非线性动力学方程的所有变量综合作用的结果, 因此, 提出的基于胞映射理论的岩体动力系统的可预测尺度是系统的平均可预测尺度; 2) 由于可预测尺度是由系统的混沌性所导致的, 且与混沌性成反比关系, 而由胞映射理论得到的可预测尺度是动力系统整体混沌性的体现, 所以可用表2中的可预测尺度值的倒数来定性表征系统的混沌性 (如表2和图5所示) 。 可见, 图5所示的岩体动力系统在发生破坏前出现了明显的混沌性下降现象, 这与岩体破坏前的“降维”
[1 ,2 ,3 ,4 ,6 ,8 ]
现象是一致的。
4结语
针对传统岩体动力系统非线性预测理论在观测资料不足的情况下无法计算系统可预测尺度的理论缺陷, 基于胞映射理论建立了一种新的岩体动力系统可预测尺度的计算方法, 成功地解决了这一难题。 由于用该方法计算得到的系统可预测尺度是建立在系统的非线性动力学方程之上的, 因此可以更深刻和科学地刻画和描述系统总体变量的综合作用与整体混沌性。 提出的方法可以计算各种岩体动力系统的可预测尺度, 具有重要的理论和应用价值, 是对岩体非线性动力学理论的进一步发展和完善。
图2 条件Ⅰ Shannon信息量 (a) 和状态概率向量的均方误差 (b) 随迭代步数的变化曲线
Fig.2 Curves of Shannon information (a) and mean squared error of state probability vector (b) according to iteration steps on condition Ⅰ
图3 条件Ⅱ Shannon信息量 (a) 和状态概率向量的均方误差 (b) 随迭代步数的变化曲线
Fig.3 Curves of Shannon information (a) and mean squared error of state probability vector (b) according to iteration steps on condition Ⅱ
图4 岩体动力系统不同时段的可预测尺度变化曲线
Fig.4 Curve of predictable time scale of rock mass dynamic system in different time segments
表2 岩体动力系统不同时段的可预测尺度和混沌性
Table 2 Predictable time scale and chaos of rock mass dynamic system in different time segments
Time segments
Predictable time scale/d
Chaos (qualitatively)
3~9
6
0.167
4~10
8
0.125
6~12
4
0.25
7~13
10
0.1
8~14
11
0.091
9~15
8
0.125
图5 矿震系统不同时段的混沌变化曲线
Fig.5 Curve of chaos of rock mass dynamic system in different time segments
参考文献
[1] 谢和平.矿山非线性岩石力学的研究与展望[J].煤炭学报, 1997 (Suppl.) :180-186. XIEHe ping.Researchandoutlookonminenonlinearrockmechanics[J].JournalofChinaCoalSociety, 1997 (Suppl.) :180-186.
[2] 秦四清, 张倬元, 王士天, 等.非线性工程地质学导引[M ].成都:西南交通大学出版社, 1993.109-110. QINSi qing, ZHANGZhuo yuan, WANGShi tian, etal.AnIntroductiontoNonlinearEngineeringGeology[M ].Chengdu:SouthwesternJiaotongUniversityPress, 1993.109-110.
[3] XieH , PariseauWG .Fractalcharacterandmechanismofrockbursts[J].IntJRockMechMinSci, 1993, 30 (4) :343-350.
[4] FengXia ting, SetoM .Fractalstructureofthetimedis tributionofmicrofracturinginrocks[J].GeophysJInt, 1999, 136:257-285.
[5] 陈子林, 周硕愚.蕴震系统前兆场的混沌吸引子及分维[J].地震学报, 1993 (4) :463-469. CHENZi lin, ZHOUShuo yu.Chaoticattractorandfractalofearthquakeprecursors[J].JournalofSeism, 1993 (4) :463-469.
[6] 安镇文, 王林瑛, 姚栋华, 等.大震孕育过程中的混沌特征[J].地震学报, 1992 (4) :393-399. ANZhen wen, WANGLin ying, YAODong hua, etal.Chaosintheprocessofgreatearthquake[J].JournalofSeism, 1992 (4) :393-399.
[7] HirataT , etal.Fractalstructureofspatialdistributionearthquake the two pointcorrelationfunction[J].Geo physJRAstronSoc, 1980, 62:303-320.
[8] 周 辉.矿震孕育过程中的混沌性及非线性预测理论研究[D].沈阳:东北大学, 2000.31-56. ZHOUHui.ChaosandNonlinearForecastingTheoryintheProcessofMineQuake[D].Shenyang:NortheasternUniversity, 2000.31-56.
[9] 周 辉, 冯夏庭, 谭云亮, 等.矿震系统的胞映射突变预测模型[J].中国有色金属学报, 2002, 12 (1) :155-160. ZHOUHui, FENGXia ting, TANYun liang, etal.Cell mapping catastrophicforecastingmodelofminequakesystem[J].TheChineseJournalofNonferrousMetals, 2002, 12 (1) :155-160.
[10] WolfA .Determininglyapunovexponentsfromtimese ries[J].Physica, 1988, 16D :285.
[11] HsuCS .Atheoryofcell to cellmappingdynamicssys tems[J].ASMEJApplMech, 1980, 47:931-93[12] HsuCS , GuttaluRS .Anuntravellingalgorithmforglobalanalysisofdynamicalsystems:anapplicationofcell to cellmappings[J].ASMEJourApplMech, 1980, 47:940-948.
[13] HsuCS .Generalizedtheoryofcell to cellmappingfornonlineardynamicalsystems[J].JourApplMech, 1981, 49:634-642.
[14] HsuCS , GuttaluRS , ZhuWH .Amethodofanalyz inggeneralizedcellmappings[J].ASMEJourApplMech, 1982, 49:885-894.
[15] HsuCS .Aprobabilistictheoryofnonlineardynamicalsystemsbasedoncellstatespaceconcept[J].ASMEJourApplMech, 1982, 49:895-902.
[16] HsuCS .Globalanalysisofdynamicsystemsusingposetsanddigraphs[J].IntJofBifurandChaos, 1995, 5 (4) :1085-1118