陈 浩， 董廷旭， 李 勇， 文星跃

（1.绵阳师范学院 资源环境工程学院，四川 绵阳621006； 2.成都理工大学 油气藏地质及开发工程国家重点实验室，四川 成都610059；3.西华师范大学 国土资源学院，四川 南充637002）

水系对于构造、气候等外部变量因素的改变非常敏感，水系演化通常记录了造山带系统新近时期的细微变动，因而，地质构造活跃区域水系发育模式及其对构造活动、气候条件的响应机制一直以来是构造地貌与河流地貌研究的热点［1］.大致呈南北向展布的岷江上游流域盆地，位于第四纪以来活动较为强烈的岷山构造带西缘和龙门山构造带北侧，虽然一些学者对岷江水系特征进行了较为深入的分析，但研究重心主要集中于岷江下切过程、河流阶地样式与区域隆升的耦合关系［2-10］，有关岷江上游流域地貌发育期次的研究目前还未见报道.近年来，地理信息系统技术的不断发展和数字高程模型（DEM）数据的应用使地貌过程研究由定性向定量化阶段转变［11］.基于DEM数据提取符合实际状况的水系网络是岷江上游流域地貌发育阶段研究的前提，最佳集水面积阈值的确定是水系网络提取过程中的关键环节，目前，基于DEM确定最佳集水面积阈值的方法主要包括流域宽度法、河网密度法、分形法和二阶导数法等［12-17］，虽然这些方法可以较准确地提取最佳河网阈值，但是仍然存在一定的主观性及不便捷性.常直杨等［18］通过均值变点法获取了青藏高原东缘杂谷……