沙子的开采-数学建模

《沙子的开采-数学建模》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、沙子的开采问题摘要本文主要研究的是估算某一地区的含沙量，并做出开采沙子的收益分析。针对问题一，我们根据题中表格中已有的数据，利用线性插值的方法对相邻两点作线性外推，估算出了缺少的数据。然后通过matlab画出地表、沙层顶部以及沙层低部的等高线。为了更精确的计算出表土层与沙层的体积，用matlab对数据进行线性插值，将50米50米的方格精确到1米1米的小方块。再利用体积公式近似计算出表土层与沙层的体积，进而得到该矿含沙量超过400万吨，表土体积小于沙子体积19%，可以对这项工程进行投资。针对问题二，构造单位面积

2、的利润函数，利用matlab软件画出利润等值线。由于贴现的缘故，施工应按照利润由大到小的顺序进行。所以施工顺序是按照等值线的分布进行的。在确定每个阶段的施工时间后，分别计算利润，进而求得总利润。关键词：沙子的开采线性插值matlab贴现率11一.问题简述沙子作为最重要建筑材料之一，一直有着稳定的市场。沙子在形成之后上面覆盖着表土，由于表土层结构松散，沙子只能露天开采即先将表土移走再挖沙子。现有一开采地，四周有几块沼泽地。蓝色隧道公司为生产建筑材料委托一支地质队调查该开采地，在开采地上画的网格点处用空心管垂直钻

3、入地下．利用进入空心管的厚度的开采地的面积可得出表土和沙子的体积，其中每个网格点中最上面一项是地表高度，中间一项是沙层顶部的高度，最下面一项是沙层底部的高度．由于地表存在沼泽地，部分网格点的数据没有采集，即所给表中的数据空白处．现有该公司提供的相关材料，需要解决以下工作：（1）通过计算表土和沙子体积决定是否该对这项工程进行投资，并画出地表、沙层顶部以及沙层低部的等高线。（2）如果可以投资，根据已给数据和条件设计开采方案，做出收益分析作并求出相应的利润。二.问题分析问题一：开采地四周的沼泽地也存在沙子，所以要补

4、充完整表格中的数据，即分别把表土高度、沙层顶部高度、沙层底部高度的数据估算出来。根据网格中已有的相关数据，取前后或上下两个连着的数据，利用线性插值法进行数据处理，然后用Matlab编程运算分别得出表土高度、沙层顶部高度、沙层底部高度的等高线图。题中给出的地表高度、沙层顶部和沙层底部高度的原始数据，相临的点相距50米。为了更精确地计算表层体积，对沙子的储量作精确的估算，用matlab软件对数据进行插值，使得相临的点彼此相距1米，每个网格是11的小方块，经过插值后的数据可更精确地反映沙层厚度分布情形。我们将所有的

5、沙土汇集在单位网格，网格点补充完整后就根据网格点中的数据分别算出每个小区域的表土厚度（表土高度-沙层顶部高度）和沙层厚度（沙层顶部高度-沙层底部高度），利用体积公式累加求出沙层、土层的总体积。11问题二：为解决开采沙子的收益分析，首先要构造利润函数。在不考虑成本中每月10万元的地区管理费的情况下，单位面积利润函数(单位:元/米2)只存在沙和土的密度与厚度这俩个变量，沙的销售单价，沙和土的开挖费和运输、处理费都为固定值。根据所得利润函数利用matlab软件求得函数的等值线。由于每个月交纳的10万元地区管理费,再

6、考虑到每月最多开挖量10万吨的限制,只作出了该函数大于等于0的范围,我们的施工设计局限于该范围。由于10万吨的限制将施工区域分成几个部分，按阶段施工，分别求出各阶段的利润，求和得出总利润。三.模型假设1.假设题目中所给数据都是准确的。2.开采地沙层与表土层都是连续的，近几年不会发生地质以及自然灾害造成地质构造的变化。3.月贴现只针对利润，并不包括其他费用。4.沙子的单价，开采费用以及其他投入价格不会发生改变。5.每月利润当月结算，所得收入在月底统一贴现。四.符号说明符号符号说明单位m1m2p表土的质量表土的密

7、度表土的体积沙子的质量沙子的密度沙子的体积单位面积利润函数开采单位面积沙土所获利润吨吨/立方米立方米吨吨/立方米立方米磅/平方米磅11五.模型的建立与求解1.1问题一：要解决表土和沙子的体积就要填充网格中所空白数据，首先对网格中的所给出的数据进行分析，利用二维线性插值法较精确地估算出表格中因沼泽地原因无法测得的数据。二维线性插值法方法如下：首先在X方向进行线性插值，得到：然后在Y方向进行线性插值，得到：这样就得到所要的结果：即得到如下数据表格：11ABCDEFGHIJKLMNO01234567822.322.

8、0421.922.022.423.223.423.222.822.523.023.223.222.621.919.119.019.319.520.020.620.419.919.3518.417.818.019.019.018.816.414.212.18.75.84.12.11.10.50.40.40.60.20.23.622.422.422.522.623.023.123.223.423.524