一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法及裝置制造方法
【專利摘要】本發(fā)明實(shí)施例提供一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法及裝置，所述方法包括：獲取三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)；針對(duì)某一深度的目標(biāo)地層，按照面元位置，對(duì)三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息進(jìn)行重新排序；根據(jù)重新排序的三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)：根據(jù)目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，獲取三維多分量地震觀測(cè)系統(tǒng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖。本發(fā)明可在地震觀測(cè)系統(tǒng)設(shè)計(jì)階段實(shí)現(xiàn)對(duì)三維多分量地震觀測(cè)系統(tǒng)的炮檢距均勻性效果的評(píng)價(jià)。
【專利說明】一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法及裝置

【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及三維多分量地震觀測(cè)系統(tǒng)分析評(píng)價(jià)技術(shù)，尤其涉及一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法及裝置。

【背景技術(shù)】
[0002]三維地震勘探是石油與天然氣勘探的主要工具。人們很早就認(rèn)識(shí)到地震波是一種彈性波，這種彈性波包含縱波和橫波信息，即多分量地震波信息。多分量地震勘探具有許多優(yōu)勢(shì)，比如可以得到更高分辨率的地震資料和更多的地下巖石物性參數(shù)等。但長(zhǎng)期以來，由于技術(shù)和成本的限制，地震勘探實(shí)際上大多只利用縱波信息，多分量地震勘探仍然面臨許多的困難。由于橫波震源的施工難度極大，如今的多分量地震一般采用縱波激發(fā)，縱波和橫波同時(shí)接收的勘探方式，其主要過程包括:(I)數(shù)據(jù)資料采集。在陸上地震勘探數(shù)據(jù)采集現(xiàn)場(chǎng)一般進(jìn)行下述三個(gè)工作:地震觀測(cè)系統(tǒng)設(shè)計(jì)、根據(jù)設(shè)計(jì)方案在野外布設(shè)震源和檢波器，地震波的激發(fā)和接收。首先在室內(nèi)進(jìn)行地震觀測(cè)系統(tǒng)設(shè)計(jì)，以確定震源點(diǎn)和檢波點(diǎn)的最佳擺放位置。然后根據(jù)設(shè)計(jì)方案，在野外布設(shè)縱波震源、縱波檢波器和橫波檢波器，陸上地震勘探的震源點(diǎn)一般采用炸藥震源，并沿地震測(cè)線等間距布置多個(gè)檢波器來接收地震波信號(hào)，現(xiàn)代地震勘探中檢波器的數(shù)量多在1000或10000個(gè)。震源點(diǎn)在爆炸后產(chǎn)生出地震波，地震波遇巖層界面反射回來被檢波器接收并傳到儀器車，儀器車將檢波器傳來的信號(hào)記錄下來，這就獲得了用以研究地下油氣埋藏情況的地震記錄。(2)地震數(shù)據(jù)處理。地震數(shù)據(jù)處理是把第一步驟采集到的地震數(shù)據(jù)資料(包括縱波地震數(shù)據(jù)和橫波地震數(shù)據(jù))輸入專用電子計(jì)算機(jī)，按不同要求用一系列功能不同的程序進(jìn)行處理運(yùn)算，把數(shù)據(jù)進(jìn)行歸類編排，突出有效的，除去無效和干擾的，最后把經(jīng)過各種處理的數(shù)據(jù)進(jìn)行疊加和偏移，最終分別得到二維或三維的縱波和橫波地震數(shù)據(jù)體文件。(3)資料解釋。資料解釋是把經(jīng)過處理的地震信息變成地質(zhì)成果的過程，包括運(yùn)用波動(dòng)理論和地質(zhì)知識(shí)，綜合地質(zhì)、鉆井、測(cè)井等各項(xiàng)資料，作出構(gòu)造解釋、地層解釋、巖性和烴類檢測(cè)解釋及綜合解釋，繪出有關(guān)成果圖件，對(duì)勘探區(qū)域作出含油氣評(píng)價(jià)，提出鉆探井位置等。
[0003]為了擺脫人工判斷三維地震觀測(cè)系統(tǒng)采樣均衡性所帶來的不確定性，中國(guó)專利申請(qǐng)?zhí)?00610114254.1提供了一種地震觀測(cè)系統(tǒng)定量分析方法，通過定量的分析三維地震觀測(cè)系統(tǒng)采樣均衡性，使得判斷結(jié)果更準(zhǔn)確。中國(guó)專利申請(qǐng)?zhí)?00610114254.1提供的定量分析方法中考慮了四個(gè)象限內(nèi)炮檢距分布的均衡性，沒有從更加精細(xì)的網(wǎng)格尺度上來分析面元的屬性，同時(shí)也沒有綜合考慮各個(gè)面元所有屬性因素(包括覆蓋次數(shù)、炮檢距、方位角等)給出該面元的整體均衡度值(需分別依據(jù)五個(gè)參數(shù)進(jìn)行判斷)，因此在應(yīng)用過程中會(huì)產(chǎn)生判斷結(jié)果不確定的問題。中國(guó)專利申請(qǐng)?zhí)?01010569364.3提供了一種三維地震觀測(cè)系統(tǒng)方案的整體均衡度定量分析方法，通過引入信息論中熵的概念分析三維地震觀測(cè)系統(tǒng)整體采樣的均衡性。但這些專利技術(shù)均基于傳統(tǒng)的縱波勘探，未考慮三維多分量地震觀測(cè)系統(tǒng)的特殊性，及其對(duì)炮檢距分布均勻性的影響。
[0004]對(duì)于三維多分量地震觀測(cè)系統(tǒng)，由于轉(zhuǎn)換波的下行波是縱波，上行波是橫波，故其射線路徑是不對(duì)稱的，進(jìn)而造成轉(zhuǎn)換波炮檢距分布的不均勻。這種炮檢距的不規(guī)則分布并非地震觀測(cè)系統(tǒng)本身的分布特征造成，而是由于轉(zhuǎn)換波的不對(duì)稱特性所帶來的假象，并可直接影響后續(xù)地震成像的效果。但一直以來缺乏適當(dāng)?shù)姆椒▽?duì)三維多分量地震觀測(cè)系統(tǒng)轉(zhuǎn)換波炮檢距分布的均勻性進(jìn)行有效的定量評(píng)估。


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

