[發(fā)明專利]一種基于雙譜運(yùn)算的重力異常分離方法

在審

• 申請?zhí)枺?/td>: 202010171594.8
• 申請日: 2020-03-12
• 公開（公告）號: CN113391364A
• 公開（公告）日: 2021-09-14
• 發(fā)明（設(shè)計）人: 吳身華
• 申請（專利權(quán)）人: 吳身華
• 主分類號: G01V7/00
• 分類號: G01V7/00;G01V7/06
• 代理公司: 暫無信息
• 代理人: 暫無信息

本發(fā)明公開了一種基于雙譜運(yùn)算的重力異常分離方法，所述把實(shí)測重力異常看作由區(qū)域異常和局部異常組成，由實(shí)測重力異常去掉區(qū)域異常后的剩余部分稱為剩余異常，習(xí)慣上看做局部異常；一：實(shí)測重力異常g與區(qū)域異常d和剩余異常j的關(guān)系為：g=d+j；二：剩余異常由實(shí)測重力異常減去區(qū)域異常得到：j=g?d；三：在重力勘測中，測點(diǎn)p（x，y）的區(qū)域異常d（x，y）。本發(fā)明具備結(jié)合已知的地質(zhì)和其它物理資料信息數(shù)據(jù)，可準(zhǔn)確的分離出計算出重力異常數(shù)據(jù)的優(yōu)點(diǎn)，解決了現(xiàn)有的重力異常分離方法計算的準(zhǔn)確性地，不能夠配合實(shí)際數(shù)據(jù)進(jìn)行運(yùn)算，分離出的數(shù)據(jù)難以與實(shí)際地質(zhì)情況相匹配的問題。

技術(shù)領(lǐng)域

本發(fā)明涉及雙譜運(yùn)算的重力異常分離技術(shù)領(lǐng)域，具體為一種基于雙譜運(yùn)算的重力異常分離方法。

背景技術(shù)

野外實(shí)測重力異常是由地表到地下深部所有密度不均勻體引起的重力效應(yīng)的疊加，信息非常豐富，重力資料解釋地質(zhì)構(gòu)造的一個實(shí)際問題就是如何確定不同密度地層的分界面深度，這必須首先從疊加重力異常中分離出單純由這個密度分界面引起的重力異常，然后用這個異常進(jìn)行反演解釋，因此，多密度分界面引起的重力異常的分離對其資料解釋應(yīng)用至關(guān)重要，目前，重力異常分離方法眾多，常用的有平均場法、高次導(dǎo)數(shù)法、解析延拓法、趨勢分析法、匹配濾波法、維納濾波法.切割法、優(yōu)選延拓法、小波分析法等。

現(xiàn)有的重力異常分離方法計算的準(zhǔn)確性地，不能夠配合實(shí)際數(shù)據(jù)進(jìn)行運(yùn)算，分離出的數(shù)據(jù)難以與實(shí)際地質(zhì)情況向匹配，為此我們提供了一種基于雙譜運(yùn)算的重力異常分離方法。

發(fā)明內(nèi)容

本發(fā)明的目的在于提供一種基于雙譜運(yùn)算的重力異常分離方法，具備結(jié)合已知的地質(zhì)和其它物理資料信息數(shù)據(jù)，可準(zhǔn)確的分離出計算出重力異常數(shù)據(jù)的優(yōu)點(diǎn)，解決了現(xiàn)有的重力異常分離方法計算的準(zhǔn)確性地，不能夠配合實(shí)際數(shù)據(jù)進(jìn)行運(yùn)算，分離出的數(shù)據(jù)難以與實(shí)際地質(zhì)情況相匹配的問題。

為實(shí)現(xiàn)上述目的，本發(fā)明提供如下技術(shù)方案：一種基于雙譜運(yùn)算的重力異常分離方法，所述把實(shí)測重力異常看作由區(qū)域異常和局部異常組成，由實(shí)測重力異常去掉區(qū)域異常后的剩余部分稱為剩余異常，習(xí)慣上看做局部異常；

一：實(shí)測重力異常g與區(qū)域異常d和剩余異常j的關(guān)系為：g=d+j；

二：剩余異常由實(shí)測重力異常減去區(qū)域異常得到：j=g-d；

三：在重力勘測中，測點(diǎn)p（x，y）的區(qū)域異常d（x，y），通常可以由測點(diǎn)坐標(biāo)x和有的n階多項(xiàng)式來擬合近似：d(x，y)≈a0+a1x+a2y+a3x³+a4xy+a5y²+...+an-1yⁿ，其中，ao，a1，...aN-1為N個待定多項(xiàng)式系數(shù)，N=1/2(n+2)(n+1)；

四：當(dāng)n=2，即二階多項(xiàng)式表示的區(qū)域異常為d(x，y)≈a0+a1x+a2y+a3x³+a4xy+a5y²；

五：重力異常頻率域變換及剩余異常分離，建立笛卡爾直角坐標(biāo)系，使z軸垂直向下，觀測場值的場源分布在z=h(h0)的平面下，則位函數(shù)f(x，y，z)是z=h平面以上空間變量x，y，z的連續(xù)函數(shù)，且有一次和二次可微導(dǎo)數(shù)，對f(x，y，z)取傅氏變換，F(xiàn)(u，v，z)=k^+]_-]f(x，y，z)exp1-2Pi(ux+vy)2dxdy，由尤拉公式e^iU=cosU+isinU得到

ReF(u，v，z)=k^+]_-]f(x，y，z)cos2Pux#cos2Pvydxdy-k^+]_-]f(x，y，z)sin2Pux#sin2Pvydxdy

ImF(u，v，z)=-k^+]_-]f(x，y，z)sin2Pux#cos2Pvydxdy-k^+]_-]f(x，y，z)cos2Pux#sin2Pvydxdy