武汉科技大学研究生数值分析论文.doc

研 究 生 课 程 论 文 封 面 课程名称： 论文题目：插值法在地质曲面拟合中的应用 学生班级； 机械工程2班 学生姓名： 任课教师： 学位类别： 评分标准及分值 选题与参阅资料 （分值） 论文内容 （分值） 论文表述 （分值） 创新性 （分值） 评分 论文评语： 总 评 分 评阅教师: 评阅时间 年 月 日 注：此表为每个学生的论文封面，请任课教师填写分项分值 插值法在地质曲面拟合中的应用 摘要：介绍了一种适合于非网格化样品点研究区域的插值方法——多重二次曲面函数插值法。这种方法运用二次基本曲面函数作为插值的基函数，数学方法简单，易于编程实现。该方法对原始数据的边界和分布情况要求不高，适合于处理分散的和非网格化的数据. 将该方法引入煤层地质曲面逼近, 效果良好。 关键词：.插值；多重二次曲面；地质曲面 1. 1曲面函数 基本二次曲面的种类较多, 在这里, 我们选用符合右手循环规则的圆锥作为二次曲面来讨论, 圆锥的顶点可以处于网格化或非网格化样品空间中的任意位置。一个圆锥的方程可以由下式表示: (1) 式( 1) 中的 z 表示在 xy 平面上的圆锥顶点的标高,tanα是圆锥曲面的斜率。类似的, 多重二次曲面系统可以由已知样品点的标高( 即一系列圆锥的顶点) 建立, 根据已知点, 可以得到多重二次曲面函数的基本方程组: (i = 1, 2, ?, n) (2) C′ = [ c 1 c 2 ? c n ] Z′ = [ z 1 z 2 ? z n ] 则可将式( 2) 展开成矩阵形式 AC= Z (3) 其中 cj 是代表圆锥斜率系数的列矢量, 圆锥的顶点表示定位在( x j, yj ) 坐标上, 并已知其标高( 或属性值) zj 的样品位置。 1.2基本原理 任何不规则连续曲面的构成都包含了由少量控制点插值而来的大量中间点。 典型的控制点是指已知( x , y) 坐标, 并且已经观测到了我们感兴趣的属性值的所在位置。这个属性值可以是基岩所在位置的标高, 也可以是不同地方的重力加速度( g)值, 或者是河流、 海洋及气流剖面在不同地方的流速, 还可以是基岩的深度或地下水水位的深度。事实上, 多重二次曲面函数法对任意包含xyz 信息的连续曲面都适用 [ 4] 。由样品点获得的 xyz 数据构了在研究区域内任何位置插入新值的基础.对于研究区域内的任何位置( x p , y p ) ( 可以是控制点, 也可以为任何位置) , 其标高值 z p 可由下式插值而得: （4） 令 则有 （5） 由于我们可以在研究区域内任意给定插值点, 所以在式( 2) 与( 4) 中唯一的未知数是列矢量 c j , 它可以通过求解齐次线性方程组的方法求出。 一旦求出 c j之后, 就可以通过式( 4) 或式( 5) 估值得到用于拟合不规则曲面的 z p . 显然, 不规则连续曲面的逼近程度取决于曲面单元面积内样品点( 即圆锥顶点) 的个数. 单元面积内样品点越多, 越逼近现实的连续曲面。 我们还可以运用偏微分的方法计算 z 在 x, y轴方向的偏导数, 从而进一步确定曲面斜率变化的情况. 对于( 4) 式用连锁法求解 z 对 x 的偏导数见下式: （6） （7） 插值的实现 通过上面的分析, 我们可以把利用多重二次曲面函数插值绘制连续不规则曲面的方法概括为以下几步: 1) 获取包含 x yz 信息的样品. 样品位置的设计可以是网格化的, 也可以是非网格化的. 在实际应用中可以通过钻孔数据得到初始值. 2) 利用样品数据解出式( 2) 或式( 3) 中的列矢量 c j . 对于 A X= B 型的方程组, 求解方法很多, 针对式( 3) 的特点, 可以用回归迭代奇异值分解法求出 c j . 3) 对方程组式( 4) 或式( 5)