水下扇沉积过程模拟方法研究.docxVIP

水下扇沉积过程模拟方法研究 地质沉陷过程的数字模拟的研究正在早期进行，但由于地质现象的复杂性、多变性、计算机操作速度、存储、图形显示性能等的影响，其研究成果与实际情况仍有很大差距。随着计算机工业的发展,已能轻而易举地解决早期不可能实现的海量数据处理和沉积过程动态显示等问题。如果能在沉积过程的概念和数学模型上有所突破,沉积过程模拟研究将成为一种有效的定量研究手段。 水下扇作为中国东部中、新生代断陷盆地内常见的沉积体,是油气聚集的重要场所。对其进行沉积过程数学模拟,有助于了解水下扇的形成过程,并能为沉积相、地层对比、岩性分布、岩石物性等方面的研究提供定量依据。 1 成岩阶段及沉积过程 水下扇属于河、湖过渡沉积类型,在断陷湖盆陡坡带,间歇性河流或洪水携带大量的陆源碎屑物质在湖盆中呈扇形沉积下来形成水下扇。 在水下扇体漫长的形成过程中,物源是间歇性的,不能按连续性的物源供给来描述碎屑物质进入湖盆后的搬运、沉积演化过程,而应考虑将水下扇体的形成过程归纳为若干个主要沉积时期。每个时期的洪水或河流均携带一定量的碎屑物质进入湖盆,进一步搬运并沉积下来。因此,可将水下扇体的形成描述为在多次物源条件下,不同粒度级物质搬运、沉积的循环过程,数学模拟的关键是在此基础上建立物质搬运模型和沉积模型。 2 物质的扩散过程 陆源碎屑物进入湖盆后,受水动力、重力、基准面、坡度等因素的影响继续向湖盆内搬运,这一过程可理解为物质浓度的扩散过程。Fick第二定律描述的物质的确定性扩散过程为 其中,k为扩散系数,c为物质浓度,x,y,z为坐标轴。若忽略z方向上的扩散,并考虑到x,y方向的不同搬运能力,则式(1)变为 式(2)中kx,ky为x,y方向的扩散系数,用其近似描述进入湖盆的物质随时间的搬运过程。 3 物质沉积层模型 3.1 顺序计算沉积量 物质的沉积量和浓度密切相关。若已知湖盆内某位置的物质浓度,按先沉积粗粒、后沉积细粒的顺序计算出不同粒级物质的沉积量,而且各种粒级物质中只有一部分沉积下来 (与沉积系统有关),剩余物质继续进行扩散。在单位面积内,如果将某粒级物质的浓度理解为到达该位置的物质厚度(密度为ρs)并假设为v,设该粒级物质的沉积系数为cs,沉积厚度为h,则有 3.2 不同沉积厚度沉积的控制作用 沉积作用受基准面(波基面)的控制,在基准面以上不能接受沉积。设湖盆某位置的水深为d,基准面深度为b,最大沉积物厚度(到达的 物质浓度与沉积系数之积)为L,实际沉积厚度为H,大致分3种情况处理基准面对沉积的控制作用: (1)若沉积底界和水的分界面高于或等于基准面时,任何物质不可能沉积,继续向其它位置扩散(搬运)。此时有d≤b,实际沉积厚度为0,剩余量为L。 (2)若湖盆某位置沉积底界与水的分界面至基准面的厚度不能容纳可能的总沉积厚度时,则沉积最多至基准面,剩余部分继续向其它位置扩散(搬运)。此时有L≥d-b,仅有部分沉积,实际沉积厚度为d-b,剩余量为L-(d-b)。 (3)若湖盆某位置沉积底界与水的分界面至基准面的厚度能够容纳可能的总沉积厚度时,则接受全部沉积。此时有L≤d-b,实际沉积厚度为L,剩余量为0(图1)。 3.3 沉降作用的稳定性 由于湖盆的整体沉降使得水深随时间而增加,进一步又可接受新的沉积,所以在模型中可通过增加各模拟单元水深值的方法描述湖盆的沉降作用。亦即给定各单元内的水深随时间的一个增量,即可描述湖盆随时间而产生的整体均匀沉降。若用一个简单的比例常数把沉降和沉积作用联系起来,各单元的沉降量等于沉积厚度与比例常数之积,则可用于描述湖底沉积地层由于压实作用而产生的形变。特殊情况下,若比例系数为1,则沉积量等于沉降量;若比例系数为0,则不出现沉降。采用类似的方法描述湖盆的上升作用时,湖盆的水深随时间而变浅。如果使沉降作用落后于沉积作用一定的时间增量,则产生滞后沉积效果。 3.4 湖盆深水下降或减少 水进和水退作用的结果造成湖盆沉积空间的增大或减少,同时湖盆水深相应发生变化。若假设湖盆整体呈均匀沉降或上升,则可以通过增加或减少不同位置网格单元的水深来描述水进和水退。 4 元长制网格单元 在一个包含水下扇沉积范围的矩形面积内,对工区进行网格划分,以求取各网格单元内各时期的物质浓度和沉积量。为模拟计算方便,一般使用单位面积的正方形网格,每个网格赋予行、列编号,各网格单元的控制面积依实际情况而定。网格单元包括界外(湖盆边界以外)、左上角、上边界、右上角、右边界、右下角、下边界,左下角、左边界、界内(湖盆边界以内)10种类型,分别用数字0～9表示。通过网格类型的设定,应该能看出工区的平面形态。 5 主要模拟参数 (1) 每个网格的原始深度 根据对湖盆的有关研究确定,反映湖盆沉积之前的原始形态。 (2) 平面初始状态 在湖盆周边确定某些网格单元为物源入口