【正文】 異常分離理論雖然它是一個(gè)反問(wèn)題 ,但是在確定區(qū)域領(lǐng)域和地區(qū)分布有著本質(zhì)的聯(lián)系 ,所以 ,作為解決方案的反問(wèn)題 ,分離的區(qū)域場(chǎng)也有多個(gè)解決嚴(yán)格數(shù)學(xué)的觀(guān)點(diǎn) ,不可能獲得的獨(dú)特的解決方案。在 1970年代蘇聯(lián)使用雙正交函數(shù)轉(zhuǎn)換公式來(lái)轉(zhuǎn)換多項(xiàng)式近似的引力場(chǎng) ,從而達(dá)到目的的分離區(qū)域場(chǎng)異常。 吉林大學(xué)學(xué)士學(xué)位論文 2 區(qū)域場(chǎng)概念方面的多項(xiàng)式擬 合曲面 ,最初做了很多探索正交多項(xiàng)式類(lèi)型但效果并不好 ,如重力區(qū)域領(lǐng)域多元多項(xiàng)式系數(shù)多項(xiàng)式擬合的地區(qū)區(qū)域異常效果 ,導(dǎo)致獲得的數(shù)學(xué)的表面可能不是一個(gè)真正的區(qū)域場(chǎng)近似。這種分離是多種處理方法的最基本的，準(zhǔn)確的分離區(qū)域異常與局部異常，可以為重磁數(shù)據(jù)的反演解釋打下堅(jiān)實(shí)基礎(chǔ)，可以讓資料解釋有較好的地質(zhì)效果。在生產(chǎn)過(guò)程中證實(shí) :實(shí)際的勘探工作里采集到的重磁資料，通常是由于各種不同的因素產(chǎn)生的疊加異常，其包含了許多有用的寶貴信息，是不同的規(guī)模和埋藏深度的地質(zhì)體的綜合反饋。重磁勘探具有儀器輕便、采集數(shù)據(jù)快、經(jīng)濟(jì)等特點(diǎn)，可以在地面、海洋、空中以及衛(wèi)星上獲取有用的觀(guān) 測(cè)數(shù)據(jù)。 致 謝 ............................................................ 34 吉林大學(xué)學(xué)士學(xué)位論文 1 緒論 近些年，隨著我國(guó)經(jīng)濟(jì)的飛速發(fā)展，礦產(chǎn)資源的短缺的問(wèn)題已成為限制國(guó)民經(jīng)濟(jì)持續(xù)快速發(fā)展的重大問(wèn)題。 第一節(jié)資料獲取工區(qū)概況 ......................... 錯(cuò)誤 !未定義書(shū)簽。 第二節(jié) matlab 編程成圖解釋分析對(duì)比 .............. 錯(cuò)誤 !未定義書(shū)簽。 第一節(jié) Nuady 非線(xiàn)性濾波法 ....................... 錯(cuò)誤 !未定義書(shū)簽。吉林大學(xué)學(xué)士學(xué)位論文 重磁異常分離： Nuady 非線(xiàn)性濾波的研究 Gravity and magic anomaly separation : study of the Nuady nonlinear filtering 吉林大學(xué)學(xué)士學(xué)位論文 目 錄 緒論 ............................................... 錯(cuò)誤 !未定義書(shū)簽。 第 二 節(jié) 插值切割法 .............................. 錯(cuò)誤 !未定義書(shū)簽。 第三節(jié)理論模型的誤差分析表格 ................... 錯(cuò)誤 !未定義書(shū)簽。 第 二 節(jié) matlab 編程成圖解釋分析對(duì)比 ............. 錯(cuò)誤 !未定義書(shū)簽。所以，資源勘探成了重中之重，地球物理勘探方法有了用武之地。對(duì)于沙漠、森林等環(huán)境非常惡劣的地域，資源勘探仍然存在很多盲區(qū)，重磁航測(cè)的發(fā)展對(duì)這些盲區(qū)有非常重要的意義。若是能從采集到的重磁數(shù)據(jù)中將不同地質(zhì)體產(chǎn)生的各種異常分離出來(lái)，就可以分別研究它的地質(zhì)意義。 在最開(kāi)始的時(shí)候 , 人們把區(qū)域場(chǎng)視為恒定的場(chǎng) , 采取徒手圓滑的技術(shù) , 再后來(lái)漸變?yōu)閿?shù)字處理的方法 , 例如移動(dòng)平均法等。 因此 ,我們建議根據(jù)一些數(shù) 據(jù)的一部分沒(méi)有區(qū)域場(chǎng)效應(yīng) ,使 用迭代方法的逐次逼近的方法或試圖消除區(qū)域異常效應(yīng)的數(shù)據(jù) ,然后尋求表面多項(xiàng)式系數(shù) ,這樣多項(xiàng)式表面提出法律障礙要克服上述需要適當(dāng)?shù)耐瑫r(shí) ,大大增強(qiáng)了活力的方法。 孤立的異常通常有很多的歪 曲 ,很難定量解釋。 naudy 非線(xiàn)性濾波法與插值切割法在異常分離中產(chǎn)生的效果是非 常好的。另一類(lèi)是頻率域?yàn)V波 。從那時(shí)起，非線(xiàn)性濾波用來(lái)去除各種地球物理數(shù)據(jù)異常，基本的原則是異常完全被移除并且這些剩余異常應(yīng)該保持不被改變。在 1991年的 2月，文百紅更加詳細(xì)的介紹了插值切割在消除磁異常在干擾中的應(yīng)用，在經(jīng)過(guò)前人的不斷研究完善，插值切割法已經(jīng)對(duì)異常分離處理有了非常明顯的效果。 基本原理： 異常被分成倆種類(lèi)型：簡(jiǎn)單異常 , 有一個(gè)最大。 C1=g(xr,y)+(g(x+r,y)g(xr,y))/2。 A=g(xr,y)。E=g(x+r,y)。I=g(x,yr)。 R=0 Q=C/2+(B+D+G+H)/8。 第 2節(jié) 插值切割法 吉林大學(xué)學(xué)士學(xué)位論文 6 進(jìn)行處理重磁實(shí)際或者理論數(shù)據(jù)的過(guò)程當(dāng)中，我們將采集到的數(shù)據(jù)看成由區(qū)域異常和局部異常疊加組成的，根據(jù)這些采集的重磁數(shù)據(jù)做出地質(zhì)結(jié)構(gòu)的解釋?zhuān)?要先將這兩種數(shù)據(jù)分離。實(shí)測(cè)數(shù)值與所求得的區(qū)域異常數(shù)值之差，就是局部異常值。 對(duì)于加權(quán)系數(shù)的確定，由下列算法求取 : 1( , ) [ ( , ) ( , ) ]2XB G x y G x r y G x r y? ? ? ? ? ?。 222() 0 c 1 。依次迭代下去，最終有 1l i m ( , ) ( , ) 0nnn R x y R x y??? ??于是有 ( , ) lim ( ,nnR x y R x y??? ） 稱(chēng)這個(gè)區(qū)域異常為切割半徑 r 的區(qū)域異常。為了盡量減少邊部損失，更充分地利用實(shí)測(cè)資料進(jìn)行地質(zhì)解釋?zhuān)枰獙?duì)實(shí)測(cè)數(shù)據(jù)先進(jìn)行擴(kuò)邊處理，以滿(mǎn)足后續(xù)處理的需要。 第 1節(jié) 參數(shù)設(shè)置 模型一的具體參數(shù)：球體的中心坐標(biāo)：（ 0， 0， 15），球體半徑： 15m, 球體的剩余密度： 長(zhǎng)方體 1的坐標(biāo)：（ 0， 0， 40），剩余密度： 長(zhǎng)方體 2的坐標(biāo)：（ 0， 0， 120），剩余密度： 長(zhǎng)方體 3的坐標(biāo)：（ 0， 0， 200），剩余密度： 0 模型二的具體參數(shù)：球體的中心坐標(biāo)：（ 0， 0， 15），球體半徑： 15m, 球體的剩余密度： 長(zhǎng)方體 1的坐標(biāo)：（ 0， 0， 40），剩余密度： 長(zhǎng)方體 2的坐標(biāo)：（ 0， 0， 120），剩余密度： 長(zhǎng)方體 3的坐標(biāo)：（ 0， 0， 200），剩余密度： 0 模型三的具體參數(shù)：球體的中心坐標(biāo)：（ 0， 0， 15），球體半徑： 15m, 球體的剩余密度： 長(zhǎng)方體 1的坐標(biāo)：（ 0， 0， 40），剩余密度： 0 長(zhǎng)方體 2的坐標(biāo)：（ 0， 0， 120），剩余密度： 0 長(zhǎng)方體 3的坐標(biāo)：（ 0， 0， 200），剩余密度： 0 第 2節(jié) matlab編程成圖解釋分析對(duì)比 吉林大學(xué)學(xué)士學(xué)位論文 9 一、模型一 區(qū)域場(chǎng)小于零的情況下由 matlab編程得出的圖像如下： 在切割半徑與窗口半徑都為 20的情況下： Nuandy： Nuandy非線(xiàn)性濾波 Interpolating Cut：一次插值切割 Fixed Number of cut：固定次數(shù)的插值切割 Regional：區(qū)域場(chǎng)異常 Many Cut：連續(xù)插值切割 6 0 0 4 0 0 2 0 0 0 200 400 60087654321x 1 03c o n t r a s t I n t e r p o l a t i n g C u tr e g i o n a lnuandym a n y C u tF i x e d N u m b e r o f c u t102030405060吉林大學(xué)學(xué)士學(xué)位論文 10 圖像解釋分析對(duì)比：圖像中異常的寬度范圍大致在 300到 300之間，在兩端異常變化的平緩，于 150和 150處變化開(kāi)始劇烈。在圖像中可以看出 Nuady非線(xiàn)性濾波的效果在切割半徑與切割窗口都為 20的情況下是比其他三種方法要好。一次插值切割的曲線(xiàn)有較小的波動(dòng)，并且曲線(xiàn)的形態(tài)與區(qū)域場(chǎng)的形態(tài)相差很大。在圖像中可以看出 Nuady非線(xiàn)性濾波的效果在切割半徑與切割窗口都為 40的情況下是比其他三種方法要好。一次插值切割的曲線(xiàn)有較大的波動(dòng)，并且曲線(xiàn)的形態(tài)與區(qū)域場(chǎng)的形態(tài)相差很大。在圖像中可以看出 Nuady非線(xiàn)性濾波的效果在切割半徑與切割窗口都為 60的情況下是比其他三種方法要好。連續(xù)插值切割的異常曲線(xiàn)形態(tài)上比較接近區(qū)域場(chǎng)的異常曲線(xiàn)形態(tài)，沒(méi)有曲線(xiàn)震蕩，而固定次數(shù)的插值切割的異常曲線(xiàn)有劇烈的震蕩，這可能是由于虛假異常造成的，而 Nuandy非線(xiàn)性濾波的異常曲線(xiàn) 形態(tài)上與區(qū)域場(chǎng)的異常曲線(xiàn)形態(tài)相差比較大，在此圖中沒(méi)有明確的看出來(lái)哪種方法更適合異常分離。而 Nuandy 非線(xiàn)性濾波的異常曲線(xiàn)形態(tài)上與區(qū)域場(chǎng)的異常曲線(xiàn)形態(tài)相差比較大。 在切割半徑與窗口半徑都為 70的情況下： 6 0 0