论文基于混沌优化的高阶段充填体可靠性分析.docx
《论文基于混沌优化的高阶段充填体可靠性分析.docx》由会员分享，可在线阅读，更多相关《论文基于混沌优化的高阶段充填体可靠性分析.docx（14页珍藏版）》请在冰豆网上搜索。
论文基于混沌优化的高阶段充填体可靠性分析
论文—基于混沌优化的高阶段充填体可靠性分析
基于混沌优化的高阶段充填体可靠性分析
作者刘志祥
指导教师李夕兵
摘要对高阶段充填体进行了力学分析～推导了分层充填力学计算公式～用可靠性理论研究了高阶段充填体稳定性。
考虑高阶段充填体可靠性分析的状态函数求导困难～提出了基于混沌优化的可靠性计算方法～为工程中复杂函数计算可靠性指标提供了一种新的方法。
高阶段充填体力学研究表明:
缩短采场长度和增大采场宽度有利于充填体稳定性和降低充填成本。
为了评价高阶段充填体稳定性～分别在实验室配制充填料浆和采场取样试验了尾砂胶结充填体强度～分析了采场原位充填体强度与实验室试验强度的差异。
研究认为:
采用实验室力学参数进行充填设计～最小设计安全系数为1.6，1.8～最小可靠性指标为1.8，2.0比较合理。
用本文方法对安庆铜矿3号高阶段采场充填体进行了可靠性分析～并评价了充填设计的可行性。
关键词:
尾砂胶结充填体;分层充填;可靠性;混沌优化
引言...................................................................................................................11高阶段充填体力学分析..............................................................................1
1.1高阶段充填体上部受力分析................................................................................1
1.2高阶段充填体力学分析.......................................................................................22基于混沌优化的充填体可靠性分析..........................................................43高阶段充填体可靠性分析..........................................................................6
3.1充填体强度试验与合理可靠性指标....................................................................6
3.2安庆铜矿3号采场充填体可靠性分析................................................................7结论...................................................................................................................8致谢...................................................................................................................8参考文献...........................................................................................................9附录.................................................................................................................10
防灾科技学院毕业设计（论文、综合实践报告）
引言
[1]高阶段采矿是一高效采矿技术，大都采用尾砂嗣后充填（矿房胶结充填、矿柱非胶结充填）。
采矿过程中，胶结充填体垂直暴露高度达60,140m，侧向暴露面积在3000,7000?
。
对高阶段充填体稳定性研究通常采用确定性方法（即定值法），该方法得到的安全系数指标体现了充填体稳定性的一个重要方面，但尾砂充填体是一种复杂的非线性
[2,3][4-7]力学介质，实践表明，有诸多因素影响胶结充填体强度，充填体力学参数存在不确定性与随机性，因此用可靠性理论研究高阶段充填体稳定性更符合客观实际。
为了评价高阶段充填体稳定性，笔者分别在实验室配制充填料浆和采场取样试验了尾砂胶结充填体强度，研究了采场原位充填体强度与实验室试验强度的差异，提出了采用实验室力学参数进行充填设计的合理安全系数与可靠性指标。
1高阶段充填体力学分析
1.1高阶段充填体上部受力分析
高阶段采场（矿房）采后用不同配比的尾砂胶结充填，充填接顶后，充填体与围岩形成力学相互作用系统。
根据自然平衡拱理论，充填体上部承受自然平衡拱内矿岩自重压应力（如图1所示）。
设矿体上下盘围岩内摩擦角为，则上下盘围岩移动角为45º+/2。
图1高阶段充填体力学分析
L1顶板岩石自然平衡拱跨度为:
L,L，2Hcot,1
（1）
HL式
（1）中:
为矿房长度;为矿房高度。
1
防灾科技学院毕业设计（论文、综合实践报告）
[8]b1自然平衡拱高度为:
b,L/（2f）,（L/2，Hcot,）/f11
（2）
f式
（2）中:
为顶板岩石的普氏系数。
0作用于充填体上部的顶压近似等于矩形岩柱ABCD的重力，其压应力为:
,,gb,,g（L/2，Hcot,）/f01（3）
g式（3）中:
为顶板岩体的体重;为重力加速度。
1.2高阶段充填体力学分析
如图2所示的高阶段胶结充填体，前面一侧全部暴露，后面一侧为非胶结尾砂
图2高阶段充填体力学分析
HLB充填体，左右两侧与围岩接触。
设胶结充填体长度为，宽度为，暴露高度为。
F0G充填体上部受力为，自重为，在水平方向上有来自非胶结尾砂一侧的侧压力
FTF12T1，与围岩接触两侧有抗剪切阻力，在滑移面上充填体产生下滑力，抗滑力为，
其中:
G,ghBL,12（4）
F,BL,00（5）
2,F,0.5,,hL121（6）
F,（G，F）sin,20（7）
2
防灾科技学院毕业设计（论文、综合实践报告）
,,T,Bh[c，0.5K（,，,）tan,]120（8）
,,,T,BL（C，,cos,tan,）/cos,（9）
h,H,Btan,1（10）
h,H,0.5Btan,2（11）
,,,11式（4）,（11）中:
为胶结充填体容重（有多种配比充填时，为各配比充填
o,,,,,,45，,/2,i,2体容重与其高度的加权平均值）;为胶结充填体滑移角，;为非胶
,,c结尾砂容重;、为胶结充填体与上下盘围岩作用的粘聚力和内摩擦角（有多种配比
c,,,,,,,,,iicc充填时，、为各配比充填体粘聚力、内摩擦角与其高度的加权平均值）;、
,,,,c为滑移经过区各配比充填体粘聚力和内摩擦角（有多种配比充填时，、为滑移经
c,ii,K过区各配比胶结充填体、与其高度的加权平均值）;为胶结充填体侧压系数（有
Ki,K多种配比充填时，为各配比充填体侧压系数与其高度的加权平均值）。
充填体稳定的条件是:
F，Fcos,,T，2T211（12）
将式（4）,（9）代入式（12），并解的不等式，有
,2BLc,,,,,,2,,,,,2（2tan）BhcK，，,,2[]sin（45）cos（45）BLghh，，，，,,,,2012021,,cos（45/2），,22,,,,,,,,,,,2tan2tanBhKBL，2tan2tanBhKBL，,,,,22
（13）
3,cm取安庆铜矿岩体参数（=4.05g/、f=12、=45.5º），用式（3）计算充填体
,01上部受力。
取配比1:
8的尾砂胶结充填体力学参数:
=1.80g/cm3、c=0.171MPa、
3,,,cm,2=38.7º、K=0.20、=1.63g/，=0.23，用式（13）计算，当矿房宽度为15m，
采场长度分别为60、80和100m时，充填体所要求的强度与暴露高度关系曲线如图3
所示;当采场长度为70m，矿房宽度分别为5、15和30m时，充填体所要求的强度与
暴露高度关系曲线如图4所示。
高阶段充填体力学分析结果表明，采场长度越长、宽度越小，所要求的充填体强度
越高。
充填体可暴露高度对采场长度尤其敏感，采场长度由60m增加至100m时，所要
3
防灾科技学院毕业设计（论文、综合实践报告）
求的充填体强度必须增加73.8%，因此降低采场长度可显著降低充填成本。
图3不同采场长度充填体强度与暴露高度关系曲线
图4不同采场宽度充填体强度与暴露高度关系曲线
2基于混沌优化的充填体可靠性分析
0高阶段采场分层充填时，根据式（3）计算充填体上部受应力，用式（13）可计
[,]jj算的各分层充填体应力，设各分层充填体强度为，定义安全储备:
F,g（X,X,X,X,X）,[,],,j12345jj（14）
XXXXX35124式（14）中:
、、、、为相互独立的随机变量，分别对应胶结充填
[,],j,c体粘结力、内摩擦角、侧压系数K、散体尾砂侧压系数及实验室充填体强度。
将式（13）用等号代入式（14），可得高阶段充填体可靠性分析的状态函数。
Fissler提出一种计算可靠性指标的迭代方法，其与标准化变量一起使用，计算较为[9],,Xxxiii方便。
若（i=1,2,…,m）为某一变量，且该变量的均值为，标准差为，对应的
xi标准化变量由下式给出:
4
防灾科技学院毕业设计（论文、综合实践报告）
X,ixix,i,xi（15）
该标准化变量具有均值为0和标准差为1的特性，此时状态方程可表示为:
F,g（x）,g（x,x,x,？
x）123m（16）
[9]Fiessler计算法的步骤见文献。
g（x）高阶段充填体可靠性分析的状态函数比较复杂，对内摩擦角变量求导困难，用Fissler法求解可靠性指标存在一定的局限性。
混沌优化（ChaosOptimization）利用混沌
[10]具有初始值敏感性、内在随机性及遍历性等特性，把混沌变量映射到待寻优的变量区间，采用混沌变量搜索，不用求解状态函数的导数，在全局寻优过程中有较高的搜索效[11]率。
为此本文采用混沌优化与Fiessler法相结合计算可靠性指标。
计算步骤如下:
F,g（X）
（1）建立状态方程;
F,g（x）
（2）根据式（15），把随机变量变成标准化变量，状态方程为;
**xxF,g（x）,0ii（3）采用混沌优化方法，得出标准化变量，使满足。
具体步骤如下:
z0,i（a）设混沌优化的标准化变量有m个，任意设定m个[0,1]区间相异的初值（i=1,2,…,m），代入Logistic迭代方程。
z,uz（1,z）（n,0,1,？
N;z,[0,1]）n，1,in,in,io,i（17）得到m个不同轨迹的混沌变量，置N为一较大的整数。
式（17）中，u为控制参量，u=4时，Logistic映射为[0,1]区间的满映射，且系统处于完全混沌状态。
（b）标准化变量的取值空间为[-1，1]，根据
x,,1，2zn，1,in，1,i（18）
将混沌变量映射到其取值空间。
zg（x）n，1,i（c）采用混沌变量（）进行迭代搜索，并计算每一步迭代的值（设其为F，
g（x）,F1第一步迭代）。
*x,xF,FF,FF,F1n，1,i111（d）如果，置，，转（c）步，继续进行迭代;如果，xn，1,i放弃，转（c）步，继续迭代。
5
防灾科技学院毕业设计（论文、综合实践报告）
FFFmin11（e）如果达到规定精度，计算结束;如果迭代到一定的次数后，保持不变，
*[12]xi把作为当前次优解，将模拟退火策略引入混沌动力学:
*xxykz,，,，（）（）,,niini，，1,1,（19）
ykyk
（1）
（1）（），,,,（20）
y（k）,,式中:
为的衰减因子;、为调节参数。
（f）与步骤c,d相同，采用混沌变量进行迭代，同时加入k循环，并逐步缩小搜索
*xF,g（x）,0i区域。
找到目标函数的最小值即为满足的全局最优解，并记下此时的值。
*xi（4）根据值，用
m*2,,（x）,i,1i（21）
计算可靠性指标。
根据以上步骤，编制了Matlab计算程序。
[9]为了检验上述算法的准确性，分析如下算例:
设系统状态方程F,g（X,X,c）,,tan,,cXX,XX,c121212，式中:
、对应、，=50kN/?
，、的均值分别为100kPa和35º，标准差分别为20kPa和5º，采用Fissler法得出系统可靠性指,,标=1.1560，用本文方法计算可靠性指标=1.1538（计算中N=400，k=150，δ=0.01，
F,miny（0）=0.95，=0.1，=0.2，=0.0001）。
以上分析可见，采用Fissler法与混沌优化相结合方法求解可靠性指标不用得出状态函数的导数，而是直接用混沌变量搜索最优解，适合于复杂函数求解可靠性指标。
该方法的缺点是迭代次数较多，但借助计算机可以很方便地求解，而且编程比较方便。
3高阶段充填体可靠性分析
3.1充填体强度试验与合理可靠性指标
[13-15]实践表明，采场原位充填体强度与实验室试验强度存在差异。
为了研究充填体强度特征，用安庆铜矿分级尾砂与325号普通硅酸盐水泥在实验室配制了1:
4、1:
8、1:
10和1:
12试块各1组（每组6块），所配制试块的浓度与采场充填浓度一致（为
6
防灾科技学院毕业设计（论文、综合实践报告）
74%）。
在标准恒温养护室养护28d后，进行了充填试块抗压强度试验，每组试块重复试验6次，并计算各组试验的均值与标准差，实验室试验结果如表1所示。
同时，在安庆铜矿3号矿房，对已充入采场尚未凝固的充填体取样，放入试模浇注，在井下养护28d后，测试其抗压强度。
由于高阶段采场取样困难，仅在有巷道进入的地点取了两组样品，配比1:
4和1:
10的充填体抗压强度试验结果见附录A中的表1。
从充填体抗压强度试验结果看，充填料浆充到采场后，在脱水过程中，由于少部分水泥浆流失及料浆离析等原因，采场充填体强度比实验室强度低26.53%,34.72%。
xxxx1122设为采场充填体强度，为实验室试验强度，、均为相互独立的正态随
Zmin机变量，定义最小安全储备:
Zxx,,min21（22）
可靠性指标β为
22ˆˆˆˆ,,,，（）/uuss2112（23）
ˆˆˆˆussxxu112122式中:
、分别为x1、x2的均值;、分别为、的标准差。
根据表1试验成果，对于1?
4和1?
10充填体，用式（23）计算，可靠性指标分
ˆˆu/u21别为:
1.78、1.53（平均1.655），安全系数的中值（）分别为:
1.36、1.54（平均1.45）。
在均值的基础上考虑有1.1,1.2的富余，安全系数中值为:
1.60,1.74，可靠性指标为:
1.82,1.99。
综合以上分析结果，在采场取样比较困难的情况下，如果以实验室试验的力学参数作为充填配比设计依据，设计的最小中值安全系数为1.6,1.8，最小可靠性指标为1.8,2.0比较合理。
3.2安庆铜矿3号采场充填体可靠性分析
安庆铜矿采用高阶段大直径深孔嗣后充填采矿法回采矿石。
采场结构参数为（矿房与矿柱相同）:
长70m（为矿体厚度）、宽15m、高100,120m。
采矿工艺为:
矿房采后在采场不同高度上用不同配比的尾砂胶结充填，矿柱采后用尾砂非胶结充填。
矿柱回采时两侧充填体垂直暴露高度100,120m，充填体稳定性直接关系到采矿的安全，有必要对高阶段充填体稳定性采用可靠性理论进行研究。
3号矿房所设计的充填配比如图5所示。
充填体力学参数如表1所示，1?
4、1?
8、1?
10、1?
12充填体及散体尾砂的侧压
7
防灾科技学院毕业设计（论文、综合实践报告）
系数均值分别为0.18、0.20、0.21、0.22、0.23，其变异系数均为0.10。
用各随机变量的
F,[,]/,jjj均值代入，可以求得各层充填体中值安全系数，采用Fissler法与混沌优化相结合方法可求得各层充填体可靠性指标（计算中N=500，k=200，δ=0.01，y（0）=
,0.95，=0.15，=0.3，Fmin=0.001），计算结果见附录A中的表2。
得出可靠性指标β后，根据标准正交分布表查出充填体失效概率，各分层充填体失效概率亦列在附表2。
从附表2可以看出，各分层充填体中值安全系数均大于1.6，可靠性指标均大于1.8，3号矿房充填体在采矿过程中不会发生失稳，该设计是可行的。
图53号采场充填配比
结论
（1）高阶段充填体力学研究表明，缩短采场长和增大采场宽度有利于充填体稳定性。
在高阶段采矿设计中，缩短采场长度可显著降低充填成本。
（2）采场充填体强度与实验室试验结果有一定的差异，在采场取样比较困难的情况下，如果以实验室试验的力学参数为充填配比设计依据，设计最小中值安全系数为1.6,1.8，最小可靠性指标为1.8,2.0比较合理。
（3）Fissler法与混沌优化相结合方法求解系统可靠性指标不用得出状态函数导数，而是直接用混沌变量搜索最优解，编程比较方便，为工程中复杂函数计算可靠性指标提供了一种新的方法。
致谢
感谢李夕兵老师和各位老师对我的精心指导和帮助;同时也感谢同学们为我所做的一切。
经过这次实习和毕业报告的编写，深感老师严谨细致、一丝不苟的作风是我工作、学习的榜样;老师们循循善诱的教导和不拘一格的思路给予我无尽的启迪，再次感谢老
8
防灾科技学院毕业设计（论文、综合实践报告）
师和同学们对我的关心和帮助。
参考文献
[1]刘志祥,李夕兵.充填体变形的混沌时序重构与神经网络预测[J].矿冶工程,2005,25
（1）:
16–19.
[2]李夕兵,刘志祥.基于重构相空间充填体变形规律的灰色预测研究[J].安全与环境学报,2004,
4（6）:
54–57.
[3]刘志祥,李夕兵.尾砂胶结充填体力学试验与损伤研究[J].金属矿山,2004（11）:
22–24.
[4]CHAPPELLB,HUGGINSG.Effectofbackfillstrengthandstiffnessonstopestability[J].AustralasianInstituteofMiningandMetallurgyPublicationSeries,1998,2:
213–217
[5]BENZAAZOUAM,BELEMT,BussiereB.Chemicalfactorsthatinfluencetheperformanceofminesulphidicpastebackfill[J].CementandConcreteResearch,2002,32（7）:
1133–1144.
[6]BENZAAZOUAM,FALLM,BELEMT.Acontributiontounderstandingthehardeningprocessofcementedpastefill[J].MineralsEngineering,2004,17
（2）:
141–152.
[7]LAMESAW,CLARKIH.Theinfluenceofmaterialcompositionandsamplegeometryonthestrengthofcementedbackfill[C]//HASSANIFP,SCOBLEMJ,YUTR,eds.InnovationsinMiningBackfillTechnology.Brookfield（USA）:
A.A.BalkemaPublishers,1989:
89–94.
[8]高磊.矿山岩石力学[M].北京:
机械工业出版社,1987.
[9]祝玉学.边坡可靠性分析[M].北京:
冶金工业出版社,1993.（ZHUYu-xue.Reliabilityanalysisin
slope[M].Beijing:
MetallurgicalIndustryPress,1993.）
[10]CHOIC,LEEJ.Chaoticlocalsearchalgorithm[J].ArtificialLife&Robotics,1998,2
（1）:
41–47.
[11]李兵,蒋慰孙.混沌优化方法及其应用[J].控制理论与应用,1997,14（4）:
613–615.
[12]邹恩,李祥飞,陈建国.混沌控制及其优化应用[M].长沙:
国防科技大学出版社,2002.
[13]OUELLETJ,BIDWELLTJ,SERVANTS.Physicalandmechanicalcharacterisationofpastebackfillbylaboratoryandin-situtesting[C]//BLOSSDM,eds.MINEFILL’98.Brisbane:
AustralasianInstituteof
MiningandMetallurgyPublication,1998:
249–253.
[14]刘志祥,李夕兵.尾砂充填体力学参数混沌优化反分析研究[J].湖南科技大学学报,2004,19（4）:
14–
17.
[15]GURTUNCARG,ADAMSDJ.Determinationoftheinsitumodulusoftherockmassbytheuseofbackfillmeasurements[J].JournaloftheSouthAfricanInstituteofMiningandMetallurgy,1991,91（3）:
81–88.
9
防灾科技学院毕业设计（论文、综合实践报告）
附录
附录A:
表1不同配比充填体力学参数测试结果汇总表
抗压强度/MPa粘结力/kPa内摩擦角/度试验试块容重均值均值配比标准差标准差-3ˆ地点/（g?
cm均值）,ˆ各试块值标准差sˆˆ,,ˆˆss
3.55,4.33,4.88,3.85,1.844.25170.512528038.441.55.891?
44.75,4.15
2.11,1.68,1.45,1.58,1.821.75170.244117129.238.75.191?
8室内1.96,1.73
升级会员 