吕梁山麓不同岩组段小流域地貌特征

岳科杉,赵金贵,王建民

(太原理工大学 矿业工程学院,太原 030024)

摘 要:选择吕梁山麓自Ar-Pt、C-O、C-P、T、Q五个岩组段10个小流域作为研究对象,借助ArcGIS软件,提取各小流域面积-纵比降关系以及河流宽度、分支比、弯曲度等地貌参数,展现了区内不同岩组段小流域地貌特征:1)在形成时间、形成环境、构造发育程度接近的情况下,吕梁山麓20 km2范围内小流域不同岩组段抗侵蚀能力由大及小依次为:太古-元古界河流>三叠系>石炭-二叠系>寒武-奥陶系>新近系-第四系黄土,抗侵蚀系数分别为2.294、2.232、1.486、1.394、1;2)小流域河流形成时间越长,水系分支比有逐渐变小的趋势;河流形成时间越短,支流间平均宽度越大;河流形成时间越久,纵比降变小,朝均衡化方向发展。

关键词:吕梁山麓;小流域;抗侵蚀能力;形成时间

地貌学是研究地球表面的形态、成因、分布特征及其发育规律的科学[1],河流是地球表面最常见的地貌形态,小流域是组成河流地貌的基本单元。在气候、构造、形成时间相同的条件下,岩性不同的小流域,遭受风化方式、侵蚀方式和侵蚀速率也不同,导致地貌形态各异,形成的灾害与采取的防治措施有别。

地貌是构造、营力和时间的函数[2],其形成的物质基础则是岩石,不同岩性分布区常常发育特定的地貌形态[3]。近半个世纪以来,国内外对流域地貌特征的研究主要针对活动造山带,如青藏高原周缘[4-5],利用基岩河道水力侵蚀模型,在保证岩性、气候相近的条件下,反映基岩抬升速率的空间差异[6-7]。在侵蚀地貌发育的造山带系统中,通过定量地貌参数表征流域地貌发育特征,反映构造活动强度在空间与时间上的差异[8-9],或分析侵蚀产沙[10-11],进而提出地质灾害防治措施[12-13]

近年来,随着地理信息技术的发展,利用GIS软件处理海量DEM数据,定量分析流域地貌特征[14-15],促进了地貌学研究的发展。在遥感地质应用中,地形因子是岩性解译的重要定量指标参数,可建立岩性与地形因子组合间的定量关联[16-17];在气候、演化时间相近的条件下,河道岩性发生变化时,河流可通过调整河道地形以改变侵蚀能力[18]。流域的地形条件是控制河流侵蚀程度的因素之一,在气候、演化时间接近的情况下,流域的地形条件与岩性密切相关。

吕梁山中段东西麓发育Ar-Pt、C-O、C-P、T、Q地层,是研究和对比不同岩组段小流域地貌特征及其差异的理想区域。现阶段,对本区的相关研究主要集中在地层学[19-20]、年代学[21-22]、夷平面地貌[23-24]、构造水文演化[25]等方面。本文选取吕梁山中段东西麓37°20′N~39°00′N间5个不同岩组段的10个小流域作为研究对象,利用ArcGIS软件,通过DEM对各小流域河流宽度、流域面积、纵比降、水系分支比、河流弯曲度等地貌参数进行了提取,归纳表征了不同岩组段小流域的地貌特征,以期从小流域地貌单元角度,探究不同岩组段抗风化、抗侵蚀能力的差异及其规律,进而为山区防灾减灾及针对性的各种工程措施提供依据。

1 技术路线与参数意义

目前广泛使用的基岩河道水力侵蚀模型经简化后成为一个统一的幂函数,但在模型简化过程中存在一系列假设,使得在岩性、坡度指数、侵蚀过程等方面存在不确定性[26]。本研究从数据统计出发,将面积-纵比降图分段线性拟合,以减少因函数简化导致的误差,具体技术路线如下:采用数字高程模型DEM,运用ArcGIS水文分析模块,结合地质图与野外踏勘,提取了太古元古界Ar-Pt(河流编号R1、R2)、寒武奥陶系C-O(河流编号R3、R4)、石炭二叠系C-P(河流编号R5、R6)、三叠系T(河流编号R7、R8)、第四系Q(河流编号R9、R10)10个20~100 km2面积范围内的侵蚀山区小流域空间分布图如图1所示。先在河道均匀选取数个分段点,自流域源头分段提取纵比降和该段的流域面积,绘制面积-纵比降关系的散点图;再分段线性拟合求取每段斜率k和截距b;最后根据面积-纵比降关系图的线性拟合分段结果,提取每段的河谷平均宽度、水系分支比、支流间平均宽度及河流弯曲度。其中每段纵比降和弯曲度是自河流源头到每个分段点一大段提取计算。

图1 小流域空间分布图
Fig.1 Location of the small watershed in the foot of the Lyuliang Mountains

1.1 流域面积与纵比降

