一種多道瞬變電磁法接收波形記錄合成方法與裝置的制造方法
【專利摘要】本申請提出一種多道瞬變電磁法接收波形記錄合成方法與裝置，涉及煤田水文地質與地球物理領域，包括：根據所需移位寄存器的級數獲得偽隨機碼發射源的編碼信號；獲得大地脈沖響應；通過低通濾波器得到記錄系統的儀器響應；獲得噪聲信號；根據所述發射源信號、所述大地脈沖響應、所述儀器響應和所述噪聲信號獲得多道瞬變電磁法的接收波形記錄。能夠預測出系統所觀測的信號，通過所預測出來的觀測信號，掌握對地質構造的先驗信息。
【專利說明】
一種多道瞬變電磁法接收波形記錄合成方法與裝置
技術領域
[0001 ]本發明涉及煤田水文地質與地球物理領域，具體涉及一種多道瞬變電磁法接收波 形記錄合成方法與裝置。
【背景技術】
[0002] 瞬變電磁法(Transient Electromagnetic，簡稱TEM)，是一種建立在電磁感應原 理基礎上的時間域人工源電磁探測方法。傳統瞬變電法(TEM)是利用回線或者接地導線向 地下發射一次場，在一次場斷電后，測量由地下介質感應產生的隨時間變化的二次場，來達 到尋找地質目標體的一種地球物理勘探方法。該方法發射波信號頻帶寬、頻譜信息豐富，一 次激發便可覆蓋探測所需的頻段，大大提高了工作效率。但是，常規瞬變電磁法信號弱，易 受干擾，探測深度一般只有500m，資料處理也主要限于單道處理，還不能很好適應開展第二 空間探礦的需求。
[0003] 發展深部探測和偏移成像精細解釋成為目前的研究熱點。經過多次試驗和研究， 英國愛丁堡大學的Wright和Ziolkowski等提出了多道瞬變電磁（Multi-channel Transient Electromagnetic method，簡稱MTEM)的概念和探測油氣目標體的相關技術。它 借鑒油氣勘探中的地震技術，采用電性源多次發射，陣列式多道接收多次覆蓋的全波場信 息，可以對數據進行類地震處理，在同等發射強度的條件下大幅度提高探測精度和深度，使 探測深度達到2000m以上。這種方法與傳統的瞬變電磁法不同，主要表現為：（1)采用接地源 形式；（2)編碼發射；（3)多道觀測，在測量感應電壓的同時，測量發送電流；（4)通過對接收 電壓與發送電流進行反卷積，得到大地脈沖響應，進行類地震資料處理. 在系統性地研究多道瞬變電磁(MTEM)深部探礦技術之前，很有必要通過合成的方式得 到Μ??Μ的接收波形記錄，以推進方法的發展。
[0004] 若將大地看成線性時不變系統，則多道瞬變電磁系統采集到的響應信號可以表達 成如下形式： ak(xs，xr，t) = i(k，xs，t)*r(xr，t)*g(xs，xr，t)+n(x r，t) 其中，81{(13 41'，1：)表示接收到的總響應，；[(1^43,1：)表示源電流，1'(11^1：)表示記錄系統 的系統響應，k代表第k次數據采集，g(Xs，Xr，t)為大地脈沖響應，n(Xr，t)為噪聲， Xs為發射 源位置，Xr為接收機位置，t為時間。
[0005] 在常規的研究和信號處理中，總是對MTEM系統采集到的響應信號首先經過去噪處 理，去除n( Xr，t)，然后通過解卷積的方式，根據觀測信號ak(Xs，Xr，t)和1仏，心3)，獲得大地 脈沖響應g(x s，Xr，t)。根據大地脈沖響應來推測地質構造。

【發明內容】

[0006] 本發明提供一種多道瞬變電磁法接收波形記錄合成方法與裝置，在地質斷面和噪 聲源已知的情況下，能夠預測出系統所觀測的信號。
[0007] 為了實現上述發明目的，本發明采取的技術方案如下： 一種多道瞬變電磁法接收波形記錄合成方法，包括： 根據所需移位寄存器的級數獲得偽隨機碼發射源的編碼信號； 獲得大地脈沖響應； 通過低通濾波器得到記錄系統的儀器響應； 獲得噪聲彳目號； 根據所述發射源信號、所述大地脈沖響應、所述儀器響應和所述噪聲信號獲得多道瞬 變電磁法的接收波形記錄。
[0008] 可選地，根據所需移位寄存器的級數獲得偽隨機碼發射源的編碼信號包括： 查表獲得偽隨機碼的特征多項式； 利用特征多項式求取序列多項式； 根據序列多項式計算出偽隨機碼的編碼； 根據偽隨機碼編碼獲得發射源信號。
[0009] 可選地，獲得大地脈沖響應包括： 在頻率域計算大地的響應； 通過頻率域向時間域的變換得到大地的階躍響應； 通過計算階躍響應的導數得到大地的脈沖響應。
[0010] 可選地，通過低通濾波器得到記錄系統的儀器響應包括： 通過計算海明窗口函數的低通濾波器響應得到記錄系統的儀器響應。
[0011] 可選地，根據所述發射源信號、所述大地脈沖響應、所述儀器響應和所述噪聲信號 獲得多道瞬變電磁法的接收波形記錄包括： 利用如下公式計算獲得多道瞬變電磁法的接收波形記錄： ak(Xs,Xr,t) = i(k,Xs,t)*r(xr,t)*g(Xs,Xr,t)+n(Xr,t) 其中，i(k,xs,t)表示發射源信號，r(xr,t)表示儀器響應，g(xs,x r,t)表示大地脈沖響 應，n(Xr，t)表示噪聲信號，k表示第k次數據采集，^為發射源位置，Xr為接收機位置，t為時 間。
[0012] 本發明還提供一種多道瞬變電磁法接收波形記錄合成裝置，包括： 發射源模塊，用于根據所需移位寄存器的級數獲得偽隨機碼發射源的編碼信號； 大地脈沖模塊，用于獲得大地脈沖響應； 儀器響應模塊，用于給通過低通濾波器得到記錄系統的儀器響應； 噪聲模塊，用于獲得噪聲信號； 合成模塊，用于根據所述發射源信號、所述大地脈沖響應、所述儀器響應和所述噪聲信 號獲得多道瞬變電磁法的接收波形記錄。
[0013] 可選地，所述發射源模塊所需移位寄存器的級數獲得偽隨機碼發射源的編碼信號 是指： 查表獲得偽隨機碼的特征多項式； 利用特征多項式求取序列多項式； 根據序列多項式計算出偽隨機碼的編碼； 根據偽隨機碼編碼獲得發射源信號。
[0014] 可選地，所述大地脈沖模塊獲得大地脈沖響應是指： 在頻率域計算大地的響應； 通過頻率域向時間域的變換得到大地的階躍響應； 通過計算階躍響應的導數得到大地的脈沖響應。
[0015] 可選地，所述儀器響應模塊通過低通濾波器得到記錄系統的儀器響應是指： 通過計算海明窗口函數的低通濾波器響應得到記錄系統的儀器響應。
[0016] 可選地，所述合成模塊根據所述發射源信號、所述大地脈沖響應、所述儀器響應和 所述噪聲信號獲得多道瞬變電磁法的接收波形記錄是指： 利用如下公式計算獲得多道瞬變電磁法波場的系統響應信號： ak(Xs,Xr,t) = i(k,Xs,t)*r(xr,t)*g(Xs,Xr,t)+n(Xr,t) 其中，i(k,xs,t)表示發射源信號，r(xr,t)表示儀器響應，g(xs,x r,t)表示大地脈沖響 應，n(Xr，t)表示噪聲信號，k表示第k次數據采集，^為發射源位置，Xr為接收機位置，t為時 間。
[0017] 本發明和現有技術相比，具有如下有益效果： 本發明的方法和裝置，在地質斷面和噪聲源已知的情況下，能夠預測出系統所觀測的 信號，通過所預測出來的觀測信號，掌握對地質構造的先驗信息。
【附圖說明】
[0018] 圖1是本發明實施例的多道瞬變電磁法接收波形記錄合成方法的流程圖； 圖2是本發明實施例的多道瞬變電磁法接收波形記錄合成裝置的結構示意圖； 圖3是本發明實施例的移位寄存器產生的部分偽隨機信號及其頻譜示意圖； 圖4是本發明實施例的有限元法得到的頻率域曲線響應與解析解的對比示意圖； 圖5是本發明實施例的時間域2D M-TEM正演模擬結果和解析解的對比示意圖； 圖6是本發明實施例的低通濾波器的特性曲線示意圖； 圖7是本發明實施例的4階偽隨機源激發下考慮儀器帶限的的接收波形記錄(信號)示 意圖； 其中圖7a為以4階偽隨機源為激發源的接收波形記錄(信號)示意圖； 圖8是本發明實施例的4階偽隨機源激發下考慮噪聲的儀器信號示意圖； 圖9是本發明實施例的不同收發距的部分接收信號示意圖。
【具體實施方式】
[0019] 為使本發明的發明目的、技術方案和有益效果更加清楚明了，下面結合附圖對本 發明的實施例進行說明，需要說明的是，在不沖突的情況下，本申請中的實施例和實施例中 的特征可以相互任意組合。
[0020] 如圖1所示，本發明實施例提供一種多道瞬變電磁法接收波形記錄合成方法，包 括： 根據所需移位寄存器的級數獲得偽隨機碼發射源的編碼信號； 獲得大地脈沖響應； 通過低通濾波器得到記錄系統的儀器響應； 獲得噪聲信號； 根據所述發射源信號、所述大地脈沖響應、所述儀器響應和所述噪聲信號獲得多道瞬 變電磁法的接收波形記錄。
[0021]根據所需移位寄存器的級數獲得偽隨機碼發射源的編碼信號包括： 查表獲得偽隨機碼的特征多項式； 利用特征多項式求取序列多項式； 根據序列多項式計算出偽隨機碼的編碼； 根據偽隨機碼編碼獲得發射源信號。
[0022]獲得大地脈沖響應包括： 在頻率域計算大地的響應； 通過頻率域向時間域的變換得到大地的階躍響應； 通過計算階躍響應的導數得到大地的脈沖響應。
[0023] 在頻率域計算大地的響應時，使用有限元法。
[0024] 通過低通濾波器得到記錄系統的儀器響應包括： 通過計算海明窗口函數的低通濾波器響應得到記錄系統的儀器響應。
[0025] 海明窗口函數的低通濾波器等同于儀器的帶限作用。
[0026] 本發明實施例可以利用隨機函數，生成不同水平的白噪聲信號，作為噪聲信號。
[0027] 根據所述發射源信號、所述大地脈沖響應、所述儀器響應和所述噪聲信號獲得多 道瞬變電磁法的接收波形記錄包括： 利用如下公式計算獲得多道瞬變電磁法波場的系統響應信號： ak(Xs,Xr,t) = i(k,Xs,t)*r(xr,t)*g(Xs,Xr,t)+n(Xr,t) 其中，i(k,xs,t)表示發射源信號，r(xr,t)表示儀器響應，g(xs,x r,t)表示大地脈沖響 應，n(Xr，t)表示噪聲信號，k表示第k次數據采集，^為發射源位置，Xr為接收機位置，t為時 間。
[0028] 如圖2所示，本發明實施例還提供一種多道瞬變電磁法波場合成裝置，包括： 發射源模塊，用于根據所需移位寄存器的級數獲得偽隨機碼發射源的編碼信號； 大地脈沖模塊，用于獲得大地脈沖響應； 儀器響應模塊，用于給通過低通濾波器得到記錄系統的儀器響應； 噪聲模塊，用于獲得噪聲信號； 合成模塊，用于根據所述發射源信號、所述大地脈沖響應、所述儀器響應和所述噪聲信 號獲得多道瞬變電磁法的接收波形記錄。
[0029] 所述所需移位寄存器的級數獲得偽隨機碼發射源的編碼信號是指： 查表獲得偽隨機碼的特征多項式； 利用特征多項式求取序列多項式； 根據序列多項式計算出偽隨機碼的編碼； 根據偽隨機碼編碼獲得發射源信號。
[0030] 所述大地脈沖模塊獲得大地脈沖響應是指： 在頻率域計算大地的響應； 通過頻率域向時間域的變換得到大地的階躍響應； 通過計算階躍響應的導數得到大地的脈沖響應。
[0031 ]所述儀器響應模塊通過低通濾波器得到記錄系統的儀器響應是指： 通過計算海明窗口函數的低通濾波器響應得到記錄系統的儀器響應。
[0032]所述合成模塊根據所述發射源信號、所述大地脈沖響應、所述儀器響應和所述噪 聲信號獲得多道瞬變電磁法的接收波形記錄是指： 利用如下公式計算獲得多道瞬變電磁法波場的系統響應信號： ak(Xs,Xr,t) = i(k,Xs,t)*r(xr,t)*g(Xs,Xr,t)+n(Xr,t) 其中，i(k,xs,t)表示發射源信號，r(xr,t)表示儀器響應，g(xs,x r,t)表示大地脈沖響 應，n(Xr，t)表示噪聲信號，k表示第k次數據采集，^為發射源位置，Xr為接收機位置，t為時 間。
[0033] 實施例一 1、發射源信號： 一個移位寄存器序列{an} =ao，ai，a2，a3,....與一個序列多項式
之間是 對應的。其中，X為位移運算符號。即akx表不將ak位移一位，而akx2表不將ak 位移兩位。
[0034] 序列多項式G(X)與特征多項式F(X)具有如下關系
本發明實施例直接利用特征多項式來求序列多項式。通過查表可查到不同階偽隨機碼 的特征多項式，從而可以計算出偽隨機碼的編碼形式。
[0035] 圖3a是用上面的方法得到的12級移位寄存器的部分偽隨機編碼信號示意圖，該信 號的碼元個數為4095個，周期為0.0625s，每個碼元的持續時間為1.52588E-05S。圖3b將(a) 的信號經FFT得到的頻譜示意圖，為對16個周期的信號進行采樣，采樣率為64k，用快速付氏 變換得到的頻譜，圖3c將圖3b的頻譜經逆FFT得到的信號示意圖，從圖3c可以看出，該信號 在8192Hz范圍內，各個頻點的頻譜幅值接近水平，與-函數的頻譜響應接近。圖3c是將頻譜 通過付氏反變換得到的信號，對比圖3a和圖3c，二者是一致的，驗證了所用的快速付氏變換 程序的正確性。
[0036] 2、大地脈沖響應： 大地脈沖響應是通過先在頻率域計算大地的響應，然后通過頻率域向時間域的變換來 得到大地的階躍響應，最后通過計算階躍響應的導數來得到大地的脈沖響應。
[0037]假定地下電性結構是二維的，將原點取在地面上，走向方向為X軸，垂直向下的方 向為z軸，位于地面與X軸垂直的方向為y軸，外加源的方向為X方向，此時只存在不隨走向而 變的X方向的電場，若電磁場的時間因子為e^ 1GJt，含源麥克斯韋方程可寫為：
其中，u = Ex，Ex是電場矢量，Jcx是外加源矢量，ω是角頻率，μ是介質的磁導率，σ是電導 率，ω為圓頻率。
[0038]其滿足的外邊界條件為：
其中，k為地層中的波數
介電常數。
[0039] 邊值問題可寫為等價的變分形式，
與邊界法線方向的夾角，積分路徑DABC與DEFC構成了外邊界，其包圍的區域為Ω。
[0040] 使用有限元法求解，將區域為Ω劃分成一系列的矩形單元，用雙線性插值函數來 近似表示小單元內的場值。最終得到線性方程組： KU = P 通過解線性方程便可得到每個單元結點上的場值。
[0041] 選用余弦變換作為從頻率域向時間域轉換的工具，為了得到Is內的時間響應結 果，使它適用于不同時間點的計算（因為不同的時間點所用的頻帶寬度不同）。
[0042]圖4為1000 hm-m的均勻半空間時，距源964m處的電場與磁場的頻率域響應曲線，其 中，圖4a為距源964m處的Ex電場曲線示意圖，圖4b距源964m處的Hy磁場；圖5為該點處的脈 沖響應曲線，其中，圖5a為電場階躍響應示意圖；圖5b為磁場階躍響應示意圖；圖5c電場脈 沖響應示意圖；圖5d磁場脈沖響應。
[0043]兩圖中離散點均為2D有限元的計算結果，實線為解析解。
[0044] 從兩圖可以看出，不管是大地的階躍響應，還是大地的脈沖響應，二者的吻合程度 非常高，說明2D有限元算法及改造的余弦變換程序是正確的。
[0045] 3、儀器響應： 儀器是有頻帶限制的，即接收器的頻帶寬度是有限的，所以合成接收信號時需考慮帶 限的影響。帶限的影響相當于一個低通濾波器，可以設計得到。
[0046]低通濾波器的設計方法有多種，不同的方法得到的濾波器的旁瓣幅值大小以及從 主瓣到旁瓣的過渡帶寬度不同。理想濾波器是過渡帶盡可能窄，旁瓣幅值盡可能小，以減小 能量的泄露，使能量盡可能地保留在主瓣內。本文選用海明窗函數法，因為海明窗可將 99.963 %的能量集中在窗譜的主瓣內，旁瓣的峰值小于主瓣峰值的1 %。阻帶最小衰減為 35dB〇
[0047]本發明實施例使用的濾波器，通帶為8192Hz，阻帶為9600Hz。理想的低通濾波器如 圖6a所示，通帶內的信號可以通過，通帶外全是阻帶，阻帶內的信號全部被濾除。理想的低 通濾波器相應的時間域信號如圖6b中的黑實線所示(實際為無限信號，為了和加窗后的信 號相比，只顯示了其中的一部分），圖6c是海明窗的形態，海明窗的譜如圖6e所示。用海明窗 對圖6b中的信號進行截斷，截斷后的信號如圖6b中的紅色離散點所示，加窗后的信號變換 到頻率域，即為圖6d所示，相對于理想低通濾波器圖6a，出現了通帶與阻帶之間的過濾帶， 阻帶也出現了小幅振蕩（因其相對于主瓣來說，能量較小，振蕩不明顯，若放大該部分，可以 看到振蕩現象）。可以看出其過渡帶較窄，阻帶能量泄漏也很小，所以設計的濾波器是一個 性能優良的濾波器。
[0048]用上述濾波器來模擬儀器的帶限影響，可以得到儀器響應。
[0049]以4階偽隨機源為激發源(7a)，和由2D M-TEM正演程序得到的大地脈沖響應(見圖 7d和圖7e)做褶積，得到的電場和磁場的理想信號分別如圖7b和圖7d所示。考慮了帶限影響 的電場和磁場信號分別如圖7c和7e所示。從這些圖可以看出，不加帶限前，信號比較平滑， 加了帶限后，信號出現了一些毛刺狀的干擾。
[0050] 4、噪聲信號： 在信號中加入不同水平的白噪聲后，信號出現了擾動，如圖8所示。當噪聲較小(噪音水 平為5%)時，對信號的干擾不明顯，但當噪聲達到10%時，對信號的干擾明顯顯現。其中圖 8a為考慮帶限的Ex接收信號的示意圖；圖8b為考慮帶限的Hy接收信號的示意圖；圖8c為加 5%噪聲的Ex接收信號示意圖；圖8d為加5%噪聲的Hy接收信號的示意圖；圖8e為加10%噪 聲的Ex接收信號的示意圖，圖8f為加10%噪聲的Ex接收信號的示意圖。
[0051 ] 實施例二 針對2015年在張北進行的試驗裝置做了數據模擬。給定電阻率是100〇hm-m的均勻半空 間模型，用一個長度是300m的發射源發射一個12階的偽隨機碼作為輸入信號，分別在偏移 距是300m和2700m的地方接收信號，模擬儀器接收的信號如圖9所示。為了和實際采集信號 作對比，使用1024Hz作為基頻，用16kHz采樣，這樣發射一個周期的源大約為4秒，圖9像是其 中的一部分，圖9a為計算的偏移距是300m時Ex部分接收信號的示意圖；圖9b為實測的偏移 距是300m時Ex部分接收信號的示意圖；圖9c為計算的偏移距是2700m時Ex部分接收信號的 示意圖；圖9d為實測的偏移距是2700m時Ex部分接收信號的示意圖。
[0052]從上面的計算結果可以看出計算的信號與實測的信號有一定的相似度，二者不完 全一致的原因是地下介質不同，計算模擬信號用的介質為均勻半空間，實側資料下面是層 狀結構，雖然如此，不同偏移距的信號相似性說明信號合成是正確的。
[0053]雖然本發明所揭示的實施方式如上，但其內容只是為了便于理解本發明的技術方 案而采用的實施方式，并非用于限定本發明。任何本發明所屬技術領域內的技術人員，在不 脫離本發明所揭示的核心技術方案的前提下，可以在實施的形式和細節上做任何修改與變 化，但本發明所限定的保護范圍，仍須以所附的權利要求書限定的范圍為準。
【主權項】
1. 一種多道瞬變電磁法接收波形記錄合成方法，其特征在于，包括： 根據所需移位寄存器的級數獲得偽隨機碼發射源的編碼信號； 獲得大地脈沖響應； 通過低通濾波器得到記錄系統的儀器響應； 獲得噪聲彳目號； 根據所述發射源信號、所述大地脈沖響應、所述儀器響應和所述噪聲信號獲得多道瞬 變電磁法的接收波形記錄。2. 如權利要求1所述的方法，其特征在于:根據所需移位寄存器的級數獲得偽隨機碼發 射源的編碼信號包括： 查表獲得偽隨機碼的特征多項式； 利用特征多項式求取序列多項式； 根據序列多項式計算出偽隨機碼的編碼； 根據偽隨機碼編碼獲得發射源信號。3. 如權利要求1所述的方法，其特征在于:獲得大地脈沖響應包括： 在頻率域計算大地的響應； 通過頻率域向時間域的變換得到大地的階躍響應； 通過計算階躍響應的導數得到大地的脈沖響應。4. 如權利要求1所述的方法，其特征在于，通過低通濾波器得到記錄系統的儀器響應包 括:通過計算海明窗口函數的低通濾波器響應得到記錄系統的儀器響應。5. 如權利要求1所述的方法，其特征在于:根據所述發射源信號、所述大地脈沖響應、所 述儀器響應和所述噪聲信號獲得多道瞬變電磁法的接收波形記錄包括： 利用如下公式計算獲得多道瞬變電磁法的接收波形記錄： ak(Xs，Xr，t) = i(k，Xs，t)*r(Xr，t)*g(Xs，Xr，t)+n(Xr，t) 其中，i (k，xs，t)表示發射源信號，r (xr，t)表示儀器響應，g (xs，xr，t)表示大地脈沖響 應，n(Xr，t)表示噪聲信號，k表示第k次數據采集，^為發射源位置，Xr為接收機位置，t為時 間。6. -種多道瞬變電磁法接收波形記錄合成裝置，其特征在于，包括： 發射源模塊，用于根據所需移位寄存器的級數獲得偽隨機碼發射源的編碼信號； 大地脈沖模塊，用于獲得大地脈沖響應； 儀器響應模塊，用于給通過低通濾波器得到記錄系統的儀器響應； 噪聲模塊，用于獲得噪聲信號； 合成模塊，用于根據所述發射源信號、所述大地脈沖響應、所述儀器響應和所述噪聲信 號獲得多道瞬變電磁法的接收波形記錄。7. 如權利要求6所述的裝置，其特征在于:所述發射源模塊所需移位寄存器的級數獲得 偽隨機碼發射源的編碼信號是指： 查表獲得偽隨機碼的特征多項式； 利用特征多項式求取序列多項式； 根據序列多項式計算出偽隨機碼的編碼； 根據偽隨機碼編碼獲得發射源信號。8. 如權利要求6所述的裝置，其特征在于:所述大地脈沖模塊獲得大地脈沖響應是指： 在頻率域計算大地的響應； 通過頻率域向時間域的變換得到大地的階躍響應； 通過計算階躍響應的導數得到大地的脈沖響應。9. 如權利要求6所述的裝置，其特征在于，所述儀器響應模塊通過低通濾波器得到記錄 系統的儀器響應是指： 通過計算海明窗口函數的低通濾波器響應得到記錄系統的儀器響應。10. 如權利要求6所述的裝置，其特征在于:所述合成模塊根據所述發射源信號、所述大 地脈沖響應、所述儀器響應和所述噪聲信號獲得多道瞬變電磁法的接收波形記錄是指： 利用如下公式計算獲得多道瞬變電磁法波場的系統響應信號： ak(Xs，Xr，t) = i(k，Xs，t)*r(Xr，t)*g(Xs，Xr，t)+n(Xr，t) 其中，i (k，xs，t)表示發射源信號，r (xr，t)表示儀器響應，g (xs，xr，t)表示大地脈沖響 應，n(Xr，t)表示噪聲信號，k表示第k次數據采集，^為發射源位置，Xr為接收機位置，t為時 間。
【文檔編號】G01V3/08GK105938205SQ201610125998
【公開日】2016年9月14日
【申請日】2016年3月4日
【發明人】王若, 薛國強, 王妙月, 底青云
【申請人】中國科學院地質與地球物理研究所