用沉積理论研究黄河的造陆贡献.docVIP

基金项目：河南省科技攻关项目(20025700002)。 用沉积理论分析黄河的造陆过程 杨道富1, 刘仕平 （1.黄河水利职业技术学院，河南 开封 475001 2.华北水利水电学院 机械系,河南 郑州 450008） 摘 要：从地层沉积物孔隙度随深度变化的关系模型入手，地层沉积厚度脱压恢复的目的有两个：①求算沉积速率，应采用“逐层剥皮脱压法”；②定量分析下伏地层沉降速率或沉降幅度，应采取“地层单元回剥法”。分析沉积速率的计算公式，选择牛顿法解渤海湾盆地沉积速率，研究了中国大陆东部边缘渤海湾盆地的沉积过程，盆地区域内沉积速率的变化存在3个高峰值，其年龄分别为：距今0~2Ma、24~32.6Ma、38~42Ma，为认识古代环境与地质演变提供了科学的依据。论述了黄河的造陆过程，第四纪以来黄河对渤海湾盆地沉积贡献最大。约距今10万年距今10万1万年间 （1） 式中：为随深度变化的孔隙度；为该种岩石在地表的孔隙度（初始孔隙度）；c为压实曲线的斜率；e为自然对数。 某一钻井柱状地层中，岩性在垂直方向上的变化是复杂的，因此有学者对上述公式进行了修正，提出综合岩性孔隙度计算公式[3]。 （2） 式中：fi表示第i种岩石在某段地层中所占的比例。 这一公式需要详细的岩性划分及组成资料。考虑实际工作中的困难，以及脱压计算过程中地层单元划分的实际意义，仍然采用式（1），对不同岩性的影响，建议采用调整初始孔隙度和压实曲线斜度c的数值进行修正，即不同的岩石具有不同的初始孔隙度和压实曲线斜率，这似乎更符合实际情况。 2 地层厚度脱压恢复分析 沉积物在地质历史时期，在深部可能经历比较复杂的变化，鉴于这些变化的复杂性和研究程度，目前“脱压恢复”在计算中，一般只考虑机械压实作用，即假定压实作用只减小沉积物的孔隙体积，压实前后岩石颗粒骨架体积不变，用数学公式表示为 （3） 式中：x1和x2分别代表地层单元顶界和底界埋深，和分别代表压实前（脱压后）地层单元顶界和底界埋深。 在正常压实情况下，式（3）可以表示为 （4） 因为x1和x2是已知的，上式左边可以计算出确定值，不妨设，脱压后的地层单元厚，那么 (5) 设脱压后（压实前）该地层单元的顶界埋深为0，即，那么，，整理后可以得出： （6） 用牛顿迭代方法解方程（6），求出脱压后的地层厚度D’。 地层沉积厚度脱压恢复的目的有两个：①求算沉积速率，应采用“逐层剥皮脱压法”。②定量分析下伏地层沉降速率或沉降幅度，应采取“地层单元回剥法”。， 3 牛顿法解渤海湾盆地沉积速率 对钻孔柱状地层进行上述全脱压校正以后，计算沉积速率就是求某一地层单元脱压后的总厚度与该地层单元底界和顶界地质年代之差的比值，求出这个比值即可求出该地层单元沉积期间的平均沉积速率[4]。用公式表示为 V=D/(Ai-1-Ai) 式中：V为沉积速率；D为某地层单元脱压后的总厚度；Ai-1为该地层单元底界年龄；Ai为该地层单元顶界的年龄。 根据石油地质钻井资料，采用牛顿法解上述方程子函数，省略具体计算过程，分析结果的真实性和使用价值，进行时段划分，结果见表1。 表1 渤海湾盆地沉积速率计算结果 (m/Ma) 4 渤海湾盆地沉积结果分析 运用上述方法对我国大陆东部边缘渤海湾盆地的沉积速率进行了定量计算。由于涉及石油地质资料的可用范围，在这里，我们忽略了具体钻井的位置和详细的地质资料细节，只分析沉积速率的时空变化趋势。 由于地区之间或钻孔之间各时段的沉积速率受到当地的沉积环境条件的影响，表现较为复杂，数值相差较大，详细的时段划分与对比，不能看出微观细节变化，但是，将渤海湾盆地沉积绘成示意图进行比较，从宏观上分析，总的趋势表现还是比较明显的，见图1。 图图1 渤海湾盆地沉积比较图(由表1的数值绘制，纵坐标m/Ma，横坐标Ma) 从渤海湾盆地60个点的平均计算结果可以看出，自始新世至今，沉积速率是逐渐增大的。距今42Ma（百万年）以前的沉积速率很小；距今38~42Ma沉积速率出现峰值，反映了盆地开始大规模的沉降和接受沉积；距今24~32.6Ma沉积速率又出现峰值；距今24Ma以来，沉积