專利名稱：一種疊前深度偏移方法
技術(shù)領(lǐng)域：
本發(fā)明涉及地震勘探領(lǐng)域，具體的是指一種疊前深度偏移方法。
背景技術(shù)：
偏移是地震數(shù)據(jù)處理的最后一個步驟，也是最重要、關(guān)鍵的一個過程，對于復(fù)雜的 地下介質(zhì)而言更是如此。由于地下介質(zhì)層位的非水平性，造成了地震疊加剖面并不能反映真實的層位形 態(tài)，因此需要對地震數(shù)據(jù)在疊加的基礎(chǔ)上進一步進行偏移處理，校正由于地層傾角存在所 產(chǎn)生的層位反射的假象時移，使地下層位的反射歸位，反映地層的真實起伏形態(tài)和產(chǎn)狀。偏移處理方法分為兩類，一類在疊加的基礎(chǔ)上進行，稱為疊后偏移；一類在疊加前 進行偏移處理，稱為疊前偏移。疊前偏移精度要高于疊后偏移，對復(fù)雜繞射波的歸位能力較 好；但疊前偏移處理需要的計算工作量很大，隨著近十多年來計算機技術(shù)的提高而得到了 廣泛的發(fā)展和應(yīng)用。其中疊前偏移處理按照偏移方法的不同，主要是方法所在的時空域的 差異，分為疊前時間偏移和疊前深度偏移。疊前深度偏移由于輸出剖面直接反映真實空間 深度域的地下介質(zhì)形態(tài)而在應(yīng)用領(lǐng)域得到更大的關(guān)注。目前業(yè)界關(guān)于多分量地震數(shù)據(jù)的疊前深度偏移研究非常少，相對成熟應(yīng)用的還是 疊后時間偏移。多分量地震數(shù)據(jù)的疊前偏移處理是目前地震勘探領(lǐng)域的攻關(guān)熱點和難點， 而且當前研究與嘗試應(yīng)用還大多局限于疊前時間偏移，也遠沒有走入生產(chǎn)實踐應(yīng)用中。而 且大部分研究還局限于在縱波偏移的基礎(chǔ)上，沿用與純縱波偏移相同的思路進行C波或等 效轉(zhuǎn)換波的偏移處理，將多分量地震信號間的空間矢量關(guān)系完全割裂開來，進行單純波的 偏移處理，既損失了波場中蘊含的豐富地下介質(zhì)信息，成像效果也較差。英國愛丁堡大學(xué)的李向陽教授所領(lǐng)導(dǎo)的研究小組(Hengchang Dai，and Xiang-Yang Li. Effect of errors in the migration velocity model of PS-converted waves on traveltime accuracy in prestack Kirchhoff time migration in weak anisotropic media.Geophysics,2008,73 (5) :S195_S205) (Xiang-Yang Li，Jianxin Yuan.Converted—wave moveout and conversion-point equations in layered VTI media :Theory and applications. Journal of Applied Geophysics，2003，34 :297_318) 初步試驗成功了多分量轉(zhuǎn)換波的疊前時間偏移技術(shù)，在一些地區(qū)的多分量地震勘探中初步 試驗和應(yīng)用，但尚需進一步攻關(guān)完善。此外，戴和李針對弱各向異性介質(zhì)轉(zhuǎn)換橫波疊前時 間偏移中速度對偏移結(jié)果的影響進行數(shù)值模擬和誤差的討論。他們分別討論了轉(zhuǎn)換橫波 速度、垂向速度比、有效速度比、各向異性參數(shù)的精度對旅行時計算誤差的貢獻，對于時間 域轉(zhuǎn)換波的成像提高具有較好的參考價值；且一些思路可以為轉(zhuǎn)換波的疊前深度偏移所借 鑒。上世紀80 年代王妙月在 Kuo (Kuo，J. T.，and Dai，T.，1984，Kirchhoff elastic wave migration for the case of non-coincident source and receiver :Geophysics, 49 (8)，1223-1238)的基礎(chǔ)上最早提出了矢量偏移的思想，并推出了彈性波Kirchhoff積分公式(秦福浩、郭亞曦、王妙月，彈性波克?；舴蚍e分偏移法，地球物理學(xué)報，1988，31(5) 577-587)。這是一種均勻各向同性介質(zhì)假設(shè)下的彈性波的多分量同時偏移方法，通過PP、 PS波的走時將二分量的數(shù)據(jù)逆時延拓到地下成像點，得到每個成像點剛剛被激發(fā)時的三個 分量上的振動形式。通過數(shù)值模擬的合成地震記錄，他們對該方法進行了測試，獲得相對較 好的結(jié)果，在這一方向進行了初步的探索，但截至目前沒有進一步的發(fā)展。該方法有以下兩個不足1)得到的是總的二個分量振動能量，沒有實現(xiàn)波場分 離；2)該方法得到的偏移結(jié)果是某一炮點的能量傳播到對于地下反射界面時激發(fā)的震動。 不同炮點對于地下同一反射點激發(fā)時的振動矢量不相同，因此不能用于多炮記錄的疊加。 而且他們的研究是一種近似的疊前深度偏移，即在偏移過程中仍然在時間域進行，只是后 來進行時-深轉(zhuǎn)換把輸出的偏移剖面轉(zhuǎn)換成深度剖面。而且對于走時的計算他們只停留在 單界面直射線追蹤這一層面，對于不同的深度成像也沒有考慮到偏移孔徑的差異，偏移速 度的獲得相對簡單，只是簡單借用疊加速度然后進行層速度的換算，因而效果并不理想。公開號為US 2009/0257308A1的發(fā)明專利保護偏移速度分析領(lǐng)域的一種技術(shù)，該 技術(shù)基于動校正誤差來不斷的更新速度模型，從而使得最終的速度模型逼近于地下介質(zhì)的 真實速度。該專利是基于剩余曲率分析方法的速度分析技術(shù)，而且該技術(shù)基于小炮檢距假 設(shè)和水平層位介質(zhì)假設(shè)，即要求剩余旅行時可以近似為雙曲線或拋物線。而當速度橫向變 化劇烈或大炮檢距情況下，該技術(shù)所基于的算法近似不合適的；而且該專利技術(shù)只針對常 規(guī)縱波數(shù)據(jù)和縱波速度進行處理分析。亓雪冬等(亓雪冬、賈光華、仝兆歧，Kirchhoff疊前時間偏移處理技術(shù)及應(yīng)用，物 探化探計算技術(shù)，2009，31 (2) 126-130)研究的是針對常規(guī)縱波資料處理的疊前時間偏移 技術(shù)及其應(yīng)用，這一技術(shù)在油氣地震勘探領(lǐng)域已是成熟應(yīng)用的、商業(yè)化的生產(chǎn)技術(shù)。 張猛等(張猛、孟祥賓、匡斌等，積分法疊前深度偏移技術(shù)在BS6地區(qū)的應(yīng)用，勘探 地球物理進展，2009，32(1) 48-55)介紹的是成熟的、商業(yè)化應(yīng)用的常規(guī)縱波地震數(shù)據(jù)的 疊前深度偏移技術(shù)及其軟件對勝利油田某區(qū)塊的地震數(shù)據(jù)進行處理試驗和應(yīng)用，并與常規(guī) 縱波的疊后偏移(時間)進行成像效果的對比，從而說明疊前深度偏移的優(yōu)勢和應(yīng)用效果。葉月明等(葉月明、李振春、仝兆歧，基于穩(wěn)定成像條件的保幅疊前深度偏移，石 油地球物理勘探，2009，44(1) 28-32)主要針對常規(guī)縱波疊前深度偏移處理中阻尼型反褶 積處理存在的問題，將分母趨于零時所引起的計算不穩(wěn)定現(xiàn)象通過高斯平滑處理，提高偏 移處理的計算穩(wěn)定性，并同時達到部分壓噪的目的。上述有關(guān)技術(shù)所針對的地震數(shù)據(jù)都是單分量、常規(guī)縱波數(shù)據(jù)，無一涉及多分量地 震數(shù)據(jù)的處理，所涉及的技術(shù)與理論方法也不適用于三分量地震數(shù)據(jù)的處理。公開號為CN 101598805A的發(fā)明專利只是針對縱波與轉(zhuǎn)換橫波傳播時間的差異， 通過速度比掃描和層位追蹤的方法，將同一層位在縱波與轉(zhuǎn)換橫波剖面的不同反射相位進 行標定、拾取，從而將各對應(yīng)相位的縱橫波速度比確定下來，進而實現(xiàn)將轉(zhuǎn)換橫波剖面按照 縱波傳播時間壓縮，使得在相同的縱波傳播時間標度上，方便縱波與轉(zhuǎn)換橫波的層位對比。 其保護的內(nèi)容并不涉及三分量地震數(shù)據(jù)的成像技術(shù)。公開號為CN 1797031A的發(fā)明專利針對常規(guī)縱波數(shù)據(jù)的疊前深度偏移處理中由 于海量數(shù)據(jù)處理耗時、計算效率低的問題，通過相同相位炮的匹配和組合，實現(xiàn)多炮合并一 炮，從而減少偏移處理的炮集數(shù)，在不損傷偏移效果的基礎(chǔ)上實現(xiàn)偏移處理效率的提高。其保護的內(nèi)容不涉及三分量地震數(shù)據(jù)的成像技術(shù)。公開號為CN 101419292A的發(fā)明專利保護的內(nèi)容是針對三分量采集的地震數(shù)據(jù)， 通過疊前預(yù)處理、靜校正、速度分析、動校正、疊加這樣一個常規(guī)的處理流程對轉(zhuǎn)換橫波進 行處理成像。該專利保護的內(nèi)容在三分量地震數(shù)據(jù)處理領(lǐng)域一般又稱為常規(guī)疊后處理方 法，僅涉及時間偏移技術(shù)。公開號為US 2004/0117123A1的發(fā)明專利保護的內(nèi)容是針對常規(guī)純縱波數(shù)據(jù)的 疊前時間偏移處理技術(shù)，而且該技術(shù)是在頻率_波數(shù)域?qū)崿F(xiàn)的。該專利是將地震數(shù)據(jù)變換 到頻率-波數(shù)域，然后根據(jù)旅行走時圖對給定空間任一點的射線參數(shù)進行相移的線性擬合 使得疊前時間偏移算法適應(yīng)速度的橫向變化劇烈的復(fù)雜介質(zhì)成像問題。公開號為US 2010/0054082A1的發(fā)明專利保護的內(nèi)容不涉及具體的偏移算法，更 與多分量疊前深度偏移無關(guān)。在任何偏移過程中均涉及波場的模擬，其并不是通過射線追 蹤走時計算完成的，而是通過聲波方程的有限差分模擬實現(xiàn)。該專利主要針對疊前逆時偏 移處理中需要存儲空間任一點的波場而造成偏移處理中內(nèi)存與硬盤空間的大量消耗問題， 且使得偏移效率低下，為此提出一種節(jié)省計算機內(nèi)存與硬盤空間的技術(shù)改進措施，即在正 演波場的模擬中只保存阻抗界面點上的波場值，其他空間點不參與計算，從而大大改進了 偏移處理的速度，降低了計算資源的消耗。疊前深度偏移是對大地層傾角和速度橫向變化劇烈的復(fù)雜構(gòu)造成像的有效工具。 但是疊前深度偏移依賴更精確的速度模型。常規(guī)的速度建模方法無法達到足夠的精度。同 時，因為疊前深度偏移對速度模型十分敏感，所以它作為一種有力的速度分析工具，發(fā)展出 了兩種偏移速度分析方法深度聚焦分析(DFA，depth-focusing analysis)和剩余曲率分 析(RCA，residual-curvature analysis)。要進行偏移速度分析就必須解決兩個問題1) 建立一個評價速度是否準確的標準；2)如何進行速度更新。深度聚焦分析的原理是根據(jù)疊加能量測量速度誤差。當偏移深度和聚焦深度相同 時，速度達到要求。剩余曲率分析是根據(jù)剩余旅行時差測量速度誤差。當不同炮檢距的成 像深度相同時，速度達到要求。因為速度分析通常是基于一個粗糙的初始迭代速度，因此一個健壯的速度分析方 法在減少迭代次數(shù)的同時還要保證速度的收斂性。傳統(tǒng)的速度分析方法，通常是基于以下 三個假設(shè)①層狀均勻介質(zhì)；②小炮檢距；③水平地層。MacKay禾口Abma(MacKay S. and Abma R. Imaging and velocity estimation with depth-focusing analysis. Geophysics. 1992,57 :1608_1622)的 DFA 方法用到了全部以上 三個近似。在RCA中，速度分析時所采用的偏移方法不同，基于的假設(shè)也不同1.共炮點或 ^Wlk^M^ Al-Yahay (Al-Yahya K. 1989. Velocity analysis by iterative profile migration, Geophysics, 1989，54 :718_729)提出的方法基于以上三個假設(shè)；Lee和 Zhang (Lee W. and Zhang L. Residual shot profile migration, Geophysics,1992,57 815-822)的方法將第三個假設(shè)拓展到小傾角。2.共炮檢距偏移=Deregowski (Deregowsiki S.M. Common-offset migration and velocity analysis.First Break,1990,8 (6) 225-234)的方法使用了前兩個假設(shè)。以上的方法都基于小炮檢距假設(shè)，在小炮檢距情況下，剩余旅行時可以近似為雙曲線或拋物線。當速度橫向變化劇烈時，使用這個近似是不合適的。因此，Liu (Liu Z. and Bleistein N. Migration velocity analysisi Theory and an iterative algorithm， Geophysics，1995，60 :142-153)在 Lafond 禾口 Levander(Lafond C.F. and Levander A. R. Migration moveout analysis and depth focusing，Geophysics，1993，58 :91_100) 的基礎(chǔ)上提出了一種新的RCA方法——基于深度聚焦的共成像點速度分析方法，這種方法 適用于任何炮檢距、傾角和速度變化。

