地下介质三维模型中界面的描述方法摘要.PDFVIP

第 1 7 卷 　第 3 期 　　　　　　　　地 　球 　物 　理 　学 　进 　展 　　　　　　　　　　Vol. 17 　No. 3 ( ) 2002 年 9 月 445～450 　　　 　　　　PROGRESS 　IN 　GEOPHYSICS 　　　　　　　 　　　　　Sept. 2002 地下介质三维模型中界面的描述方法 李迎春 　徐果明 　徐 　涛 ( 中国科学技术大学地球和空间科学系,合肥 230026) [摘 　要] 　介质模型可分为离散模型和连续模型. 离散模型 ,使用规则网格或不规则网格将空间 离散化 ,网格内介质参数不变. 对连续模型 ,可用解析函数描述 ,或者由插值方法得到介质参数的 空间分布. 连续模型的描述 ,要素之一是界面的描述方法. 本文借助计算机图形学和计算机辅助几 何设计中发展起来的曲线和曲面技术 ,使用参数形式的 Coons 曲面、Bezier 曲面、B 样条曲面和三角 形面片 ,描述三维地质构造. 结合地震勘探采集系统设计的实际需要 ,给出一种三角形面片的模型 描述方案. [ 关键词] 　三维模型 ;Coons 曲面;Bezier 曲面;B 样条曲面 ;三角形面片 ( ) [ 中图分类号] 　P315 　　　　[文献标识码] 　A 　　　　[文章编号] 2002 0 　引 　言 [1 ] 地震模拟技术 主要有积分变换法、波动方程法和射线追踪法. 积分变换法使用的模型 是横向均匀介质模型或者球对称介质模型. 人们已经对一维介质模型做了比较充分的研究. 考虑体波、面波和自由振荡 ,一维非均匀地球模型的理论地震图的计算方法已经建立[2 ] ,例 如利用 Cagniard 方法和 Fourier 变换方法研究体波和面波 ,用 Cagniard 方法和 Fuchs 方法研究 有很多均匀平面层组成的介质. 波动方程法是将波动方程离散化 ,利用数值方法求解方程 , 主要有有限差分方法和有限元方法. 波动方程法使用的是网格模型 ,计算量主要取决于网格 的多少. 考虑一个三维地质模型[3 ] ,大小为 3000 m ×3000 m ×3000 m ,网格间距为 10 m ,网格 7 节点数达到 2. 7 ×10 . 要使计算结果精度高 ,则要求更小的网格和更小的时间步长 ,此方法 的使用主要受计算机内存和速度的影响. 比较而言,射线追踪法的计算速度较快 ,波动方程 法能得到方程的完全解 ,如能得到更精确的振幅值和波形. 波动方程法和射线追踪法可适用 于任意复杂的介质. 射线追踪法是建立在射线级数展开基础上的高频近似方法 ,要求介质性 质在一个波长范围内变化不大. 上面几种方法所用的介质模型可分为离散模型和连续模型. 对离散模型 ,使用规则网格 或不规则网格将空间离散化 ,网格内介质参数不变. 连续模型的物性参数和它的一阶、二阶 导数在模型单元内连续变化. 连续模型的界面为物性参数或者它的一阶、二阶导数的不连续 面. 构造连续模型的一种方法是通过解析函数得到物性参数的分布 ,或者由空间的一些离散 点的数值 ,利用线性或者样条函数插值得到模型. Aki[4 ] ( 1977) 等利用地震台阵记录的地震 [5 ] ( ) 波到时反演地球内部的三维速度结构 ,用的是分区均匀的离散速度模型. 刘福田 1989 等 在球面坐标系中 ,用三维空间的六面体网格节点的速度 ,线性插值得到速度的空间函数 ,速 度界面取为球面 ,用此模型进行了层析成像方面的研究. 射线追踪[6 —8 ] 有很多种算法 ,所用 [ 收稿日期] ; 　[ 修回日期] ( ) [ 基金来源] 　国家自然科学基金项目 4997