專利名稱:各向異性三維疊前時(shí)間偏移方法
技術(shù)領(lǐng)域:
本發(fā)明屬于油氣田和礦產(chǎn)資源地震勘探中反射地震資料處理技術(shù)領(lǐng)域,涉及地震 資料處理過程中的疊前偏移成像技術(shù)范疇,是三維地震資料偏移成像的各向異性三維疊前 時(shí)間偏移方法。
背景技術(shù):
地震勘探中反射地震資料處理過程中,疊前偏移成像是關(guān)鍵的環(huán)節(jié),而疊前時(shí)間 偏移是疊前偏移成像中的一種重要方法。疊前時(shí)間偏移方法可對一類斷層較為復(fù)雜但速度 橫向變化不是很劇烈的地質(zhì)構(gòu)造較好成像。與疊前深度偏移方法相比,除具有較高的計(jì)算 效率外,其主要的優(yōu)點(diǎn)是只需使用疊加(均方根)速度;這樣可簡單地通過速度掃描等方式 得到恰當(dāng)?shù)乃俣饶P停乇芰耸褂茂B前深度偏移方法面臨的一個(gè)主要困難速度建模。因 此,疊前時(shí)間偏移方法已成為地震勘探領(lǐng)域廣泛應(yīng)用的關(guān)鍵技術(shù)。影響疊前時(shí)間偏移成像效果的因素包括偏移速度、地震波走時(shí)計(jì)算、偏移計(jì)算時(shí) 的偏移孔徑、計(jì)算偏移幅值的權(quán)系數(shù)、偏移算法實(shí)現(xiàn)流程。走時(shí)計(jì)算與偏移速度共同決定了 反射波能否正確歸位,偏移孔徑及其應(yīng)用方式?jīng)Q定了偏移噪音和偏移算法的計(jì)算量,權(quán)系 數(shù)決定了成像幅值能否正確反應(yīng)地下界面的物性參數(shù)變化,偏移算法實(shí)現(xiàn)流程對偏移的計(jì) 算效率和存儲需求有重要影響。對偏移方法而言,成像效果、計(jì)算效率和存儲需求是評價(jià)偏 移方法的重要指標(biāo)。實(shí)驗(yàn)室測試及場地?cái)?shù)據(jù)表明,許多沉積巖具有各向異性;在偏移算法中忽略介質(zhì) 的各向異性,將導(dǎo)致錯(cuò)誤的偏移歸位和丟失(或摸糊)某些地質(zhì)構(gòu)造。實(shí)際地層的速度各 向異性,可用橫向各向同性較好地近似。已有一些疊前時(shí)間偏移方法在計(jì)算地震波走時(shí)考 慮各向異性,但尚沒有準(zhǔn)確考慮各向異性對地震波幅值的影響,以及如何在各向異性情況 下實(shí)現(xiàn)保幅成像。既然在計(jì)算地震波走時(shí)與幅值時(shí)需考慮介質(zhì)的各向異性,必須研究如何在疊前時(shí) 間偏移方法中引入各向異性參數(shù),以及如何估算與偏移算法對應(yīng)的各向異性參數(shù)。類似于 各向同性介質(zhì)中使用疊加速度進(jìn)行偏移計(jì)算,采用疊加(等效)各向異性參數(shù)是一個(gè)主流 的選擇,這一方法使得各向異性參數(shù)建模變得更容易。三維觀測已成為海上和陸上地震勘探的主要觀測方式,三維觀測和三維偏移可對 陡傾角構(gòu)造正確成像,避免了二維觀測和偏移時(shí)側(cè)反射波難以正確歸位的問題。隨著采集 和計(jì)算設(shè)備的發(fā)展,三維采集和三維偏移已變的更容易實(shí)現(xiàn),因此必須發(fā)展相應(yīng)的各向異 性三維疊前時(shí)間偏移技術(shù)?,F(xiàn)行的疊前時(shí)間偏移方法是以雙平方根方程為基礎(chǔ)的,其成像條件實(shí)際上使用了 疊前深度偏移的相關(guān)成像條件,因此成像時(shí)沒有補(bǔ)償?shù)卣鸩▊鞑サ膸缀螖U(kuò)散效應(yīng)。為獲得 保幅的共反射點(diǎn)(CRP)道集,以服務(wù)于疊前反演等利用疊前信息進(jìn)行油氣和流體檢測的相 應(yīng)技術(shù),需發(fā)展保幅,即幅值能否正確反應(yīng)地下界面物性參數(shù)變化的疊前時(shí)間偏移方法。偏移孔徑對疊前時(shí)間偏移是重要的,較好的偏移孔徑對壓制偏移噪音和降低偏移計(jì)算量有更重要作用。小的偏移孔徑可減少偏移計(jì)算量,但存在著不能對陡傾角構(gòu)造正確 成像的風(fēng)險(xiǎn);過大的孔徑又帶來了偏移噪音和較大的計(jì)算量。偏移孔徑的選取與地震波的 入射角度有關(guān),地球介質(zhì)各向異性改變了地震波的傳播路徑,相應(yīng)地改變了入射角度,因此 在確定偏移孔徑時(shí)需考慮介質(zhì)的各向異性。對三維地震資料而言,由于采集信號在沿測線 和垂直測線方向有不同的空間采樣率,地震資料存在方位角的變化,偏移孔徑的選取較二 維情況變得更復(fù)雜。
發(fā)明內(nèi)容
本發(fā)明的目的是提供一種各向異性介質(zhì)中的三維疊前時(shí)間偏移方法,它通過在 偏移過程中自主決定偏移速度和各向異性參數(shù),正確考慮了地球介質(zhì)各向異性對地震波走 時(shí)和幅值的影響;通過給出保幅成像的權(quán)系數(shù),實(shí)現(xiàn)了保幅偏移;通過給出考慮各向異性 的三維時(shí)變偏移孔徑和對偏移孔徑邊緣區(qū)的幅值施加衰減,實(shí)現(xiàn)了陡傾角構(gòu)造成像和壓制 偏移噪音;通過在并行計(jì)算的各個(gè)進(jìn)程之間合理分配疊前地震資料,實(shí)現(xiàn)了高加速比的并 行計(jì)算。本發(fā)明采用的技術(shù)方案是各向異性三維疊前時(shí)間偏移方法,具體步驟包括(1)用多條拖纜或多條測線記錄人工震源激發(fā)的反射地震信號,記錄到磁帶上;( 從磁帶上讀取地震信號,對疊前地震資料進(jìn)行常規(guī)的壓制噪音處理。針對部分 共中心點(diǎn),抽取共中心點(diǎn)道集,對抽取的道集作常規(guī)的NMO動校正速度拾取,對所得結(jié)果做 橫向平均,作為初始偏移速度;(3)將疊前地震資料按偏移距大小排序,對地震資料做偏移孔徑相關(guān)的相等計(jì)算 量分組,將不同組地震資料存放到集群計(jì)算機(jī)的不同計(jì)算節(jié)點(diǎn)上;(4)利用初始偏移速度,對已存放到集群計(jì)算機(jī)各個(gè)計(jì)算節(jié)點(diǎn)上的地震資料,應(yīng)用 時(shí)變偏移孔徑的三維各向同性疊前時(shí)間偏移方法,進(jìn)行并行的各向同性偏移計(jì)算,收集各 計(jì)算節(jié)點(diǎn)的偏移結(jié)果,形成共反射點(diǎn)道集;(5)依據(jù)形成的共反射點(diǎn)道集,確定各向異性偏移使用的偏移速度和各向異性參 數(shù);(6)根據(jù)地下構(gòu)造在不同時(shí)間深度的最大傾角、偏移速度、各向異性參數(shù)和地震道 的偏移距和方位角,確定時(shí)變的三維偏移孔徑和對應(yīng)于偏移孔徑的、用于壓制偏移噪音的 衰減系數(shù);(7)采用基于單程波算子和穩(wěn)相點(diǎn)原理的查表法求得各向異性介質(zhì)中的地震波走 時(shí)和保幅成像的權(quán)系數(shù);(8)在每個(gè)計(jì)算節(jié)點(diǎn),對地震資料的所有地震道循環(huán)。對每一地震道,由偏移孔 徑確定有效的成像區(qū)域和區(qū)域中各成像點(diǎn)偏移計(jì)算的起始成像時(shí)間,利用保幅成像的權(quán)系 數(shù),從起始成像時(shí)間開始計(jì)算各成像點(diǎn)的偏移幅值。對偏移孔徑邊緣區(qū)域的成像點(diǎn),進(jìn)一步 用衰減系數(shù)來衰減偏移幅值。將偏移幅值累加到存放偏移結(jié)果的數(shù)組中對應(yīng)的偏移距上;(9)收集各計(jì)算節(jié)點(diǎn)的偏移結(jié)果,對形成的全部共反射點(diǎn)道集做常規(guī)的各向同性 剩余動校正,將出現(xiàn)明顯拉伸和噪音部分的對應(yīng)數(shù)值置為零,將不同偏移距的偏移結(jié)果疊 加,形成偏移疊加剖面;(10)通過顯示軟件將偏移疊加剖面數(shù)值轉(zhuǎn)換為地下反射構(gòu)造的剖面圖像,剖面圖像將指示地下構(gòu)造的形態(tài)、斷裂部位、斷距大小和地層沉積樣式及地層的波阻抗特征,用于 確定地下生、儲油構(gòu)造和識別油氣儲層。所述的將疊前地震資料按偏移距大小排序,對地震資料做偏移孔徑相關(guān)的相等計(jì) 算量分組,將不同組地震資料存放到集群計(jì)算機(jī)的不同計(jì)算節(jié)點(diǎn)上是這樣實(shí)現(xiàn)的設(shè)地下 構(gòu)造在沿測線和垂直測線方向的兩個(gè)最大傾角為θ y,已得到的橫向均勻的初始偏移
速度為^wii(O,計(jì)算兩個(gè)長度值a(T) = Κ““工yf+ COS2(2ΘΧ) + Λ2 sin2(2ΘΧ) — 2cos(2代)[\ + λ2 sin2(2ΘΧ)]'/2b(T) = Vms(T)-Ttan0y式中T是用單程旅行時(shí)表達(dá)的時(shí)間深度,單位是秒,h是地震道的半偏移距,單位 是米,;1 = /2/(巧皿(7>0是無量綱參數(shù)。時(shí)間深度T處的成像區(qū)域可看作以a (T)和b(T)為 長、短半軸的橢圓面。令偏移計(jì)算的起始時(shí)間深度為Ttl,最大時(shí)間深度為T2,而大于T1偏移 孔徑將不隨時(shí)間深度變化,則偏移距為的地震道的相對偏移計(jì)算量可近似為G(h) = O(T1)-b(Tx ){Τ2-ΤΧ) (1/ΔΓ) + X α(Τ0 + jAT) ■ b(T0 + JAT)
M式中η = (Τ「Τ。)/ΔΤ是一整數(shù),Δ T是時(shí)間深度方向上成像的采樣間距,B(T1) 和b(I\)等可由前一個(gè)公式計(jì)算得到。設(shè)地震資料的最大半偏移距為hmax,最小半偏移距為 hmin,偏移距采樣間距為Ah,可用整數(shù)i= (h_hmin)/Ah+l作為偏移距的索引;統(tǒng)計(jì)有相同i 的地震道,可得到偏移距滿足hmin+ (i-1) Δ h彡h < hmin+i Δ h的地震道道數(shù)Hii (i = 1,1),其 中1 = (hmax_hmin)Mh+l。記所有Hii中最大的為m」,可得無量綱的權(quán)系數(shù)P i = G(hmin+(i_l) Ah)/G(hmin+(j-l) ΔΙι)。如集群計(jì)算機(jī)共有k個(gè)核,即希望同時(shí)有k個(gè)進(jìn)程并行計(jì)算,則計(jì)算 叫=(IjAA)/^,可按如下方式確定各進(jìn)程中包含的地震資料的偏移距索引對進(jìn)程1,求 滿足下式的叫,
/I1M1+1J PiTYil <w0 PlTYil
1=1 (=1若PiTni >2]P1TH1 一m0
/=1 /=1取叫=ηι+1,索引數(shù)為i = 1,II1的地震道將分配給進(jìn)程1 ;對進(jìn)程2,求滿足下 式的n2,
”2n2 +1J P1Ifil <m0< P1TH1
I=M1 +1J=M1+!若
n2M2 +1W0- ^ PiIni > J PiTni — W0取 =n2+l,索引數(shù)為i =叫+1,n2的地震道將分配給進(jìn)程2 ;如此類推,完成全CN 102141633 A
部k個(gè)進(jìn)程的地震資料分配,這樣可使得各進(jìn)程的偏移計(jì)算量近似相同。一般集群計(jì)算機(jī) 的一個(gè)計(jì)算節(jié)點(diǎn)可同時(shí)進(jìn)行幾個(gè)進(jìn)程的計(jì)算,可將相關(guān)進(jìn)程的地震資料一起存放到這個(gè)計(jì) 算節(jié)點(diǎn)上。所述的依據(jù)形成的共反射點(diǎn)道集,確定各向異性偏移使用的偏移速度和各向異性 參數(shù)是這樣實(shí)現(xiàn)的對共反射點(diǎn)道集,利用初始偏移速度做反動校,再做動校正得到新的速 度,對這一速度做空間平滑處理,作為各向異性偏移的初始速度V·。對反動校后的共反射 點(diǎn)道集,再次利用下式進(jìn)行各向異性動校正
權(quán)利要求
1.一種各向異性三維疊前時(shí)間偏移方法,其特征在于采用以下步驟A)用拖纜或測 線記錄人工震源激發(fā)的反射地震信號,記錄到磁帶上;B)從磁帶上讀取地震信號,對疊前 地震資料進(jìn)行常規(guī)的壓制噪音處理,針對部分共中心點(diǎn),抽取共中心點(diǎn)道集,對抽取的道集 作常規(guī)的NMO動校正速度拾取,對所得結(jié)果做橫向平均,作為初始偏移速度;C)將疊前地震 資料按偏移距大小排序,對地震資料做偏移孔徑相關(guān)的相等計(jì)算量分組,將不同組地震資 料存放到集群計(jì)算機(jī)的不同計(jì)算節(jié)點(diǎn)上;D)利用初始偏移速度,對已存放到集群計(jì)算機(jī)各 個(gè)計(jì)算節(jié)點(diǎn)上的地震資料,應(yīng)用時(shí)變偏移孔徑的三維各向同性疊前時(shí)間偏移方法,進(jìn)行并 行的各向同性偏移計(jì)算,收集各計(jì)算節(jié)點(diǎn)的偏移結(jié)果,形成共反射點(diǎn)道集;E)依據(jù)形成的 共反射點(diǎn)道集,確定各向異性偏移使用的偏移速度和各向異性參數(shù);F)根據(jù)地下構(gòu)造在不 同時(shí)間深度的最大傾角、偏移速度、各向異性參數(shù)和地震道的偏移距和方位角,確定時(shí)變的 三維偏移孔徑和對應(yīng)于偏移孔徑的、用于壓制偏移噪音的衰減系數(shù);G)采用基于單程波算 子和穩(wěn)相點(diǎn)原理的查表法求得各向異性介質(zhì)中的地震波走時(shí)和保幅成像的權(quán)系數(shù);H)在 每個(gè)計(jì)算節(jié)點(diǎn),對地震資料的所有地震道循環(huán),對每一地震道,由偏移孔徑確定有效的成像 區(qū)域和區(qū)域中各成像點(diǎn)偏移計(jì)算的起始成像時(shí)間,利用保幅成像的權(quán)系數(shù),從起始成像時(shí) 間開始計(jì)算各成像點(diǎn)的偏移幅值,對偏移孔徑邊緣區(qū)域的成像點(diǎn),進(jìn)一步用衰減系數(shù)來衰 減偏移幅值,將偏移幅值累加到存放偏移結(jié)果的數(shù)組中對應(yīng)的偏移距上;I)收集各計(jì)算節(jié) 點(diǎn)的偏移結(jié)果,對形成的全部共反射點(diǎn)道集做常規(guī)的各向同性剩余動校正,將出現(xiàn)明顯拉 伸和噪音部分的對應(yīng)數(shù)值置為零,將不同偏移距的偏移結(jié)果疊加,形成偏移疊加剖面;J) 通過顯示軟件將偏移疊加剖面數(shù)值轉(zhuǎn)換為地下反射構(gòu)造的剖面圖像,剖面圖像將指示地下 構(gòu)造的形態(tài)、斷裂部位、斷距大小和地層沉積樣式及地層的波阻抗特征,用于確定地下生、 儲油構(gòu)造和識別油氣儲層。
2.根據(jù)權(quán)力要求1所述的一種各向異性三維疊前時(shí)間偏移方法,其特征在于在C步 驟中,所述的將疊前地震資料按偏移距大小排序,對地震資料做偏移孔徑相關(guān)的相等計(jì)算 量分組,將不同組地震資料存放到集群計(jì)算機(jī)的不同計(jì)算節(jié)點(diǎn)上是這樣實(shí)現(xiàn)的設(shè)地下構(gòu) 造在沿測線和垂直測線方向的兩個(gè)最大傾角為θ y,已得到的橫向均勻的初始偏移速 度為^mi(T),計(jì)算兩個(gè)長度值a(T) 二λ + ^\2θχ) + λ28 η\2θ^-2οοΒ(2θχ)[\ + λ2 sm2(10x)Tsin(2《)vKO = Vms(T)Ttaney式中T是用單程旅行時(shí)表達(dá)的時(shí)間深度,單位是秒,h是地震道的半偏移距,單位是米, 夂二 VdiOTO是無量綱參數(shù),時(shí)間深度T處的成像區(qū)域可看作以a(T)和b(T)為長、短半 軸的橢圓面,令偏移計(jì)算的起始時(shí)間深度為Ttl,最大時(shí)間深度為T2,而大于T1偏移孔徑將不 隨時(shí)間深度變化,則偏移距為池的地震道的相對偏移計(jì)算量可近似為GQi) = α{Τλ)-b(T,)(T2 -η)(1/ΔΓ) +1 α(Τ0 + jAT)·b(T0 + jAT)7=1式中n= (T1-Ttl)/Δ T是一整數(shù),ΔΤ是時(shí)間深度方向上成像的采樣間距,a(I\)和b(I\) 等可由前一個(gè)公式計(jì)算得到,設(shè)地震資料的最大半偏移距為hmax,最小半偏移距為hmin,偏移 距采樣間距為Ah,可用整數(shù)i = (h_hmin)/Ah+l作為偏移距的索引;統(tǒng)計(jì)有相同i的地震道,可得到偏移距滿足hmin+ (i-1) Δ h彡h < hmin+i Δ h的地震道道數(shù)Hii (i = 1,1),其中1 = (hmax-hmin)/Ah+l,記所有Hii中最大的為Hij,可得無量綱的權(quán)系數(shù)P丄=G(hmin+(i_l) Ah)/ G(hmin+(j-l) Ah);(ι \ /如集群計(jì)算機(jī)共有k個(gè)核,即希望同時(shí)有k個(gè)進(jìn)程并行計(jì)算,則計(jì)算Y4PlTni k,ν '=I )!可按如下方式確定各進(jìn)程中包含的地震資料的偏移距索引對進(jìn)程1,求滿足下式的Ii1,Iil+1YiPlIni <m0<J^ PiInli=l i=\若ηλW1+!講0-Σ Aw, > Σ Aw,-%取叫二叫+丨,索引數(shù)為i =丨,叫的地震道將分配給進(jìn)程丨;對進(jìn)程?,求滿足下式的 ,"2 2+1Σ p>mi ^wO ^ Σ P'm> /=^1+1 /=^1+1若n2H1 +1^ο-Σ p<m'> Σ p<mi~moJ=^1+1 = | +1取 =n2+l,索引數(shù)為i = ni+l,n2的地震道將分配給進(jìn)程2 ;如此類推,完成全部k 個(gè)進(jìn)程的地震資料分配,這樣可使得各進(jìn)程的偏移計(jì)算量近似相同;一般集群計(jì)算機(jī)的一 個(gè)計(jì)算節(jié)點(diǎn)可同時(shí)進(jìn)行幾個(gè)進(jìn)程的計(jì)算,可將相關(guān)進(jìn)程的地震資料一起存放到這個(gè)計(jì)算節(jié);卜.ο
3.根據(jù)權(quán)力要求1所述的一種各向異性三維疊前時(shí)間偏移方法,其特征在于在E步 驟中,所述的依據(jù)形成的共反射點(diǎn)道集,確定各向異性偏移使用的偏移速度和各向異性參 數(shù)是這樣實(shí)現(xiàn)的對共反射點(diǎn)道集,利用初始偏移速度做反動校,再做動校正得到新的速 度,對這一速度做空間平滑處理,作為各向異性偏移的初始速度V·,對反動校后的共反射 點(diǎn)道集,再次利用下式進(jìn)行各向異性動校正一 二 ^綱么扣例+ 1啊F^ff4(1+ η)4(1+ η) V(1 + 2η)^οH(2h) = (2Tf + ~(2/ ),(1 + 2T1F丄。式中h是地震道的半偏移距,單位是米,T是用單程旅行時(shí)表達(dá)的時(shí)間深度,單位是秒, n是無量綱的各向異性參數(shù),可根據(jù)疊加能量最大和n的空間聯(lián)續(xù)性,S卩n的橫向變化不 應(yīng)過大,確定各空間點(diǎn)的各向異性參數(shù)η ;再分別取速度為0.95V·和1.05V·進(jìn)行上式各 向異性動校正,(0.95V·,η)、(Vnmo, n)和(1.05V·,n)三組參數(shù)中使疊加能量最大、同 相軸最平直的那組即是該點(diǎn)最終的偏移速度和各向異性參數(shù),對全部各向異性參數(shù)還需進(jìn) 行空間平滑處理。
4.根據(jù)權(quán)力要求1所述的一種各向異性三維疊前時(shí)間偏移方法,其特征在于在F步 驟中,所述的根據(jù)地下構(gòu)造在不同時(shí)間深度的最大傾角、偏移速度、各向異性參數(shù)和地震道的偏移距和方位角,確定時(shí)變的三維偏移孔徑和對應(yīng)于偏移孔徑的、用于壓制偏移噪音的 衰減系數(shù)是這樣實(shí)現(xiàn)的用一組整數(shù)("K,咚)(/ = 1,幻描述三維偏移孔徑,其中k是有效 成像區(qū)域中總的離散點(diǎn)數(shù),令ΔΧ和Ay是偏移成像結(jié)果在兩個(gè)水平坐標(biāo)方向的空間采樣, ΔΤ是時(shí)間深度方向的采樣,則是成像區(qū)域中某離散點(diǎn)與地震道中心點(diǎn)的沿兩 個(gè)水平坐標(biāo)方向的距離,而是該點(diǎn)偏移計(jì)算的起始成像時(shí)間;令地震道中心點(diǎn)坐標(biāo)為 和<Δ^,引入時(shí)變的偏移孔徑,就是在該地震道的偏移計(jì)算中,對成像區(qū)域中水平坐標(biāo) 為和如的點(diǎn)(i = 1,k),僅從Γ = 開始進(jìn)行偏移計(jì)算;在幾個(gè)關(guān)鍵深度,定義地質(zhì)構(gòu)造在沿測線和垂直測線方向的最大傾角為《和《,構(gòu)建如 下不等式方程組
5.根據(jù)權(quán)力要求1所述的一種各向異性三維疊前時(shí)間偏移方法,其特征在于在G步 驟中,所述的采用基于單程波算子和穩(wěn)相點(diǎn)原理的查表法求得各向異性介質(zhì)中的地震波走 時(shí)和保幅成像的權(quán)系數(shù)是這樣實(shí)現(xiàn)的定義
全文摘要
各向異性三維疊前時(shí)間偏移方法,應(yīng)用于地震勘探中反射地震資料處理,是針對三維采集地震資料的疊前偏移成像方法。該方法考慮地球介質(zhì)速度各向異性對地震波傳播的走時(shí)和幅值的影響,能在偏移過程中自主決定偏移速度和各向異性參數(shù),因而能得到歸位準(zhǔn)確、保幅的偏移圖像。該方法根據(jù)地下構(gòu)造時(shí)變的傾角確定各向異性介質(zhì)中的三維時(shí)變偏移孔徑,且能在偏移過程中同時(shí)壓制偏移噪音。該方法通過在集群計(jì)算機(jī)的計(jì)算節(jié)點(diǎn)間合理分配地震資料,實(shí)現(xiàn)高效的并行計(jì)算。該方法的核心是應(yīng)用深度偏移的單程波算子和穩(wěn)相點(diǎn)原理求得各向異性介質(zhì)中地震波的走時(shí)、幅值、成像權(quán)系數(shù)以及入射角度。該方法對油氣、礦產(chǎn)資源勘探有重要應(yīng)用價(jià)值。
文檔編號G01V1/28GK102141633SQ201010597160
公開日2011年8月3日 申請日期2010年12月10日 優(yōu)先權(quán)日2010年12月10日
發(fā)明者井涌泉, 張劍鋒, 張江杰, 張 浩 申請人:中國科學(xué)院地質(zhì)與地球物理研究所