發(fā)明內(nèi)容
本發(fā)明所要解決的技術(shù)問題是針對多分量地震數(shù)據(jù)，提供一種基于矢量波場的三 分量疊前深度偏移方法，不需要進行波場分離，并實現(xiàn)多分量偏移，實現(xiàn)轉(zhuǎn)換波的準確歸 位。本發(fā)明解決上述技術(shù)問題的技術(shù)方案如下一種疊前深度偏移方法，包括以下步 驟步驟1、對地震數(shù)據(jù)進行偏移速度分析；步驟2、通過射線追蹤計算走時、弧長、出射角、入射角；步驟3、進行偏移孔徑計算；步驟4、利用KirchhofT積分法矢量偏移公式進行疊前深度偏移。本發(fā)明的有益效果是本方法將偏移速度分析、偏移孔徑選取、射線追蹤、 Kirchhoff積分公式結(jié)合于一體，不需要進行波場分離，并實現(xiàn)多分量偏移，實現(xiàn)轉(zhuǎn)換波的 準確歸位。在上述技術(shù)方案的基礎(chǔ)上，本發(fā)明還可以做如下改進。進一步，步驟1中所述偏移速度分析采用基于深度聚焦的共成像點速度分析方 法。采用上述進一步方案的有益效果是，基于深度聚焦的共成像點速度分析方法是在 剩余曲率分析方法上修改而成，其收斂速度快，適用于各種速度模型和炮檢距。進一步，所述步驟1包括如下步驟步驟al、初始地質(zhì)模型，并進入步驟a2 ；步驟a2、建立初始速度模型，并進入步驟a3 ；步驟a3、進行基于目標層位底界面的PSDM(pre_stack depth migration，疊前深 度偏移)，并進入步驟a4;步驟a4、進行偏移速度分析，并進入步驟a5 ；步驟a5、判斷是否建立了完整的速度模型，如果為“是”則進入步驟2，如果為“否” 則進入步驟a6 ；步驟a6、進行基于目標層位底界面的PSDM，并進入步驟a7 ；步驟a7、拾取目標層底界面，并進入步驟a8 ；步驟a8、更新速度模型，將目標層移至下一層位，并進入步驟a3。其中，步驟a5中所述完整的速度模型為從淺到深各層位齊全的以速度、密度、層 位埋深作為參數(shù)的地震地質(zhì)模型。進一步，步驟2中，采用直射線追蹤法或者采用基于費馬原理的最小走時射線追蹤法。采用上述進一步方案的有益效果是，采用直射線追蹤法計算速度快，采用基于費 馬原理的最小走時射線追蹤法可以避免射線盲區(qū)。進一步，步驟4中采用適用于非均勻各向異性介質(zhì)的彈性波Kirchhoff積分法矢 量偏移公式進行疊前深度偏移。采用上述進一步方案的有益效果是，適用于非均勻各向異性介質(zhì)的彈性波 Kirchhoff積分法矢量偏移公式，具有波場分離的特性，可以做到偏移的同時分離PP波 和PS波；偏移結(jié)果可分別輸出PP波剖面和PS波剖面，利于分別分析PP波和PS波的 AVO(Amplitude Versus Offset，振幅隨偏移距的變化)特征；使用了彈性波動方程的格林 函數(shù)，提高算法的保幅性；適用于非均勻各向異性介質(zhì)。本發(fā)明在王妙月等人(秦福浩、郭亞曦、王妙月，彈性波克?；舴蚍e分偏移法，地 球物理學(xué)報，1988，31 (5) 577-587)研究的基礎(chǔ)上推出新的偏移公式，并針對其所存在的 兩個不足進行改進。在偏移過程中直接輸出PP波和PS波，在進行多炮疊加的時候?qū)崿F(xiàn)了 相干加強的效果，得到PP和PS波的剖面。


圖1為本發(fā)明疊前深度偏移方法的原理圖一；圖2為本發(fā)明疊前深度偏移方法的原理圖二 ；圖3為本發(fā)明疊前深度偏移方法的原理圖三；圖4為本發(fā)明疊前深度偏移方法的流程圖；圖5為本發(fā)明中偏移速度分析的實現(xiàn)流程；圖6A為采用基于共CCP (Common Converted Point，共轉(zhuǎn)換點)道集的疊后時間偏 移方法對大慶一條2D3C地震數(shù)據(jù)進行偏移之后的PP波效果圖；圖6B為采用基于共CCP道集的疊后時間偏移方法對大慶一條2D3C地震數(shù)據(jù)進行 偏移之后的PS波效果圖；圖7A為采用本發(fā)明方法對大慶的一條2D3C地震數(shù)據(jù)進行的偏移處理之后的PP 波效果圖；圖7B為采用本發(fā)明方法對大慶的一條2D3C地震數(shù)據(jù)進行的偏移處理之后的PS 波效果圖；圖8為本發(fā)明疊前深度偏移方法的模擬數(shù)據(jù)模型；圖9A為圖8的數(shù)據(jù)模型的地震數(shù)據(jù)Z分量記錄圖；圖9B為圖8的數(shù)據(jù)模型的地震數(shù)據(jù)X分量記錄圖；圖10為圖8的數(shù)據(jù)模型中采用本發(fā)明的疊前深度偏移方法得到的縱波偏移深度 剖面結(jié)果；圖11為圖8的數(shù)據(jù)模型中采用本發(fā)明的疊前深度偏移方法得到的轉(zhuǎn)換橫波偏移 深度剖面結(jié)果；圖12為圖8的數(shù)據(jù)模型中采用聲波方程偏移并隨傾角開孔徑得到的縱波偏移深 度剖面結(jié)果；圖13為圖8的數(shù)據(jù)模型中采用聲波方程偏移并隨傾角開孔徑得到的轉(zhuǎn)換橫波偏移深度剖面結(jié)果；圖14為圖8的數(shù)據(jù)模型中采用聲波方程偏移并固定孔徑得到的縱波偏移深度剖 面結(jié)果；圖15為圖8的數(shù)據(jù)模型中采用聲波方程偏移并固定孔徑得到的轉(zhuǎn)換橫波偏移深 度剖面結(jié)果。
具體實施例方式以下結(jié)合附圖對本發(fā)明的原理和特征進行描述，所舉實例只用于解釋本發(fā)明，并 非用于限定本發(fā)明的范圍。方法原理如圖1所示，均勻各向同性介質(zhì)的三維矢量波場延拓方程為
'[(1 - ) Sm3Pn + Iv2PJnP, ] X iim (x', t + tP )/(r Vp ) +(^3 + s^L - 2rJnr3) χ um (x',t + ts )/(rVs ) Hv )Tim(x\t + tP)/(prVP2) +(《 -V,)T3m(^t + ts)/(prVs2)
1 OO OO
>dx{dxf2 其中風(fēng)
'fl^
(1)
》表示質(zhì)點位移，單位為m/s ；,t)表示質(zhì)點速度，單位m ;T(x'
表示應(yīng)力，單位為N/m2 ；m, η = 1,2,3分別表示x，y，ζ坐標，單位m ；f = f/r為單位矢徑， fm, fs分別表示射線出射角的單位向量在m、n坐標上的投影；t = t(x|x')表示震源到 地下繞射點的走時，單位為s ;tp、ts分別代表繞射點到接收點的PP波、PS波走時，單位為s r表示繞射點到接收點的距離，單位為m ；Vp、Vs分別表示縱、橫波速度，單位為m/s ；ν = Vs/ Vp表示縱橫波速度比。 由于自由表面邊界條件，地表應(yīng)力為零，方程(1)的后兩項可以舍去，有
'[(1 ~ 2v2)Sm,rn + 2vXrnr, ] χ um (x',t + tP)/(rVPj
1 OO QQ
>dx\dx\
(2)將此方程展開，將三個分量上的能量分別投影到PP分量(3)式和PS分量(4)式 上，如圖2所示Upp = U1Cos θ sin Φ +u2sin θ sin Φ +u3cos Φ (3)Ups = U1Cos θ cos Φ +u2sin θ cos Φ +u3sin Φ (4)將⑵式展開后帶入⑶式、⑷式，得
1 CO OO
uPP(x,y,z,t) = — f f 4π
2v2frIv2Pr
—f^ ■ ux(x',y',0,t + tp) +. uy(x',y',0,t + tp)
rVprVp
\-Iv2+Iv2Kfz、
+--— · u2 (x', y, o,/+/P)
rVP
dx'dy
\J 5 /fv
8
權(quán)利要求
一種疊前深度偏移方法，包括以下步驟步驟1、對地震數(shù)據(jù)進行偏移速度分析；步驟2、通過射線追蹤計算走時、弧長、出射角、入射角；步驟3、進行偏移孔徑計算；步驟4、利用Kirchhoff積分法矢量偏移公式進行疊前深度偏移。
2.根據(jù)權(quán)利要求1所述的疊前深度偏移方法，其特征在于步驟1中所述偏移速度分 析采用基于深度聚焦的共成像點速度分析方法。
3.根據(jù)權(quán)利要求2所述的疊前深度偏移方法，其特征在于，所述步驟1包括如下步驟 步驟al、初始地質(zhì)模型，并進入步驟a2 ；步驟a2、建立初始速度模型，并進入步驟a3 ；步驟a3、進行基于目標層位底界面的疊前深度偏移，并進入步驟a4 ；步驟a4、進行偏移速度分析，并進入步驟a5 ；步驟a5、判斷是否建立了完整的速度模型，如果為“是”則進入步驟2，如果為“否”則 進入步驟a6 ；步驟a6、進行基于目標層位底界面的疊前深度偏移，并進入步驟a7 ；步驟a7、拾取目標層底界面，并進入步驟a8 ；步驟a8、更新速度模型，將目標層移至下一層位，并進入步驟a3。
4.根據(jù)權(quán)利要求1所述的疊前深度偏移方法，其特征在于步驟2中，采用直射線追蹤 法或者采用基于費馬原理的最小走時射線追蹤法。
5.根據(jù)權(quán)利要求1所述的疊前深度偏移方法，其特征在于步驟4中采用適用于非均 勻各向異性介質(zhì)的彈性波KirchhofT積分法矢量偏移公式進行疊前深度偏移。
全文摘要
本發(fā)明涉及地震勘探領(lǐng)域，具體的是指一種疊前深度偏移方法。該方法包括對地震數(shù)據(jù)進行偏移速度分析；通過射線追蹤計算走時、弧長、出射角、入射角；進行偏移孔徑計算；利用Kirchhoff積分法矢量偏移公式進行疊前深度偏移。本方法將偏移速度分析、偏移孔徑選取、射線追蹤、Kirchhoff積分公式結(jié)合于一體，不需要進行波場分離，并實現(xiàn)多分量同時偏移，實現(xiàn)轉(zhuǎn)換波的準確歸位。從實際勘探數(shù)據(jù)來看，采用本發(fā)明的疊前深度偏移方法，不僅地表成像效果好，成像分辨率高，而且轉(zhuǎn)換波成像剖面上的深層反射同相軸的連續(xù)性提高，構(gòu)造更加清晰。
文檔編號G01V1/30GK101937100SQ20101025532
公開日2011年1月5日 申請日期2010年8月17日 優(yōu)先權(quán)日2010年8月17日
發(fā)明者王赟, 蘆俊 申請人:中國科學(xué)院地質(zhì)與地球物理研究所