《岩土工程学报论文模板.docx》由会员分享，可在线阅读，更多相关《岩土工程学报论文模板.docx（6页珍藏版）》请在课桌文档上搜索。
1、岩土工程学报论文模板刘志祥，李夕兵，张义平(中南大学，湖南长沙410083)摘要：对高阶段充填体进行了力学分析，推导了分层充填力学计算公式，用可靠性理论研究了高阶段充填体稳固性。考虑高阶段充填体可靠性分析的状态函数求导困难，提出了基于混沌优化的可靠性计算方法，为工程中复杂函数计算可靠性指标提供了一种新的方法。高阶段充填体力学研究说明：缩短采场长度与增大采场宽度有利于充填体稳固性与降低充填成本。为了评价高阶段充填体稳固性，分别在实验室配制充填料浆与采场取样试验了尾砂胶结充填体强度，分析了采场原位充填体强度与实验室试验强度的差异。研究认为：使用实验室力学参数进行充填设计，最小设计安全系数为L61.
2、8,最小可靠性指标为L82.0比较合理。用本文方法对安庆铜矿3号高阶段采场充填体进行了可靠性分析，并评价了充填设计的可行性。关键词：尾砂胶结充填体；分层充填；可靠性；混沌优化中图分类号：TD853.34文献标识码：A文章编号：1000-4548(2006)03-0348-05作者简介：刘志祥(1967-),男，湖南宁乡人，中南大学博士后，从事采矿与岩石力学研究。Reliabilityanalysisofhighlevelbackfillbasedonchaoticoptimization1.IUZhi-xiang,LIXi-bing,ZHANGYi-ping(CentralSouthUnive
3、rsity.Changsha410083.China)Abstract:Tliroughthemechanicalanalysisofhighbackfill,thecalculationformulaeofthestratifiedbackfillpracticewerededucedandreliabilityanalysisoftheirstabilitieswerestudied.Inreliabilityanalysisofhighbackfill,becausethedifferentialofstatusfunctionwasdifficulttobeobtained,acalc
4、ulatingmethodofreliabilitybasedonchaoticoptimizationwasproposed,whichwasanewmethodofreliabilityanalysisforcomplexstatusfunctioninengineeringpractice.Themechanicalresearchesshowedthatreducinglengthorincreasingwidthofstopewerepropitioustoimprovestabilityofbackfillandtolowerfillingcost.Inordertoevaluat
5、ethestabilityofbackfill,aseriesofstrengthexperimentsofcementedtailingsbackfillbothintheIaboratorjfandunderfillingstopeconditionsweredone,andtheirdifferencesofstrengthwereanalyzed,aswellasaconclusionwasdrawnthatifthemechanicalparametersinlabwereusedasthebasisoffillingdesign,theminimumsafetyfactorshou
6、ldbeL6to1.8,andtheminimumindexofreliabilityshouldbe1.8to2.0.Usingthepresentmethod,thereliableindexesofbackfillinstopeNo.3inAnqingcoppermine,Anhuiprovince,werecalculated,andthefeasibilityoffillingdesignwereevaluated.Keywords:consolidatedtailingsbackfill;stratifiedfilling;reliability;chaoticoptimizati
7、on0引言高阶段采矿是一高效采矿技术，大都使用尾砂嗣后充填（矿房胶结充填、矿柱非胶结充填）。采矿过程中，胶结充填体垂直暴露高度达60140m,侧向暴露面积在30007000m2o对高阶段充填体稳固性研窕通常使用确定性方法（即定值法），该方法得到的安全系数指标表达了充填体稳固性的一个重要方面，但尾砂充填体是一种复杂的非线性力学介质叫实践说明“刀，有诸多因素影响胶结充填体强度，充填体力学参数存在不确定性与随机性，因此用可靠性理论研究高阶段充填体稳固性更符合客观实际。为了评价高阶段充填体稳固性，笔者分别在实验室配制充填料浆与采场取样试验了尾砂胶结充填体强度，研究了采场原位充填体强度与实验室试验强度的
8、差异，提出了使用实验室力学参数进行充填设计的合理安全系数与可靠性指标。1高阶段充填体力学分析1.1高阶段充填体上部受力分析高阶段采场（矿房）采后用不一致配比的尾砂胶结充填，充填接顶后，充填体与围岩形成力学相互作基金项目：XXXXXXXXXXXXXXXX(编号XXXXXXXX)收稿日期：XXXX-XX-XX用系统。根据自然平衡拱理论，充填体上部承受自然平衡拱内矿岩自重压应力(如图1所示)。设矿体上下盘围岩内摩擦角为夕，则上下盘围岩移动角。为45+2o图1高阶段充填体力学分析Fig. 1 Mechanicalanalysisofhighback11ll顶板岩石自然平衡拱跨度Ly为Z1=L+27/c
9、ot0(1)式中L为矿房长度；H为矿房高度。自然平衡拱高度也为U=LJQf)=(Ll2+Hcot/于,(2)式中/为顶板岩石的普氏系数。作用于充填体上部的顶压近似等于矩形岩柱ABCO的重力，其压应力5)为a=rsh=g(L2+Hcot0)f，(3)式中脑顶板岩体的体重；g为重力加速度。1.2高阶段充填体力学分析如图2所示的高阶段胶结充填体，前面侧全部暴露，后面一侧为非胶结尾砂充填体，左右两侧与围岩接触。设胶结充填体长度为Lf宽度为8,暴露高度为从图2高阶段充填体力学分析Fig. 2 Mechanicalanalysisofhighbackfill充填体上部受力为6，自重为G,在水平方向上有来自
10、非胶结尾砂侧的侧压力K，与围岩接触两侧有抗剪切阻力(，在滑移面上充填体产生下滑力F2,抗滑力为T,其中：G=gh2BL,(4)4=BL,，(5)Fi=0.512L,(6)6=(G+Ej)Sina,(7)T=Bl2c+0.5X(+0)tan,(8)T=BLe+costan夕)/COSa,(9)hy=H-Btana,(10)h2=H-0.5Btanzo(11)式中y：为胶结充填体容重(有多种配比充填时，为各配比充填体容重九与其高度的加权平均值)；a为胶结充填体滑移角，1=45+0/2；%为非胶结尾砂容重；。、”为胶结充填体与上下盘围岩作用的粘聚力与内摩擦角(有多种配比充填时，c、”为各配比充填体粘
11、聚力G、内摩擦角%与其高度的加权平均值)；c、为滑移通过区各配比充填体粘聚力与内摩擦角(有多种配比充填时，/、“为滑移通过区各配比胶结充填体q、%与其高度的加权平均值)；K为胶结充填体侧压系数(有多种配比充填时，K为各配比充填体侧压系数Ki与其高度的加权平均值)。充填体稳固的条件是5+4CoSaT+27JO(12)将式(4)(9)代入式(12),并解b的不等式，有2BLl+0sin(45+)+月CoS(45+)Z：2Bh2Ktan+2BLtan2B%(2c+K0tan,)+-cos(45+2)、o2Bh2K,tan,+2BLtan取安庆铜矿岩体参数(y=4.05gCm3、/=12、-45.5)
12、,用式(3)计算充填体上部受力b。取配比1：8的尾砂胶结充填体力学参数：l=1.80g/cn?、c=0.171MPa、e=38.7。、K=O.20、6=1.63g/Cm，=0.23,用式(13)计算，当矿房宽建为15m,采场长度分别为60、80与100m时，充填体所要求的强度与暴露高度关系曲线如图3所示;当采场长度为70m,矿房宽度分别为5、15与30m时，充填体所要求的强度与暴露高度关系曲线如图4所示。Udwa甯60 7() 0 90 I(X) IlO 120充填体暴褥四度/m图3不一致采场长度充填体强度与暴露高度关系曲线Fig.3Curvesbetweenrequiredstrength&
13、exposureheightofbackfillinstopesofdifferentlengths高阶段充填体力学分析结果说明,采场长度越长、宽度越小，所要求的充填体强度越高。充填体可暴露高度对采场长度特别敏感，采场长度由60m增加至100m时，所要求的充填体强度务必增加73.8%,因此降低采场长度可显著降低充填成本。图4不一致采场宽度充填体强度与暴露高度关系曲线Fig.4Curvesbetweenrequiredstrength&exposureheightofbackfillinstopesofdifferentwidths2基于混沌优化的充填体可靠性分析高阶段采场分层充填时，根据式(3
14、)计算充填体上部受应力气,用式(13)可计算的各分层充填体应力%,设各分层充填体强度为%,定义安全储备.=g(X1,X2,X3,X4,X5)=J-.,(14)式中：X、X2.X3、X4、X5为相互独立的随机变量，分别对应胶结充填体粘结力C、内摩擦角、侧压系数K、散体尾砂侧压系数及实验室充填体强度.。将式(13)用等号代入式(14),可得高阶段充填体可靠性分析的状态函数。FiSSIer提出一种计算可靠性指标的迭代方法，其与标准化变量一起使用，计算较为方便。若X#=1,2,用)为某一变量，且该变量的均值为j,标准差为,对应的标准化变量无由下式给出：该标准化变量具有均值为。与标准差为1的特性,如今状
15、态方程可表示为=g(x)=g(xvx2.x3,)。(16)Fiessler计算法的步骤见文献9。高阶段充填体可靠性分析的状态函数比较复杂，g(x)对内摩擦角变量求导困难，用FiSSlei法求解可靠性指标存在一定的局限性。混沌优化(ChaOSOptimization)利用混沌具有初始值敏感性、内在随机性及遍历性等特性I，把混沌变量映射到待寻优的变量区间，使用混沌变量搜索，不用求解状态函数的导数，在全局寻优过程中有较高的搜索效率WL为此本文使用混沌优化与FieSSler法相结合计算可靠性指标。计算步骤如下：(1)建立状态方程尸=g(X);(2)根据式(15),把随机变量变成标准化变量，状态方程为产
16、=g(x)；(3)使用混沌优化方法,得出标准化变量片，使x：满足尸=g(x)=00具体步骤如下：(a)设混沌优化的标准化变量有m个，任意设定m个0,1区间相异的初值Z(M(i=1,2,第7),代入Logistic迭代方程。ZgU=UZni(I-Zni)(=Oj,N;ZOj0,l),(17)得到m个不一致轨迹的混沌变量，置N为一较大的整数。式(17)中，为操纵参量，=4时，LOgiStiC映射为。1区间的满映射，且系统处于完全混沌状态。(b)标准化变量的取值空间为T,1,根据+u=T+2zn/，(18)将混沌变量映射到其取值空间。(C)使用混沌变量(z,ii)进行迭代搜索，并计算每一步迭代的g(
17、x)值(设其为F,第一步迭代g(x)=)。(d)假如/耳，置耳=尸，K=/.),转(C)步，继续进行迭代；假如尸E,放弃乙+以，转(C)步，继续迭代。(e)假如耳达到规定精度%m,计算结束;假如迭代到一定的次数后，耳保持不变，把阳作为当前次优解，将模拟退火策略引入混沌动力学【：Xft+u=X；+yW(-ct+zni),(19)义攵+1)=(1-Ry(Q,(20)式中6为Nk)的衰减因子；。、为调节参数。与步骤cd相同，使用混沌变量进行迭代，同时加入左循环，并逐步缩小搜索区域。找到目标函数的最小值即为满足F=g(x)=0的全局最优解，并记下如今的匕值。(4)根据：值，用TE2(2D计算可靠性指标
18、o根据以上步骤，编制了MaUab计算程序。为了检验上述算法的准确性，分析如下算例：设系统状态方程F=g(X,X2,c)=6an-c,式中：X、X2对应。、,c=50kNm2,XI、X2的均值分别为100kPa与35,标准差分别为20kPa与50,使用FiSSIer法得出系统可靠性指标P=1.1560,用本文方法计算可靠性指标夕=1.1538(计算中N=400,k=150,5=0.01,y(0)=0.95,a=0.1,77=0.2,Fmin=0.0001)以上分析可见，使用FiSSler法与混沌优化相结合方法求解可靠性指标不用得出状态函数的导数，而是直接用混沌变量搜索最优解，适合于复杂函数求解可
19、靠性指标。该方法的缺点是迭代次数较多，但借助计算机能够很方便地求解，而且编程比较方便。表1不一致配比充填体力学参数Table1Mechanicalparametersofbackfillwithdifferenttailing-cementratios试验地点配比试块容重/(gcm3)抗压强度/MPa粘结力/kPa内摩擦角/度各试块值均值标准差S均值a标准差S均值标准差S1:41.843.55,4.33,4.88,3.85,4.75,4.154.25170.512528038.441.55.89室内1:81.822.11,1.68,1.45,1.58,1.96,1.731.75170.2441
20、17129.238.75.19试验1：101.791.05,1.62,1.15,1.21,1.52,1.081.27150.239814330.136.95.251：121.770.66,0.95,0.82,0.85,0.75,0.830.81(X)0.097810512.733.24.01充填1:41.832.82,2.95,2.69,2.76,3.35.3.75,3.38,3.293.1240.374采场1：101.780.78,0.65,0.78,0.87,0.68,0.95,1.15,0.780.8300.1617号桥7号房5弓桥5号房-.1 尾砂胶结充填f.,*:?足砂非胶结充填 ;
21、hw1:12 20 m1:10 26m1:8 8m3高阶段充填体可靠性分析3.1 充填体强度试验与合理可靠性指标实践说明，采场原位充填体强度与实验室试验强度存在差异1。为了研究充填体强度特征，用安庆铜矿分级尾砂与325号普通硅酸盐水泥在实验室配制了1：4、1：8、1：10与1：12试块各1组（每组6块），所配制试块的浓度与采场充填浓度一致（为74%）在标准恒温养护室养护28d后，进行了充填试块抗压强度试验，每组试块重复试验6次，并计算各组试验的均值与标准差，实验室试验结果如表1所示。同时，在安庆铜矿3号矿房，对已充入采场尚未凝固的充填体取样，放入试模浇注，在井下养护28d后，测试其抗压强度。由
22、于高阶段采场取样困难，仅在有巷道进入的地点取了两组样品，配比1：4与1：10的充填体抗压强度试验结果亦列在表1。从表1可看出：充填料浆充到采场后，在脱水过程中，由于少部分水泥浆流失及料浆离析等原因，采场充填体强度比实验室强度低26.53%34.72%。设R为采场充填体强度，X2为实验室试验强度，总、M均为相互独立的正态随机变量，定义最小安全储备Zmin:Zmin=七一，（22）可靠性指标S为=（）+：，Q3）式中M1X心分别为即、X2的均值；S1a冬分别为为、M的标准差。根据表1试验成果，关于1：4与1：10充填体,用式（23）计算，可靠性指标分别为：1.78、1.53（平均1.655）,安全
23、系数的中值（u2iil）分别为：1.36、1.54（平均1.45）在均值的基础上考虑有1.11.2的富余，安全系数中值为：1.601.74,可靠性指标为:1.821.99.综合以上分析结果，在采场取样比较困难的情况下，假如以实验室试验的力学参数作为充填配比设计根据，设计的最小中值安全系数为161.8,最小可靠性指标为L82.0比较合理。3.2 安庆铜矿3号采场充填体可靠性分析安庆铜矿使用高阶段大直径深孔嗣后充填采矿法回采矿石。采场结构参数为（矿房与矿柱相同）：长70m（为矿体厚度）、宽15m、高100120m采矿工艺为：矿房采后在采场不一致高度上用不一致配比的尾砂胶结充填，矿柱采后用尾砂非胶结
24、充填。矿柱回采时两侧充填体垂直暴露高度100120m,充填体稳固性直接关系到采矿的安全，有必要对高阶段充填体稳固性使用可靠性理论进行研究。3号矿房所设计的充填配比如图5所示。充填体力学参数如表1所示，1：4、1：8、1：10、1：12充填体及散体尾砂的侧压系数均值分别为0.18、0.20、0.21、0.22、0.23,其变异系数均为0.10。用各随机变量的均值代入I=%/?，能够求得各层充填体中值安全系数，使用FiSSler法与混沌优化相结合方法可求得各层充填体可靠性指标（计算中N=500,k=200,S=O.01,y（0）=0.95,=0.15,77=0.3,Fmin=0.001）,计算结果
25、如表2所示。得出可靠性指标后，根据标准正交分布表查出充填体失效概率，各分层充填体失效概率亦列在表2o从表2能够看出，各分层充填体中值安全系数均大于1.6,可靠性指标均大于1.8,3号矿房充填体在采矿过程中不可能发生失稳，该设计是可行的。配比与岛度1:44m1:88m1:1016m三1:46m*1:814m-图53号采场充填配比Fig.5Cement-tailingratiosinstopeNo.3表23号矿房配比设计与可靠性计算Table2Thedesignofcement-tailingratiosandreliabilitycalculationinstopeNo.3序号分层高/m高程/m
26、配比水泥用量中值安全系数可苑性指标我失效概率/%1441：41092.021.810.502124-161：81820.06.898.13032016-361：122100.02.193.760.00941636-521：101985.52.787.8605652-581：41638.08.549.56061458-721：82123.33.365.93072672-981：103226.42.212.231.298898-1061：81213.32.785.060.000396106-1121：41638.06.438.560合计11216835.74结语(1)高阶段充填体力学研究说明，缩短
27、采场长与增大采场宽度有利于充填体稳固性。在高阶段采矿设计中，缩短采场长度可显著降低充填成本。(2)采场充填体强度与实验室试验结果有一定的差异，在采场取样比较困难的情况下，假如以实验室试验的力学参数为充填配比设计根据，设计最小中值安全系数为1.6-1.8,最小可靠性指标为1.82.0比较合理。(3 Fissler法与混沌优化相结合方法求解系统可靠性指标不用得出状态函数导数，而是直接用混沌变量搜索最优解，编程比较方便，为工程中复杂函数计算可靠性指标提供了一种新的方法。参考文献：|1|刘志祥,李夕兵.充填体变形的混沌时序重构与神经网络预测J.矿冶工程，2005,25(1):16-19.(LlUZhi
28、-xiang,LIXi-bing.Reconstructionofchaotictimeseriesforbackfilldeformationandpredictionwithneuralnetwork!J.MiningandMetallurgicalEngineering,2005,25(1):16-19.)2李夕兵,刘志祥.基于重构相空间充填体变形规律的灰色预测研究UL安全与环境学报，2004,4(6):54-57.(LIXi-bing,LIUZhi-xiang.Researchongreypredictionofdeformationlawsinbackfillbasedonphase
29、spacereconstructionJ.JournalofSafetyandEnvironment,2004,4(6):54-57.)3刘志祥,李夕兵.尾砂胶结充填体力学试验与损伤研窕国.金属矿山，2004(11):22-24.(LIUZhi-xiang,LIXi-bing.MechanicstestanddamagestudyoncementedtailingfillingJ.MetalMine,2004(11):22-24.)(4 CHAPPELLB,HUGGINSGEffectofbackfillstrengthandstiffnessonstopestabilityJ.Austral
30、asianInstituteOfMiningandMetallurgyPublicationSeries,1998,2:213-217151BENZAAZOUAM.BELEMT,BUSSiereB.ChemicalfactorsthatinfluencetheperformanceofmineSulphidicpastebackfill(JJ.CementandConcreteResearch,2002,32(7):1133-1144.6 BENZAAZOUAM.FALLM,BELEMT.Acontributiontounderstandingthehardeningprocessofceme
31、ntedpastefillJ.MineralsEngineering,24,17(2):141-152.7 1.AMESAW,CLARKIH.TheinfluenceofmaterialcompositionandsamplegeometryonthestrengthofcementedbackfillC/HASSANIFP.SCOBLEMJ,YUTR.eds.InnovationsinMiningBackHlITechnology.BrookfieId(USA):A.A.BalkemaPublishers,1989:89-94.8高磊.矿山岩石力学M.北京：机械工业出版社,1987.(GAO
32、Lei.RockmechanicsinmineM.Beijing:MechanicalIndustryPress,1987.)9祝玉学.边坡可靠性分析M.北京：冶金工业出版社.1993.(ZHUYu-xue.ReliabilityanalysisinslopeM.Beijing:MetallurgicalIndustryPress,1993.)10 CHOIC,LEEJ.ChaoticlocalsearchalgorithmJ.ArtificialLife&Robotics,1998,2(1):41-47.11李兵,蒋慰孙.混沌优化方法及其应用J.操纵理论与应用，1997,14(4):613-
33、615.(LIBing,JIANGWei-sun.ChaosoptimizationmethodanditsapplicationJ.ControlTheor),andApplications,1997,14(4):613-615.)(12邹恩,李祥飞,陈建国.混沌操纵及其优化应用M.长沙:国防科技大学出版社，2002.(ZOUEn,LIXiang-fei.CHENJian-guo.ChaoticcontrolanditsapplicationinoptimizationM.Changsha:NationalDefenseUniversityofScienceandTechnologyPres
34、s,2002.)13 OUELLETJ,BIDWELLTJ,SERVANTS.PhysicalandmechanicalcharacterisationofpastebackfillbyIaboratorjfandin-situtesiingC/BLOSSDM,eds.MINEFILL,98.Brisbane:AustralasianInstituteofMiningandMetallurgyPublication.1998:249-253.14刘志祥,李夕兵.尾砂充填体力学参数混沌优化反分析研究J.湖南科技大学学报,2004,19(4):14-17.(LIUZhi-xiang,LIXi-bi
35、ng.ResearchonbackanalysisformechanicalparametersoftailingsbackIIwithchaoticoptimizationJ.JournalofHunanUniversityofScience&Technology(NaturalScience),2004,19(4):14-17.)15 GURTUNCArC.ADAMSDJ.DeterminationoftheinsitumodulusoftherockmassbytheuseofbackfillmeasurementsJ.JournaloftheSouthAfricanInstituteofMiningandMetallurgy,1991,91(3):81-88.
