資源描述：
《長(zhǎng)安大學(xué)重磁報(bào)告二》由會(huì)員上傳分享，免費(fèi)在線(xiàn)閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。
1、重磁實(shí)驗(yàn)報(bào)告（二）位場(chǎng)邊緣識(shí)別姓名:學(xué)號(hào):專(zhuān)業(yè)：勘查技術(shù)與工程指導(dǎo)教師：魯寶亮、王萬(wàn)銀完成日期：2013年12月10日基本原理:1、重力異常計(jì)算公式.2、磁力異常計(jì)算公式四、五、目錄輸入/輸出數(shù)據(jù)格式設(shè)計(jì):1、輸入數(shù)據(jù)文件名的格式設(shè)計(jì)2、輸出數(shù)據(jù)文件名的格式設(shè)計(jì)3、重要變量的名稱(chēng)總體設(shè)計(jì):測(cè)試結(jié)果:11233345結(jié)論及建議:錯(cuò)誤!未定義書(shū)簽。錯(cuò)誤!未定義書(shū)簽。附錄：源程序代碼Efff<-z）物^如今叫J松一彳埔一力2+（卜z）2f2式中：V為地質(zhì)體的剩余質(zhì)量對(duì)A點(diǎn)的單位質(zhì)量產(chǎn)生的引力位；V為地質(zhì)體的體積。我們還
2、可以推導(dǎo)出計(jì)算重力異常垂向梯度或重力垂向梯度異常的基本公式為：=v=r>Tfff2((—z)2—(f—x)2—(7—y)2dzJJJ[(^-%)24-(Z7-J)24-(^-z)2]5/2d^d/jdC計(jì)算重力異常水平梯度或重力水平梯度異常的基本公式為:(fz)(f—X)[0-又)2+("-y)2+((-z)2]5/2(^-z)(Z7-y)d^dTjdC[(^-x)2+(7-y)2+(^-z)2]5/2d^drjd^一階導(dǎo)數(shù)類(lèi)邊緣識(shí)別計(jì)算公式：垂向?qū)?shù)：VDR（x,z）=3Ag（；V°?，Z）dz總水平導(dǎo)數(shù)：TH
3、DR(x,y,z)=.^S^y.z)}2+z3A^(x,y,z)dx解析信號(hào)振幅：ASM=JTHDR2+VDR2(2)(x,y，z)/at)/J“0，)’0，z0)Lz//zc直立長(zhǎng)方體示意圖均勻磁化直立長(zhǎng)方體磁場(chǎng)各分量表達(dá)式:H(?(x,>',z)^{-Mvarctan(^-x)(r]-y)(^-x)2+r(^-z)+(^-z)2aIbIc(1)Mvlnk+(r_z)]+M2lnk+(7_y)l}廣+!Z°12-V°~iZ°~2H“y(x，y,z)=^{Mx[r+(^-z)]-Mvarctan(7—y)2
4、+r(f—z)+(f—z)2aabX~2K-三
5、+MJn[r+(^-x)]}?22—(2)Z“(x，y，z)=^{Mxln
6、r+(Z7—y)l+Mvln[r+(^-x)
7、4兀Mzarctan(￡zM}abx.+一a2?v々abv—一。2^2e?+2e"2(3)把(1)、(2)、(3)式代入磁異常AT計(jì)算公式屮得：AT(%,>z)=^M{k^[r+^-x)]+女2?Ink+(77-y)]+&?ln[r+((-z)J+k4^arctan+眾6?arctan(^-x)2+r(^-z)+(^-z)2>"I7I+k
8、5?arctan(7-y)2+r(^-z)+(^-z)2ab?、勹y^ClbX''-2門(mén)1廠(chǎng)((—z)其中r=[(
9、-x)2+(7-y)2+((-2)2]1/2,么77,(為異常體的坐標(biāo)，Lo、Mo、為地磁場(chǎng)的方向余弦；貸為總磁化強(qiáng)度，M。、My、Ma為其在x、Y、Z坐標(biāo)方向的分量。設(shè)M為H的模，I為磁化強(qiáng)度傾角，D為磁化強(qiáng)度磁偏角，且滿(mǎn)足關(guān)系式Ma==McosIcosDMy=McosIsinD二、輸入/輸出數(shù)據(jù)格式設(shè)計(jì):1、輸入數(shù)據(jù)文件名格式設(shè)計(jì)：（1）輸入數(shù)據(jù)及文件名對(duì)應(yīng)如下：inputgravity_
10、fi1ename:anomaly,grdsquare_num:0i叩ut_gravity_filename對(duì)應(yīng)anomaly,grd文件，存放重力異常數(shù)據(jù)；square_num賦值力0表不?將區(qū)域擴(kuò)力256*256的方形區(qū)域。（2）觀(guān)測(cè)面上的重力異常存放在“anomaly,grd”中，坐標(biāo)單位為m，格式如下：2011000.0000001000.000000fieldmaxDSAA201-1000.000000-1000.000000fieldmin可知，網(wǎng)格化數(shù)據(jù)按GRD格式保存，有201條線(xiàn)、每條線(xiàn)201個(gè)點(diǎn)
11、，點(diǎn)（線(xiàn)）坐標(biāo)從一1000m?1000m，點(diǎn)（線(xiàn)）距為10m，fieldmin,fieldmax分別表示網(wǎng)格化數(shù)據(jù)的最小值與最大值。2、輸出數(shù)據(jù)文件名格式設(shè)計(jì)：(1)輸出數(shù)據(jù)及文件名對(duì)應(yīng)如下：verticalderivativefilename:verticalderivative.grdtotal_horizontal_derivative_filename:total_horizontal_derivative.grdsignalamplitudefi1ename:signalamplitude.grd(2)
12、對(duì)給定平面上的重力異常，計(jì)算其垂向一階導(dǎo)數(shù)存放在vertical_derivative.grd文件中；總水平導(dǎo)數(shù)存放在total_horizontal_derivative.grd文件中；解析信號(hào)振幅存放在signal_amplitude.grd文件中。坐標(biāo)單位為m，格式如下：DSAA2011000.0000001000.000000fieldmax201-1000.00
