六方各向异性介质本构方程课件.ppt

《六方各向异性介质本构方程课件.ppt》由会员分享，可在线阅读，更多相关《六方各向异性介质本构方程课件.ppt（48页珍藏版）》请在三一办公上搜索。

1、第二章 六方各向异性介质本构方程,各向异性介质的物性是通过物性矩阵来刻画的。沉积型地层中普遍存在的基本类型是六方各向异性介质，六方各向异性中主要有 TI、EDA 和 PTL 三种形式。六方各向异性介质中存在一个对称主轴，本章讨论的六方各向异性介质，其主轴皆为 Z 轴；与主轴正交的平面为各向同性面。在各向同性面上，介质的速度是同性即相同的。对称主轴空间分布的任意性可以引出方位各向异性介质。坐标旋转是各向异性介质研究中的一个特殊的问题。坐标旋转的核心问题是构造 6 阶邦德变换矩阵。,本章主要内容,六方介质三种典型的物性矩阵应力应变和物性矩阵的坐标变换方位六方各向异性介质,2.1 六方介质三种典型的

2、物性矩阵,六方各向异性介质本构方程和物性矩阵:1、横向各向同性介质的物性矩阵2、扩张型裂隙各向异性介质的物性矩阵3、薄互层各向异性介质的物性矩阵,2.1.1 六方各向异性介质本构方程 和物性矩阵,均匀弹性六方各向异性介质本构方程,取单双角标关系：,说明：1.独立的分量共5个2.Z轴为主轴的柱状对称介质，是六方介质的普遍特性3.规定Z轴垂直向下,物性矩阵,是3阶零矩阵,和 分别是 左上和右下三阶子矩阵,的分块矩阵形式为,2.1.2 横向各向同性介质的物性矩阵,（1）VTI介质,另外，是3阶零矩阵，上式Z轴为对称主轴.因此，也称其为VTI介质(Vertical Transversely Isatm

3、pical)，即对称轴是垂向的横向各向同性介质，简称为垂向横向各向同性介质;或垂向六方各向异性介质。分别是与介质水平性质和垂向性质有关的拉梅系数。,式中，是介质的密度；是qP波垂直于对称轴即在各向同性面传播的速度；是qP波平行于对称轴传播的速度；是与qP波有关的过渡速度；是qS波垂直于对称轴传播的速度；是qS波平行于对称轴传播的速度。为了与各向同性介质中的纵波和横波区别，这里用qP和qS分别表示视纵波和视横波,也称为准纵波和准横波。上述五个速度称为本构方程的初态速度。且有,用这5个参数可以引出5个速度,（2）HTI介质 这种介质的对称轴是水平的，现在取对称轴主轴为X轴，这时介质的物性矩阵 中，

4、是3阶零矩阵，另外，,（3）VTI介质+HTI介质组成单斜各向异性介质 1）VTI介质与XOZ倾斜HTI介质组成的单斜介质,弹性矩阵中有13个弹性参数。初态时，若介质A即HTI介质对称主轴是水平X轴，其各向同性面是YOZ，该对称面含有VTI介质的对称主轴；因此，这时的介质具有一个垂向的XOZ对称面。当XOZ面绕Y轴旋转角度时，则介质的对称主轴已不再是水平，且其对称面内不再包含VTI介质的对称主轴，从而使介质表现为单斜各向异性。,2）各向同性与微褶皱表面的垂直裂隙组合 这种模型是Schoenberg等（1989，1995）根据线性滑动理论建立的。裂隙的法向与X轴平行，具有一个垂直对称面。微褶皱的

5、含义是指裂隙表面在小于地震波波长尺度上观测是不规则的。这种不规则性可以视为一个锯齿状的剖面，裂隙表面的小褶皱之间相互有一定的错动，从侧面看呈锯齿状。裂隙表面的不规则性是相对于垂向而言的，这种不规则性使得裂隙受到外力的作用时，在法向和切向上位移之间产生相互的联系，从而使介质表现为单斜各向异性。,弹性矩阵中有12弹性参数。,介质的物性矩阵如下：,3）各向同性与两组非正交的垂直裂隙组合 在各向同性介质背景下，存在两组非正交的HTI介质（据Grechka等，2000）。弹性矩阵中有13个弹性参数，矩阵如下,六方介质三种典型的物性矩阵,六方各向异性介质本构方程和物性矩阵横向各向同性介质的物性矩阵扩张型裂

6、隙各向异性介质的物性矩阵薄互层各向异性介质的物性矩阵,（1）流体型EDA介质,这种模型介质裂隙内充填物质是流体。其物性矩阵 中，子矩阵,2.1.3 扩张型裂隙各向异性介质的物性矩阵,其中，和 是均匀弹性各向同性介质中纵波和横波的传播相速度。是裂隙参数；是裂隙密度，它表示半径为b的球体内所含裂隙的数目；是单位球体内裂隙的数目。,是3阶零矩阵,（2）气体型EDA介质,是3阶零矩阵,六方介质三种典型的物性矩阵,六方各向异性介质本构方程和物性矩阵横向各向同性介质的物性矩阵扩张型裂隙各向异性介质的物性矩阵薄互层各向异性介质的物性矩阵,薄互层各向异性介质即PTL（Periodic Thin-Layers）

