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

研 究 生 课 程 论 文 封 面 课程名称： 数值分析 论文题目：插值法在地质曲面拟合中的应用 学生班级； 机械工程2班 学生姓名： wtfgw 任课教师： ryhersgy 学位类别： 学位课（2学分，32学时） 评分标准及分值 选题与参阅资料 （分值 ） 论文内容 （分值 ） 论文表述 （分值 ） 创新性 （分值 ） 评分 论文评语： 总 评 分 评阅教师: 评阅时间 年 月 日 注：此表为每个学生的论文封面，请任课教师填写分项分值 插值法在地质曲面拟合中的应用 摘要：介绍了一种适合于非网格化样品点研究区域的插值方法——多重二次曲面函数插值法。这种方法运用二次基本曲面函数作为插值的基函数，数学方法简单，易于编程实现。该方法对原始数据的边界和分布情况要求不高，适合于处理分散的和非网格化的数据. 将该方法引入煤层地质曲面逼近, 效果良好。 关键词：.插值；多重二次曲面；地质曲面 目前国内常用的插值方法有双三次样条函数插值法、 曲面样条函数插值法和曲面磨光插值法等。这几种方法各有优缺点: 双三次样条函数法虽然从理论上可以作为插值的方法, 但是由于该方法必须满足已知点的平面坐标按规则点阵排列, 同时需知四条边界的一阶导数和四个角点的二阶导数,地质问题的实际情况往往无法满足。 中国矿业大学余志伟于1987年首次将曲面样条函数引入数学地中, 用来计算地质曲面, 取得了良好的效果。 这种方法的优点是: 已知点平面坐标不必按规则点阵排列, 采用自然边界条件而不需边界导数的信息, 就可以得到任意阶可微的光滑曲面 [ 1, 2] . 但是这种方法数学原理比较复杂, 已知样品点多时, 数值实现也比较困难. 曲面磨光函数插值法在一定程度上符合地质曲面的实际情况, 它能够在曲面磨光图上显示断层、 挠曲、 向斜、 背斜、 构造盆地和穹隆等各种地质构造信息[ 3] 。但是这种方法必须预先对已知点作网格化处理, 增加了工作量。 本文所描述的多重二次曲面函数法( Multiquadric Function) 是一种新的地质曲面插值方法, 这种运用二次曲面基本函数插值的方法最早是哈迪在1971年提出的, 现在已经应用到了地形图绘制、 流体力学计算、 在河流剖面上绘制流速图、雨量区域化分布图绘制、 矿产资源勘探等多个研究领域, 并取得了良好的效果[ 4～8] 。但是, 这种方法在国内的应用却较少, 本文将该方法引入煤层地质曲面逼近, 取得了良好的效果。 数学原理 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 ) ( 可以是控制点, 也可以为任何位置) , 其标