1D应变速率反演的A MATLAB程序
Hai-Bin Song,Lin Chen, Jiong Zhang, Chang Yu Zhao , Chong-Zhi Dong
1.导言
过去30多年来,裂谷和被动大陆边缘的演化一直受到普遍关注。构造上活动的裂谷、古裂谷和被动大陆边缘在不同类型的沉积盆地中形成了一类与伸展有关的重要盆地组合(Bally and Snelson, 1980)。伸展盆地覆盖了全球相当大的面积,并且蕴藏了重要的沉积矿产和能源资源。大量的重要含油气省与裂谷和被动大陆边缘有关(White et al., 2003;Ziegle and Cloetingh,2004)。为了进一步认识裂谷盆地的演化,一些学者已提出了几种模式(McKenzie, 1978; Royden and Keen,1980; Hellinger and Sclater, 1983; Rowley and Sahagian, 1986)。在模拟沉降作用时,伸展作用一般被假定为瞬时的(McKenzie, 1978),或者是以一定的速率在有限的时期内发生(Jarvis and McKenzie, 1980)。White(1993, 1994)提出了一种利用沉降资料反演的方法,该方法现已应用于许多不同的盆地(Newman and White, 1999; Xie et al., 2006)。
在本文中,我们介绍的研究成果,主要是针对已有模型的重要进展方面,为研究者研究裂谷盆地的演化提供有用的工具。该项工作的目的是提供一个开放的MATLAB程序,以用于根据沉降资料反演应变速率和估算伸展系数。
图1 岩石圈随时间伸展的物理模型示意图。
a=岩石圈厚度;T1=软流圈温度;u(x)=水平速度; a(z)=垂直速度。
2.物理模型
我们所考虑的2D物理模型如图1所示(其细节参见Jarvis and McKenzie (1980)资料)。地壳和岩石圈地幔以水平速度u(x)进行伸展,与此同时,软流圈地幔向上流动穿过平面z=0,进而取代向外流出的岩石圈物质。上部表面即z=a,这里的a是岩石圈厚度,上部表面温度为T=0℃;平面z=0处的温度保持为T=T1,也就是软流圈的温度。这种边界条件显然是大致的,因为岩石圈的顶面伴随着伸展作用的进行一定会发生沉降,并且也会发生沉积。变形作用限定为纯剪应变方式,并且垂直速度在z=0处为G(t)a,在这里G(t)为垂直应变速率,它是时间t 的函数。垂直速度要求在z=a处消失,因此可表达为: