一種時移重力增量的去噪方法和裝置的制造方法
【專利摘要】本發(fā)明提供了一種時移重力增量的去噪方法和裝置，其中，該方法包括：通過對工區(qū)內(nèi)所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和時移重力增量進行加權(quán)處理，得到去噪后的當(dāng)前測點的時移重力增量。在本發(fā)明實施例中，采用加權(quán)處理的方法有效解決了現(xiàn)有技術(shù)中進行去噪處理時沒有考慮進行時移微重力多次觀測得到的數(shù)據(jù)與單次觀測多次讀數(shù)所得到的數(shù)據(jù)之間的聯(lián)系的問題，去除了時移重力增量的觀測誤差，即觀測噪聲，使得最終得到的數(shù)據(jù)與實際的地質(zhì)情況吻合度更高。
【專利說明】
_種時移重力増量的去噪方法和裝置
技術(shù)領(lǐng)域
[0001]本發(fā)明涉及石油勘探技術(shù)領(lǐng)域，特別涉及一種時移重力增量的去噪方法和裝置。
【背景技術(shù)】
[0002]在油氣生產(chǎn)過程中，地下的油氣物質(zhì)會重新分布，從而導(dǎo)致地下的重力場產(chǎn)生局部的變化。針對地下重力場變化所進行的測量稱為時移微重力測量，對地球微重力場產(chǎn)生的局部變化進行重復(fù)性測量所得到的數(shù)據(jù)稱為時移微重力數(shù)據(jù)。因此，進行時移微重力測量時需要對同一工區(qū)進行多次觀測，對相同測點進行多次觀測多次讀數(shù)，從而工區(qū)內(nèi)的測點均存在時移微重力增量。
[0003]目前，對微重力勘探數(shù)據(jù)進行去噪處理，S卩，對勘探數(shù)據(jù)中的個別突跳點進行平滑濾波所采用的方法主要有:向上延拓的去噪處理方法以及低通濾波的去噪處理方法。
[0004]然而上述兩種方法并沒有考慮進行時移微重力多次觀測得到的數(shù)據(jù)與單次觀測多次讀數(shù)所得到的數(shù)據(jù)之間的聯(lián)系，因而并不適用于時移微重力增量的去噪處理情況。
[0005]針對上述問題，目前尚未提出有效的解決方案。

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