流域面积是接受降雨补给的承载,流域面积越大,接受到的降雨通量也就越大,河流遭受侵蚀的能力也就越大。

纵比降是指在某一河段上,河床落差与其河流长度之比。河流的纵比降反映的是盆山分异强度,纵比降越大,盆山分异越强烈。计算公式如下:

IZ/L.

(1)

式中:I为纵比降,ΔZ为河床落差,L为对应的河流长度。

1.2 水系分支比

采用Strahler分类系统[27]对各小流域的水系进行分级,无支流河道定义为一级水系,两个一级水系交汇为二级水系,一个一级水系与一个二级水系汇合后仍为二级,两个二级水系汇合构成三级水系,依此类推。水系分支比是流域盆地内所有每一级别水系总数与下一级别水系总数比值的平均值。计算公式如下:

(2)

式中:Rb为河流的分支比;n为流域内水系的最高级别数;Ni为第i级水系所包括的河段的总条数。

水系分支比表征了流域内水系发育的成熟度,分支比越大,表明水系发育得越成熟;相反,分支比越小,表明水系仍未完全发育。水系分支比可以代表水系延展繁殖能力。

1.3 河流弯曲度

河流弯曲度是指河段的实际长度与该河段首末直线长度之比,称为该河段的河流弯曲度。可用下式表示:

Ka=L/l.

(3)

式中:Ka为河流弯曲度;L为河段实际长度;l为河段的直线长度。

河流的弯曲度可以代表小流域内构造发育强度。

1.4 支流间平均宽度

支流间平均宽度定义为在干流两侧所有相邻一级支流间距离的平均值。支流间平均宽度可以代表水系繁殖能力。

1.5 河谷平均宽度

河谷的平均宽度可反映河谷形成时间的长短,河谷越宽,表明形成的时间越长。当然,不同岩性段河谷宽度可代表抗侵蚀能力不同的岩性段河谷形成的时间。

2 不同岩组段小流域地貌特征

基于DEM提取了太古元古界Ar-Pt河流(河流编号R1、R2)、寒武-奥陶系C-O河流(河流编号R3、R4)、石炭-二叠系C-P河流(河流编号R5、R6)、三叠系T河流(河流编号R7、R8)、第四系Q河流(河流编号R9、R10),先在河道均匀选取数个分段点,自流域源头分段提取纵比降和该段的流域面积(表1、3、5、7、9),绘制面积-纵比降关系图如图2-6所示,再分段线性拟合求取每段斜率k和截距b.最后根据面积-纵比降关系图的线性拟合分段结果,提取每段的河谷平均宽度、水系分支比等地貌参数如表2、4、6、8、10所示。具体不同岩组段小流域地貌特征如下。

表1 河流R1、R2流域面积和纵比降值
Table 1 Watershed area and vertical ratio of R1 and R2

R1流域面积/km22.152 04.365 66.992 78.686 810.351 412.658 213.518 816.835 818.217 920.601 5纵比降0.163 60.132 30.111 90.096 10.093 80.087 40.082 80.079 30.077 40.072 4R2流域面积/km21.015 43.777 66.557 111.069 014.149 016.030 017.442 019.025 122.140 0-纵比降0.258 90.206 30.157 50.135 50.125 10.115 00.110 50.100 60.093 2-

表2 河流R1、R2各段地貌参数值
Table 2 Geomorphic parameter values of each segment of river R1 and R2

河流分段纵比降最高海拔/m最低海拔/m河谷平均宽度/m分支比支流间平均宽度/m河流曲度分段函数R1第一段0.096 12 2671 7132 2004.55101.28y=-0.010 0x+0.181 7(x=0-8.68)第二段0.082 82 2671 5942 2005.14151.08y=-0.002 6x+0.120 4(x=8.68-13.51)第三段0.072 42 2671 4882 5005.14101.01y=-0.001 4x+0.1031(x=13.51-20.60)R2第一段0.206 32 3661 6888004.58051.11y=-0.018 3x+0.276 8(x=0-3.78)第二段0.100 62 3661 0641 30010.04901.15y=-0.004 4x+0.186 2(x=3.78-19.02)第三段0.093 22 3669611 0009.14901.16y=-0.002 3x+0.146 1(x=19.02-22.14)

图2 河流R1、R2流域面积-纵比降关系分段图
Fig.2 Sectional graph of relationship between basin area and vertical gradient of R1 and R2

2.1 Ar-Pt小流域地貌特征

R1流域位于吕梁山中段关帝山东麓交城县庞泉沟,为文峪河末端支流,流域面积20.6 km2,R2流域位于云中山东麓原平市,为滹沱河一级支流,流域面积22.14 km2,主要由元古界的吕梁山花岗岩组成,间夹太古界片麻岩。二者的宽度不同,R1平均宽度为R2的2.23倍,根据华北山地地貌演化的阶段性[28],以及N2j红黏土在小流域内的分布状况,可判断R1为第三纪或第三纪之前更古老的河流,R2形成于汾渭地堑裂解之后的第三纪末。二者的局地侵蚀基准不同,R1为1 488 m,R2为961 m,进而导致二者流域面积-纵比降图不同。R1河流的分支比与支流间平均宽度上下游近似相等,表明此小流域接近平衡态,繁殖能力差;R2河流的分支比自上游至下游从4.5增长到9~10,支流间平均宽度上游大,下游小,表明此小流域生命力旺盛,繁殖能力强。二者河流曲度接近,反映二者构造发育强度接近。

