重磁異常分離nuady非線(xiàn)性濾波的研究畢業(yè)論文-在線(xiàn)瀏覽
2024-07-30 05:15本頁(yè)面
【正文】 更加詳細(xì)的介紹了插值切割在消除磁異常在干擾中的應(yīng)用，在經(jīng)過(guò)前人的不斷研究完善，插值切割法已經(jīng)對(duì)異常分離處理有了非常明顯的效果。從那時(shí)起，非線(xiàn)性濾波用來(lái)去除各種地球物理數(shù)據(jù)異常，基本的原則是異常完全被移除并且這些剩余異常應(yīng)該保持不被改變?；驹恚寒惓１环殖蓚z種類(lèi)型：簡(jiǎn)單異常, 有一個(gè)最大。對(duì)于正常空間數(shù)據(jù)的1/2差的定義：S=C(B+D+G+H)/4。 C1=g(xr,y)+(g(x+r,y)g(xr,y))/2。 G1=g(x,yr)+(g(x,y+r)g(x,yr))*3/4。 A=g(xr,y)。C=g(x,y)。E=g(x+r,y)。G=g(x,y+r)。I=g(x,yr)。 異常形狀取決于比率R或者 T/S,R=T/S。R=0 Q=C/2+(B+D+G+H)/8。 0=R=4 Q=(B+D+G+H)/3(A+E+F+I)/12。第2節(jié) 插值切割法 進(jìn)行處理重磁實(shí)際或者理論數(shù)據(jù)的過(guò)程當(dāng)中，我們將采集到的數(shù)據(jù)看成由區(qū)域異常和局部異常疊加組成的，根據(jù)這些采集的重磁數(shù)據(jù)做出地質(zhì)結(jié)構(gòu)的解釋?zhuān)鸵葘⑦@兩種數(shù)據(jù)分離。 插值切割法是在多次切割方法的基礎(chǔ)上成長(zhǎng)起來(lái)的一種新的空間域分離區(qū)域與局部異常的新方法[16】。實(shí)測(cè)數(shù)值與所求得的區(qū)域異常數(shù)值之差，就是局部異常值?；驹恚涸O(shè)測(cè)區(qū)D內(nèi)重磁異常G(x，y)由區(qū)域異常側(cè)x，y和局部異常L(x，y)組成，令A(yù)(x，y)代表與點(diǎn)(x，y)相距r的某4點(diǎn)重(磁)異常的平均值，即:(1)式中，r為切割半徑。對(duì)于加權(quán)系數(shù)的確定，由下列算法求取:。取用上述方法得到的區(qū)域異常稱(chēng)為第1次切割的區(qū)域異常，用Rl(x，y)表示。依次迭代下去，最終有于是有稱(chēng)這個(gè)區(qū)域異常為切割半徑r的區(qū)域異常。一般切割半徑r越大，切割出的局部異常所反映的地質(zhì)體深度和規(guī)模就越大。為了盡量減少邊部損失，更充分地利用實(shí)測(cè)資料進(jìn)行地質(zhì)解釋?zhuān)枰獙?duì)實(shí)測(cè)數(shù)據(jù)先進(jìn)行擴(kuò)邊處理，以滿(mǎn)足后續(xù)處理的需要。另一個(gè)假設(shè)是先沉積的地層中由于地質(zhì)構(gòu)造有巖漿侵入并在一定的條件下形成了火成巖礦體，先沉積的地層的密度比周?chē)蟪练e的圍巖的密度要大，形成一個(gè)正的剩余密度。第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ì)比一、模型一 區(qū)域場(chǎng)小于零的情況下由matlab編程得出的圖像如下：在切割半徑與窗口半徑都為20的情況下：Nuandy：Nuandy非線(xiàn)性濾波Interpolating Cut：一次插值切割 Fixed Number of cut：固定次數(shù)的插值切割Regional：區(qū)域場(chǎng)異常 Many Cut：連續(xù)插值切割 圖像解釋分析對(duì)比：圖像中異常的寬度范圍大致在300到300之間，在兩端異常變化的平緩，于150和150處變化開(kāi)始劇烈。連續(xù)插值切割與固定次數(shù)的插值切割在曲線(xiàn)形態(tài)上也和區(qū)域場(chǎng)的形態(tài)相差比較大，波動(dòng)比較小，比一次插值切割的效果要好一點(diǎn)。在圖像中可以看出Nuady非線(xiàn)性濾波的效果在切割半徑與切割窗口都為20的情況下是比其他三種方法要好。粗淺的來(lái)看連續(xù)差值切割的效果是插值切割法中效果最好的。一次插值切割的曲線(xiàn)有較小的波動(dòng)，并且曲線(xiàn)的形態(tài)與區(qū)域場(chǎng)的形態(tài)相差很大。連續(xù)的插值切割效果比較好，曲線(xiàn)比較圓滑，但是形態(tài)上仍然是和區(qū)域場(chǎng)的形態(tài)差別較大。在圖像中可以看出Nuady非線(xiàn)性濾波的效果在切割半徑與切割窗口都為40的情況下是比其他三種方法要好。并且Nuady非線(xiàn)性濾波和連續(xù)差值切割的曲線(xiàn)在不斷的接近區(qū)域場(chǎng)的曲線(xiàn)。一次插值切割的曲線(xiàn)有較大的波動(dòng)，并且曲線(xiàn)的形態(tài)與區(qū)域場(chǎng)的形態(tài)相差很大。連續(xù)的插值切割效果比較好，曲線(xiàn)比較圓滑，但是形態(tài)上仍然是和區(qū)域場(chǎng)的形態(tài)差別較大。在圖像中可以看出Nuady非線(xiàn)性濾波的效果在切割半徑與切割窗口都為60的情況下是比其他三種方法要好。 模型二 區(qū)域場(chǎng)大于零的情況下由matlab編程得出的圖像如下：在切割半徑與窗口半徑都為20的情況下：Cut：一次插值切割 cishu：固定次數(shù)的插值切割Regional：區(qū)域場(chǎng)異常 Nuandy：Nuandy非線(xiàn)性濾波duoci：連續(xù)插值切割圖像解釋對(duì)比分析：圖像中異常的寬度范圍大致在300到300之間，在兩端異常變化的平緩，于150和150處變化開(kāi)始劇烈。連續(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)哪種方法更適合異常分離。一次插值切割的異常曲線(xiàn)與區(qū)域場(chǎng)的異常曲線(xiàn)形態(tài)上相差仍然比較大，但是有向相同的形態(tài)靠近的趨勢(shì)。而Nuandy非線(xiàn)性濾波的異常曲線(xiàn)形態(tài)上與區(qū)域場(chǎng)的異常曲線(xiàn)形態(tài)相差比較大。在切割半徑與窗口半徑都為50的情況下：圖像解釋分析對(duì)比：圖像中異常的寬度范圍大致在300到300之間，在兩端異常變化
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1