[0005]本發(fā)明實(shí)施例提供一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法及裝置，以對(duì)三維多分量地震觀測(cè)系統(tǒng)轉(zhuǎn)換波炮檢距分布的均勻性進(jìn)行有效的定量評(píng)估提供了一種技術(shù)方案。
[0006]一方面，本發(fā)明實(shí)施例提供了一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法，所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法包括:
[0007]獲取三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)；
[0008]針對(duì)某一深度的目標(biāo)地層，按照面元位置，對(duì)所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息進(jìn)行重新排序；
[0009]根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù):
[0010]根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，獲取所述三維多分量地震觀測(cè)系統(tǒng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖。
[0011]另一方面，本發(fā)明實(shí)施例提供了一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取裝置，所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取裝置包括:
[0012]獲取單元，用于獲取三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)；
[0013]排序單元，用于針對(duì)某一深度的目標(biāo)地層，按照面元位置，對(duì)所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息進(jìn)行重新排序；
[0014]非均勻性系數(shù)計(jì)算單元，用于根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù):
[0015]非均勻性系數(shù)分布圖繪制單元，用于根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，獲取所述三維多分量地震觀測(cè)系統(tǒng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖。
[0016]上述技術(shù)方案具有如下有益效果:因?yàn)椴捎盟鋈S多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法包括:獲取三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)；針對(duì)某一深度的目標(biāo)地層，按照面元位置，對(duì)所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息進(jìn)行重新排序；根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù):根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，獲取所述三維多分量地震觀測(cè)系統(tǒng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖的技術(shù)手段，所以達(dá)到了如下的技術(shù)效果:克服了現(xiàn)有炮檢距均勻性定量分析方法無法被應(yīng)用于多分量地震觀測(cè)系統(tǒng)轉(zhuǎn)換波偏移噪聲分析的缺陷，提出了一種三維多分量地震觀測(cè)系統(tǒng)轉(zhuǎn)換波炮檢距均勻性的分析方法。利用本發(fā)明所述技術(shù)方案，可將觀測(cè)系統(tǒng)炮檢距均勻性定量分析算法推廣至多分量地震勘探的轉(zhuǎn)換波分析，在地震觀測(cè)系統(tǒng)設(shè)計(jì)階段實(shí)現(xiàn)對(duì)多分量三維地震觀測(cè)系統(tǒng)的炮檢距均勻性效果的評(píng)價(jià)，有力地促進(jìn)了多分量地震勘探的實(shí)際應(yīng)用。

【專利附圖】

【附圖說明】
[0017]為了更清楚地說明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見地，下面描述中的附圖僅僅是本發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以根據(jù)這些附圖獲得其他的附圖。
[0018]圖1為本發(fā)明實(shí)施例一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法流程圖；
[0019]圖2為本發(fā)明實(shí)施例一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取裝置結(jié)構(gòu)示意圖；
[0020]圖3為本發(fā)明實(shí)施例非均勻性系數(shù)計(jì)算單元結(jié)構(gòu)示意圖；
[0021]圖4為本發(fā)明應(yīng)用實(shí)例三維多分量地震觀測(cè)系統(tǒng)轉(zhuǎn)換波傳播過程的路徑示意圖；
[0022]圖5為本發(fā)明應(yīng)用實(shí)例某三維多分量地震觀測(cè)系統(tǒng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖。

【具體實(shí)施方式】
[0023]下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例。基于本發(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。
[0024]如圖1所示，為本發(fā)明實(shí)施例一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法流程圖，所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法包括:
[0025]101、獲取三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)；
[0026]102、針對(duì)某一深度的目標(biāo)地層，按照面元位置，對(duì)所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息進(jìn)行重新排序；
[0027]103、根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù):
[0028]104、根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，獲取所述三維多分量地震觀測(cè)系統(tǒng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖。
[0029]優(yōu)選的，所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)包括:震源點(diǎn)的三維空間坐標(biāo)、檢波點(diǎn)三維空間坐標(biāo)、震源點(diǎn)與檢波點(diǎn)之間的對(duì)應(yīng)關(guān)系；其中，所述震源點(diǎn)的三維空間坐標(biāo)包括:高程；所述檢波點(diǎn)三維空間坐標(biāo)包括:高程；所述面元為棋盤狀水平網(wǎng)格。
[0030]優(yōu)選的，所述根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，包括:根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距；根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距，計(jì)算所述轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)。
[0031]優(yōu)選的，所述根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距，包括:
[0032]利用如下公式計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距:
[0033]R=--^j(xo-xs) +(凡-凡)+(Z0-Zi) +yj(xr -xo) + (yr - V0) + (zr -Z0)
[0034]其中，震源點(diǎn)為(xs, ys, Zs),目標(biāo)點(diǎn)為(X。，y。，Z。)，檢波點(diǎn)為(x" yr, zr)。
[0035]優(yōu)選的，所述根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距，計(jì)算所述轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，包括:
[0036]假設(shè)一個(gè)面元所對(duì)應(yīng)的炮檢距為R1, R2,……，Rn，則利用如下公式計(jì)算該面元炮檢距分布的非均勻系數(shù):
[0037](.如丨Λ,...Α)=ΣΣ 17?

2=1Λ/

1' J
[0038]C代表了該面元處炮檢距分布的非均勻性，C數(shù)值越小則炮檢距分布越均勻。
[0039]如圖2所示，為本發(fā)明實(shí)施例一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取裝置結(jié)構(gòu)示意圖，所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取裝置包括:
[0040]獲取單元21，用于獲取三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)；
[0041]排序單元22，用于針對(duì)某一深度的目標(biāo)地層，按照面元位置，對(duì)所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息進(jìn)行重新排序；
[0042]非均勻性系數(shù)計(jì)算單元23，用于根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù):
[0043]非均勻性系數(shù)分布圖繪制單元24，用于根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，獲取所述三維多分量地震觀測(cè)系統(tǒng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖。
[0044]優(yōu)選的，所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)包括:震源點(diǎn)的三維空間坐標(biāo)、檢波點(diǎn)三維空間坐標(biāo)、震源點(diǎn)與檢波點(diǎn)之間的對(duì)應(yīng)關(guān)系；其中，所述震源點(diǎn)的三維空間坐標(biāo)包括:高程；所述檢波點(diǎn)三維空間坐標(biāo)包括:高程；所述面元為棋盤狀水平網(wǎng)格。
[0045]優(yōu)選的，如圖3所示，為本發(fā)明實(shí)施例非均勻性系數(shù)計(jì)算單元結(jié)構(gòu)示意圖，非均勻性系數(shù)計(jì)算單元23包括:
[0046]轉(zhuǎn)換波炮檢距計(jì)算模塊231，用于根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距；
[0047]非均勻性系數(shù)計(jì)算模塊232，用于根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距，計(jì)算所述轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)。
[0048]優(yōu)選的，所述轉(zhuǎn)換波炮檢距計(jì)算模塊231，進(jìn)一步具體用于利用如下公式計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距:
[0049]R = ^(xo- Xs )2 + (V0 - Vs )2 + (z0 - Zs )2 +」(x,.-X0)2+ (v;.— V0 )2 + (z, - Z0 )2








