数学建模辽宁工程技术大学作业.docVIP

辽宁工程技术大学 数 学 建 模 课 程 成 绩 评 定 表 学 期 08-09学年1学期 姓 名 姜德海王智吴洪博 专 业 交通土建 班 级 土建07-8班 课程名称 数学建模 论文题目 储量计算模型分析 评 定 标 准 评定指标 分值 得分 知识创新性 20 理论正确性 20 内容难易性 15 结合实际性 10 知识掌握程度 15 书写规范性 10 工作量 10 总成绩 100 评语: 任课教师 时 间 09年 9 月 19 日 备 注 储量计算模型分析 1 摘要： 矿产储量是指具有一定地质研究与控制程度的已查明的矿产资源。本文从矿产资源开发利用和保护以及储量管理的现状分析入手，采用的储量计算理论，应用广泛的——MATLAB，结合的储量问题加以研究，技术进行系统开发研究。根据地质勘查工作获得的矿床资料矿体形态和矿石质量分布特征通过计算以确定的量。通过采集与处理测量数据及其他相关数据，进行矿山地质模型建模，实现矿山矿产储量的图形生成，任意个区域的储量计算以及矿体在该区域的倾角、报表的生成，动态了解资源的动用情况。 这对于指导矿山生产，提高资源的回采率和提高经济效益等有重要的意义。对于矿政主管部门监督税收、审批办矿，合理配置资源，强化储量的动态管理也都有重要的现实意义。该矿物的储量。 表2-1 地区勘探数据表 编号 1 2 3 4 5 6 7 8 9 10 x坐标(KM) 1 1 1 1 1 2 2 2 2 2 y坐标(KM) 1 2 3 4 5 1 2 3 4 5 矿物体厚度H(M) 13.72 25.80 8.47 25.27 22.32 15.47 21.33 14.49 24.83 26.19 编号 11 12 13 14 15 16 17 18 19 20 x坐标(KM) 3 3 3 3 3 4 4 4 4 4 y坐标(KM) 1 2 3 4 5 1 2 3 4 5 矿物体厚度H(M) 23.28 26.48 29.14 12.04 14.58 19.95 23.73 15.35 18.01 16.29 3 模型假设 3.1 假设该地区短时间内地理环境、地形地貌保持不变。 3.2 假设该地区的矿产资源不受气候条件的影响。 3.3 假设该地区近期不受自然灾害及地质灾害的影响。 3.4 假设调查期间技术人员和建模经费得到保证。 4 问题分析与建立模型 由表1可以知矿物体的厚度H是坐标x和坐标y的二元函数H=H(x,y)，根据二重积 分的知识有所求该矿物的储量就是二重积分的值。由于函数H(x,y)没有给出具体表达式，而只给出一些离散点上的函数值，下面采用数值积分的方法来近似计算该矿物的储量。 在计算方法中有很多计算定积分的公式，根据本题特点，这里采用计算定积分的复化 梯形公式来计算二重积分。为得到计算二重积分的复化梯形公式，注意到计算定积分的复 化梯形公式为： 这里节点是等距节点，是相邻节点间的距离。 对二重积分令，则有 于是二重积分计算可以由g(x)的复化梯形公式： 计算，而其中的也是用复化梯形公式计算的，于是可以利用逐次复化梯形求积公式近似计算出来。通过具体推 导，可以得到计算二重积分的复化梯形公式： （4-1） 其中是下面矩阵A的元素： 5 模型求解 表2-1的数据是网格化的（数据点位置见图5-1） 图5-1数据点位置示意图 由上图可得：h x= h =1,n=3，m=4， 相应的矩阵： 利用公式(4-1) 编程计算有： In[1]:= da={{13.72/4,15.47/2,23.28/2,19.95/4},{25.8/2, 21.33，26.48,23.73/2}, {8.47/2, 14.49, 29.14,15.35/2}, {25.27/2, 24.83, 12.04,18.01/2}, {22.32/4,26.19/2,14.58/2,16.29/4}}; daa=4da; Sum[daa[[i,j]],{i,1,5},{j,1,4}]/4 Out[1]= 244.455 单位换算即得出该矿物的储量约为0.244455(10 9(立方米)。 6 模型检验 本次建模实验的目的对储量建立科学的价值分级和评估模型从而科学地评价其储量资产的价值采用综合评价方法研究了储量资源在勘探开发中由于储量资源受地质条件、田所处地理位置、交通运输条件和社会经济条件的不同影响使各储量资源形成级差成本和级差收益传统法与法计算结果相对误差大小不能说明计算结果的正确性和真实性，就是说，两种方法储量计算值相近