基于MATLAB岩体结构面三维模拟与表征单元体研究.docVIP

基于MATLAB岩体结构面三维模拟与表征单元体研究 　　摘 要：岩体中的结构面对岩体等效物理力学参数有重要影响，当岩体体积大于表征单元体体积（REV）时，岩体等效参数趋于稳定。以新疆提依尔金矿为工程背景，使用Matlab软件，采用蒙特卡洛方法进行岩体结构面网络三维模拟，建立了与实际岩体具有相同统计特征的结构面三维模型；通过定义REV指标，综合确定岩体表征单元体为13m×13m×13m，约为最大节理迹长的3.7倍，为后续研究提供了基础数据。 　　关键词：Matlab 蒙特卡洛方法 结构面三维模拟 表征单元体 　　中图分类号：TU452 文献标识码：A 文章编号：1672-3791（2015）07（a）-0093-02 　　研究岩体工程问题前需要合理确定岩体的等效力学参数，然而由于岩体存在复杂的结构面（断层、节理、裂隙等），其物理力学性质表现出不均匀、不连续性。分析大尺度岩体的力学行为时，应当注意一般情况下岩体的物理力学性质随岩体体积增大逐渐趋于稳定，即岩体存在表征单元体积（REV）。目前通常可采用数值分析的方法确定岩体的等效力学参数，国内外不少学者对此进行了深入研究[1-2]。使用数值分析方法确定岩体等效参数时，需要先明确岩体的REV尺度，只有当岩体模型体积大于岩体的REV时，才能保证数值分析法得到的力学参数与岩体实际的力学参数具有等效性。 　　岩体内结构面分布具有随机性，可用概率分布描述[3]。该文以新疆提依尔金矿露天采场岩质边坡为工程背景，根据现场结构面测量数据统计分析所得的结构面产状概率分布及特征参数，利用Matlab软件，运用蒙特卡洛方法建立了与边坡岩体具有相同分布特征的结构面三维模型，然后以此为基础，研究了边坡岩体的REV大小，为后续边坡稳定性研究提供了依据。 　　1 岩体结构面三维模拟 　　该文的研究对象为新疆提依尔金矿露天采场南部边坡岩体。该部分边坡以凝灰岩为主，结构面发育，主要表现为节理，属于Ⅳ级结构面。该矿最终边坡高度228 m，为合理确定边坡岩体的等效参数以便进行下一步研究，需先确定边坡岩体的表征单元体体积（REV），该文通过建立岩体结构面三维模型来研究岩体的REV大小。 　　1.1 结构面三维模拟原理 　　提依尔金矿南部边坡岩体结构面主要表现为节理，因此，该文针对岩体节理建立结构面三维模型。 　　该文使用Baecher模型模 拟岩体中的节理，该模型被广泛用于岩石力学相关分析[4]。该模型假设节理面形状为薄圆盘，在三维坐标系中，可由6个参数确定：中心点坐标（x，y，z），圆盘直径D，产状（倾向、倾角），其中需重点说明的是节理面的圆盘直径和模型内节理面数目的确定。 　　（1）节理面圆盘直径D。 　　圆盘直径D的概率分布可用负指数分布或者对数正态分布来描述，Villaesusa（1991）研究得出节理面直径呈负指数分布时，根据实测节理迹长数据估算节理面直径概率分布参数的公式，可以在工程中直接使用[3]。 　　如果节理面直径的概率分布fD（D）为负指数分布，则有μ=σ=E（D）=πμ′/8；E（D2）=π2μ′2/32。其中，μ是节理面直径的均值；μ′是节理面迹线长度的均值。 　　假设节理面出露迹线即为节理面圆盘的弦长，平均迹长μ′即为平均弦长，则有 　　，其中，r和D表示节理面半径与直径。 　　假设节理面直径D服从分布fD（D），由μ′=0.25πD可得fD（D）=0.25π?f（0.25πD）。其中，f （x）为节理面迹长的概率分布密度函数。该次研究中，实测节理迹长数据均服从负指数分布，因此假设节理面直径D也服从负指数分布，然后根据公式μ=σ=E（D）=πμ′/8，求取分布参数。 　　（2）节理面数目。 　　伍法权等[5]指出，在相对均质区域内，节理面的体密度相对恒定，与节理面方向无关。因此，模拟空间中节理面数目n可以表示为结构面体积密度λv与模拟空间体积V的乘积，即n=λvV。 　　节理面体积密度λv=4λL/[πE（D2）]，其中，λL为节理面线密度[3]。 　　使用节理面间距测量数据的均值μ〃求出λL=1/μ〃后，再通过（1）中的公式求出E（D2），就可得到节理体积密度。 　　如果岩体中有m组节理面，则总体密度λv总=（λv1+λv2+……+λvm）/2π，其中λv1，λv2，……，λvm分别为第1组至第m组节理的体积密度。 　　1.2 蒙特卡洛法建立结构面三维模型 　　在现场勘察中，共得到208条节理测量数据。各组结构面几何参数的概率分布见表1，其中倾向和倾角数据服从正态分布，迹长数据服从负指数分布。 　　根据以上数据，该文使用蒙特卡洛法建立岩体结构面（节理面）三维模型的步骤为：（1） 假设生成节理三维模型的空间为立方体，坐标原点为立方体中心。（2）如前