表3 河流R3、R4流域面积和纵比降值
Table 3 Watershed area and vertical ratio of R3 and R4

R3流域面积/km21.791 45.846 411.177 415.477 123.177 326.851 034.075 641.010 944.369 849.752 854.178 0纵比降0.128 80.104 20.075 60.066 80.064 90.062 60.051 90.047 90.041 70.040 000.038 1R4流域面积/km23.824 411.499 719.018 831.517 365.218 086.064 896.217 0109.770 8---纵比降0.071 10.062 80.055 70.050 10.043 50.040 10.035 40.035 7---

图3 河流R3、R4流域面积-纵比降关系分段图
Fig.3 Sectional graph of relationship between basin area and vertical gradient of R3 and R4

表4 河流R3、R4各段地貌参数值
Table 4 Geomorphic parameter values of each segment of river R3 and R4

河流分段纵比降最高海拔/m最低海拔/m河谷平均宽度/m分支比支流间平均宽度/m河流曲度分段函数R3第一段0.066 82 0741 5813 0003.85101.04y=-0.004 6x+0.133 6(x=0-15.48)第二段0.047 92 0741 4603 7004.24751.02y=-0.000 8x+0.081 8(x=15.48-41.01)第三段0.038 12 0741 3631 6005.55201.05y=-0.000 6x+0.074 0(x=41.01-54.17)R4第一段0.055 71 9901 5101 7004.35251.30y=-0.001 0x+0.748 8(x=0-19.01)第二段0.043 51 9901 4051 2004.24951.22y=-0.000 2x+0.059 4(x=19.01-65.21)第三段0.035 71 9901 2171 3004.54901.27y=-0.000 2x+0.055 8(x=65.21-109.77)

表5 河流R5、R6流域面积和纵比降值
Table 5 Watershed area and vertical ratio of R5 and R6

R5流域面积/km23.401 87.418 611.092 716.406 219.615 926.248 038.789 842.781 252.719 755.798 459.171 462.127 3纵比降0.119 80.102 10.090 70.082 40.076 00.070 40.064 30.061 40.059 00.055 70.054 10.051 2R6流域面积/km21.300 12.132 13.126 54.373 88.418 111.369 613.828 216.866 517.691 818.563 4--纵比降0.193 30.143 60.123 20.098 20.095 40.081 80.074 70.074 30.071 20.066 0--

图4 河流R5、R6流域面积-纵比降关系分段图
Fig.4 Sectional graph of relationship between basin area and vertical gradient of R5 and R6

表6 河流R5、R6各段地貌参数值
Table 6 Geomorphic parameter values of each segment of river R5 and R6

河流分段纵比降最高海拔/m最低海拔/m河谷平均宽度/m分支比支流间平均宽度/m河流曲度分段函数R5第一段0.076 01 8341 1402 5004.55001.14y=-0.002 5x+0.124 1(x=0-19.61)第二段0.059 01 8349203 0004.45101.17y=-0.000 5x+0.084 8(x=19.61-52.71)第三段0.051 21 8347882 7005.24151.25y=-0.000 7x+0.099 8(x=52.71-62.13)R6第一段0.098 21 5951 0501 5954.15051.24y=-0.029 2x+0.219 5(x=0-4.37)第二段0.074 31 5958501 6005.24151.17y=-0.002 2x+0.109 4(x=4.37-16.87)第三段0.066 01 5957941 5905.15201.26y=-0.004 8x+0.156 6(x=16.87-18.56)

表7 河流R7、R8流域面积和纵比降值
Table 7 Watershed area and vertical ratio of R7 and R8

R7流域面积/km20.970 72.221 03.160 65.560 47.525 28.653 710.297 312.186 613.918 4-纵比降0.088 40.081 80.070 20.063 80.058 60.056 70.058 10.060 80.059 6-R8流域面积/km23.509 29.618 017.250 225.204 329.888 236.109 140.093 343.655 547.796 653.981 1纵比降0.050 70.038 80.031 50.029 90.027 60.026 80.023 50.022 80.023 20.023 6

图5 河流R7、R8流域面积-纵比降关系分段图
Fig.5 Sectional graph of relationship between basin area and vertical gradient of R7 and R8

表8 河流R7、R8各段地貌参数值
Table 8 Geomorphic parameter values of each segment of river R7 and R8

