基于Delaunay三角剖分的重力正反演研究

来源 :成都理工大学 | 被引量 : 2次 | 上传用户:AsiaITt
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
地球重力场信息在地球构造研究、国土资源调查、矿产勘探以及国民经济等领域中发挥着重大作用。重力的正演问题与反演问题是解释重力资料的核心研究内容,其中前者主要研究不同形状大小、空间位置和物性参数等场源体或地质体所引起的重力异常特征及分布情况,其计算过程为正演计算;后者主要根据重力异常的分布来计算地质体或场源体的形状大小、空间位置、物性参数,其计算过程为反演计算。目标地质体的重力建模建立了地质体的剩余密度值与地面重力异常值的关系,既可利用重力正演公式计算地质体所引起的重力异常,也可以应用于求解重力反演问题。因此,重力建模是重力正反演问题的关键环节,建模方法的优劣则直接影响到重力正反演的计算精度。首先,研究Delaunay三角剖分的基本原理和方法,通过研究不同的三角剖分算法,分析其各自的优势,为构建更好的三角网模型提供依据;重点研究在二维平面内,根据所给定的节点位置以及数量,实现对平面内任意形状的地质模型进行有效Delaunay三角剖分,为后文地质体模型构建提供帮助。其次,在研究常规则地质体重力异常的基础上,重点研究了截面为任意形状的变密度二度体的重力问题,针对计算截面为任意形状的非均匀密度二度体的重力问题,本文提出基于Delaunay三角剖分的非规则重力建模方法,即利用二维Delaunay三角剖分方法,将其截面分割为若干三角形,即将二度体分解成为若干三角单元的组合,建立剩余密度值与三角形坐标的线性方程组,推导基于三角单元的重力计算公式;通过求取三角单元的重力异常,并将其累加,最后得出截面为任意形状的水平无限延伸的二度体的重力计算公式,以变密度的水平圆柱体为例,对比该模型的理论重力异常值,其结果相同,论证该方法的有效性。再次,分别对截面为圆形、矩形、梯形、工形、倾斜台阶的非均匀二度体模型进行Delaunay三角剖分,计算各自模型理论重力异常值,再将理论异常与反演质量计算异常进行拟合,达到最佳拟合时输出反演结果,并对其进行分析和解释,验证该方法的有效性及可行性。最后以西藏那曲地区依拉山一条实测重力剖面为例,用该方法进行了重力正演计算并反演,结合地质填图资料以及工区岩矿石标本的物性测试结果,对反演结果进行了相应分析,具有实用价值。
其他文献
新一轮课程改革尽管已经开展了一段时间,在各地也引起了广泛的重视,但在实际教学中,仍然有不少教师沉迷在应试教育的模式之中,他们在教学中仍然是把英语作为一种考试的工具,忽视英语的人文性,学生成为背单词、默写课文的工具,不能把英语学习与学生的日常生活相联系,使学生感到英语枯燥无味。然而,语言是在生活中产生的,语言即生活,生活原本是鲜活的、丰富多彩的、美丽的,所以英语也应该是美丽的,英语课堂也应该是美丽的