重力勘探報告finish

1、地球物理與空間信息學(xué)院2010級本科生重力勘探大作業(yè)班 級：作 者:學(xué) 號:指導(dǎo)老師：、原始數(shù)據(jù)繪圖二、布格校正首先在excel里對數(shù)據(jù)進(jìn)行布格校正，本次處理忽略了地形校正，從而減少了工作量。先 后進(jìn)行了自由空間校正、中間層校正和緯度校正，最后將每個點的校正值與相對重力值g相加得到了各點的布格重力異常值應(yīng)用公式為緯度校正：合 g 如£.814sin2 那 X(mGal , X:Km )中間層校正：S gB= 2 n Gp =0。0419 ch( h:m ,cr :t/m3 mGal)自由空間校正：gf=0.3086h(h : m,mGal)布格重力異常：gB = g 測 + gf+

2、S g申 S g $三、畫出布格異常等值線圖然后將獲得的布格重力異常數(shù)據(jù)導(dǎo)入surfer軟件中，直接可以生成平面等值線圖如下：15.554.5k?,53£.51.515001DCOH0.5-0.5:-3 5b.s1-5500-o-p65Cd三維圖如下從該圖中可以看出區(qū)域場大致為平行的直線，初步判斷為單斜異常， 正好與老師所給的背斜翼部的背景吻合。平面等值線向異常降低的方向扭曲，說明地下為正異常，而且在中間位置形成了圈閉，證明地下異常體的水平梯度大于單斜異常的水平梯度。這時根據(jù)形狀初步四、分離區(qū)域異常和局部異常，對局部異常進(jìn)行處理1,運(yùn)用趨勢分析法選用n階多項式擬合區(qū)域異常，之后求解該

3、方程組用以描述多項式函數(shù)從而得到區(qū)域背景場 的特征擬合成的一次方程：各 g (x,y)=0.00200226x 0.0019722y繪圖如下用布格異常-區(qū)域場，消去區(qū)域場，得到局部異常，如下圖:O'J0.971.90I從上圖中已經(jīng)可以分辨出地下異常體的大致形狀為一個球體。接著應(yīng)用頻率域濾波法分離出區(qū)域異常。由于區(qū)域異常為低頻，則使用低通濾波器，進(jìn)行低通濾波后得到下圖在對數(shù)據(jù)進(jìn)行一階導(dǎo)數(shù)濾波后，可以明顯的看到一個極大值和一個極小值 ，但這對于判 斷地下異常體的形狀及位置雖有一定的幫助， 但是不能直觀的顯示出來。之后選用5點差分 法和9節(jié)點平均場法得到下面兩張圖.7.5.31 011.76

4、Z.2.2LIN1 .1 1B:M_ii:;_i三 MM其中異常導(dǎo)數(shù)轉(zhuǎn)換的作用是突出埋深較淺范圍較小的地質(zhì)體的重力異常特征，壓制深部的區(qū)域性重力異常根據(jù)以上的濾波，可以判斷出地下異常體就是一埋深不大的高密度球體。無、反演，計算埋深和質(zhì)量在數(shù)據(jù)中選擇兩條剖面，由局部重力異常等值線圖與橫縱向剖面圖綜合分析，該異常極可能是由球狀地質(zhì)異常體引起。剖面極大值對應(yīng)坐標(biāo)即為該球體的質(zhì)量中心在地面的投影位置。3.532.521.510.500-5從圖中可以看出，異常體位于最大重力異常的正下方，根據(jù)公式：ID1 X1/2X1/22：n 2廠1對局部重力異常的數(shù)據(jù)進(jìn)行提取，X=1250 時，有極大值=2.8740

5、74127mGal取n=2計算得到球體的埋深：D=1.305X1/2=0.6524 (X'1/2 -X1/2) 156.576mM= g m ax/G X D 2 10569139t六、總結(jié)本次大作業(yè)的總體過程是先計算布格重力異常， 畫出等值線圖， 再進(jìn)行分離區(qū)域異常和 局部異常，對局部異常進(jìn)行處理處理，分析出地下異常體的分布情況,最后選用合適的方法進(jìn)行反演計算 .做完之后再回顧這一過程就顯得十分清楚了，這其中許多部分都是根據(jù)經(jīng)驗 來判斷的，相信在以后的工作中我們會更加熟練的掌握并運(yùn)用這些技巧。本次實踐讓我收獲不小，原來聽老師講各種濾波處理以及反演運(yùn)算都只是公式上的運(yùn) 算，這次實際使用 surfer 時還沒有非常熟練，相信在以后的學(xué)習(xí)和工作中會熟練掌握 .本次大作業(yè)中存在的問題主要是軟件使用不熟練，