河流分段纵比降最高海拔/m最低海拔/m河谷平均宽度/m分支比支流间平均宽度/m河流曲度分段函数R7第一段0.070 21 4911 2521 0108.15151.08y=-0.008 1x+0.097 4(x=0-3.16)第二段0.056 71 4911 1301 3109.33601.06y=-0.002 5x+0.078(x=3.16-8.65)第三段0.060 81 4918601 0206.21 1001.13y=0.001 2x+0.046 5(x=8.65-12.18)第四段0.059 61 4917729102.07051.14y=-0.000 7x+0.069 4(x=12.18-13.91)R8第一段0.031 51 8521 4912 0507.14801.25y=-0.001 3x+0.054 3(x=0-17.25)第二段0.022 81 8521 3603 2004.35101.10y=-0.000 3x+0.038 0(x=17.25-43.65)第三段0.023 61 8521 2682 7005.84501.14y=-0.000 06x+0.019 8(x=43.65-53.98)

表9 河流R9、R10流域面积和纵比降值
Table 9 Watershed area and vertical ratio of R9 and R10

R9流域面积/km20.727 61.660 52.668 46.320 67.324 89.939 716.434 517.481 7纵比降0.100 80.058 20.051 70.041 70.038 60.034 30.030 10.022 4R10流域面积/km22.406 47.496 412.047 717.789 822.620 849.156 159.001 779.437 285.274 6纵比降0.046 30.026 20.017 40.014 50.012 50.010 80.009 60.009 10.008 7

图6 河流R9、R10流域面积-纵比降关系分段图
Fig.6 Sectional graph of relationship between basin area and vertical gradient of R9 and R10

表10 河流R9、R10各段地貌参数值
Table 10 Geomorphic parameter values of each segment of river R9 and R10

河流分段纵比降最高海拔/m最低海拔/m河谷平均宽度/m分支比支流间平均宽度/m河流曲度分段函数R9第一段0.058 21 2141 010 9004.16151.01y=-0.045 7x+0.134 1(x=0-1.66)第二段0.038 61 2149308054.36201.01y=-0.003 2x+0.062 1(x=1.66-7.32)第三段0.030 11 2148904204.76601.02y=-0.000 8x+0.044 3(x=7.32-16.43)第四段0.022 41 2148736104.06851.03y=-0.007 3x+0.151 4(x=16.43-17.48)R10第一段0.017 41 6811 4622 2103.55201.36y=-0.003 0x+0.052 1(x=0-12.04)第二段0.012 51 6811 4502 5303.54401.40y=-0.0000 4x+0.023 0(x=12.04-22.62)第三段0.008 71 6811 4043 1203.55751.82y=-0.000 06x+0.0137 5(x=22.62-85.27)

2.2 C-O小流域地貌特征

R3位于吕梁山东麓与晋中断陷交界部位,流域面积54.17 km2,R4位于汾河峡谷北部,流域面积109.77 km2,主要由寒武奥陶系灰岩组成。二者的宽度不同,R3上中段宽度达3 000~3 700 m,R4平均宽度与R3下段宽度近一致,表明R3是在更古老的河流基础上形成,R4则是汾河峡谷带动下形成,形成于汾渭地堑裂解之后的第三纪末。二者的局地侵蚀基准差异较小,R3为1 363 m,R4为1 217 m,进而导致二者流域面积-纵比降图类似。R3河流的分支比自上而下逐渐递增,支流间平均宽度上下游近似相等,表明此小流域上中游接近平衡态,繁殖能力差;R4河流的分支比与支流间平均宽度自上游至下游接近一致,表明本区域灰岩地区繁殖能力极限值在4附近。R4河流曲度是R3的1.2倍以上,表明R4构造发育强度大于R3.

2.3 C-P小流域地貌特征

R5、R6均位于太原西山煤田与晋中断陷盆地交界部位,主要由C-P煤系地层组成,侵蚀基准均为晋中断陷盆地。R5平均宽度达2 733 m,R6平均宽度为1 595 m,表明R5形成时间早于R6,均形成于汾渭地堑裂解之后的第三纪末。二者的局地侵蚀基准、分支比、支流间平均宽度近一致,表明本区域煤系繁殖能力极限值在5附近。二者河流曲度也近一致,表明构造发育强度近一致。

2.4 T小流域地貌特征

R7位于太原西山与晋中断陷交界部位,R8位于太原西山内部。二者的宽度不同,R7平均宽度达1 000 m左右,R8平均宽度在2 050 m以上,是R7的2倍以上,表明R8是在更古老的河流基础上形成,R7则是形成于汾渭地堑裂解之后的第三纪末。二者的局地侵蚀基准差异较大,R7为772 m,R8为1 268 m,进而导致二者流域面积-纵比降图不同。R7支流间平均宽度振荡较大,可能与边山构造发育有关,同样也反映在河流曲度的振荡上。R7的纵比降是R8的2.5倍,反映了R7处于河流侵蚀的旺盛期。

2.5 Q小流域地貌特征

