ACTASEDIMENTOLOGICASINICA
沉积学报
Vol.36 No3Jun.2018
文章编号:1000⁃0550(2018)03⁃0427⁃19DOI:10.14027/j.issn.1000⁃0550.2018.045
基于地貌形态学交融的现代曲流河道迁移构型表征
2.长江大学地球科学学院,武汉 430100
1.长江大学油气资源与勘探技术教育部重点实验室,武汉 430100
林志鹏1,2,单敬福1,2,陈乐1,2,孙千军3,王义武4,汪菲5
5.大庆油田采气分公司,黑龙江大庆 163000
4.北京大学地球与空间科学学院,北京 100871
3.中国地质大学能源学院,北京 100083
摘 要 曲流河道地貌迁移规律是重建与恢复古河道储层构型与历史演变的基础,是预测未来迁移演化趋势的前提,目前正成为国际河流沉积演化研究的热点。通过地貌形态学与河流沉积学的学科交融,以额尔齐斯河和诺威特纳河两条典型天然曲流河道为研究对象,结合GoogleEarth和ACMEMapper技术获取高分辨率历史卫星图像,提出现代曲流河道地貌迁移构型表征新方法与思路。在前人基础上,综合采用28种刻画要素对曲流河道平面构型进行层次解剖,更加精准定义弯度指数(S)与曲率(C),学分析。通过对12条优选典型河曲的构型刻画,表征6种迁移结构:对称扩张、上游旋转扩张、下游旋转扩张、对称收敛、上游旋转收敛和下游旋转收敛结构,归纳出9种迁移模式,并结合地貌过程学定量表征思路,探讨其对现今地下储层构型方法的重要指导意义,以期推进河流地貌形态学与河流沉积学的共同发展。关键词 曲流河道;地貌形态学;迁移模式;构型表征;表征参数
第一作者简介 林志鹏,男,1993年出生,硕士,河流沉积学、河流地貌形态学、储层构型与不整合结构,E⁃mail:linzp93@yangtzeu.通信作者 单敬福,男,副教授,E⁃mail:shanjingfu2003@163.com中图分类号 P512.2 文献标志码 Aedu.cn
并首次运用3种新构型表征参数:上游偏转角(Δθ)、下游偏转角(Δθ′)和扩张系数(KM),实现曲流河道迁移构型表征的定量过程
0 引言
Leopoldetal.[16]对河道水动力学特征、曲流结构、模式[17]以及河流地貌学特征[18]等做过详细研究,并取得一定的成果;Wolmanetal.[19]通过对地貌过程分析认为河流周期性的流量变化是控制河流形态的重要因素,并结合实例追踪研究河流演化过程;Wil⁃liams[20]通过对水流通量的研究提出满岸流量的定义和算法;Brice[21]提出曲流环的4种迁移模式,包括对称、不对称、复合对称和复合不对称模式;Hooker[22⁃24]对曲流河道地貌形态学做了大量研究,总结并提出了河弯变形的7种基本类型和3种复合类型;Miall[25⁃28]对河流体系的研究得出最好的线索还是结合野外露头的结论,从构型要素到沉积演化的不断深入重建历史过程。Willisetal.[13]通过三维模拟河道迁移模式,进一步提出曲流环的平移、旋转、扩张和复杂模式。来则开始尝试通过结合地貌形态学的方法来恢复古Ielpietal.[15],Ghinassietal.[29]和Wuetal.[30]近年
推演重建曲流河道构型及其沉积演化过程需重
点把握地貌形态迁移结构和演化规律,现阶段这项基础性理论研究相对不足。对现代曲流河地貌迁移规律的认知是深入解剖地下河流储层构型的钥匙,也是且河流相砂体内部结构十分复杂,传统方法难以精确恢复古河道历史演化过程的一种重要科学方法[1⁃5],式[6⁃10]。通过曲流河道地貌迁移规律的构型表征,有
揭示储层内部结构单元分布和预测油气的运移模
助于指导古河道砂体地下空间构型的分布特关重要的研究手段之一。尽管对古代和现代河流沉分析依然是极具挑战的任务[13,15]。国内外许多经典积体系已存大量研究,但曲流河动态地貌的过程演化征[11⁃14]。结合地貌形态学的理论和方法则是其中至
文章从河道的规模估算到水动力学及数值模拟等方面尝试对河流三维动态迁移结构进行深层次解析,如
河道特征,实现动态地貌过程演化分析。相比国外,
收稿日期:2017⁃05⁃19;收修改稿日期:2017⁃08⁃22
基金项目:国家自然科学基金项目(41372125);国家科技重大专项(2016ZX05063002⁃006)[Foundation:NationalNaturalSciencesFoundationof
China,No.41372125;NationalScienceandTechnologyMajorProject,No.2016ZX05063002⁃006]
428
沉 积 学 报 第36卷
国内则更多是关注地下储层结构,对曲流河形态的研究主要集中在侧积体的概念性划分。裘怿楠[31]对河流形态对河型的控制做了详细说明;薛培华[14]将曲流点坝沉积建造单元分为侧积面、侧积层、侧积体,结合拒马河现代沉积剖面对侧积体发育与演化过程进行解剖;尹燕义等[32]通过对现代饮马河点坝的三维解剖,对其内部侧积体沉积特征进行了详细研究;岳大力等[33]提出地下点坝识别及内部结构精细解剖的方法,并将其应用于识别古代点坝沉积;周银邦等[34]而得到多种演化模型以指导和预测同等地质条件下的河道形态演化。Schuurmanetal.[5]通过建立地形数值模型实验来观测曲流河点坝的形态变化和迁移状况。由此进一步证明,揭示河道过程演化离不开河流地貌形态学与沉积学科的交融。
(基于GoogleEarth的高精度通用绘图应用程序)获
鉴于此,本文通过GoogleEarth和ACMEMapper
取高分辨率历史卫星图像,结合两条天然曲流河道的迁移结构,采用新的精细构型表征参数法,实现河道通过现代露头与构型研究分析了点坝侧积体结构形成的控制因素“再现古代沉积演化过程的意义成因”的研究;,单敬福等[35]探讨了曲流河道现代沉积的规律对
结合对点坝内部侧积体;林志鹏等[3,36]通过结合曲流河地貌动力学理论,对曲流河道迁移构型结构提出了新的量化表征思路。综合看来,对于曲流河迁移构型的研究,国外起步虽早,但总体较为分散,国内则稍显欠缺,大部分研究尚未有效结合地貌形态学理
论,在没有弄清河道迁移演化规律的基础上进行的储层表征思路,需要得到纠正。
对现代曲流河道迁移结构的精细表征是解剖河流沉积过程和恢复古河道的重要方法和参考依据。曲流河道地貌迁移演化规律的不同,会影响储层单元的岩性和储集物性,河道平移、旋转、扩张等变化都会引起河道形态的改变,进而影响储层构型[23,37⁃38]年来诸多学者着手河流地貌形态学理论,取得了显著。近的成果,如Ielpietal.[15]通过对英国约克郡斯卡伯勒北部曲流河构型和动态地貌学的分析,结合现代露头的平面和剖面特征,提出了适用于该流域的一套综合河道迁移结构模型;Wuetal.[30,39]通过对复合点坝露头地貌迁移形态的重建分析了河流沉积演化过程,并建立了一种半定量的河道迁移模式;Ghinassietal.[29,40]通过研究曲流河道的迁移模式,探讨了重建古河流系统和开发河流相油气藏的重要方法;吴胜和等[41]提出“层次约束、模式拟合和多维互动”的构型分析与建模思路,论述地下古河道储层结构的演化方法;单敬福等[42]通过探讨“初末流线包络线法”还原古河流迁移形态,重建河道演化过程。近年来也有大量成果通过数值模拟技术加以研究,如Asahietal.[43]通过数值模拟计算曲流河繁复的演化过程,并探讨河道演化过程的可预测性。Schuurman[44]结合河流动力学数值模拟技术解释河流地貌演化过程Rousseauetal.;
[45]通过给定曲流河初始的水流参数和地形条件,同时改变相关河流参数进行数值模拟,从
地貌迁移构型表征,得出其迁移规律与模式。拟主要解决以下四个问题:1)提出一套完善的曲流河地貌迁移形态构型表征方案;2)结合不同曲流河道迁移特征进行历史演化过程形态学分析,并对其结构加以详细刻画;3)针对不同迁移结构,应用表征参数归纳一套完整的曲流河道迁移模型;4)探讨迁移结构对地下储层构型的重要指导意义。
1 地质背景
天然河道的选择需要着重考虑人类活动因素,且
不同地区不同曲流河具有的特征差异性也需要加以研究,为更好展示不同地域河流的特征,本文选取2条来自不同国家、人口分布分散、自然环境和河道结构保存较为完整、易于观察的典型曲流河作为研究对象,每条河流各选取50条河曲进行统计分析,并选取其中12段河曲进行精细表征。这2条河流分别为额尔齐斯河和诺威特纳河。
额尔齐斯河(IrtyshRiver)是鄂毕河(ObRiver)最大的支流之一,全长约4248km,源起中国新疆维吾尔自治区阿尔泰山南坡,自东南向西北流出中国,经哈萨克斯坦向北流入俄罗斯,在汉特—曼西斯克汇入鄂毕河60°56′N,69°19′E,。研究河段南起托博尔斯克北部km,如图1中IR:I-II,,坐标为较好直线距离为河段长44.51,易于观察。
22.99km。该河段曲流河结构保存诺威特纳河(NowitnaRiver)位于美国阿拉斯加州国家野生动物保护区,为育空河(YukonRiver)一支流,属于典型的曲流河,自北部卡斯科奎姆山流经诺威特纳国家野生动物保护区一直流向塔纳诺西南576.27地区汇入育空河。全长约455.45km,流域面积达到
或有连绵起伏的山川丘陵km2,拥有众多的湖泊和沼泽,周围地形平坦,km因素的干扰被指定为野生河流,河流结构保存十分完好,也因此现阶段很少受到人为,河岸植被茂盛;其中362。研究河段北起
第3期 林志鹏等:基于地貌形态学交融的现代曲流河道迁移构型表征
429
图1 额尔齐斯河与诺威特纳河研究区位置图
A.额尔齐斯河,IR:I-II为研究河段;B.诺威特纳河,NR:I-II为研究河段;坐标信息源于GoogleEarth和ACMEMapper
Fig.1 LocationmapoftheIrtyshRiverandNowitnaRiver
205°45′E,南抵坐标为64°36′N,205°33′E的位置,如图1中NR:I-II,河段长97.89km,直线距离为
诺威特纳河与育空河分流河口处,坐标为64°53N,曲率圆、点坝(上游坝、中游坝、下游坝)和深泓线。
曲流环(Meanderloop),又称河弯,是河道蜿蜒过程中的某段环状弯曲,如图2中(4)和(5)所限定的范围就是一个单一曲流环,是曲流河道的基本构成单元。一条曲流河本质上总是可以分为数段连续的曲流环。中心线(Centerline)是河道所有横截面外岸与内岸边缘中点的连线,是抽象概念线条,其长度表示为LC(图2);线上每点都能代表河道某个横截面的
35.53km。该河段河道迁移现象显著,特征明显,有利
于对原生态曲流河的迁移结构特征等进行深入研究。
2 迁移构型表征
尽管现代冲积平原上大量曲流河的形态特征已
被识别,但在沉积学领域对河道的刻画依然停留在沉积微相的解释中[46⁃48],结合地貌形态学加以精细表征的体现较少,故此结合地貌形态学与河流沉积学理论,系统完善曲流河道平面结构的构型表征。综合前为静态要素、动态要素和表征参数三大类。
人的分析[11,13,21,49],将曲流河道迁移构型要素具体分2.1 静态要素
平面几何中心,在对河道进行线性研究和形态学研究时,可近似将中心线作为河道本身,此概念的使用使得对河道流线的研究和迁移轨迹有了更加便捷和明确的表示,不用考虑外岸与内岸距离的差异。曲流环转折端(Bendinflection)是两相邻曲流环的过渡位
静态要素指对曲流河形态结构进行定性描述的地貌沉积单元和抽象概念流线,主要包括河道内要素:河道、深洼、浅滩、深泓线、中心线、曲流环转折端和曲率顶点;河道外要素:曲流带、点坝(上游坝、中游坝、下游坝)、串沟河道、外岸(凹岸)、内岸(凸岸)、曲流环(河弯)、废弃河道、河漫滩、曲率圆和曲流环顶点等。其中对迁移结构起到重要表征作用的主要是曲流环(河弯)、中心线、曲流环转折端、曲率顶点、
2);转折端的出现表明河道的迁移方向发生改变,是识别不同曲流环的重要标志。曲率顶点(Curvatureapex)是曲流环外岸顶点与内岸顶点连线的中点,也是曲流环轴向线与中心线的交点(图2);其迁移方向也可指示曲流环的迁移方向。曲率圆(Curvaturecir⁃cle)是表示单一曲流环弯曲程度的几何圆(图2),为相邻两个曲流环转折端连线的中点,以连接该中点和曲率顶点之间的线段作为直径所的圆就是曲率圆,其圆心一般落在点坝中游坝附近,圆周越大,表示曲流
置,形态学上即是曲流带轴向线与中心线的交点(图
430
沉 积 学 报 第36卷
环弯曲程度越大;圆周越小,表示曲流环弯曲程度越小。点坝(Pointbar)又称心滩,是河道迁移摆动和沉积物堆积的重要沉积单元,发育于河道内岸,是曲流环的主体(图2),据其形态位置又可定性分为上游坝(Upstreambar)、中游坝(Centralbar)和下游坝
种类各有差异。深泓线(Thalweg),也称最深谷底线,是河道内海拔最低点的连线(图2);是抽象概念线条,所代表的流域流速相对较高,谷底遭受冲蚀程度最大,总是偏向外岸一侧。
(Downstreambar),其水流强度和沉积物颗粒大小与
图2 曲流河道迁移结构精细表征
像,摄于2016年8月23日;C理想化曲流河模型刻画其迁移结构,共包括17个静态要素和11个动态要素
A美国诺威特纳河NR:I-II河段河曲15、16、17现代露头(露头照片源于OliverKurmis),坐标见表1;B龙溪河GoogleEarth卫星图
Fig.2 Planformelementsofmigrationarchitectureofthemeanderingriver
第3期 林志鹏等:基于地貌形态学交融的现代曲流河道迁移构型表征
431
2.2 动态要素
动态要素指从静态要素中提取出来的动态参数,是河道结构的量化反映,其数值变化一定程度上可指示河道演化过程。包括标量要素:曲流带宽、曲流带长、单一曲流带宽、单一曲流带长、单一曲流环宽、单一曲流环长、中心线长和曲率半径;矢量要素:曲流带轴向、曲流环轴向、上游偏转角和下游偏转角(表1、图2)。
曲流带宽(WM)指完整多期复合河道的迁移摆动宽距离;其值随河流规模不同而不同,范围从几百米长(LM)指完整多期复合河道的跨域长度,即一条河示所研究的某单一曲流带长,所以没有标注LM;用轴到上万米,本研究区WM为8700~22000m。曲流带
度,即单一河道轴向线的长度,如图2中(3)长度所示。其值变化灵活,根据所研究对象规模不同而改变,本研究区LSM为3300~23000m。单一曲流环宽2中(4)宽度所示。其值一般为数百米到上千米,本相邻两曲流环转折端之间的直线距离,如图2中(5)研究区WML为700~3800m。单一曲流环长(LML)指(WML)指曲流环顶点到曲流带轴向的垂直距离,如图
长度所示。其值一般为数十米到上千米,本研究区LML为400~7600m。中心线长(LC)即是河道中心线的长度,其值可反映研究河段的实际长度,根据河流11000m。曲率半径(R)即曲率圆的半径,其值为曲率顶点到相邻曲流环转折端中点距离的一半,如图3中(6)长度所示。其值一般为数十米到上百米,本研带蜿蜒摆动的中心轴,是各曲流环转折端的连线,其方向指向河流流动方向,当研究某单一曲流环时,其长度近似等于LML长度。本研究区AXMB长度为400~7600m,方向为50°~340°。曲流环轴向(AXML)指究区R为300~2000m。曲流带轴向(AXMB)指曲流规模其值变化范围较大,本研究区LC为1600~
宽度,如图2中(1)宽度所示,包含整个复合河道最
流从开始到结束曲流带轴向线的长度,图3中只能表向线代替直线长度融入了河道地貌形态的考虑,其变化范围区间很大,可从数千米到上万米,本研究区两带宽(WSM)指单一河道的迁移摆动宽度,如图2中一般为数百米到上万米,本研究区WSM为1800~条河流LM分别为302741m和55458m。单一曲流
(2)宽度所示,一般大于两倍单一曲流环宽度,其值12000m。单一曲流带长(LSM)指单一河道的跨域长
曲流环流动的中心轴,是曲流环顶点与对应两转折端中点的连线,其方向指向外岸扩张方向,本研究区AXML方向为44°~310°。上游偏转角(θU)指曲率顶点
表1 研究区12段河曲结构表征要素统计表
Table1 Theplanformarchitectureelementsandparametersofthe12keystudyingmeandersinthesturdyreaches
河流河曲编号纬度经度WM/mWSM/mWML/mLML/mLC/mR/m|AXMB|/mAXMB/°AXML/°θU/°SCΔθ/°Δθ′/°KMθD/°LSM/mLM/m
158°52′N16520.912263.03821.3919.730274111094.810166.92155.3919.7344.344.012968°47′E
268°50′E3027419765.02541.6
额尔齐斯河1168°52′E3027412467.31758.1729.78409.43551.94626.33551.9211.85055311.91.300.0014-52.435
1869°17′E3027413225.52148.5226.994.41855
3169°48′E3027417394.32642.61318.149.73028
3869°52′E3027413309.8
1864°45′N9299.9
2964°41′N10283.34111.81022.0465.86970.91078.12487.31078.1292.32146174.22.310.0021-251.1625
58°56′N59°20′N59°34′N60°14′N60°40′N14479.611813.22472.97305.9103.72013
13068.9
14855.112230.87637.923525.611994.37637.9
10522.017068.82282.76999.22282.7316.83.070.00082-2
22592.713075.53145.018441.56315.5336.53415
诺威特纳河3064°42′N10327.03941.52269.6866.57339.3
3664°40′N12014.76217.5700.2684.81873.8
4164°39′N10044.37918.2715.93201.4
5064°35′N8770.7
205°40′E205°34′E205°32′E205°27′E205°30′E205°39′E55458.01774.7654.8727.93660.91608.9300.1237.23234727.9
55458.0
55458.0
55458.0
55458.0
55458.03336.5592.0425.62217.4
1186.1
1389.6
2472.8136.92.95
1311.9
4578.91186.1278.62311311.23.86
1884.5355.3279.816.8262.750.0028-40.96430684.8
1221.3456.6
2300.01221.3109.858.54322
1547.8296.9130.929.31822425.6
3145.0110.12.01
143.72.21
11.050.0005-30.223
1.570.0005-371.7837
1.880.002221-21
3.640.0034-40.724
0.00077-7
0.000819-19
0.0033-21.212
0.001212-12
0.890.871.200.681.34
432
沉 积 学 报 第36卷
图3 对称扩张结构迁移构型表征
摄于1984年12月,②摄于1994年12月,③摄于2004年12月,④摄于2014年12月;IR2为对IR1中河曲11迁移结构要素拍摄的卫星图像,坐标为64°45′N,205°40′E,选取2张,其中①摄于2011年3月,②摄于2012年6月;NR2为对NR1中河曲18迁移结构要素的平面刻画;NR3为河曲18的迁移结构表征参数统计数据
IR1为额尔齐斯河研究区内河曲11近30年来通过GoogleEarth拍摄的卫星图像,坐标为59°20′N,68°52′E,选取4张,其中①
的平面刻画;IR3为河曲11的迁移结构表征参数统计数据。NR1为诺威特纳河研究区内河曲18近2年来通过GoogleEarth
Fig.3 Planformmigrationarchitectureofthesymmetricalexpansionstructure
和对应上游曲流转折端连线与曲流环轴向的夹角,如
图2中(9)角度所示;下游偏转角(θD)指曲率顶点和对应下游曲流转折端连线与曲流环轴向的夹角,如图
构型表征参数:弯度指数、曲率、顺流偏转角差、逆流偏转角差和扩张系数,其中顺流偏转角差、逆流偏转角差和扩张系数为本文首次提出的表征参数;同时弯度指数与曲率的定义也更为精确。
弯度指数(S)是指研究河道范围内,中心线长度与曲流带轴向线长的比值。前人对弯度指数的定义是河道长度与河谷长度的比值,但这一定义存在的问题是,河道与河谷长度区分不清晰,是外岸线、内岸线、深泓线、中心线长度,还是河道起点和终点的直线段长度,容易引起混淆。采用中心线长度代表河道长度可统一内、外岸线的问题,曲流带轴向线长度代表河谷长度,其一不会引起概念模糊;其二,轴向线长度比
2中(10)角度所示,根据曲流环形态的不同,夹角范差值是判断迁移结构的重要依据。本研究区θU为9°2.3 构型表征参数~50°,θD为12°~55°。
围介于0°~90°之间,且上游偏转角和下游偏转角的
静态要素和动态要素可刻画河道形态结构,但地貌迁移特征的量化表征与规律性模式,还需要对结构参数做特征分析,得到河道结构动态演化的变量,即表征参数。依据结构要素,提炼出5个曲流河道迁移
第3期 林志鹏等:基于地貌形态学交融的现代曲流河道迁移构型表征
433
河道起点和终点的直线段长度能更真实地反映地貌迁移情况,因为它考虑了河道随地形蜿蜒变化的因素,而非单一简化地当做一条直线处理。用公式表示为: S=LC/AXMB(1)
化表征,本文结合形态结构表征参数,对研究区曲流河道迁移演化规律进行定量研究。因为河流对堤岸的作用总是侵蚀凹岸,堆积凸岸,当点坝侧向堆积至一定程度时易发生截弯取直,随之重复新一轮侧积侵蚀过程,所以曲流河道整体迁移构型宏观上总呈现出重复性不断向外扩展的趋势。具体根据迁移前后河道地貌形态演化特征,并结合表征参数的定量刻画,可将曲流河道迁移构型具体分为扩张结构和收敛结构两大类。
曲率(C)指曲率半径的倒数,表示河弯的弯曲程
度,曲率越小,表示河弯程度越大;反之则越小。已有
通过点坝弧长直接除点坝直径计算曲率的方法[50],因为点坝形态不规则,提取的直径和弧长很不明确,误差较大。选择曲率半径,通过转折端来确定河曲开始迁移的起点,再结合曲率顶点,控制整个曲流环形态,最终以曲率顶点和对应两转折端3点确定的曲率圆则能更有效反映曲率的真实情况 。顺流偏转角差C=1/R用公式表示为:(Δθ)指上游偏转角与下游偏转角(2)
的差值,用以反映河弯对称性。差值越接近于0,表明曲流环对称性越高,差值为正且越大,表明曲率顶点更靠近上游转折端,呈逆流旋转趋势;差值为负且越小,表明曲率顶点更靠近下游转折端,呈顺流旋转趋势 。Δθ用公式表示为=θ:U-θ逆流偏转角差D
(Δθ′),与顺流偏转角差对应相
(3)
反,指下游偏转角与上游偏转角的差值,差值越接近于0,表明曲流环对称性越高,差值为正且越大,表明曲率顶点更靠近下游转折端,呈顺流旋转趋势;差值为负且越小,表明曲率顶点更靠近上游转折端,呈逆流旋转趋势 。用公式表示为:扩张系Δθ′=θ数D-θ(KC
M)指单一曲流环长与曲率直(4)
(2径
所以基本能控制曲流环形态R)的比值,因为它们都由转折端和曲率顶点控制,曲流环若不断往外扩,张,受转折端控制,曲率直径会逐渐变大,而曲流环长则改变幅度较小;当曲流环扩大到一定程度开始收敛,走向截弯趋势时,受转折端控制,曲率直径一般会大于曲流环长。因此,他们可以对河道的扩张收敛情况提供基本的量化反映。比值大于1,则说明处于扩张期,比值越大,扩张程度越小;比值小于1,则说明处于收敛期 K,比值越小,收敛程度越大。用公式表示为:M=LML/2R(5)研究河段的基本结构表征要素和参数数据可见
表1。
3 地貌形态迁移构型特征
对点坝侧积形态和河道演化的定性描述缺乏量
3.1 扩张结构
曲流河扩张结构(ExpansionStructure,ES)指河道向外不断侧积侵蚀迁移,K环和点坝向上游和下游摆动的幅度与M通常大于1。根据曲流Δθ值的差异,可细分为对称扩张结构、上游旋转(逆流旋转)扩张3.1.1 结构和下游旋转对称扩张结构
(顺流旋转)扩张结构。
1,ΔExpansionθ曲流环不断向外岸侧积侵蚀靠近Structure,0°,此为曲流河对称扩张结构,对称性好SES)。如图3中(IR1Symmetrical
,KM大于1984表明从
移不明显年到,2014通过年IR2,额尔齐斯河研究区内河曲迁移结构表征和IR3表征参数的11迁计算分析,可见该曲流环S为1.3,C为0.0014,反映弯曲程度较低;Δθ为-5°,表明曲流环对称性较好,呈弱顺流旋转趋势;K且扩张程度较小。综合表明河曲M为2.43,表明河弯处于扩张期,
11为SES。NR1中18诺威特纳河几乎没有发生迁移2011年到,2012通过年的卫星图像显示河曲NR3C为表征参数的计算分析0.0033,反映弯曲程度较大,可见该曲流环NR2迁移结构表征和;Δθ为-2°,S为表明曲2.21,流环对称性较好,呈弱顺流旋转趋势;K河弯处于扩张期,且扩张程度较大。综合表明河曲M为1.21,表明183.1.2 为SES。
曲流环不断向外岸侧积侵蚀上游旋转(逆流旋转)扩张结构
,且曲率顶点更靠近
上游转折端,呈逆流旋转趋势,KM大于1,Δθ大于0°,此为曲流河上游旋转Rotation表明从1984Expansion年到2014Structure,URES)。(逆流旋转)扩张结构年,额尔齐斯河研究区内河曲
如图(Upstream4中IR4
38数的计算分析迁移不明显,,可见该曲流环通过IR5迁移结构表征和S为2.01,CIR6为0.000表征参8,
反映弯曲程度较大;Δθ为19°,表明曲率顶点更靠近上游转折端,呈逆流旋转趋势;K张程度较大。综M为1.20,表明河弯处于扩张期,且扩合表明河曲38为
434
沉 积 学 报 第36卷
图4 上游旋转(逆流旋转)扩张结构迁移构型表征
IR4为额尔齐斯河研究区内河曲38近30年来通过GoogleEarth拍摄的卫星图像,坐标为60°40′N,69°52′E,选取4张,其中①摄于1984年12月,②摄于1994年12月,③摄于2004年12月,④摄于2014年12月;IR5为对IR4中河曲38迁移结构要素的平面刻画;IR6为河曲38的迁移结构表征参数统计数据。NR4为诺威特纳河研究区内河曲41近2年来通过GoogleEarth拍摄的卫星图像,坐标为64°39′N,205°30′E,选取2张,其中①摄于2011年3月,②摄于2013年3月;NR5为对NR4中河曲41迁移结构要素的平面刻画;NR6为河曲41的迁移结构表征参数统计数据
Fig.4 Planformmigrationarchitectureoftheupstreamrotationexpansionstructure
URES。NR4中诺威特纳河2011年到2013年的卫星图像显示河曲41迁移不明显,通过NR5迁移结构表1.88,C为0.0028,反映弯曲程度较大;Δθ为21°,表为1.34,表明河弯处于扩张期,且扩张程度较大。综3.1.3 下游旋转(顺流旋转)扩张结构合表明河曲41为URES。
明曲率顶点更靠近上游转折端,呈逆流旋转趋势;KM征和NR6表征参数的计算分析,可见该曲流环S为
内河曲18迁移缓慢,通过IR8迁移结构表征和IR9表征参数的计算分析,可见该曲流环S为1.57,C为0.0005,反映弯曲程度较大;Δθ为-37°,表明曲率顶表明河弯处于扩张期,且扩张程度较大。综合表明河曲18为DRES。NR7中诺威特纳河2011年到
点更靠近下游转折端,呈顺流旋转趋势;KM为1.78,
2012年的卫星图像显示河曲29河道迁移缓慢,通过NR8迁移结构表征和NR9表征参数的计算分
曲流环不断向外岸侧积侵蚀,且曲率顶点更靠近
下游转折端,呈顺流旋转趋势,KM大于1,Δθ小于0°,streamRotationExpansionStructure,DRES)。如图5此为曲流河下游旋转(顺流旋转)扩张结构(Down⁃
析,可见该曲流环S为2.31,C为0.0021,反映弯曲程度较大;Δθ为-25°,表明曲率顶点更靠近下游转张期,且扩张程度较大。综合表明河曲29为折端,呈顺流旋转趋势;KM为1.16,表明河弯处于扩
中IR7表明从1984年到2013年,额尔齐斯河研究区
DRES。
第3期 林志鹏等:基于地貌形态学交融的现代曲流河道迁移构型表征
435
图5 下游旋转(顺流旋转)扩张结构迁移构型表征
1984年12月,②摄于1994年12月,③摄于2004年12月,④摄于2013年12月;IR8为对IR7中河曲18迁移结构要素的平面刻画;标为64°41′N,205°34′E,选取2张,其中①摄于2011年3月,②摄于2012年3月;NR8为对NR7中河曲29迁移结构要素的平面刻画;NR9为河曲29的迁移结构表征参数统计数据
IR9为河曲18的迁移结构表征参数统计数据。NR7为诺威特纳河研究区内河曲29近2年来通过GoogleEarth拍摄的卫星图像,坐IR7为额尔齐斯河研究区内河曲18近30年来通过GoogleEarth拍摄的卫星图像,坐标为64°45′N,154°20′E,选取4张,其中①摄于
Fig.5 Planformmigrationarchitectureofthedownstreamrotationexpansionstructure
3.2 收敛结构
曲流河收敛结构(ConstrictionStructure,CS)指河道在向外扩张过程到达一定阶段,开始发生向曲流带轴向截弯侵蚀的过程,表现为河道扩张趋势明显减缓,并在转折端附近侵蚀程度加深,最终截弯废弃。这种结构截弯KM通常小于1。根据曲流环和点坝向上游和下游摆动幅度与Δθ值的差异,可细分为对称收敛结构、上游旋转(逆流旋转)收敛结构和下游旋3.2.1 对称收敛结构转(顺流旋转)收敛结构。
曲流环向外扩张过程减缓,在转折端附近开始侵
河对称收敛结构(SymmetricalConstrictionStructure,SCS)。如图6中IR10表明从1984年到2015年,额尔齐斯河研究区内河曲31河道迁移较缓,通过IR11
迁移结构表征和IR12表征参数的计算分析,可见该曲流环S为3.07,C为0.0008,反映弯曲程度较大;Δθ为2°,形态近似对称,表明曲流环对称性较好,呈微弱逆流旋转趋势;KM为0.87,表明河弯处于收敛NR10中诺威特纳河2011年到2013年的卫星图像显示河曲36河道迁移不显著,通过NR11迁移结构表征和NR12表征参数的计算分析,可见该曲流环S为2.75,C为0.0028,反映弯曲程度较大;Δθ为-4°,形
期,且收敛程度较小。综合表明河曲31为SCS。
蚀截弯,KM小于1,Δθ靠近0°且形态对称,此为曲流
436
沉 积 学 报 第36卷
图6 对称收敛结构迁移构型表征
1984年12月,②摄于1994年12月,③摄于2004年12月,④摄于2015年7月;IR11为对IR10中河曲31迁移结构要素的平面刻画;IR12为河曲31的迁移结构表征参数统计数据。NR10为诺威特纳河研究区内河曲36近2年来通过GoogleEarth拍摄的卫星图像,坐标为64°40′N,205°27′E,选取2张,其中①摄于2011年12月,②摄于2013年3月;NR11为对NR10中河曲36迁移结构要素的平面刻画;NR12为河曲36的迁移结构表征参数统计数据
IR10为额尔齐斯河研究区内河曲31近30年来通过GoogleEarth拍摄的卫星图像,坐标为60°14′N,69°48′E,选取4张,其中①摄于
Fig.6 Planformmigrationarchitectureofthesymmetricalconstrictionstructure
态近似对称,表明曲流环对称性较好,呈微弱顺流旋转趋势;KM为0.96,表明河弯处于收敛期,且收敛程3.2.2 上游旋转(逆流旋转)收敛结构度较小。综合表明河曲36为SCS。
程度较大;Δθ为7°,曲率顶点更靠近上游转折端,呈逆流旋转趋势;KM为0.89,表明河弯处于收敛期,且收敛程度较大。综合表明河曲2为URCS。NR13中30河道迁移不明显,通过NR14迁移结构表征和C为0.0012,反映弯曲程度较大;Δθ为12°,曲率顶表明河弯处于收敛期,且收敛程度较大。综合表明河诺威特纳河2011年到2013年的卫星图像显示河曲NR15表征参数的计算分析,可见该曲流环S为3.86,
曲流环向外扩张过程减缓,在转折端附近开始侵
蚀截弯,且曲率顶点更靠近上游转折端,呈逆流旋转趋势,KM小于1,Δθ大于0°,此为曲流河上游旋转(逆2016年,额尔齐斯河研究区内河曲2河道迁移较缓,通过IR14迁移结构表征和IR15表征参数的计算分析,可见该曲流环S为2.95,C为0.0007,反映弯曲Structure,URCS)。如图7中IR13表明从1984年到流旋转)收敛结构(UpstreamRotationConstriction
点更靠近上游转折端,呈逆流旋转趋势;KM为0.68,3.2.3 下游旋转(顺流旋转)收敛结构曲30为URCS。
曲流环向外扩张过程减缓,在转折端附近开始侵
第3期 林志鹏等:基于地貌形态学交融的现代曲流河道迁移构型表征
437
图7 上游旋转(逆流旋转)收敛结构迁移构型表征
1984年12月,②摄于1994年12月,③摄于2005年8月,④摄于2016年8月;IR14为对IR13中河曲2迁移结构要素的平面刻画;IR15为河曲2的迁移结构表征参数统计数据。NR13为诺威特纳河研究区内河曲30近2年来通过GoogleEarth拍摄的卫星图像,坐标为64°42′N,205°32′E,选取2张,其中①摄于2011年12月,②摄于2013年3月;NR14为对NR13中河曲30迁移结构要素的平面刻画;NR15为河曲30的迁移结构表征参数统计数据
IR13为额尔齐斯河研究区内河曲2近30年来通过GoogleEarth拍摄的卫星图像,坐标为58°56′N,68°50′E,选取4张,其中①摄于
Fig.7 Planformmigrationarchitectureoftheupstreamrotationconstrictionstructure
蚀截弯,且曲率顶点更靠近下游转折端,呈顺流旋转(顺流旋转)收敛结构(DownstreamRotationConstric⁃趋势,KM大于1,Δθ小于0°,此为曲流河下游旋转tionStructure,DRCS)。如图8中IR16表明从1984
结构表征和NR18表征参数的计算分析,可见该曲流环S为3.64,C为0.0034,反映弯曲程度较大;Δθ为-4°,曲率顶点更靠近下游转折端,呈顺流旋转趋势;综合表明河曲50为DRCS。
年到2016年,额尔齐斯河研究区内河曲1河道迁移过程较缓,通过IR17迁移结构表征和IR18表征参数的计算分析,可见该曲流环S为11.05,C为0.0005,反映弯曲程度很大;Δθ为-3°,曲率顶点更靠近下游转折端,呈顺流旋转趋势;KM为0.22,表明河弯处于收敛期,且收敛程度较大。综合表明河曲1为DRCS。NR16中诺威特纳河2010年到2013年的卫星图像显示河曲50河道迁移较缓,通过NR17迁移
KM为0.72,表明河弯处于收敛期,且收敛程度较大。
4 曲流河道地貌迁移构型模式
通过对迁移构型表征的刻画与分析,可以对曲流
河道迁移模式的普适性规律进行探讨。由于河道中心线可近似代替主河道,故本文采用中心线代表河道结构以简化表示曲流河迁移过程和建立迁移构型模式的目的;结合上述构型表征的特征,建立如下曲流
438
沉 积 学 报 第36卷
图8 下游旋转(顺流旋转)收敛结构迁移构型表征
1984年12月,②摄于1994年12月,③摄于2005年8月,④摄于2016年8月;IR17为对IR16中河曲1迁移结构要素的平面刻画;IR18为河曲1的迁移结构表征参数统计数据。NR16为诺威特纳河研究区内河曲50近2年来通过GoogleEarth拍摄的卫星图像,坐标为64°35′N,205°39′E,选取2张,其中①摄于2000年12月,②摄于2013年3月;NR17为对NR16中河曲50迁移结构要素的平面刻画;NR18为河曲50的迁移结构表征参数统计数据
IR16为额尔齐斯河研究区内河曲1近30年来通过GoogleEarth拍摄的卫星图像,坐标为58°52′N,68°47′E,选取4张,其中①摄于
3种基本复合模式(图9)。4.1 扩张迁移模式4.1.1 对称扩张迁移模式
河道理想化条件下迁移构型模式的6种常规模式和
Fig.8 Planformmigrationarchitectureofthedownstreamrotationconstrictionstructure
曲线迁移,LC逐渐增大,|AXMB|保持相对稳定,缓慢大,θD逐渐减小,Δθ大于0°且逐渐增大,Δθ′小于0°
减小,S逐渐增大;R逐渐增加,C逐渐减小;θU逐渐增且逐渐减小;LML保持相对稳定,缓慢减小,KM始终大
移,LC逐渐增大,|AXMB|保持相对稳定,缓慢减小,S
曲流河道不断侧向迁移,曲率顶点近似呈直线迁于1;如图9中UREM,为曲流河上游旋转(逆流旋转)扩张迁移模式(UpstreamRotationExpansionMi⁃4.1.3 下游旋转(顺流旋转)扩张迁移模式gration,UREM),见图4中IR4和NR4。
逐渐增大;R逐渐增加,C逐渐减小;θU和θD基本保持稳定不变,Δθ和Δθ′基本维持在0°左右;LML保持相对稳定,缓慢减小,KM始终大于1;如图9中SEM,此4.1.2 上游旋转(逆流旋转)扩张迁移模式Migration,SEM),见图3中IR1和NR1。
为曲流河对称扩张迁移模式(SymmetricalExpansion
曲线迁移,LC逐渐增大,|AXMB|保持相对稳定,缓慢减小,S逐渐增大;R逐渐增加,C逐渐减小;θU逐渐减小,θD逐渐增大,Δθ小于0°且逐渐减小,Δθ′大于0°且逐渐增大;LML保持相对稳定,缓慢减小,KM始终大
曲流河道不断侧向迁移,曲率顶点向下游方向呈
曲流河道不断侧向迁移,曲率顶点向上游方向呈
第3期 林志鹏等:基于地貌形态学交融的现代曲流河道迁移构型表征
439
图9 理想条件下曲流河道地貌迁移模式
SEM为对称扩张迁移模式,UREM为上游旋转扩张迁移模式,DREM为下游旋转扩张迁移模式,SCM为对称收敛迁移模式,URCM为上游DRE⁃CM下游旋转扩张—收敛迁移模式。A为河道迁移过程简化示意图,B为迁移结构表征参数值随河道迁移的变化趋势示意图旋转收敛迁移模式,DRCM为下游旋转收敛迁移模式,SE⁃CM为对称扩张—收敛迁移模式,URE⁃CM为上游旋转扩张—收敛迁移模式,
于1;如图9DREM,此为曲流河下游旋转(顺流旋gration,DREM),见图5中IR7和NR7。
转)扩张迁移模式(DownstreamRotationExpansionMi⁃
Fig.9 Planformmigrationmodelofthemeanderingchannelundertheidealcondition
4.2.3 下游旋转(顺流旋转)收敛迁移模式
和NR13。
4.2 收敛迁移模式
曲流河道向外扩张过程减缓,在转折端附近开始
4.2.1 对称收敛迁移模式
曲流河道向外扩张过程减缓,在转折端附近开始
增大,|AXMΒ|逐渐减小,S逐渐增大;R逐渐增加,C
侵蚀截弯,曲率顶点向下游方向呈曲线迁移,LC逐渐
侵蚀截弯,曲率顶点近似呈直线迁移,LC逐渐增大,小;θU和θD基本保持稳定不变,Δθ和Δθ′基本维持在此为曲流河对称收敛迁移模式(SymmetricalConstric⁃|AXMB|逐渐减小,S逐渐增大;R逐渐增加,C逐渐减
逐渐减小;θU先减小后增大再减小,θD先增大后减小,KM始终小于1;如图9DRCM,此为曲流河下游旋转(顺流旋转)收敛迁移模式(DownstreamRotationCon⁃strictionMigration,DRCM),见图8中IR16和NR16。Δθ小于且靠近0°,Δθ′大于且靠近0°;LML逐渐减小,
0°左右;LML逐渐减小,KM始终小于1;如图9中SCM,4.2.2 上游旋转(逆流旋转)收敛迁移模式tionMigration,SCM),见图6中IR10和NR10。
4.3 复合迁移模式
河道的迁移演化伴随着地貌、水动力、沉积环境等改变而不断变化,在实际迁移过程中更多是上述6种模式的复合,其中基本的复合模式主要可分为对称扩张—收敛迁移模式(SymmetricalExpansion⁃Con⁃strictionMigration,SE⁃CM)、上游旋转扩张—收敛迁gration,URE⁃CM)、下游旋转扩张—收敛迁移模式(DownstreamRotationExpansion⁃ConstrictionMigra⁃
曲流河道向外扩张过程减缓,在转折端附近开
始侵蚀截弯,曲率顶点向上游方向呈曲线迁移,LC
加,C逐渐减小;θU先增大后减小,θD先减小后增大逐渐减小,KM始终小于1;如图9中URCM,此为曲流河上游旋转(逆流旋转)收敛迁移模式(UpstreamRo⁃tationConstrictionMigration,URCM),见图7中IR13再减小,Δθ大于且靠近0°,Δθ′小于且靠近0°;LML
逐渐增大,|AXMΒ|逐渐减小,S逐渐增大;R逐渐增
移模式(UpstreamRotationExpansion⁃ConstrictionMi⁃
tion,DRE⁃CM)。理论上其他更为复杂的迁移构型都可厘定为这些模式。
440
沉 积 学 报 第36卷
5 讨论
堤岸抗蚀性强度、水流强度、坡度和气候等,在此多变量影响下,河流宏观迁移规律在陆相环境中的呈现方式就尤为重要;什么环境下发育哪种迁移模式,什么因素决定了某种迁移过程,迁移结构对沉积环境具有什么指示作用?而且这种因素对于曲流河和辫状河是具有普适性的,是否能对曲辫转换和共存研究提供新思路?从曲流河迁移构型衍生运用到整个河流沉积迁移体系中去。
以上这些思考都有助于深入解剖现代沉积过程。进一步考虑不同沉积体系下河道迁移构型模式的差在整个分流体系(DFS)中[51⁃52],近物源区河道迁移改道十分迅速,且受阻于局限环境,是否迁移结构较异(图10),这种差异性能否反馈沉积环境的特殊性?
以地貌形态学基础,揭示地质要素历史演化规
律,进而指导认识古代沉积的过程演绎法正成为国内外地球科学研究发展的新趋势,并逐渐渗透到各地质学科。本文研究旨在融入地貌形态学的理论,深入分析和探讨曲流河道迁移演化规律,针对结构要素提出一种新的具体可行的表征参数来定量表征曲流河道迁移构型。这种基于地貌形态学与河流沉积学交融的迁移构型表征法对现代和古代河流沉积的预测与5.1 现代河流迁移构型表征与预测
河道的迁移演化过程受多因素控制,包括流量、重建意义最主要体现在两方面:
图10 不同沉积体系下河流迁移构型模式分布图
河流露头照片源于现代4K超高清分辨率图像截取,A、B摄于冰岛,C摄于美国黄石国家公园,D摄于阿根廷巴塔哥尼亚高原,E摄于智利孔吉利奥国家公园,F摄于美国密西西比河三角洲,G摄于墨西哥尤卡坦半岛
Fig.10 Distributionpatternsofrivermigrationunderdifferentdepositionalsystems
第3期 林志鹏等:基于地貌形态学交融的现代曲流河道迁移构型表征
441
为单一,以SE⁃CM为主;平原区河岸抗侵蚀性相对较差,是否迁移结构多样化,以URE⁃CM和DRE⁃CM复合为主;入海口地区,受控于三角洲河流与潮汐相互作用,迁移结构是否具有特殊性;深海河流体系中,如墨西哥海底的水下河流(图10G),受海底地貌和洋流作用,河道演化是否遵循陆相河流迁移构型模式,5.2 河道迁移模式对地下古河道内部构型演化过程重建的意义
地下古河道储层构型研究虽已蔚然成风,但正如前文所言,由于对地貌迁移规律认知不清,在进行河道内部砂体解剖时,容易采用错误的理论模型,形而2013年发表的地下点坝内部结构解剖与河道演化特上学。此处仅以一例加以说明,如图11中,A为亟待深究。
篇文章的成果图,然而,通过对其河流历史演化的过程还原,如图B;可明显看出矛盾所在,第0期河道期次对应缺失,在第1期河道到第11期河道的演化过程中,河流结构明显与天然河流的演化结构相违背,自然条件下,河流不可能出现像第1期河道这样的迁移结构。由此可见,此处对点坝构型表征的结果整个不符合自然规律,而这还是在基于该区较密的井网条件下得出的结果。无疑将引起质疑,是否整套研究方法与思路都存在问题?如果结合本文的研究思路与方法,结合地貌形态学理论,考虑河道的迁移构型模式,理论上应该得出图C和D结果,如果文章还对古流向或者表征参数有测量,那么还可以对古河道的演化进行量化表征与精确重建。这也同时反映了基于
征分析成果,可参见原文献图7[53]。图A为原文整
Fig.11 Reservoircharacterizationbasedonthepatternofmeanderingchannelmigration
图11 基于曲流河道地貌迁移构型模式的储层表征
442
流储层构型的重要指导、重建与预测意义。
沉 积 学 报 第36卷
地貌形态学交融的曲流河道迁移构型表征方法对河
基础理论架构的不足很可能导致方法甚至最终结果的错误,并进一步误导深入研究,2017年最新发表的成果中还依然存在该案列呈现的问题,所以本文的研究十分必要。张昌民等[54]通过回顾20年来河流沉积学的研究也指出国内河流相研究重视应用而轻视基础理论研究的问题。今后的研究也应当积极重视基础研究,贯彻从实践到理论,理论到实践的思提出曲流河道理想化条件下迁移构型的6种常规模式和3种基本复合模式,分别是对称扩张迁移模式(SEM)、上游旋转(逆流旋转)扩张迁移模式(DREM)、对称收敛迁移模式(SCM)、上游旋转(逆流旋转)收敛迁移模式(URCM)、下游旋转(顺流旋转)(SE⁃CM)、上游旋转扩张—收敛迁移模式(URE⁃CM)和下游旋转扩张—收敛迁移模式(DRE⁃CM)。收敛迁移模式(DRCM)、对称扩张—收敛迁移模式(UREM)、下游旋转(顺流旋转)扩张迁移模式
想路线,将河流沉积学与地貌学联系起来,相互印证指导,以期进一步推进河流研究领域各学科间的交融 与共同发展 然而,本文也具有一定局限性。
家典型天然河流地貌构型进行识别与研究,虽然针对两个国,但对之世界繁不甚数的河流而言,只是冰山一角,且迁移构型模式的表征方法过于理想化,在实际运用中必须加以考虑具体因素具体分析。并且,由此衍生出的以下四点问题将在今后工作中加以深入研究:1)不同迁移模式之间的内在联系,2)迁移构型与垂向沉积结构4)之间的对应关系,3)河道迁移模式形成的主控因素,
本研究所得出的构型表征法与迁移模式对于河流的曲流河迁移结构中平移结构存在的争议。因此,普遍应用,仍需在今后更多实践工作基础上加以深入佐证。
6 结论
Google通过对额尔齐斯河和诺威特纳河典型河曲采用
合地貌形态学与河流沉积学Earth和ACMEMapper,对曲流河道迁移构型表现代卫星图像分析,结征结构主要得出如下几点认识提。(1)曲流河地貌形态结构表征可分细为刻画河道迁移结构是分析其演化规律的前
:
17个静态要素、11个动态要素,深入定义了弯度指数与曲率,提出并运用了3种新构型表征参数:顺流偏转角、逆流偏转角和扩张系数(2)5个表征参数。
S并据此综合提炼出、C、Δθ、Δθ′和KM可对曲流
河迁移结构做定量表征,6种迁移结构,分别是对称扩张结构(SES)、上游旋转(逆流旋转)扩张结构(URES)、下游旋转(顺流旋转)扩张结构(DRES)、对称收敛结构(SCS)、上游旋转(逆流旋转)收敛结构(URCS)和下游旋转(顺流旋转)收敛结构(DRCS)。
(3)通过对研究区12段河曲的精细结构表征,
艾,而跨学科结合地貌形态学理论综合表征河流演化(4)河流沉积学与地下储层构型的研究方兴未
结构的发展势在必行,在现代河流迁移构型表征预测和地下古河道过程重建的讨论中,进一步明确将来发展要结合地貌形态学规律,综合多因素动态分析河流演化过程,以辩证思想真正实现将今论古、以古鉴今的科学研究。
参考文献(References)
[1] BlumfromquaternaryM,MartinanalogsJ,MillikenandK,experimentsetal.Paleovalley[J].systems:insights
[2] views,Kasvi2013,116:128⁃169.
Earth⁃ScienceRe⁃logicalE,andtemporallychangesLaamanenintensiveinaL,sandyLotsariADCPmeanderE,etal.Flowpatternsandmorpho⁃
measurementbendduringapproach[aflood:J].spatially
Water,
[3] 2017,LinZP,9(2):Shan106.
planformJF,ChenL.Geomorphologyprocessesof[4] caLin(EnglishmigrationEdition),onmeandering2017,91:rivers[J].134⁃135.
ActaGeographicachannel
Sini⁃anismZP,ofunconformityShanJF,ChenbeltL.structure[MeticulousJ].depictionEarthScienceandgeneticResearch,
mech⁃
[5] 2017,Schuurman6(2):19.
disturbancesF,inKleinhanslargesand⁃bedMG,braidedMiddelkooprivers[J].H.Networkresponseto
[6] namics,EarthSurfaceDy⁃胡光义流相高精度概念模型的表征方法,2016,范廷恩4(1):,陈飞,25⁃45.
等.从储层构型到[J].地质学报“地震构型相,2017,91(2):
”:一种河
465⁃478.ervoirodofaarchitecture[HuGuangyi,FanTing’En,ChenFei,etal.Fromres⁃high⁃resolutiontoseismicarchitecturefacies:characteristicGeologicameth⁃ca,2017,91(2):465⁃478.]
fluvialfaciesmodel[J].ActaSini⁃
[7] 胡浩藏,2016,28(.基于砂体结构的剩余油挖潜调整措施研究4):113⁃120.[HuHao.Adjustment[J].measures岩性油气
ofre⁃mainingoiltappingbasedonsandbodystructure[J].Lithologic
[8] 穆龙新Reservoirs,2016,28(4):113⁃120.]
103⁃108..油藏描述的阶段性及特点[9] scription[J].[Mu[J].石油学报,2000,21(5):
吴胜和,王仲林ActaLongxin.Stagesandcharacteristicofreservoirde⁃.陆相储层流动单元研究的新思路PetroleiSinica,2000,21(5):103⁃108.]
[J].沉积学
第3期 林志鹏等:基于地貌形态学交融的现代曲流河道迁移构型表征
报,1999,17(2):252⁃257.[WuShenghe,WangZhonglin.AnewlogicaSinica,1999,17(2):252⁃257.]
methodofnon⁃marinereservoirflowunitstudy[J].ActaSedimento⁃
JohnWiley,1997:265⁃280.
443
[25] MiallAD.Architectural⁃elementanalysis:anewmethodoffacies
1985,22(4):261⁃308.
analysisappliedtofluvialdeposits[J].Earth⁃ScienceReviews,
[10] 尹太举,张昌民,赵红静,等.依据高分辨率层序地层学进行剩
Taiju,ZhangChangmin,ZhaoHongjing,etal.Remainingoildis⁃
余油分布预测[J].石油勘探与开发,2001,28(4):79⁃82.[Yin
[26] MiallAD.Thegeologyoffluvialdeposits:sedimentaryfacies,ba⁃
Springer,1996.
[11] DebM,DasD,UddinM.Evaluationofmeanderingcharacteristics
82.]
[J].PetroleumExplorationandDevelopment,2001,28(4):79⁃
tributionpredictionbasedonhigh⁃resolutionsequencestratigraphy
sinanalysis,andpetroleumgeology[M].BerlinHeidelberg:
[27] MiallAD.Fluvialdepositionalsystems[M].Cham:Springer,[28] MiallAD.Faciesmodels[M]//MIALLAD.Stratigraphy:A2014.
Protection,usingRS&2012,GISof4(3):ManuRiver[J].163⁃171.
JournalofWaterResourceand
隋新光油学院.,2006.曲流河道砂体内部建筑结构研究[SuiXinguang.Astudyoninternal[D].大庆architecture:大庆石
of
Institute,channelsandWillisBJ,2006.]
inmeanderingriver[D].Daqing:DaqingPetroleumdeposits[J].TangJournalH.ofThree⁃dimensionalSedimentaryResearch,connectivity2010,ofpoint⁃bar
440⁃454.
80(5):
薛培华社,1991..河流点坝相储层模式概论[XuePeihua.Anintroduction[M].to北京reservoir:石油工业出版
modelsof
pointIelpibarfacies[M].Beijing:PetroleumIndustryPress,1991.]
andmorphodynamicsA,GhinassiM.ofPlanformanexhumedarchitecture,stratigraphicsignature
by(7):Formation,1923⁃1960.
Yorkshire,UK)[J].JurassicSedimentology,meanderplain2014,(Scal⁃
61LeopoldchannelsLandB,someMaddockphysiographicTJr.ThehydraulicgeometryProfessionalofstream
252[M].WashingtonDC:U.S.implications:GovernmentPrintingOffice,paper
Leopold1953.
andering,LandB,Wolmanstraight:MProfessionalG.Riverchannelpaper282⁃B[M].patterns:braided,Washington
me⁃
DC:LeopoldU.S.GovernmentPrintingOffice,1957.
insinuousLB,orBagnoldRA,WolmanMG,etal.Flowresistance
WashingtonDirregularC:U.S.Governmentchannels:ProfessionalPrintingOffice,paper1960.282⁃D[M].
WolmanomorphicMprocesses[G,MillerJ].JP.TheMagnitudeJournalandofGeology,frequency1960,offorces68in(1):
ge⁃
54⁃74.
WilliamsResearch,G1978,P.Bank⁃full14(6):discharge1141⁃1154.
ofrivers[J].WaterResourcesBriceJC.Evolutionofmeanderloops[J].GSABulletin,1974,Hooke85(4):JofM.581⁃586.
resultsanalyses[ChangesJ].inriverProgressmeanders:inPhysicalareviewGeography,oftechniques1984,and
(4):473⁃508.8
HookechannelsJM.intheProcessesUK[M]ofchannelplanformchangeonmeandering
Channels.Chichester:John//GurnellWileyandA,Sons,PettsG1995:E.Changing87⁃115.
River
Hookepattern[JM]M.The//GregorydistributionKJ.andRivernatureChannelofchangesChanges.inriverChichester:channel
[29] GhinassiModernSynthesis.M,NemecCham:W,AldinucciSpringer,M,et2016:al.161⁃214.Plan⁃formevolutionoutcropof
[30] sections[J].ancientmeanderingWurisingCandL,UllanSedimentology,riversreconstructedfallingMfloodS,Lustages:J,et2014,evidenceal.Formation61(4):from952⁃977.
longitudinalfrombarofpointmorphology,barsthrough
sedi⁃(6):menttransport1458⁃1473.
andbedshearstress[J].Sedimentology,2016,63
[31] 裘亦楠1985,12(2):72⁃74..河流沉积学中的河型分alsedimentology[J].[QiuPetroleumYinan.River类[ExplorationtypeJ].classification石油勘探与andDevelopment,
in开fluvi⁃发,
[32] 1985,尹燕义12(2):72⁃74.]
[J].,王国娟,祁小明.曲流河点坝储集层侧积体类型研究
Guojuan,石油勘探与开发ofthemeanderingQiXiaoming.,1998,25(2):37⁃40.riverpointAstudy[YinYanyi,Wangbaronthelateralaccretionbodytype[33] tion岳大力andDevelopment,1998,25(2):reservoirs[J].37⁃40.]
PetroleumExplora⁃方法[J].,吴胜和石油,学刘建民报,2007,.曲流河点坝地下储层构型精细解剖
28(4):99⁃103.[YueDali,WuShenghe,tectureLiuJianmin.Anaccuratemethodforanatomizingarchi⁃ActaPetroleiofsubsurfaceSinica,reservoir2007,28(4):inpoint99⁃103.]
barofmeanderingriver[J].
[34] 周银邦析及识别方法,吴胜和[J].,岳大力中国石油大学学报,等.点坝内部侧积层倾角控制因素分
(自然科学版),2009,33(2):7⁃11.trollingbeddingsfactor[ZhouanalysisYinbang,angleinandWuShenghe,YueDali,etal.Con⁃pointidentificationbar[J].JournalmethodofofChinalateralUniversity
accretion[35] ofshale单敬福Petroleum,2009,33(2):7⁃11.]
分析:以吉林油田杨大城子油层第,张吉,赵忠军,等.地下曲流河点坝砂体沉积演化过程
23小层为例[J].石油学报,2015,36(7):809⁃819.[ShanJingfu,ZhangJi,ZhaoZhongjun,undergroundetlayerofmeanderingal.AnalysisYangdachengziriverofsedimentaryoilpointreservoirbar:andevolutionprocessforinAcaseJilinstudyoilfield[fromJ].No.23PetroleithinSinica,2015,36(7):809⁃819.]
Acta
[36] Linrivers[J].ZP.ChannelJournalofplanformGeologymigration&Geophysics,architecture2017,of6(S3):meandering[37] Gilvear63.
planformD,bia[J].changeWinterbottomandmeanderS,SichingabulaH.Characterofchannel
421⁃436.
EarthSurfaceProcessesdevelopment:andLandforms,Luangwa2000,River,25(Zam⁃4):[12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] 444
沉 积 学 报
sion[J].EarthSurfaceProcessesandLandforms,1980,5(2):
PetroleumIndustryPress,2013.]
第36卷
[38] HookeJM.Magnitudeanddistributionofratesofriverbankero⁃
143⁃157.
[47] 何幼斌,王文广.沉积岩与沉积相[M].北京:石油工业出版
社,2008.[HeYoubin,WangWenguang.Sedimentaryrocksandsedimentaryfacies[M].Beijing:PetroleumIndustryPress,
[39] WuCL,BhattacharyaJP,UllahMS.Paleohydrologyand3Dfa⁃
delta,South⁃CentralUtah,U.S.A.[J].JournalofSedimentaryRe⁃
ciesarchitectureofancientpointbars,ferronsandstone,notom
[40] GhinassiM,IelpiA,AldinucciM,etal.Downstream⁃migrating
2016,334:66⁃96.
search,2015,85(4):399⁃418.
[48] 朱筱敏.沉积岩石学[M].4版.北京:石油工业出版社,2008.
leumIndustryPress,2008.]
2008.]
[ZhuXiaomin.Sedimentarypetrology[M].4thed.Beijing:Petro⁃
fluvialpointbarsintherockrecord[J].SedimentaryGeology,[49] BluckBJ.Sedimentationinthemeanderingriverendrick[J].[50] 石书缘,胡素云,冯文杰,等.基于GoogleEarth软件建立曲流
ScottishJournalofGeology,1971,7(2):93⁃138.
[41] 吴胜和,岳大力,刘建民,等.地下古河道储层构型的层次建模
研究[J].中国科学D辑:地球科学,2008,38(增刊1):111⁃121.eling[WuofsubsurfaceShenghe,palaeochannelYuDali,LiuJianmin,etal.Hierarchy[J].mod⁃
Sci⁃单敬福111⁃121.]
enceinChinaSeriesD:EarthSciences,reservoir2008,architecture38(Suppl.1):
建方法,[J].李占东中国矿业大学学报,葛雪,等.一种古沉积期曲流河道演化过程重
,2015,44(5):843⁃852.[ShanJingfu,evolutionaryLiZhandong,processofGemeanderingXue,etal.channelAmethodinpaleo⁃sedimentary
ofreconstructing43] Asahi2015,period[44(5):J].Journal843⁃852.]
ofChinaUniversityofMining&Technology,
meanderingK,Shimizuself⁃evolvingY,NelsonJ,banksetal.[NumericalJ].JournalsimulationofGeophysical
ofriver
Research:SchuurmanEarthwithSurface,2013,118(4):2208⁃2229.
ingF.Barandversity,riversusingphysics⁃basedchannelmodeling[D].evolutioninmeanderingUtrecht:Utrechtandbraid⁃
Uni⁃Rousseau2015.
tedflowfieldsYY,BironandbathymetriesPM,vanDeWielMJ.Sensitivityofsimula⁃
model[J].inmeanderingEarthSurfacechannelsProcessestoand
the
冯增昭Landforms,choiceofamorphodynamic[FengZengzhao..中国沉积学2016,41(9):Sedimentology[M].1169⁃1184.
2版.of北京China:石油工业出版社[M].2nded.Beijing:
,2013.
河地质知识库[J].沉积学报,2012,30(5):869⁃878.[Shi
knowledgeShuyuan,HudatabaseSuyun,basedFengongoogleWenjie,earthetsoftware[J].al.Buildinggeological
[51] mentologicaHartleyfluvialAJ,Sinica,Weissmann2012,G30(5):S,Nichols869⁃878.]
ActaSedi⁃Gvelopment[systems:167⁃183.
J].Journalcharacteristics,ofSedimentarydistribution,J,etal.Research,andLargecontrolsdistributive
2010,80(on2):
de⁃[52] WeissmannmodernGS,HartleyAJ,NicholsGJ,etal.Fluvialformin
[J].张昌民Geology,continental2010,sedimentary38(1):39⁃42.
basins:distributivefluvialsystems
[53] 储层内部建筑结构分析,尹太举,喻辰,等:.以大庆油田萨北地区为例基于过程的分流平原高弯河道砂体
[J].沉积学报,2013,31(4):653⁃662.[ZhangChangmin,YinTaiju,Yu
Chen,sandstoneetal.[54] Acta张昌民SedimentologicainReservoirthedeltaarchitecturalanalysisofmeanderingchannelSinica,plainbased2013,on31(4):thedepositional653⁃662.]
process[J].沉积学,报张尚锋,2004,,李少华22(2):,等183⁃192..中国河流沉积学研究[ZhangChangmin,20年Zhang[J].
entologyShangfeng,fromLi1983Shaohua,to2003etal.[JAdvances].ActainSedimentologicaChinesefluvialSinica,sedim⁃
2004,22(2):183⁃192.]
[42] [[44] [45] [46] 第3期 林志鹏等:基于地貌形态学交融的现代曲流河道迁移构型表征
445
GeomorphologyProcessesofChannelPlanformMigrationonModernMe⁃anderingRivers
1.KeyLaboratoryofExplorationTechnologiesforOilandGasResources,MinistryofEducation,YangtzeUniversity,Wuhan430100,China3.SchoolofEnergyResources,ChinaUniversityofGeosciences,Beijing100083,China4.SchoolofEarthandSpaceSciences,PekingUniversity,Beijing100871,China
LINZhiPeng1,2,SHANJingFu1,2,CHENLe1,2,SUNQianJun3,WANGYiWu4,WANGFei5
2.SchoolofGeosciences,YangtzeUniversity,Wuhan430100,China
5.GasProductionCompanyofDaqingOilfieldLimitedCompany,Daqing,Heilongjiang163000,China
Abstract:Geomorphologicalmigrationstructuresarethebasistoreconstructandrestoresedimentaryevolutionofthemorphologymigrationstructuresofmeanderingrivers,high⁃resolutionhistoricalsatelliteimagesareacquiredprimarilythroughGoogleEarthandACMEMapper.ThroughthefinefeaturecharacterizationfortheIrtyshRiverandtheNowit⁃naRiver,whicharebothwellinpreservationconditionsofnaturalstructures,28kindsofstructureelementsareuti⁃lizedtodepictthemeanderingchannels.Moreover,5kindsofcharacterizationparametersareproposedtomakequan⁃thedifferenceofcounter⁃currentdeflectionangle(Δθ′),andexpansioncoefficient(KM)arefirstlybroughtforwardS,C,Δθ,Δθ′,KMfactorsareadoptedtosummarizefor6kindsofplanformmigrationstructures:symmetricalexpan⁃titativedescriptiononthechannelstructure,extraordinarily,thedifferenceofalong⁃currentdeflectionangle(Δθ),frompreviousoriginaldefinitionandmoremeticulous.Byfineanatomyofthestructurefor12typicalmeanderloops,sionstructure,upstreamrotationexpansionstructure;downstreamrotationexpansionstructure;symmetricalconstric⁃tivecharacterization,theevolutionofthemeanderingriverhasacertaindegreeofuniversallawandguidingroleinor⁃paleochannel,aswellasthepremisetoforecastmigrationtrendofrivers.Inordertofigureouttheregularityofgeo⁃
andappliedforintroduction.Inaddition,theconceptionofsinuosityindex(S)andcurvature(C)arealsodifferent
tionstructure;upstreamrotationconstrictionstructure;downstreamrotationconstrictionstructure.Ultimately,9mean⁃dertopromotethelandscapemorphologyandfluvialsedimentologyofcommondevelopment.ters
deringchannelmigrationpatternsareconcludedanddiscussed.Combinedwiththeideaofgeomorphologicalquantita⁃Keywords:meanderingchannel;Geomorphology;migrationpatterns;structureelements;characterizationparame⁃
因篇幅问题不能全部显示,请点此查看更多更全内容