，
[0050]其中，震源點(diǎn)為(xs, ys, zs),目標(biāo)點(diǎn)為(X。，Y。, z。)，檢波點(diǎn)為(X" yr, zr)。
[0051]優(yōu)選的，所述非均勻性系數(shù)計(jì)算模塊232，進(jìn)一步具體用于假設(shè)一個(gè)面元所對(duì)應(yīng)的炮檢距為R1, R2,……，Rn，則利用如下公式計(jì)算該面元炮檢距分布的非均勻系數(shù):
_嘩各…⑷韻我
[0053]C代表了該面元處炮檢距分布的非均勻性，C數(shù)值越小則炮檢距分布越均勻。
[0054]本發(fā)明實(shí)施例上述技術(shù)方案具有如下有益效果:因?yàn)椴捎盟鋈S多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法包括:獲取三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)；針對(duì)某一深度的目標(biāo)地層，按照面元位置，對(duì)所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息進(jìn)行重新排序；根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù):根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，獲取所述三維多分量地震觀測(cè)系統(tǒng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖的技術(shù)手段，所以達(dá)到了如下的技術(shù)效果:克服了現(xiàn)有炮檢距均勻性定量分析方法無法被應(yīng)用于多分量地震觀測(cè)系統(tǒng)轉(zhuǎn)換波偏移噪聲分析的缺陷，提出了一種三維多分量地震觀測(cè)系統(tǒng)轉(zhuǎn)換波炮檢距均勻性的分析方法。利用本發(fā)明所述技術(shù)方案，可將觀測(cè)系統(tǒng)炮檢距均勻性定量分析算法推廣至多分量地震勘探的轉(zhuǎn)換波分析，在地震觀測(cè)系統(tǒng)設(shè)計(jì)階段實(shí)現(xiàn)對(duì)多分量三維地震觀測(cè)系統(tǒng)的炮檢距均勻性效果的評(píng)價(jià)，有力地促進(jìn)了多分量地震勘探的實(shí)際應(yīng)用。
[0055]以下結(jié)合應(yīng)用實(shí)例對(duì)本發(fā)明實(shí)施例進(jìn)行詳細(xì)說明:
[0056]如圖4所示，為本發(fā)明應(yīng)用實(shí)例三維多分量地震觀測(cè)系統(tǒng)轉(zhuǎn)換波傳播過程的路徑示意圖。由于轉(zhuǎn)換波的下行波是縱波，上行波是橫波，故其射線路徑是不對(duì)稱的。其反射點(diǎn)位于x/(l+vs/vp)處(X為震源點(diǎn)到檢波點(diǎn)的距離，Vs和Vp分別為橫波和縱波的速度)，與縱波的共反射點(diǎn)位于x/2處有較大的差別，且對(duì)于不同深度、不同速度比的地層，其反射點(diǎn)位置各不相同?？傮w而言，反射點(diǎn)比較靠近接收點(diǎn)的一邊，且隨著地層由深至淺更加靠近接收點(diǎn)。因此，對(duì)于轉(zhuǎn)換波地震觀測(cè)系統(tǒng)設(shè)計(jì)而言，應(yīng)該首先確定目標(biāo)地層的深度，再以此深度為基礎(chǔ)進(jìn)行后續(xù)分析工作。
[0057]依據(jù)地下地震速度模型與三維地震觀測(cè)系統(tǒng)設(shè)計(jì)方案，直接計(jì)算三維地震觀測(cè)系統(tǒng)的預(yù)期偏移噪聲。由于采用離散退化方法求解波動(dòng)方程，該方法對(duì)于地下強(qiáng)對(duì)比度介質(zhì)條件與廣角成像范圍也具有很高的計(jì)算精度。由于輸入地震速度模型為三維網(wǎng)格模型，通過合理設(shè)計(jì)網(wǎng)格間距，該方法可適用于任意復(fù)雜的介質(zhì)模型。本發(fā)明應(yīng)用實(shí)例的技術(shù)方案如下:
[0058]I)計(jì)算多分量地震轉(zhuǎn)換波的射線路徑
[0059]地震波總是從震源點(diǎn)(xs，ys，zs)開始下行,在目標(biāo)點(diǎn)hi)發(fā)生等角度反射，然后上行到達(dá)檢波點(diǎn)因此，可以針對(duì)某一深度的目標(biāo)地層，根據(jù)入射波與反射波角度相同的原理得到目標(biāo)點(diǎn)Z0)的位置，進(jìn)而計(jì)算多分量地震轉(zhuǎn)換波的射線路徑。
[0060]2)計(jì)算多分量地震轉(zhuǎn)換波的炮檢距
[0061]將震源點(diǎn)(xs, ys, Zs)與目標(biāo)點(diǎn)(X。，y。，Z0)的距離與目標(biāo)點(diǎn)(x。，y。，z。)和檢波點(diǎn)(xr, yr, zr)的距離的和作為多分量地震轉(zhuǎn)換波的炮檢距
[0062]R = -xJ+ (y0-ys)2 + (-0-zs)2 + V(xr-χο)2 + (>v-凡)2+(-,.-:0f ⑴
[0063]其中，震源點(diǎn)為(xs, ys, zs),目標(biāo)點(diǎn)為(x。，y0, z。)，檢波點(diǎn)為(x^ yr, zr)。
[0064]3)計(jì)算一個(gè)面元所對(duì)應(yīng)的所有炮檢對(duì)(一個(gè)震源點(diǎn)和其對(duì)應(yīng)的一個(gè)檢波點(diǎn))。面元為棋盤狀水平網(wǎng)格，比如200*200個(gè)，每個(gè)網(wǎng)格大小可為12.5米*12.5米或25米*25米等。當(dāng)震源點(diǎn)和檢波點(diǎn)的水平坐標(biāo)中點(diǎn)位于面元內(nèi)時(shí)，我們認(rèn)為該震源點(diǎn)和檢波點(diǎn)對(duì)信息屬于該面元。計(jì)算同一個(gè)面元內(nèi)所有炮檢對(duì)的炮檢距差的倒數(shù)和。假設(shè)一個(gè)面元所對(duì)應(yīng)的炮檢距為Rl，R2,......, Rn,該面元炮檢距非均勻系數(shù)為:
_5] ⑵
[0066]C代表了該面元處炮檢距分布的非均勻性，C數(shù)值越小則炮檢距分布越均勻。
[0067]I)獲取三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)
[0068]三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)包括檢波點(diǎn)的三維空間坐標(biāo)(含高程)、震源點(diǎn)三維空間坐標(biāo)(含高程)，和震源點(diǎn)與檢波點(diǎn)之間的對(duì)應(yīng)關(guān)系。一個(gè)地震觀測(cè)系統(tǒng)設(shè)計(jì)方案包含大量的震源點(diǎn)和檢波點(diǎn)的三維坐標(biāo)信息。這里的震源點(diǎn)指地震勘探普遍采用的人工炸藥震源。檢波點(diǎn)指地震勘探中在地表布設(shè)的檢波器，一般至少有上千個(gè)，用于接收震源發(fā)出經(jīng)地下反射后回到地表的地震波。
[0069]2)針對(duì)某一深度的目標(biāo)地層，按照面元位置，對(duì)多分量觀測(cè)系統(tǒng)數(shù)據(jù)的震源點(diǎn)與檢波點(diǎn)位置信息進(jìn)行重新排序。
[0070]面元為棋盤狀水平網(wǎng)格，比如200*200個(gè)，每個(gè)網(wǎng)格大小可為12.5米*12.5米或25米*25米等。炮檢對(duì)為地震觀測(cè)系統(tǒng)的最小數(shù)據(jù)單位，包含震源點(diǎn)三維空間坐標(biāo)信息和其所對(duì)應(yīng)的一個(gè)檢波點(diǎn)三維空間坐標(biāo)信息。每個(gè)炮檢對(duì)需要用6個(gè)參數(shù)進(jìn)行描述:震源點(diǎn)橫坐標(biāo)xs、震源點(diǎn)縱坐標(biāo)ys、震源點(diǎn)高程zs、檢波點(diǎn)橫坐標(biāo)xr、檢波點(diǎn)縱坐標(biāo)yr和檢波點(diǎn)高程zr。如果一個(gè)炮檢對(duì)所對(duì)應(yīng)的目標(biāo)點(diǎn)(nz。)位于某個(gè)面元位置，則認(rèn)為該炮檢對(duì)的信息屬于該面元位置。值得注意的是，由于轉(zhuǎn)換波的下行波是縱波，上行波是橫波，其射線路徑不對(duì)稱，故目標(biāo)點(diǎn)的水平位置一般不位于炮檢對(duì)的中點(diǎn)。重新排序前，觀測(cè)系統(tǒng)數(shù)據(jù)一般按照震源點(diǎn)位置排序，即按照震源點(diǎn)的信息和其所對(duì)應(yīng)的所有檢波點(diǎn)的信息排序。基于面元位置的重新排列過程，將每個(gè)面元位置所對(duì)應(yīng)的觀測(cè)系統(tǒng)信息整理在一起，為下一步的基于單個(gè)面元的分析計(jì)算進(jìn)行數(shù)據(jù)的準(zhǔn)備。
[0071]3)計(jì)算一個(gè)面元所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布非均勻性系數(shù):
[0072]為了剝離地下介質(zhì)因素的影響，獨(dú)立的分析復(fù)雜地表對(duì)面元振幅能量均衡性的影響，我們對(duì)地下介質(zhì)情況進(jìn)行簡(jiǎn)化，假設(shè)其為均勻介質(zhì)，即介質(zhì)速度為一恒定值，而且僅考慮水平目標(biāo)層位情況。
[0073]地震轉(zhuǎn)換波傳播示意圖如圖4所示。利用公式(I)計(jì)算一個(gè)面元所對(duì)應(yīng)所有轉(zhuǎn)換波的炮檢距，然后利用公式(2)計(jì)算該面元處轉(zhuǎn)換波炮檢距分布的非均勻性。
[0074]4)分別對(duì)某一目的層位的每個(gè)面元，分別計(jì)算各個(gè)面元位置處轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，并繪制轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖:
[0075]如圖5所示，為本發(fā)明應(yīng)用實(shí)例某三維多分量地震觀測(cè)系統(tǒng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖，橫縱坐標(biāo)分別表示X和Y坐標(biāo)，單位米。圖5中的顏色深度值代表每個(gè)面元的波炮檢距分布的非均勻性系數(shù)。每個(gè)面元的值越小，該面元的炮檢距均勻性越好。該技術(shù)可以實(shí)現(xiàn)多分量地震勘探的炮檢距分布均勻性的定量評(píng)估，克服了以往定性方法的缺點(diǎn)。應(yīng)用本發(fā)明可對(duì)轉(zhuǎn)換波炮檢距分布進(jìn)行定量評(píng)價(jià)，得到的轉(zhuǎn)換波炮檢距分布圖。利用該技術(shù)可實(shí)現(xiàn)對(duì)三維多分量地震觀測(cè)系統(tǒng)轉(zhuǎn)換波炮檢距分布的定量分析。
[0076]本發(fā)明實(shí)施例及應(yīng)用實(shí)例從地震勘探采集的源頭實(shí)現(xiàn)了對(duì)轉(zhuǎn)換波所引起的地震波炮檢距均勻性的定量評(píng)估，對(duì)提高多分量地震采集資料保真度提供了保障，進(jìn)而為多分量地震偏移成像、油氣儲(chǔ)層預(yù)測(cè)、油藏描述的可靠性奠定了基礎(chǔ)，具有重要應(yīng)用價(jià)值。
[0077]本領(lǐng)域技術(shù)人員還可以了解到本發(fā)明實(shí)施例列出的各種說明性邏輯塊(illustrative logical block),單元,和步驟可以通過電子硬件、電腦軟件,或兩者的結(jié)合進(jìn)行實(shí)現(xiàn)。為清楚展示硬件和軟件的可替換性(interchangeability),上述的各種說明性部件(illustrative components),單元和步驟已經(jīng)通用地描述了它們的功能。這樣的功能是通過硬件還是軟件來實(shí)現(xiàn)取決于特定的應(yīng)用和整個(gè)系統(tǒng)的設(shè)計(jì)要求。本領(lǐng)域技術(shù)人員可以對(duì)于每種特定的應(yīng)用，可以使用各種方法實(shí)現(xiàn)所述的功能，但這種實(shí)現(xiàn)不應(yīng)被理解為超出本發(fā)明實(shí)施例保護(hù)的范圍。
[0078]本發(fā)明實(shí)施例中所描述的各種說明性的邏輯塊，或單元都可以通過通用處理器，數(shù)字信號(hào)處理器，專用集成電路(ASIC)，現(xiàn)場(chǎng)可編程門陣列或其它可編程邏輯裝置，離散門或晶體管邏輯，離散硬件部件，或上述任何組合的設(shè)計(jì)來實(shí)現(xiàn)或操作所描述的功能。通用處理器可以為微處理器，可選地，該通用處理器也可以為任何傳統(tǒng)的處理器、控制器、微控制器或狀態(tài)機(jī)。處理器也可以通過計(jì)算裝置的組合來實(shí)現(xiàn)，例如數(shù)字信號(hào)處理器和微處理器，多個(gè)微處理器，一個(gè)或多個(gè)微處理器聯(lián)合一個(gè)數(shù)字信號(hào)處理器核，或任何其它類似的配置來實(shí)現(xiàn)。
[0079]本發(fā)明實(shí)施例中所描述的方法或算法的步驟可以直接嵌入硬件、處理器執(zhí)行的軟件模塊、或者這兩者的結(jié)合。軟件模塊可以存儲(chǔ)于RAM存儲(chǔ)器、閃存、ROM存儲(chǔ)器、EPROM存儲(chǔ)器、EEPROM存儲(chǔ)器、寄存器、硬盤、可移動(dòng)磁盤、⑶-ROM或本領(lǐng)域中其它任意形式的存儲(chǔ)媒介中。示例性地，存儲(chǔ)媒介可以與處理器連接，以使得處理器可以從存儲(chǔ)媒介中讀取信息，并可以向存儲(chǔ)媒介存寫信息?？蛇x地，存儲(chǔ)媒介還可以集成到處理器中。處理器和存儲(chǔ)媒介可以設(shè)置于ASIC中，ASIC可以設(shè)置于用戶終端中?？蛇x地，處理器和存儲(chǔ)媒介也可以設(shè)置于用戶終端中的不同的部件中。
[0080]在一個(gè)或多個(gè)示例性的設(shè)計(jì)中，本發(fā)明實(shí)施例所描述的上述功能可以在硬件、軟件、固件或這三者的任意組合來實(shí)現(xiàn)。如果在軟件中實(shí)現(xiàn)，這些功能可以存儲(chǔ)與電腦可讀的媒介上，或以一個(gè)或多個(gè)指令或代碼形式傳輸于電腦可讀的媒介上。電腦可讀媒介包括電腦存儲(chǔ)媒介和便于使得讓電腦程序從一個(gè)地方轉(zhuǎn)移到其它地方的通信媒介。存儲(chǔ)媒介可以是任何通用或特殊電腦可以接入訪問的可用媒體。例如，這樣的電腦可讀媒體可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盤存儲(chǔ)、磁盤存儲(chǔ)或其它磁性存儲(chǔ)裝置，或其它任何可以用于承載或存儲(chǔ)以指令或數(shù)據(jù)結(jié)構(gòu)和其它可被通用或特殊電腦、或通用或特殊處理器讀取形式的程序代碼的媒介。此外，任何連接都可以被適當(dāng)?shù)囟x為電腦可讀媒介，例如，如果軟件是從一個(gè)網(wǎng)站站點(diǎn)、服務(wù)器或其它遠(yuǎn)程資源通過一個(gè)同軸電纜、光纖電纜、雙絞線、數(shù)字用戶線(DSL)或以例如紅外、無線和微波等無線方式傳輸?shù)囊脖话谒x的電腦可讀媒介中。所述的碟片(disk)和磁盤(disc)包括壓縮磁盤、鐳射盤、光盤、DVD、軟盤和藍(lán)光光盤，磁盤通常以磁性復(fù)制數(shù)據(jù)，而碟片通常以激光進(jìn)行光學(xué)復(fù)制數(shù)據(jù)。上述的組合也可以包含在電腦可讀媒介中。
[0081]以上所述的【具體實(shí)施方式】，對(duì)本發(fā)明的目的、技術(shù)方案和有益效果進(jìn)行了進(jìn)一步詳細(xì)說明，所應(yīng)理解的是，以上所述僅為本發(fā)明的【具體實(shí)施方式】而已，并不用于限定本發(fā)明的保護(hù)范圍，凡在本發(fā)明的精神和原則之內(nèi)，所做的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。
【權(quán)利要求】
1.一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法，其特征在于，所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法包括: 獲取三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)； 針對(duì)某一深度的目標(biāo)地層，按照面元位置，對(duì)所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息進(jìn)行重新排序； 根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù): 根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，獲取所述三維多分量地震觀測(cè)系統(tǒng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖。
2.如權(quán)利要求1所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法，其特征在于，所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)包括:震源點(diǎn)的三維空間坐標(biāo)、檢波點(diǎn)三維空間坐標(biāo)、震源點(diǎn)與檢波點(diǎn)之間的對(duì)應(yīng)關(guān)系；其中，所述震源點(diǎn)的三維空間坐標(biāo)包括:高程；所述檢波點(diǎn)三維空間坐標(biāo)包括:高程；所述面元為棋盤狀水平網(wǎng)格。
3.如權(quán)利要求1所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法，其特征在于，所述根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，包括: 根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距； 根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距，計(jì)算所述轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)。
4.如權(quán)利要求3所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法，其特征在于，所述根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距，包括: 利用如下公式計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距:
R = yj{xo - λ?)+ {y0 _ λ ) + (zo ~zs) +」(χ,.-χ0) + (yr ~?ο) +(λ-ζο) 其中，震源點(diǎn)為(Xs, ys, Zs),目標(biāo)點(diǎn)為(X。，y。，z。)，檢波點(diǎn)為(Xr, yr, zr)。
5.如權(quán)利要求4所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取方法，其特征在于，所述根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距，計(jì)算所述轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，包括: 假設(shè)一個(gè)面元所對(duì)應(yīng)的炮檢距為R1, R2,……，Rn，則利用如下公式計(jì)算該面元炮檢距分布的非均勻系數(shù):