R9位于系舟山南麓阳曲断陷盆地北侧,R10位于吕梁山西侧岢岚县境内。二者的宽度不同,R10平均宽度是R9的3.7倍,表明R10形成时间早于R9,现场调查显示,R10切割过最老地层为N2j,R9切割过最老地层为Q2l下部。二者的局地侵蚀基准差异较大,R9为890 m,R10为1 404 m,进而导致二者流域面积-纵比降图不同。二者分支比不同,R9为4左右,R10为3.5左右,反映二者繁殖能力不同,由于切割岩性差异,R9支流间平均宽度略大于R10.R9较R10而言,流域面积较小,干流河谷宽度较窄,纵比降较大,支流间平均宽度较大,却有较大的水系分支比。

3 讨论

流域地貌是气候、构造、时间共同作用于不同地质体的结果,只有在小流域的流域规模、气候环境、构造条件、形成时间长短类似的情况下,才能通过对比河谷平均宽度、水系分支比、支流间平均宽度、河流弯曲度、河流纵比降等形态学参数,探究吕梁山麓不同岩性小流域抗侵蚀能力的差异。

经拟合结果与流域实际状况综合对比表明,面积-纵比降关系图的线性分段点多位于两种位置:其一是在较大的支流汇入点附近,如R1、R3、R6、R7、R10的第一个线性分段点,R2的所有线性分段点。较大支流的汇入导致流域面积大幅增大,使纵比降原本随面积的下降的趋势减缓;其二是在新构造运动形成的河流裂点附近,如R7和R8的第二个线性分段点。侵蚀性河流中的裂点可以看做是下游新河道和老河道之间的移动分界点[29]。裂点上游河流纵比降随面积增大而递减,裂点下游纵比降却随面积增大有升高趋势。

因本次研究提取的小流域均位于吕梁山中段37°20′N~39°00′N之间,故可认为各小流域气候环境基本一致。为了控制各流域面积规模接近,选择20 km2附近的面积-纵比降图线性分段点的上游流域为对比研究区域。由于面积-纵比降关系图的线性分段点大多位于较大支流汇入点及河流裂点附近,选取分段点的上游流域为研究对象,既确保了研究对比区域的统计参数不受较大支流汇入的影响,又由于裂点上游为溯源侵蚀未及的老河道,从而减弱了构造运动对统计参数的影响。

河谷的平均宽度可反映河谷形成的时间,若以河谷宽度与时间成正比推断,在相同岩性、流域面积大小(20 km2范围左右)接近的情况下,各河流形成的时间关系如下:tR1=2.23tR2tR3=1.76tR4tR5=1.57tR6tR8=1.93tR7tR10=3.47tR9,野外调查表明,R2、R4、R6、R7、R10小流域中均残留有N2j红黏土,以此判断,其形成的时间均为新近纪末。

对比形成于新近纪末的R2、R4、R6、R7、R10小流域可知,河谷的平均宽度大小为:R10(2 370 m)>R4(1 700 m)>R6(1 595 m)>R7(1 062 m)>R2(1 033 m),若以R10的河谷平均宽度为参照分子,其余为分母,构建抗侵蚀系数E,则ER10=1、ER4=1.394、ER6=1.486、ER7=2.232、ER2=2.294.

水系分支比的大小代表了水系发育程度的高低[30]。河流的分支比越大,则河网对流域的填充程度越大,河网分布范围越大[28]。水系分支比可以代表水系的延展繁殖能力。Horton通过分析大量自然水系数据,得出自然水系平均分支比为3.5~5.0[31].同一时间段内,5种岩组20 km2流域的水系分支比提取结果显示,BR2=7.87、BR4=4.30、BR6=4.80、BR7=6.40、BR9=4.28,更古老河流不同岩性段分支比BR1=4.90、BR3=3.80、BR5=4.50、BR8=7.10、BR10=3.50,总体上年老河流比对应的相同岩性段年轻河流分支比小,反映了小流域河流形成时间越长,分支比有逐渐变小的趋势。

支流间平均宽度同样是构造与岩性的函数,20 km2流域范围支流间平均宽度WR2=595 m、WR4=525 m、WR6=480 m、WR7=670 m、WR9=645 m,更古老河流不同岩性段支流间平均宽度WR1=445 m、WR3=510 m、WR5=500 m、WR8=480 m、WR10=480 m,总体上年老河流比对应的相同岩性段的年轻河流的支流间平均宽度小。对比结果表明,河谷形成时间越短,支流间平均宽度越大。

河流曲度是河流平面地貌形态的重要表征指标[32]。河流往往就位于构造薄弱部位,河流曲度可以代表小流域内构造发育的强度,同一时间段内河流弯曲度为CR2=1.16、CR4=1.30、CR6=1.26、CR7=1.14、CR9=1.03;更古老河流不同岩性段弯曲度为CR1=1.01、CR3=1.04、CR5=1.14、CR8=1.25、CR10=1.40.