[0006]本發(fā)明實施例提供了一種時移重力增量的去噪方法和裝置，以解決現(xiàn)有技術(shù)中采用常規(guī)的去噪處理方法對時移重力數(shù)據(jù)進行去噪處理時，無法考慮時移微重力多次觀測得到的數(shù)據(jù)與單次觀測多次讀數(shù)所得到的數(shù)據(jù)之間的聯(lián)系的問題。
[0007]本發(fā)明實施例提供了一種時移重力數(shù)據(jù)的去噪方法，包括:獲取工區(qū)內(nèi)各個測點的時移重力增量，其中，所述時移重力增量為測點在第一時間段內(nèi)多次測量的重力平均值和在第二時間段內(nèi)多次測量的重力平均值之間的差值;獲取所述各個測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差;對所述各個測點中每個測點采用以下操作得到所述各個測點中每個測點去噪后的時移重力增量:選擇當(dāng)前測點作為中心點，將以所述中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和時移重力增量進行加權(quán)處理，得到去噪后的所述當(dāng)前測點的時移重力增量。
[0008]在一個實施例中，在得到所述各個測點中每個測點去噪后的時移重力增量之后，還包括:將所述各個測點中每個測點去噪后的時移重力增量作為未去噪的時移重力增量，對所述各個測點中每個測點重復(fù)執(zhí)行以下去噪操作直至達到預(yù)設(shè)的重復(fù)次數(shù):選擇當(dāng)前測點作為中心點，將以所述中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和未去噪的時移重力增量進行加權(quán)處理，得到去噪后的所述當(dāng)前測點的時移重力增量。
[0009]在一個實施例中，所述預(yù)設(shè)的重復(fù)次數(shù)不小于10次。
[0010]在一個實施例中，按照以下公式計算得到去噪后的所述當(dāng)前測點的時移重力增量:
[0011]當(dāng)A〈(D/2)和B〈(D/2)時，M = N;
[0012]當(dāng)(D/2)〈A〈(3D/4)和(D/2)〈B〈(3D/4)時，M= NXa+PXb;
[0013]當(dāng)(0/2)〈厶〈(30/4)和8>(30/4)、(0/2)〈8〈(30/4)和厶>(30/4)滿足二者之一時，]?=N X c+P X d ；
[0014]當(dāng)A> (3D/4)且 B〉(3D/4)時，M = NXe+PXf;
[0015]其中，A表示以所述預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差，B表示在第二時間段內(nèi)多次測量的重力方差，D表示以所述預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差中的最大值，M表示去噪后的所述當(dāng)前測點的時移重力增量，N表示所述時移重力增量，P表示以所述預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點的時移重力增量的平均值，a、b、c、d、e、f均表示在O %至100 %之間的加權(quán)系數(shù)。
[0016]在一個實施例中，獲取工區(qū)內(nèi)各個測點的時移重力增量，包括:對所述時移動重力增量進行網(wǎng)格化處理;獲取所述各個測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差，包括:對所述各個測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差進行網(wǎng)格化處理。
[0017]在一個實施例中，所述網(wǎng)格化處理方法包括:克里金插值法。
[0018]在一個實施例中，所述預(yù)設(shè)去噪?yún)^(qū)間包括:nXn的矩陣，其中，所述η表示所述預(yù)設(shè)去噪?yún)^(qū)間中的測點個數(shù)。
[0019]本發(fā)明實施例還提供了一種時移重力增量的去噪裝置，包括:平均值獲取模塊，用于獲取工區(qū)內(nèi)各個測點的時移重力增量，其中，所述時移重力增量為測點在第一時間段內(nèi)多次測量的重力平均值和在第二時間段內(nèi)多次測量的重力平均值之間的差值;方差獲取模塊，用于獲取所述各個測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差;加權(quán)處理模塊，用于對所述各個測點中每個測點采用以下操作得到所述各個測點中每個測點去噪后的時移重力增量:選擇當(dāng)前測點作為中心點，將以所述中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和時移重力增量進行加權(quán)處理，得到去噪后的所述當(dāng)前測點的時移重力增量。
[0020]在一個實施例中，所述時移重力增量的去噪裝置還包括:重復(fù)模塊，用于在得到所述各個測點中每個測點去噪后的時移重力增量之后，將所述各個測點中每個測點去噪后的時移重力增量作為未去噪的時移重力增量，對所述各個測點中每個測點重復(fù)執(zhí)行以下去噪操作直至達到預(yù)設(shè)的重復(fù)次數(shù):單次重復(fù)模塊，用于選擇當(dāng)前測點作為中心點，將以所述中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和未去噪的時移重力增量進行加權(quán)處理，得到去噪后的所述當(dāng)前測點的時移重力增量。
[0021 ]在一個實施例中，所述預(yù)設(shè)的重復(fù)次數(shù)不小于10次。
[0022]在上述實施例中，在對時移重力增量進行去噪的過程中，通過對所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和時移重力增量進行加權(quán)處理，得到去噪后的當(dāng)前測點的時移重力增量，即采用加權(quán)處理的方法有效解決了現(xiàn)有技術(shù)中進行去噪處理時，沒有考慮進行時移微重力多次觀測得到的數(shù)據(jù)與單次觀測多次讀數(shù)所得到的數(shù)據(jù)之間的聯(lián)系的問題，去除了時移重力增量的觀測誤差，即觀測噪聲，使得最終得到的數(shù)據(jù)與實際的地質(zhì)情況吻合度更高。
【附圖說明】
[0023]此處所說明的附圖用來提供對本發(fā)明的進一步理解，構(gòu)成本申請的一部分，并不構(gòu)成對本發(fā)明的限定。在附圖中:
[0024]圖1是本發(fā)明實施例的一種時移重力增量的去噪方法的流程圖；
[0025]圖2是本發(fā)明實施例的一種時移重力增量的去噪裝置的一種結(jié)構(gòu)框圖。
【具體實施方式】
[0026]為使本發(fā)明的目的、技術(shù)方案和優(yōu)點更加清楚明白，下面結(jié)合實施方式和附圖，對本發(fā)明做進一步詳細說明。在此，本發(fā)明的示意性實施方式及其說明用于解釋本發(fā)明，但并不作為對本發(fā)明的限定。
[0027]針對采用常規(guī)的去噪處理方法對時移重力數(shù)據(jù)進行去噪處理時所出現(xiàn)的問題，發(fā)明人提出了一種時移重力增量的去噪方法，以解決采用向上延拓的去噪處理方法以及低通濾波的去噪處理方法進行去噪處理時，并沒有考慮進行時移微重力多次觀測得到的數(shù)據(jù)與單次觀測多次讀數(shù)所得到的數(shù)據(jù)之間的聯(lián)系的問題。具體地，該一種時移重力增量的去噪方法可以如圖1所示，包括以下步驟:
[0028]步驟101:獲取工區(qū)內(nèi)各個測點的時移重力增量，其中，所述時移重力增量為測點在第一時間段內(nèi)多次測量的重力平均值和在第二時間段內(nèi)多次測量的重力平均值之間的差值；
[0029]步驟102:獲取所述各個測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差；
[0030]其中，本實施例中和重力相關(guān)的數(shù)據(jù)均特指微重力。微重力測量具體指的是，微弱重力異常的測量技術(shù)與方法。微重力也稱為微伽級重力，毫伽是重力測量的常用單位，微伽是毫伽的千分之一。
[0031]在實際應(yīng)用的過程中進行重力測量時，需要使用多臺儀器對同一個測點進行多次測量。其中，重力平均值指的是多臺儀器多次重力測量結(jié)果的平均值。同樣地，重力方差指的是多臺儀器多次重力測量結(jié)果的方差值。
[0032]其中，上述第一時間段和第二時間段為兩個獨立的時間段，二者之間沒有重疊。
[0033]分別對上述觀測得到的各個測點的時移重力增量、第一重力方差以及第二重力方差進行網(wǎng)格化處理。在本實施例中，可以采用克里金插值法分別對上述觀測得到的各個測點的時移重力增量、第一重力方差以及第二重力方差進行網(wǎng)格化處理，從而可以得到時移重力增量網(wǎng)格文件、第一重力方差網(wǎng)格文件以及第二重力方差網(wǎng)格文件。
[0034]步驟103:對所述各個測點中每個測點采用以下操作得到所述各個測點中每個測點去噪后的時移重力增量:選擇當(dāng)前測點作為中心點，將以所述中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和時移重力增量進行加權(quán)處理，得到去噪后的當(dāng)前測點的時移重力增量。
[0035]進一步地，在得到各個測點中每個測點去噪后的時移重力增量之后，可以將各個測點中每個測點去噪后的時移重力增量作為未去噪的時移重力增量，對各個測點中每個測點重復(fù)執(zhí)行以下去噪操作直至達到預(yù)設(shè)的重復(fù)次數(shù):選擇當(dāng)前測點作為中心點，將以該中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和未去噪的時移重力增量進行加權(quán)處理，得到去噪后的當(dāng)前測點的時移重力增量。
[0036]其中，上述預(yù)設(shè)去噪?yún)^(qū)間可以包括但不限于:nXn的矩陣，其中，所述η表示預(yù)設(shè)去噪?yún)^(qū)間中的測點個數(shù);同時，上述預(yù)設(shè)的重復(fù)次數(shù)可以設(shè)定為不小于10次。
[0037]在本實施例中，可以按照以下公式計算得到去噪后的當(dāng)前測點的時移重力增量:
[0038]當(dāng)A〈(D/2)和B〈(D/2)時，M = N;
[0039 ]當(dāng)(D/2) <A〈( 3D/4)和(D/2) <B〈( 3D/4)時，M = N X a+P X b;
[0040]當(dāng)(0/2)〈厶〈(30/4)和8>(30/4)、(0/2)〈8〈(30/4)和厶>(30/4)滿足二者之一時，]? =N X c+P X d ；
[0041 ]當(dāng) A> (3D/4)且 B〉(3D/4)時，M = NXe+PXf;
[0042]其中，A表示以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差，B表示在第二時間段內(nèi)多次測量的重力方差，D表示以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差中的最大值，M表示去噪后的當(dāng)前測點的時移重力增量，N表示時移重力增量，P表示以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點的時移重力增量的平均值，a、b、c、d、e、f均表示在O %至100 %之間的加權(quán)系數(shù)。
[0043]例如:a可以設(shè)置為10%，此時，相應(yīng)的，b可以設(shè)置為90% ;c可以設(shè)置為21%，此時，相應(yīng)的，d可以設(shè)置為79 % ; e可以設(shè)置為41 %，此時，相應(yīng)的，f可以設(shè)置為59 %。值得注意的是，上述對a、b、c、d、e、f的取值只是為了更好的說明本申請，也可以取別的數(shù)值，本申請對此不作限定。但是，8、13、0、(1、6、;1^需要滿足以下要求:3+&=100% ,c+d= 100%，e+f =100%。
[0044]下面以對某一測點進行時移重力增量去噪處理為例說明本實施例的去噪過程，可以包括:首先，確定進行去噪處理的當(dāng)前測點并確定去噪處理的重復(fù)次數(shù);再以該測點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和未去噪的時移重力增量進行加權(quán)處理，得到去噪后的當(dāng)前測點的時移重力增量;然后，將得到的去噪后的當(dāng)前測點的時移重力增量作為未去噪的時移重力增量，重復(fù)上述去噪處理的過程，直至達到上述確定的去噪處理的重復(fù)次數(shù);最后，得到去噪后的當(dāng)前測點的時移重力增量。對所測定的時移重力增量網(wǎng)格文件、第一重力方差網(wǎng)格文件以及第二重力方差網(wǎng)格文件中的所有測點均執(zhí)行上述去噪過程，可以得到去噪后的時移重力增量網(wǎng)格文件。
[0045]下面以某地區(qū)的時移重力增量網(wǎng)格文件、第一重力方差網(wǎng)格文件以及第二重力方差網(wǎng)格文件進行去噪處理為例對本實施例所提供的一種時移重力增量的去噪方法進行說明，然而值得注意的是，該具體實施例僅是為了更好地說明本發(fā)明，并不構(gòu)成對本發(fā)明的不當(dāng)限定。
[0046]在本實施例中，可以選擇5X5的矩陣，S卩n = 5，作為進行時移重力增量去噪運算的區(qū)間。當(dāng)然，也可以選擇3X3的矩陣作為進行時移重力增量去噪運算的區(qū)間，本發(fā)明對此不作限定。
[0047]進一步地，設(shè)置時移重力增量去噪處理的重復(fù)次數(shù)為50次。在本實施例中，選擇50次作為去噪處理的重復(fù)次數(shù)是因為進行50去噪后所得到的時移重力增量網(wǎng)格文件的去噪效果最佳。同樣地，進行時移重力去噪處理時也可以進行60次去噪處理，本發(fā)明對此不作限定。
[0048]相應(yīng)地，進行去噪處理時，對選定的以待去噪測點區(qū)間，可以按照以下規(guī)則進行加權(quán)處理:
[0049]當(dāng)A〈4微伽并且B〈4微伽時，M= N;
[0050]當(dāng)4微伽<A〈6微伽并且4微伽<B〈6微伽時，M = NX80%+PX20% ；
[0051 ]當(dāng)4微伽<A〈6微伽并且B>6微伽、或者4微伽<B〈6微伽并且A>6微伽，這兩個條件滿足二者之一時，M=NX60%+PX40% ；
[0052]當(dāng)A>6微伽并且B>6微伽時，M= NX30%+PX70% ；
[0053]其中，A表示以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差，B表示在第二時間段內(nèi)多次測量的重力方差，D表示以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差中的最大值，M表示去噪后的當(dāng)前測點的時移重力增量，N表示時移重力增量，P表示以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點的時移重力增量的平均值。
[0054]對某地區(qū)的時移重力增量網(wǎng)格文件、第一重力方差網(wǎng)格文件以及第二重力方差網(wǎng)格文件進行去噪處理的過程中。首先，可以確定進行去噪處理的當(dāng)前測點;再以該測點為中心，以5X5的矩陣區(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和未去噪的時移重力增量進行加權(quán)處理，得到去噪后的當(dāng)前測點的時移重力增量;然后，將得到的去噪后的當(dāng)前測點的時移重力增量作為未去噪的時移重力增量，重復(fù)上述去噪處理的過程，直至重復(fù)次數(shù)達到50次后，可以得到去噪后的當(dāng)前測點的時移重力增量。對某地區(qū)所測定的時移重力增量網(wǎng)格文件、第一重力方差網(wǎng)格文件以及第二重力方差網(wǎng)格文件中的所有測點均執(zhí)行上述去噪過程，可以得到去噪后的時移重力增量網(wǎng)格文件。
[0055]基于同一發(fā)明構(gòu)思，本發(fā)明實施例中還提供了一種時移重力增量的去噪裝置，如下面的實施例所述。由于時移重力增量的去噪裝置解決問題的原理與時移重力增量的去噪方法相似，因此時移重力增量的去噪裝置的實施可以參見時移重力增量的去噪方法的實施，重復(fù)之處不再贅述。以下所使用的，術(shù)語“單元”或者“模塊”可以實現(xiàn)預(yù)定功能的軟件和/或硬件的組合。盡管以下實施例所描述的裝置較佳地以軟件來實現(xiàn)，但是硬件，或者軟件和硬件的組合的實現(xiàn)也是可能并被構(gòu)想的。圖2是本發(fā)明實施例的時移重力增量的去噪裝置的一種結(jié)構(gòu)框圖，如圖2所示，包括:平均值獲取模塊201、方差獲取模塊202和加權(quán)處理模塊203，下面對該結(jié)構(gòu)進行說明。
[0056]平均值獲取模塊201，用于獲取工區(qū)內(nèi)各個測點的時移重力增量，其中，所述時移重力增量為測點在第一時間段內(nèi)多次測量的重力平均值和在第二時間段內(nèi)多次測量的重力平均值之間的差值；
[0057]方差獲取模塊202，用于獲取所述各個測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差；
[0058]加權(quán)處理模塊203，用于對所述各個測點中每個測點采用以下操作得到所述各個測點中每個測點去噪后的時移重力增量:選擇當(dāng)前測點作為中心點，將以所述中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和時移重力增量進行加權(quán)處理，得到去噪后的所述當(dāng)前測點的時移重力增量。
[0059]在一個實施例中，所述裝置還可以包括:重復(fù)模塊，用于在得到所述各個測點中每個測點去噪后的時移重力增量之后，將所述各個測點中每個測點去噪后的時移重力增量作為未去噪的時移重力增量，對所述各個測點中每個測點重復(fù)執(zhí)行以下去噪操作直至達到預(yù)設(shè)的重復(fù)次數(shù):選擇當(dāng)前測點作為中心點，將以所述中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和未去噪的時移重力增量進行加權(quán)處理，得到去噪后的所述當(dāng)前測點的時移重力增量。
[0060]在一個實施例中，所述預(yù)設(shè)的重復(fù)次數(shù)不小于10次。
[0061]在一個實施例中，所述加權(quán)處理模塊可以按照以下公式計算得到去噪后的所述當(dāng)前測點的時移重力增量:
[0062]當(dāng)A〈(D/2)和B〈(D/2)時，M = N;
[0063 ]當(dāng)(D/2) <A〈( 3D/4)和(D/2) <B〈( 3D/4)時，M = N X a+P X b;
[0064]當(dāng)(0/2)〈厶〈(30/4)和8>(30/4)、(0/2)〈8〈(30/4)和厶>(30/4)滿足二者之一時，]?=N X c+P X d ；
[0065]當(dāng)A>(3D/4)且 B>(3D/4)時，M = NXe+PXf;
[0066]其中，A表示以所述預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差，B表示在第二時間段內(nèi)多次測量的重力方差，D表示以所述預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差中的最大值，M表示去噪后的所述當(dāng)前測點的時移重力增量，N表示所述時移重力增量，P表示以所述預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點的時移重力增量的平均值，a、b、c、d、e、f均表示在O %至100 %之間的加權(quán)系數(shù)。
[0067]在一個實施例中，所述平均值獲取模塊用于對所述時移動重力增量進行網(wǎng)格化處理;所述方差獲取模塊用于對所述各個測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差進行網(wǎng)格化處理。
[0068]在一個實施例中，所述網(wǎng)格化處理方法包括:克里金插值法。
[0069]在一個實施例中，所述預(yù)設(shè)去噪?yún)^(qū)間包括:nXn的矩陣，其中，所述η表示所述預(yù)設(shè)去噪?yún)^(qū)間中的測點個數(shù)。
[0070]從以上的描述中，可以看出，本發(fā)明實施例實現(xiàn)了如下技術(shù)效果:在對時移重力增量進行去噪的過程中，通過對所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和時移重力增量進行加權(quán)處理，得到去噪后的當(dāng)前測點的時移重力增量，即采用加權(quán)處理的方法有效解決了現(xiàn)有技術(shù)中進行去噪處理時沒有考慮進行時移微重力多次觀測得到的數(shù)據(jù)與單次觀測多次讀數(shù)所得到的數(shù)據(jù)之間的聯(lián)系的問題，去除了時移重力增量的觀測誤差，即觀測噪聲，使得最終得到的數(shù)據(jù)與實際的地質(zhì)情況吻合度更尚O
[0071]顯然，本領(lǐng)域的技術(shù)人員應(yīng)該明白，上述的本發(fā)明實施例的各模塊或各步驟可以用通用的計算裝置來實現(xiàn)，它們可以集中在單個的計算裝置上，或者分布在多個計算裝置所組成的網(wǎng)絡(luò)上，可選地，它們可以用計算裝置可執(zhí)行的程序代碼來實現(xiàn)，從而，可以將它們存儲在存儲裝置中由計算裝置來執(zhí)行，并且在某些情況下，可以以不同于此處的順序執(zhí)行所示出或描述的步驟，或者將它們分別制作成各個集成電路模塊，或者將它們中的多個模塊或步驟制作成單個集成電路模塊來實現(xiàn)。這樣，本發(fā)明實施例不限制于任何特定的硬件和軟件結(jié)合。
[0072]以上所述僅為本發(fā)明的優(yōu)選實施例而已，并不用于限制本發(fā)明，對于本領(lǐng)域的技術(shù)人員來說，本發(fā)明實施例可以有各種更改和變化。凡在本發(fā)明的精神和原則之內(nèi)，所作的任何修改、等同替換、改進等，均應(yīng)包含在本發(fā)明的保護范圍之內(nèi)。
【主權(quán)項】
1.一種時移重力增量的去噪方法，其特征在于，包括: 獲取工區(qū)內(nèi)各個測點的時移重力增量，其中，所述時移重力增量為測點在第一時間段內(nèi)多次測量的重力平均值和在第二時間段內(nèi)多次測量的重力平均值之間的差值； 獲取所述各個測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差； 對所述各個測點中每個測點采用以下操作得到所述各個測點中每個測點去噪后的時移重力增量:選擇當(dāng)前測點作為中心點，將以所述中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和時移重力增量進行加權(quán)處理，得到去噪后的當(dāng)前測點的時移重力增量。2.如權(quán)利要求1所述的方法，其特征在于，在得到所述各個測點中每個測點去噪后的時移重力增量之后，所述方法還包括: 將所述各個測點中每個測點去噪后的時移重力增量作為未去噪的時移重力增量，對所述各個測點中每個測點重復(fù)執(zhí)行以下去噪操作直至達到預(yù)設(shè)的重復(fù)次數(shù):選擇當(dāng)前測點作為中心點，將以所述中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和未去噪的時移重力增量進行加權(quán)處理，得到去噪后的當(dāng)前測點的時移重力增量。3.如權(quán)利要求2所述的方法，其特征在于，所述預(yù)設(shè)的重復(fù)次數(shù)不小于10次。4.如權(quán)利要求1所述的方法，其特征在于，按照以下公式計算得到去噪后的所述當(dāng)前測點的時移重力增量: 當(dāng) A〈(D/2)和 B〈(D/2)時，M=N; 當(dāng)(D/2) <A< (3D/4)和(D/2) <B< (3D/4)時，M=N X a+P X b; 當(dāng)(0/2)〈厶〈(30/4)和8>(30/4)、(0/2)〈8〈(30/4)和厶>(30/4)滿足二者之一時，]?=~\0+PXd； 當(dāng)八>(30/4)且8>(30/4)時，]\1=~\6+?\匕 其中，A表示以所述預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差，B表示在第二時間段內(nèi)多次測量的重力方差，D表示以所述預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差中的最大值，M表示去噪后的所述當(dāng)前測點的時移重力增量，N表示所述時移重力增量，P表示以所述預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點的時移重力增量的平均值，a、b、c、d、e、f均表示在O %至100 %之間的加權(quán)系數(shù)。5.如權(quán)利要求1所述的方法，其特征在于，獲取工區(qū)內(nèi)各個測點的時移重力增量，包括:對所述時移動重力增量進行網(wǎng)格化處理； 獲取所述各個測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差，包括:對所述各個測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差進行網(wǎng)格化處理。6.如權(quán)利要求5所述的方法，其特征在于，所述網(wǎng)格化處理方法包括:克里金插值法。7.如權(quán)利要求1所述的方法，其特征在于，所述預(yù)設(shè)去噪?yún)^(qū)間包括:nXn的矩陣，其中，所述η表示所述預(yù)設(shè)去噪?yún)^(qū)間中的測點個數(shù)。8.一種時移重力增量的去噪裝置，其特征在于，包括: 平均值獲取模塊，用于獲取工區(qū)內(nèi)各個測點的時移重力增量，其中，所述時移重力增量為測點在第一時間段內(nèi)多次測量的重力平均值和在第二時間段內(nèi)多次測量的重力平均值之間的差值； 方差獲取模塊，用于獲取所述各個測點在第一時間段內(nèi)多次測量的重力方差以及在第二時間段內(nèi)多次測量的重力方差； 加權(quán)處理模塊，用于對所述各個測點中每個測點采用以下操作得到所述各個測點中每個測點去噪后的時移重力增量:選擇當(dāng)前測點作為中心點，將以所述中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和時移重力增量進行加權(quán)處理，得到去噪后的所述當(dāng)前測點的時移重力增量。9.如權(quán)利要求8所述的裝置，其特征在于，所述裝置還包括: 重復(fù)模塊，用于在得到所述各個測點中每個測點去噪后的時移重力增量之后，將所述各個測點中每個測點去噪后的時移重力增量作為未去噪的時移重力增量，對所述各個測點中每個測點重復(fù)執(zhí)行以下去噪操作直至達到預(yù)設(shè)的重復(fù)次數(shù):選擇當(dāng)前測點作為中心點，將以所述中心點為中心，以預(yù)設(shè)去噪?yún)^(qū)間作為邊界所限定的范圍內(nèi)的所有測點在第一時間段內(nèi)多次測量的重力方差、在第二時間段內(nèi)多次測量的重力方差和未去噪的時移重力增量進行加權(quán)處理，得到去噪后的所述當(dāng)前測點的時移重力增量。10.如權(quán)利要求9所述的裝置，其特征在于，所述預(yù)設(shè)的重復(fù)次數(shù)不小于10次。
【文檔編號】G01V7/06GK105974488SQ201610517205
【公開日】2016年9月28日
【申請日】2016年7月4日
【發(fā)明人】劉云祥, 趙文舉, 鄭勝賢
【申請人】中國石油天然氣集團公司, 中國石油集團東方地球物理勘探有限責(zé)任公司