LITHOSPHERIC THERMAL-RHEOLOGICAL STRUCTURE AND DEEP GEODYNAMICS IN THE YANGTZE RIVER ECONOMIC BELT
-
摘要: 通过Crust 2.0模型构建有限元三维数值模型，以地表观测温度、深部反演温度和地表热流作为约束，计算了长江经济带地区岩石圈温度结构；在温度结构的基础上，通过GPS观测数据得到的地表应变率和选取代表性岩石物性，计算了长江经济带岩石圈流变强度和等效粘滞性系数。结果表明：长江经济带地区岩石圈表现出明显的横向不均匀性，其中四川盆地表现为低温、高强度和高粘滞性的特征，相同深度，四川盆地核心区比周围地块的温度低100~300℃，强度和粘滞性分别比周缘高1~2个数量级；在温度、强度和粘滞性的过渡地带构造活动性比较强烈。结合深部地球物理观测综合分析认为，长江经济带东、西部分别受到太平洋板块深俯冲和印度-欧亚大陆板块碰撞的影响，其深部影响范围分别达到四川盆地的东、西边界，深部动力学过程可能导致了岩石圈的横向不均匀性。Abstract: Constructing the finite three-dimensional numerical model through Crust 2.0, under constraints of the surface temperature, deep inversion temperature and surface heat flow, we have calculated the thermal structure of the continental lithosphere beneath Yangtze River Economic Belt and adjacent regions. On the basis of 3D thermal structure, deriving surface strain rates from GPS observation, and selecting representative rock physical properties, we have calculated the strength and effective viscosity of the continental lithosphere. The results show that there is obvious lateral heterogeneity in the Yangtze River Economic Belt. The Sichuan basin has low temperature, high strength and high viscosity. In the same depth, temperature beneath the Sichuan basin is lower than ambient region for 100~300℃. Strength and viscosity is lower for 1~2 orders of magnitude. The transition zone of temperature, strength and viscosity has strong active tectonics. Combining with deep geophysical observation, the study proposed that the Pacific Plate subduction and India-Eurasia Plate collision played important roles in the deep dynamic process respectively in the eastern and western parts of the Yangtze River Economic Belt. The effective sphere reach the eastern and western boundaries of Sichuan basin, deep dynamic process may lead to the horizontal heterpgeneity of the lithosphere.
-
0. 引言
影响滑坡稳定的因素较多, 且各因素存在复杂性及不确定性, 目前对其研究大多从地形地貌、岩土体性质、岩土体结构面、地下水、地震、工程影响等方面进行分析评价[1~2]。在一定的滑坡工程地质条件下, 有些因素对滑坡失稳影响较大, 而有些则较小, 这就是说, 影响因素对滑坡失稳的贡献是不同的。
滑坡稳定影响因素的敏感性分析就是定量分析影响滑坡稳定性的各因素与滑坡稳定系数之间的相关性, 即分析各因素的变化对于滑坡稳定系数的影响。分析滑坡影响因素的变化与滑坡稳定性的相关关系以及滑坡在不同外界条件下的稳定性, 不但可以找出滑坡失稳的主导因素, 而且根据敏感性分析得出的主要影响因素, 在滑坡治理及优化设计中就可以有针对性的采取相应的整治措施, 使滑坡治理达到安全、经济和有效的目的。因此, 在滑坡及滑坡治理工程中, 影响因素的敏感性分析显得非常重要。
1. 王良寺滑坡概况及其工程地质条件[3]
1.1 滑坡工程概况
王良寺滑坡位于延安市宝塔区桥沟镇东十里铺村, 延河河谷北岸山坡上, 西距延安市中心延河大桥约5.5km, 东临延安飞机场, 南约1km为210国道。
该滑波在延安市宝塔区政府修建的经济适用房场区———又称窑洞小区。窑洞小区投资2亿多元, 共分三期, 其中第一期已交付使用, 在建的第二期工程顺坡台阶式展布, 共六排四列, 多为跃层式别墅房, 少者两层半, 多者四层半, 占地面积100余亩。据现场勘察和走访调查, 第二期工程场地位于一老滑坡体(即王良寺滑坡)上, 其滑坡体体积超过100 × 104m3, 按滑坡规模和物质组成划分, 该滑坡属于大型黄土滑坡。在平场切坡开挖前是稳定的高陡滑坡, 由于第二期工程的大量开挖活动的影响, 改变了滑坡的应力分布, 使滑坡失稳破坏。
1.2 滑坡工程地质条件
1.2.1 地形地貌
滑坡区地貌单元属黄土梁峁斜坡区, 总体由西北向东南延河河谷倾斜, 平均坡降约20.5 %。
1.2.2 地层岩性
根据工程地质测绘及钻探资料, 滑坡岩土体自上而下共分为两个岩土层, 对其特征分述如下:
① 层新黄土(Q3eol):厚15~25m, 滑坡体的主要组成部分, 土质较均, 岩性为粉质粘土, 褐黄~黄褐色, 硬塑~坚硬, 分布有棕褐色古土壤条带, 多含有钙质结核, 直径3~8cm, 局部呈富集状, 孔隙不发育。
② 层老黄土(Q1-2eol):厚20~25m, 多呈黄褐色, 土质均一, 可塑~硬塑, 针状孔隙发育, 含较多钙质结核, 夹杂有棕褐色古土壤层。
1.2.3 水文地质条件
据钻探及井探资料显示, 滑坡区为贫水区, 地下水位埋深较大, 在滑体前缘的井探中发现地下水(据地面8m)。由此可见, 坡体大部分在地下水位以上。地下水主要接受大气降水补给, 以径流方式向延河排泄。
2. 滑坡稳定敏感性分析的计算方法[3]
敏感性分析是系统分析中分析系统稳定的一种方法。对于滑坡系统, 敏感性分析可转化为以滑坡稳定系数FS为考察对象的单指标多因素的显著性分析。本文分析拟采用正交分析法。
单指标多因素的显著性分析可采用线性模型如下:
(1) 式中:β0为常数项; βi为自变量Xi的回归系数; e为随机误差, 服从标准正态分布。
如果在模型中令某些因素的主效应或交互效应为零, 而其余效应的最小二乘估计不受影响, 即与在不假定上述效应为零时所得的估计一致。这保证对每个效应的估计不受到其它效应的影响。则设计矩阵X必须满足如下条件:
(2) 式中:S11、S22、…、Srr都是方阵, 每一块相应于一组效应。
对于某个因素变量Xi对指标Y的显著性次序分析, 不要求做定量结论, 只要求辨明自变量Xi对因变量Y的显著性影响次序。因此无须求解(1)式中的回归系数, 只需按式(2)对各影响因素进行试验设计。此时, 正交试验可满足模型要求。
设A、B、…表示不同的因素; r为各因素水平数; Ai表示因素A的第i水平(i =1, 2, L, r); Xij表示因素j的第i水平的值(i =1, 2, L, r; j =A, B, L)。
在Xij下进行试验得到因素j第i水平的试验结果指标Yij, Yij是服从正态分布的随机变量。在Xij下做了n次试验得到n个试验结果, 分别为Yijk (k =1, 2, L, n)。有计算参数如下:
(3) 式中:Kij为因素j在i水平下的统计参数; n为因素j在i水平下的试验次数; Yijk为因素j在i水平下第k个试验结果指标值。
评价因素显著性的参数为极差Rj, 其计算公式如下:
(4) 极差越大说明该因素的水平改变对试验结果影响也越大。极差最大的因素也就是最主要的因素, 极差较小的因素为较次要的因素, 依此类推。
利用正交分析法可以解决以下几个问题:
(1) 分清各因素对指标影响的主次顺序, 即分清哪些是主要因素, 哪些是次要因素;
(2) 找出优化的设计方案, 即考察的每个因素各取什么水平才能达到试验指标的要求;
(3) 分析因素与指标间的关系, 即当因素变化时指标是怎样变化的, 找出指标随因素的变化规律和趋势。
3. 滑坡稳定敏感性分析
3.1 正交试验设计
王良寺滑坡稳定系数的计算, 采用圆弧条分法中的简化Bishop法[4]。在对其进行正交试验设计时, 考虑的因素主要有粘聚力(c)、内摩擦角(φ)、重度(γ)、地下水位(hw)、坡高(h)等5种参数变化。各影响因素分析均只考虑土体参数的自相关性, 不考虑其的互相关性影响。为了减少由于水平次序引起的系统误差, 各因素水平的次序应随机排列。参数取值范围和按抽签方式确定的因素水平次序如表 1所示。
表 1 各参数取值范围Table 1. The range of each parameter
假设各因素间无交互作用, 可选5因素正交表安排试验。对于4水平5因素正交试验, 最少试验次数为16次, 记为L16 (45)。
3.2 计算结果分析
选用数理统计L16 (45)正交表[5], 按表 1确定的试验方案, 对王良寺滑坡进行以稳定系数为指标的多因素显著性计算分析。正交试验因素水平概化值及计算结果如表 2所示。
表 2 正交试验计算结果Table 2. Results of orthogonality test
对于因素γ, 它的第一水平(γ=16)下的4次计算(表 2中试验号1, 2, 3, 4)求得的稳定系数之和K1j=3.955;第二水平(γ=18)下的4次计算(试验号5, 6, 7, 8)求得的稳定系数之和K2j =3.825;依此类推, 如表 3。
表 3 不同重度值的稳定性系数Table 3. Stability coefficient of different weights
对于因素c, 它的第一水平(c =57)下的4次计算(表 2中试验号1, 5, 9, 13)求得的稳定系数之和K1 j =4.417;第二水平(c =43)下的4次计算(试验号2, 6, 10, 14)求得的稳定系数之和K2j =3.891;依此类推, 如表 4。
表 4 不同粘聚力的稳定性系数Table 4. Stability coefficient of different cohesive forces
用同样的方法计算φ、h、hw相应的稳定系数之和K1j、K2j、K3j、K4j, 如表 5、表 6、表 7。
表 5 不同内摩擦角的稳定性系数Table 5. Stability coefficient of different internal friction angles
表 6 不同坡高的稳定性系数Table 6. Stability coefficient of different slope heights
表 7 不同地下水位的稳定性系数Table 7. Stability coefficient of different groundwater levels
各参数的极差分析如表 8所示。
表 8 各参数极差分析Table 8. Range analysis of parameters
正交试验设计中各因素不同水平之间的搭配是均衡的, 例如根据因素γ的4种不同水平分成的4组计算的每一组计算中, c、φ等其他因素的不同水平皆出现一次。因此γ的K1j、K2j、K3j、K4j之间的差异可以认为主要是γ取不同水平造成的。同理, 哪个因素的极差R大, 就可以认为该因素水平不同对滑坡稳定性产生的影响就大, 影响最大的因素即最敏感因素。从表 8中的极差值可以看出, φ是最敏感的因素, 其他依次为c、hw、γ、h。
4. 滑坡治理措施探讨
王良寺滑坡为一大型黄土老滑坡因不适当开挖引发复活, 滑坡后部宽约130m, 前部宽约270m, 南北长约350m。坡体中部及前部为拟建的经济适用窑小区, 8排24幢房屋, 建成后将有千余人入住, 依据《建筑边坡工程技术规范》 (GB05330-2002), 其安全等级为Ⅱ级。
根据滑坡稳定影响因子敏感性分析结果, 本工程采用上部截排地表水、滑体后部推力段依势适当削坡、下部抗滑桩支挡三种工程相结合的综合治理措施。抗滑桩位置的选择既考虑坡体的整体稳定又能防止坡体的局部剪出。
4.1 截、排水工程
上部截水沟在滑坡后壁后侧, 依滑坡后壁走势布置, 内截面倒梯形, 长度约400m, 排水沟纵向坡率为3 ‰, 水沟材料采用浆砌片石, 等级为M10。
4.2 削坡、护坡工程
削坡工程充分利用黄土局部直立性能好的特点, 分级削坡, 每级坡高为10m, 坡比为1: 1.5, 共5级, 两级削坡间设3m宽平台。部分削坡土方回填于下部抗滑桩后, 既减少土方运距又可适当压脚。护坡工程采用2m ×2m尺寸的菱形网格格构, 格体采用M10浆砌片石构筑, 截面尺寸为300mm×400mm。
4.3 支挡工程
支挡工程采用半悬臂式抗滑桩。设计中采用的土体力学参数为:滑面以上土体重度γ= 16kN/m3, 粘聚力c =31kPa, 内摩擦角φ=23°; 滑体后部滑动带土体重度γ=18kN/m3, 粘聚力c =25kPa, 内摩擦角φ=18°; 滑面以下土体重度为γ=18kN/m3, 粘聚力c =31kPa, 内摩擦角φ=23°。取边坡稳定安全系数为1.25, 采用圆弧滑动法计算, 滑面以上适当考虑桩前土的抗力作用, 桩后剩余下滑力和桩前土抗力均按矩形分布, 滑面下按弹性地基梁计算, 地基水平抗力系数m=2000kN/m4。设计抗滑桩截面尺寸为2.5m×3.5m, 中至中间距为6m。
5. 结论
采用概率统计原理和计算机电算技术的正交分析法对滑坡稳定影响因子进行敏感性分析, 其工作量小、易于完成, 具有一定的通用性, 对工程技术人员来说是一种较实用的分析方法。实例分析表明, 运用正交分析法得到的结论与一般工程经验相符, 因此该方法具有一定的可靠性。
本文在对滑坡稳定性影响因子进行敏感性分析时, 没有考虑参数的相关性对其稳定性的影响, 在今后的工作可以考虑改进正交分析法以进行因素间有交互作用时的分析。实际上岩土体参数c、φ通常呈负相关性, 考虑参数的相关性可以使稳定性分析模型更接近于工程实际, 提高对滑坡稳定系数的计算精度, 以期达到安全、可靠, 经济、合理的治理效果。
-
图 2 长江经济带地区岩石圈不同深度温度场
Figure 2. Temperature distribution at different depths of the lithosphere in the Yangtze River Economic Belt
图 3 长江经济带地区各个剖面岩石圈温度、强度和等效粘滞性系数分布(剖面位置见图 1)
Figure 3. Distributions of temperature, strength and viscosity of the lithosphere at different profiles in the Yangtze River Economic Belt
图 4 长江经济带地区岩石圈不同深度的强度分布
Figure 4. Lithospheric strength distribution at different depths in the Yangtze River Economic Belt
图 5 长江经济带地区岩石圈不同深度的等效粘滞性系数分布
Figure 5. Viscosity distribution at different lithospheric depths in the Yangtze River Economic Belt
-
[1] 舒良树.华南构造演化的基本特征[J].地质通报, 2012, 31(7):1035~1053 http://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD201207004.htmSHU Liang-shu. An analysis of principal features of tectonic evolution in South China Block[J]. Geological Bulletin of China, 2012, 31(7):1035~1053. http://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD201207004.htm [2] 张国伟, 郭安林, 王岳军, 等.中国华南大陆构造与问题[J].中国科学:地球科学, 2013, 43(10):1553~1582 http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201310003.htmZHANG Guo-wei, GUO An-lin, WANG Yue-jun, et al. Tectonics of South China continent and its implications[J]. Science China:Earth Sciences, 2013, 43(10):1553~1582. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201310003.htm [3] 张岳桥, 徐先兵, 贾东, 等.华南早中生代从印支期碰撞构造体系向燕山期俯冲构造体系转换的形变记录[J].地学前缘, 2009, 16(1):234~247 http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200901033.htmZHANG Yue-qiao, XU Xian-bing, JIA Dong, et al. Deformation record of the change from Indosinian collision-related tectonic system to Yanshanian subduction-related tectonic systemin South China during the Early Mesozoic[J]. Earth Sciences Frontiers, 2009, 16(1):234~247. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200901033.htm [4] Gan W, Zhang P, Shen Z-K, et al. Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements[J]. Journal of Geophysical Research, 2007, 112:B08416. doi: 10.1029/2005JB004120/references [5] Kreemer C, Blewitt G, Klein E C. A geodetic plate motion and Global Strain Rate Model[J]. Geochemistry, Geophysics, Geosystems, 2014, 15(10):3849~3889. doi: 10.1002/2014GC005407 [6] Deng Q, Zhang P, Ran Y, et al. Basic characteristics of active tectonics of China[J]. Science in China:Earth Sciences, 2003, 46(4):356~372. doi: 10.1360%2F03yd9032.pdf [7] Goes S, Govers R, Vacher P. Shallow mantle temperatures under Europe from P and S wave tomography[J]. Journal of Geophysical Research, 2000, 105(B5):11153~11169. doi: 10.1029/1999JB900300 [8] Goes S, Lee S v d. Thermal structure of the North American uppermost mantle inferred from seismic tomography[J]. Journal of Geophysical Research, 2002, 107(B3):2050. doi: 10.1029/2000JB000049 [9] An M, Shi Y. Lithospheric thickness of the Chinese continent[J]. Physics of the Earth and Planetary Interiors, 2006, 159(3/4):257~266. http://www.sciencedirect.com/science/article/pii/S0031920106002494 [10] An M, Shi Y. Three-dimensional thermal structure of the Chinese continental crust and upper mantle[J]. Science in China:Series D, 2007, 50(10):1441~1451. doi: 10.1007/s11430-007-0071-3 [11] Hu S, He L, Wang J. Heat flow in the continental area of China:a new data set[J]. Earth and Planetary Science Letters, 2000, 179:407~419. doi: 10.1016/S0012~821X(00)00126~6 [12] Sun Y, Dong S, Zhang H, et al. 3D thermal structure of the continental lithosphere beneath China and adjacent regions[J]. Journal of Asian Earth Sciences, 2013, 62:697~704. doi: 10.1016/j.jseaes.2012.11.020 [13] Kohlstedt D L, Evans B, Mackwell S J. Strength of the lithosphere:Constraints imposed by laboratory experiments[J]. Journal of Geophysical Research, 1995, 100(B9):17587~17602. doi: 10.1029/95JB01460 [14] 臧绍先, 李昶, 魏荣强.岩石圈流变机制的确定及影响岩石圈流变强度的因素[J].地球物理学进展, 2002, 17(1):50~60 http://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ200201007.htmZANG Shao-xian, LI Chang, WEI Rong-qiang. The determination of the deformational mechanism of the lithosphere and the factors influenced the lithospheric rheological strength[J]. Progress in Geophysics, 2002, 17(1):50~60. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ200201007.htm [15] Zang S X, Wei R Q, Liu Y G. Three-dimensional rheological structure of the lithosphere in the Ordos block and its adjacent area[J]. Geophysical Journal International, 2005, 163(1):339~356. doi: 10.1111/gji.2005.163.issue-1 [16] Zang S X, Wei R Q, Ning J Y. Effect of brittle fracture on the rheological structure of the lithosphere and its application in the Ordos[J]. Tectonophysics, 2007, 429:267~285. doi: 10.1016/j.tecto.2006.10.006 [17] 魏荣强, 臧绍先.岩石破裂强度的温度和应变率效应及其对岩石圈流变结构的影响[J].地球物理学报, 2006, 49(6):1730~1737 http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200606019.htmWEI Rong-qiang, ZANG Shao-xian. Effects of temperature and strain rate on the fracture strength of rock and their influences on the rheological structure of the lithosphere[J]. Chinese Journal of Geophysics, 2006, 49(6):1730~1737. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200606019.htm [18] Shimada M. Lithosphere strength inferred from fracture strength of rocks at high confining pressures and temperatures[J]. Tectonophysics, 1993, 217:55~64. doi: 10.1016/0040-1951(93)90202-U [19] Kirby S H. Rheology of the lithosphere[J]. Reviews of Geophysics and Space Physics, 1983, 21(6):1458~1487. doi: 10.1029/RG021i006p01458 [20] Ranalli G, Murphy D C. Rheological stratification of the lithosphere[J]. Tectonophysics, 1987, 132:281~295. doi: 10.1016/0040-1951(87)90348-9 [21] Weertman J. The creep strength of the Earth's mantle[J]. Reviews of Geophysics and Space Physics, 1970, 8(1):145~168. doi: 10.1029/RG008i001p00145 [22] Kirby S H, Kronenberg A K. Rheology of the lithosphere:Selected topics[J]. Review of Geophysics, 1987, 25(6):1219~1244. doi: 10.1029/RG025i006p01219 [23] Zhu S, Shi Y. Estimation of GPS strain rate and its error analysis in the Chinese continent[J]. Journal of Asian Earth Sciences, 2010, 40(1):351~362. http://www.escience.cn/system/download/71645 [24] 孙玉军, 董树文, 范桃园, 等.中国大陆及邻区岩石圈三维流变结构[J].地球物理学报, 2013, 56(9):2936~2946 http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201309008.htmSUN Yu-jun, DONG Shu-wen, FAN Tao-yuan, et al. 3D rheological structure of the continental lithosphere beneath China and adjacent regions[J]. Chinese Journal of Geophysics, 2013, 56(5):546~558. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201309008.htm [25] Huang J, Zhao D. High-resolution mantle tomography of China and surrounding regions[J]. Journal of Geophysical Research, 2006, 111(B9):B09305. http://www.oalib.com/references/15772398 [26] Wei W, Xu J, Zhao D, et al. East Asia mantle tomography:New insight into plate subduction and intraplate volcanism[J]. Journal of Asian Earth Sciences, 2012, 60:88~103. doi: 10.1016/j.jseaes.2012.08.001 期刊类型引用(6)
1. 王建强，梁杰，陈建文，赵青芳，张银国，董贺平. 中国近海基岩油气藏分布特征及其控制因素. 科学技术与工程. 2022(03): 863-872 . 
百度学术2. 胡贺伟，李慧勇，肖述光，李俊伟，丁熠然. 沙垒田凸起西段走滑断裂发育特征及其对油气的控制作用. 海洋地质前沿. 2022(03): 36-44 . 百度学术3. 关蕴文，周锋，蒲仁海，范昌育. 西湖凹陷平北缓坡带盖层特征及封堵性评价. 海洋地质前沿. 2022(10): 34-41 . 百度学术4. 蔡益栋，高国森，刘大锰，邱峰. 鄂尔多斯盆地东缘临兴中区煤系气富集地质条件及成藏模式. 天然气工业. 2022(11): 25-36 . 百度学术5. 康海亮，林畅松，牛成民. 渤海西部沙东南构造带东营组古地貌特征及对沉积的控制作用. 地质力学学报. 2021(01): 19-30 . 
本站查看6. 吴迅达，廖晋，孙文钊，刘平，李春雷. 莺歌海盆地天然气运聚成藏条件与分布富集规律. 地质力学学报. 2021(06): 963-974 . 本站查看其他类型引用(6)
-
下载:
下载:
百度学术
计量
- 文章访问数: 319
- HTML全文浏览量: 150
- PDF下载量: 27
- 被引次数: 12