沟谷纵比降和流域面积的关系有明显的规律性,流域面积越大的沟谷,纵比降就相对较小,反之亦然[33]。沟谷的流域面积较小时,说明沟谷正处于河谷地貌发育的旺盛阶段,流水侵蚀作用强烈,沟床纵比降较大[34]。当然,流域面积和纵比降的关系也受流域发育所处的岩性段和构造环境的影响。纵比降反映下蚀强度,同样也反映河流所处的演化阶段。同一小流域不同分段纵比降不同,反映不同河段下蚀强度不同,同样也反映在分段函数的斜率上。同一时间段内年轻河流纵比降为VgR2=0.093 2、VgR4=0.055 7、VgR6=0.066 0、VgR7=0.059 6、VgR9=0.022 4,更古老河流纵比降为VgR1=0.072 4、VgR3=0.066 8、VgR5=0.076 0、VgR8=0.031 5、VgR10=0.012 5,总体上年老河流比对应的相同岩性段年轻河流要小。对比结果表明,小流域形成时间越久,纵比降变小,小流域朝均衡化方向发展。

4 结论

1)在形成时间、形成环境、构造发育程度接近的情况下,吕梁山麓20 km2范围内小流域不同岩组段抗侵蚀能力由大及小依次为:太古-元古界河流>三叠系>石炭-二叠系>寒武-奥陶系>新近系-第四系黄土,抗侵蚀系数分别为2.294、2.232、1.486、1.394、1.

2)小流域河流形成时间越长,分支比有逐渐变小的趋势;河谷形成时间越短,支流间平均宽度越大。小流域形成时间越久,纵比降变小,朝均衡化方向发展。

参考文献:

[1] 黄骁力,丁浒,那嘉明,等.地貌发育演化研究的空代时理论与方法[J].地理学报,2017,72(1):94-104.

HUANG X L,DING H,NA J M,et al.Theories and methods of space-for-time substitution in geomorphology[J].Acta Geographica Sinica,2017,72(1):94-104.

[2] DAVIS W M.The geographical cycle[J].Geographical Journal,1973,14(5):481-504.

[3] 焦树林,梁虹.喀斯特地区流域地貌与岩性的统计关系探讨——以贵州省为例[J].中国岩溶,2002,21(2):95-100.

JIAO S L,LIANG H.The research on the relationship between drainage basins’landforms and its lithologic features in karst region-a case study in Guizhou province[J].Carsologica Sinica,2002,21(2):95-100.

[4] HU X F,PAN B T,KIRBY E,et al.Spatial diferences in rock uplift rates inferred from channel steepness indices along the northern flank of the Qilian Mountain,northeast Tibetan Plateau[J].Chinese Science Bulletin,2010,55(27/28):3205-3214.

[5] WANG J,HU Z B,PAN B T,et al.Spatial distribution pattern of channel steepness index as evidence for differential rock uplift along the eastern Altun Shan on the northern Tibetan Plateau[J].Global and Planetary Change,2019,181:102979.1-102979.10.

[6] 赵洪壮,李有利,杨景春,等.天山北麓河流纵剖面与基岩侵蚀模型特征分析[J].地理学报,2009(5):563-570.

ZHAO H Z,LI Y L,YANG J C,et al.The longitudinal profiles of the ten rivers in North Tianshan Mountains and their tectonic significance[J].Acta Geographica Sinica,2009(5):563-570.

[7] 李小强,王军,熊仁伟,等.六盘山地区河道陡峭指数对隆升速率差异的响应[J].第四纪研究,2016,36(2):443-452.

LI X Q,WANG J,XIONG R W,et al.The response of the change channel steepness index to the difference of uplift rate in the Liupanshan mountain area[J].Quaternary Sciences,2016,36(2):443-452.

[8] 徐伟,刘志成,袁兆德,等.华山山前河流地貌参数及其构造意义[J].地震地质,2017,39(6):1316-1335.

XU W,LIU Z C,YUAN Z D,et al.River geomorphic parameters of the Huashan piedmont and their tectonic implications[J].Seismology and Geology,2017,39(6):1316-1335.

[9] LEE C S,TSAI L L.A quantitative analysis for geomorphic indices of longitudinal river profile:a case study of the Choushui River,Central Taiwan[J].Environmental Earth Sciences,2010,59(7):1549-1558.

[10] 廖义善,蔡强国,秦奋,等.基于DEM的黄土丘陵沟壑区流域地貌特征及侵蚀产沙尺度研究[J].水土保持学报,2008,22(1):1-6.

LIAO Y S,CAI Q G,QIN F,et al.Geomorphological feature and scale effect on sediment yield and soil erosion in Hilly Loess Areas,North China[J].Journal of Soil and Water Conservation,2008,22(1):1-6.

[11] 李晨瑞,李发源,马锦,等.黄河中游流域地貌形态特征研究[J].地理与地理信息科学,2017,33(4):107-112.

LI C R,LI F Y,MA J,et al.The study of watershed topography characteristics in the middle reaches of the Yellow River[J].Geography and Geo-Information Science,2017,33(4):107-112.

