資源描述：
《含沙量的測(cè)定建模分析》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫(kù)。
1、含沙量的測(cè)定一、問(wèn)題的重述與分析1、原題某地區(qū)有優(yōu)質(zhì)細(xì)沙埋在地下，某公司欲在該地開(kāi)采鉆探資料如下表：在每個(gè)格點(diǎn)上有3個(gè)數(shù)字，分別代表覆蓋表層的標(biāo)高，沙層頂部的標(biāo)高，沙層底部的標(biāo)高，陰影部分表示沼澤地地帶，沒(méi)有鉆探資料，試估計(jì)整個(gè)矩形區(qū)域內(nèi)的含沙量。ABCDEFGH22.4****22.523.023.2020.0****18.417.818.05.8****0.50.40.423.023.123.223.423.524.024.024.0119.920.020.019.819.920.019.819.66.03.21.
2、61.01.11.00.80.923.123.223.423.423.524.224.124.1219.819.719.420.020.120.320.320.52.21.40.60.50.3-0.2-0.10.02、分析我們假設(shè)沙層是連續(xù)的。首先，通過(guò)插值法，得到沼澤地帶的大致沙層數(shù)據(jù)。然后，把整個(gè)沙層切割成多個(gè)小方塊，體積求和即為最終結(jié)果。二、算法和模型1、插值法計(jì)算沙層高度分別對(duì)沼澤地帶覆蓋表層的標(biāo)高，沙層頂部的標(biāo)高，沙層底部的標(biāo)高進(jìn)行插值估計(jì)，得到以下結(jié)果：BCDE22.4222.4422.4622.4819.
3、680019.360019.040018.72004.74003.68002.62001.5600進(jìn)一步，通過(guò)excel的減法運(yùn)算可以得到沙層高度：12345678114.214.9415.6816.4217.1617.917.417.6213.916.818.418.818.8191918.7317.618.318.819.519.820.520.420.5我們令步長(zhǎng)為0.1，通過(guò)曲面插值得到更精確的沙層高度分布，如圖：2、分割求沙層體積原本數(shù)據(jù)點(diǎn)有701*201個(gè)，現(xiàn)取700*200個(gè)小長(zhǎng)方塊，每個(gè)長(zhǎng)方塊底為1*1
4、=1平方米，高為該數(shù)據(jù)點(diǎn)與(x+0.01,y)、(x,y+0.01)、(x+0.01,y+0.01)高度的平均值，疊加即最后結(jié)果。求出最后結(jié)果：251.7482。又原題中一個(gè)采樣格長(zhǎng)度為50m，故結(jié)果乘2500為最終結(jié)果：629370.5立方米。三、附錄>>x1=[1678];>>y1=[22.422.52323.2];>>x2=[2345];>>y2=interp1(x1,y1,x2)y2=22.420022.440022.460022.4800>>y1=[2018.417.818];>>y2=interp1(x1,
5、y1,x2)y2=19.680019.360019.040018.7200>>y1=[5.80.50.40.4];>>y2=interp1(x1,y1,x2)y2=4.74003.68002.62001.5600用EXCEL求出沙層高度z導(dǎo)入matlab>>x=[111222333444555666777888];>>y=[123123123123123123123123];>>z=[14.213.917.614.9416.818.315.6818.418.816.4218.819.517.1618.819.817.9
6、1920.517.41920.417.618.720.5];>>cx=1:0.01:8;>>cy=1:0.01:3;>>cz=griddata(x,y,z,cx,cy','cubic');whilei<3;whilej<8;v1=0.01*0.01*(cz(i,j)+cz(i,j+0.01)+cz(i+0.01,j)+cz(i+0.01,j+0.01))/4;v=v1+v;j=j+0.01;endi=i+0.01;j=1;end>>vv=251.7482
