六方各向异性介质本构方程课件_第1页
六方各向异性介质本构方程课件_第2页
六方各向异性介质本构方程课件_第3页
六方各向异性介质本构方程课件_第4页
六方各向异性介质本构方程课件_第5页
已阅读5页,还剩43页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

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 Isatmpical),即对称轴是

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

4、, (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)介质,这种介质模型首先由帕斯特玛(P

7、ostman,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轴的旋转变换直角坐标系下物性矩阵为D,有3阶变换矩阵是:绕Y轴的旋转变换直角坐标系下物性矩阵为D,有3阶变换矩阵是:绕X轴的旋转变换直角坐标系下物性矩阵

10、为D,有3阶变换矩阵是: 一般情况下,变换前后物性矩阵中的元素表达形式不同了,原先为零的矩阵元素,有一些元素现在不是零了。对各向同性介质而言,变换前后的物性矩阵是相等的;由此可以知道,研究各向同性介质中地震波传播规律时,与坐标系的选择没有关系。对各向异性介质,变换前后的物性矩阵是不相同的;研究各向异性介质中地震波传播规律时,与坐标系的选择有密切关系;波的传播特点与波传播方向有关系。方位六方各向异性介质直角坐标系中绕三个坐标轴的旋转变换直角坐标系任意方向的旋转变换 由空间解析几何知识可以知道,凡是直角坐标系的任意转动都能通过绕不同坐标轴的相继旋转来完成,下面说明其中的两种实现方式。第一种方式。第1步,绕原坐标系的 Z 轴逆时针旋转 角度。第2步,绕第一步变换后得到的 Y 轴逆时针旋转 角。第3步,再绕第二步后得到的 X 轴逆时针旋转 角度。三次坐标

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论