[12] 赵岩,孟兴民,郑娇玉,等.地貌学在泥石流研究中的应用与理论初探[J].灾害学,2017,32(1):43-49.

ZHAO Y,MENG X M,ZHENG J Y,et al.Application of geomorphological theory in study of debris flow and exploration of its applied theory[J].Journal of Catastrophology,2017,32(1):43-49.

[13] 倪化勇.海螺沟景区典型泥石流流域地貌特征及灾害防治[J].水土保持研究,2010,17(1):154-158.

NI H Y.Geomorphologic characteristics of typical debris-flow basins in Hailuogou scenic spot and disaster prevention[J].Research of Soil and Water Conservation,2010,17(1):154-158.

[14] 程三友,王红梅,李英杰.渭河水系流域地貌特征及其成因分析[J].地理与地理信息科学,2011,27(3):45-49.

CHENG S Y,WANG H M,LI Y J.Geomorphology characteristics of the Wei River Basin and its formation reasons[J].Geography and Geo-Information Science,2011,27(3):45-49.

[15] 梁欧博,任俊杰,吕延武.涪江流域河流地貌特征对虎牙断裂带活动性的响应[J].地震地质,2018,40(1):48-62.

LIANG O B,REN J J,LYU Y W.The response of fluvial geomorphologic characteristics of the Fujiang drainge basin to activity of the Huya fault zone[J].Seismology and Geology,2018,40(1):48-62.

[16] 陈霄燕,潘军,邢立新,等.桂林-阳朔地区DEM地形特征与岩性相关性分析及分类研究[J].地球信息科学学报,2019,21(12):1867-1876.

CHEN X Y,PAN J,XING L X,et al.Correlation analysis and classification of DEM topographic features and lithology in GuilinYangshuo,China[J].Journal of Geo-information Science,2019,21(12):1867-1876.

[17] 姜莎莎,李培军.基于ASTER图像和地形因子的岩性单元分类——以新疆木垒地区为例[J].地球信息科学学报,2011,13(6):825-832.

JIANG S S,LI P J.Lithologic unit mapping using ASTER data and topographic variables:a case study of Mulei area of Xinjiang[J].Journal of Geo-information Science,2011,13(6):825-832.

[18] ALLEN G H,BARNES J B,PAVELSKY T M,et al.Lithologic and tectonic controls on bedrock channel form at the northwest Himalayan front[J].Journal of Geophysical Research:Earth Surface,2013,118(3):1806-1825.

[19] 李建星,岳乐平,徐永,等.吕梁山地区新近纪晚期地层研究进展[J].地层学杂志,2009,33(2):177-182.

LI J X,YUE L P,XU Y,et al.New advances in the study of late neogene in the Lvliang mountain stratigraphy[J].Journal of Stratigraphy,2009,33(2):177-182.

[20] 岳乐平,邓涛,张云翔,等.保德阶层型剖面磁性地层学研究[J].地层学杂志,2004(1):48-51.

YUE L P,DENG T,ZHANG Y X,et al.Magnetostratigraphy of stratotype of the Baode stage section[J].Journal of Stratigraphy,2004(1):48-51.

[21] 朱大岗,孟宪刚,邵兆刚,等.山西保德-静乐地区新近纪地层时代讨论[J].地质通报,2008(4):510-516.

ZHU D G,MENG X G,SHAO Z G,et al.Ages of Neogene strata in the Baode-Jingle area,China[J].Geological Bulletin of China,2008(4):510-516.

[22] 李建星,刘池洋,岳乐平,等.吕梁山新生代隆升的裂变径迹证据及其隆升机制探讨[J].中国地质,2015,42(4):960-972.

LI J X,LIU C Y,YUE L P,et al.Apatite fission track evidence for the Cenozoic uplift of the Lvliang Mountains and a discussion on the uplift mechanism[J].Geology in China,2015,42(4):960-972.

[23] 吴忱.华北地貌环境及其形成演化[M].北京:科学出版社,2008.

[24] 熊建国,李有利,张培震.夷平面研究新进展[J].地球科学进展,2020,35(4):378-388.

XIONG J G,LI Y L,ZHANG P Z.New advances in planation surface research[J].Advances in Earth Science,2020,35(4):378-388.

[25] 赵金贵.西山煤田岩溶陷落柱形态学特征及构造水文演化[D].太原:太原理工大学,2004.

[26] 胡小飞,潘保田,李琼.基岩河道水力侵蚀模型原理及其最新研究进展[J].兰州大学学报(自然科学版),2014,50(6):824-831.

HU X F,PAN B T,LI Q.Principles of in research the stream power erosion model and its latest progress[J].Journal of Lanzhou University(Natural Sciences),2014,50(6):824-831.

[27] 王林,何仲太,马保起.岱海流域地貌演化及其对断裂活动性的指示意义[J].第四纪研究,2008,28(2):310-318.