i=l j=lJH ~ C代表了該面元處炮檢距分布的非均勻性，C數(shù)值越小則炮檢距分布越均勻。
6.一種三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取裝置，其特征在于，所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取裝置包括: 獲取單元，用于獲取三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)； 排序單元，用于針對(duì)某一深度的目標(biāo)地層，按照面元位置，對(duì)所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息進(jìn)行重新排序； 非均勻性系數(shù)計(jì)算單元，用于根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù): 非均勻性系數(shù)分布圖繪制單元，用于根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)，獲取所述三維多分量地震觀測(cè)系統(tǒng)的轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)分布圖。
7.如權(quán)利要求6所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取裝置，其特征在于，所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)包括:震源點(diǎn)的三維空間坐標(biāo)、檢波點(diǎn)三維空間坐標(biāo)、震源點(diǎn)與檢波點(diǎn)之間的對(duì)應(yīng)關(guān)系；其中，所述震源點(diǎn)的三維空間坐標(biāo)包括:高程；所述檢波點(diǎn)三維空間坐標(biāo)包括:高程；所述面元為棋盤狀水平網(wǎng)格。
8.如權(quán)利要求6所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取裝置，其特征在于，所述非均勻性系數(shù)計(jì)算單元包括: 轉(zhuǎn)換波炮檢距計(jì)算模塊，用于根據(jù)重新排序的所述三維多分量地震觀測(cè)系統(tǒng)數(shù)據(jù)中的震源點(diǎn)與檢波點(diǎn)位置信息，計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距；非均勻性系數(shù)計(jì)算模塊，用于根據(jù)所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距，計(jì)算所述轉(zhuǎn)換波炮檢距分布的非均勻性系數(shù)。
9.如權(quán)利要求8所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取裝置，其特征在于， 所述轉(zhuǎn)換波炮檢距計(jì)算模塊，進(jìn)一步具體用于利用如下公式計(jì)算所述目標(biāo)地層每一個(gè)面元位置處所對(duì)應(yīng)的轉(zhuǎn)換波炮檢距:
R = yl(x?-^)2 +(yr,-)\)2 +(z?- ) +」(χ, -)2 + Gv-y,,? + (-,-?








， 其中，震源點(diǎn)為(xs, ys, zs),目標(biāo)點(diǎn)為(χ。，y。，z。)，檢波點(diǎn)為(Xr, yr, zr)。
10.如權(quán)利要求9所述三維多分量地震觀測(cè)系統(tǒng)炮檢距分布獲取裝置，其特征在于， 所述非均勻性系數(shù)計(jì)算模塊，進(jìn)一步具體用于假設(shè)一個(gè)面元所對(duì)應(yīng)的炮檢距為R1,R2,……，Rn，則利用如下公式計(jì)算該面元炮檢距分布的非均勻系數(shù): c(^，…
. 1 1， C代表了該面元處炮檢距分布的非均勻性，C數(shù)值越小則炮檢距分布越均勻。
【文檔編號(hào)】G01V1/30GK104267434SQ201410466484
【公開日】2015年1月7日 申請(qǐng)日期:2014年9月12日 優(yōu)先權(quán)日:2014年9月12日
【發(fā)明者】魏偉, 符力耘, 管西竹 申請(qǐng)人:中國(guó)科學(xué)院地質(zhì)與地球物理研究所