7、介质，这种介质模型首先由帕斯特玛（Postman，1955）提出，怀特（White，1983）指出其物性矩阵是,是3阶零矩阵,2.1.4 薄互层各向异性介质的物性矩阵,2.2 应力应变和物性矩阵的坐标变换,三阶和六阶变换矩阵邦德矩阵及其变换六阶矩阵变换的结构和构造方法应变矩阵的六阶变换矩阵物性矩阵的坐标变换,令变换前老坐标系下的应力矩阵是 或，变换后新坐标系下的应力矩阵是 或。按照第一篇应力和应变矩阵单双角标关系，有,总是存在下面的3阶正交变换矩阵,2.2.1 三阶和六阶变换矩阵邦德矩阵及其变换,使得：,展开,于是有,同理,然后把应力分量按单角标形式整理成,或,等价,应力应变和物性矩阵的坐标变

8、换,三阶和六阶变换矩阵邦德矩阵及其变换六阶矩阵变换的结构和构造方法应变矩阵的六阶变换矩阵物性矩阵的坐标变换,把 对应的6阶M变换矩阵划分成 四个3阶子块矩阵。,2.2.2 六阶变换矩阵的结构和构造方法,（1）六阶变换矩阵M的结构和构造方法,子矩阵的构造,子矩阵的构造,子矩阵的构造,子矩阵的构造,是正交矩阵,（2）六阶变换矩阵 的结构和构造方法,采用构造M的方法可以得到6阶变换矩阵，它满足,把 式代入上式，得到,同理得到,于是得到关系,由 派生而来，方法和M的构造等同。,三阶和六阶变换矩阵邦德矩阵及其变换六阶矩阵变换的结构和构造方法应变矩阵的六阶变换矩阵物性矩阵的坐标变换,应力应变和物性矩阵的坐

9、标变换,科西方程 有两种表示形式，做应力矩阵变换时也有同样的问题,第一种情况：应力矩阵的非对角线元素表达式系数为“1/2”时，即,应变的变换满足,只要把矩阵 中元素 的角标对换就可以得到,2.2.3 应变矩阵的六阶变换矩阵,第二种情况：应变矩阵的非对角线元素表达式中系数为“1”时，即,应变的变换满足,应力应变和物性矩阵的坐标变换,三阶和六阶变换矩阵邦德矩阵及其变换六阶矩阵变换的结构和构造方法应变矩阵的六阶变换矩阵物性矩阵的坐标变换,2.2.4 物性矩阵的坐标变换,2.3 方位六方各向异性介质,直角坐标系中绕三个坐标轴的旋转变换直角坐标系任意方向的旋转变换,绕Z轴的旋转变换,直角坐标系下物性矩阵

10、为D，有3阶变换矩阵是：,绕Y轴的旋转变换,直角坐标系下物性矩阵为D，有3阶变换矩阵是：,绕X轴的旋转变换,直角坐标系下物性矩阵为D，有3阶变换矩阵是：,一般情况下，变换前后物性矩阵中的元素表达形式不同了，原先为零的矩阵元素，有一些元素现在不是零了。对各向同性介质而言，变换前后的物性矩阵是相等的；由此可以知道，研究各向同性介质中地震波传播规律时，与坐标系的选择没有关系。对各向异性介质，变换前后的物性矩阵是不相同的；研究各向异性介质中地震波传播规律时，与坐标系的选择有密切关系；波的传播特点与波传播方向有关系。,方位六方各向异性介质,直角坐标系中绕三个坐标轴的旋转变换直角坐标系任意方向的旋转变换,

11、由空间解析几何知识可以知道，凡是直角坐标系的任意转动都能通过绕不同坐标轴的相继旋转来完成，下面说明其中的两种实现方式。,第一种方式。第1步，绕原坐标系的 Z 轴逆时针旋转 角度。第2步，绕第一步变换后得到的 Y 轴逆时针旋转 角。第3步，再绕第二步后得到的 X 轴逆时针旋转 角度。三次坐标系的转动依次使用了、三个角度。坐标系转动需要 3 个 6 阶变换矩阵。,第二种方式。第1步，绕原坐标系的 Z 轴逆时针旋转 角度。第2步，绕第一步变换后得到的 Y 轴逆时针旋转 角度。第3步，再绕第二步后得到的 Z 轴逆时针旋转 角度。坐标系的所有转动能用、三个角度来表示。由此坐标系转动需要两个 6 阶变换矩阵。,