WANG L,HE Z T,MA B Q.Geomorphic evolution and its implication for the fault activity in the Daihai drainage basin[J].Quaternary Sciences,2008,28(2):310-318.

[28] 王玉成,董有福,和娴.基于分形的小流域水系演化试验研究[J].人民黄河,2016(7):17-20.

WANG Y C,DONG Y F,HE X.Experimental research on the evolution of drainage networks in a small watershed based on rractal dimension[J].Yellow River,2016(7):17-20.

[29] 王乃瑞,韩志勇,李徐生,等.河流纵剖面陡峭指数对庐山构造抬升的指示[J].地理学报,2015,70(9):1516-1525.

WANG N R,HAN Z Y,LI X S,et al.Tectonic uplift of Mt.Lushan indicated by the steepness indices of the river longitudinal profiles[J].Acta Geographica Sinica,2015,70(9):1516-1525.

[30] 黄帅堂,赖爱京,徐衍刚,等.基于GIS平台的石羊河流域地貌演化及构造活动性研究[J].内陆地震,2017,31(4):398-404.

HUANG S T,LAI A J,XU Y G,et al.Geomorphic evolution and its tectonic activity study of Shiyanghe River Basin based on GIS[J].Inland Earthquake,2017,31(4):398-404.

[31] 任泓颖,何林芝,程济帆,等.田湾河流域地貌和水系对泥石流分布的影响[J].资源环境与工程,2019,33(2):212-216.

HE H Y,HE L Z,CHENG J F,et al.Effect of landform and water system on debris flow distribution in Tianwan River Basin[J].Resources Environment &Engineering,2019,33(2):212-216.

[32] 赵军,单福征,杨凯,等.平原河网地区河流曲度及城市化响应[J].水科学进展,2011,22(5):631-637.

ZHAO J,SHAN F Z,YANG K,et al.River sinuosity in tidal plain and its response to rapid urbanization[J].Advances in Water Science,2011,22(5):631-637.

[33] 黄成,张友谊,眭静,等.地形因子对沟谷泥石流发育的影响-以都汶高速“7.10”群发泥石流为例[J].人民长江,2019,50(8):115-119.

HUANG C,ZHANG Y Y,SUI J,et al.Influence of topographic factors on valley debris flow development:case of “7.10”cluster debris flows in Duwen Expressway[J].Yangtze River,2019,50(8):115-119.

[34] 朱渊,余斌,亓星,等.地形条件对泥石流发育的影响-以岷江流域上游为例[J].吉林大学学报(地球科学版),2014,44(1):268-277.

ZHU Y,YU B,QI X,et al.Topographical factors in the formation of gully type debris flows in the upper reaches of Minjiang River[J].Journal of Jilin University(Earth Science Edition),2014,44(1):268-277.

Geomorphologic Characteristics of Small Watersheds in Different Formations from the Foot of the Lyuliang Mountains

YUE Keshan,ZHAO Jingui,WANG Jianmin

(CollogeofMiningEngineering,TaiyuanUniversityofTechnology,Taiyuan030024,China)

Abstract:In order to illustrate geomorphic characteristics in different sorts of formations,10 small watersheds from Ar-Pt,C-O,C-P,T,and Q formations in the foot of Lyuliang Mountains were selected as the research objectives and the relationship between the area of these watersheds and its vertical gradient,as well as the geomorphic parameters including the width of the rivers,branch ratio,and curvature,were obtained by using ArcGIS software.We got such conclusions:(1)In 20 km2research area,the erosion resistance of these formations is Archean Proterozoic >Triassic >Carboniferous Permian>Cambrian Ordovician>Neogene Quaternary loess and its resistance coefficient is 2.294,2.232,1.486,1.394,and 1,respectively,if the formation time,paleoenvironment,and structural development do not show a clear difference.(2)The longer the river formation time is,the smaller the branch ratio of the river system is;The shorter the river formation time is,the larger the average width between tributaries is;The longer the river formation time is,the smaller the vertical gradient is,and the development trend is towards equalization.

Keywords:the foot of the Lyuliang Mountains;small watershed;erosion resistance;time of formation

中图分类号:P931

文献标识码:A

DOI:10.16355/j.cnki.issn1007-9432tyut.2022.02.014

文章编号:1007-9432(2022)02-0299-09

引文格式:岳科杉,赵金贵,王建民.吕梁山麓不同岩组段小流域地貌特征[J].太原理工大学学报,2022,53(2):299-307.

YUE Keshan,ZHAO Jingui,WANG Jianmin.Geomorphologic Characteristics of Small Watersheds in Different Formations from the Foot of the Lyuliang Mountains[J].Journal of Taiyuan University of Technology,2022,53(2):299-307.

收稿日期:2020-07-27

第一作者:岳科杉(1995-),硕士,(E-mail)jerryyuea@163.com

通信作者:赵金贵(1974-),博士,副教授,主要从事构造地貌学研究,(E-mail)22971865@qq.com

(编辑:贾丽红)

Baidu
map