天然氣水合物相關的Slipstream海底滑坡體速度結構模型反演

藍坤 朱賀 何濤,? 梁前勇 吳學敏 董一飛 張毅

天然氣水合物相關的Slipstream海底滑坡體速度結構模型反演

藍坤1,2朱賀1,2何濤1,2,?梁前勇3吳學敏3董一飛3張毅4

1.造山帶與地殼演化教育部重點實驗室, 北京大學地球與空間科學學院, 北京 100871; 2.北京天然氣水合物國際研究中心, 北京 100871; 3.中國地質調查局廣州海洋地質調查局, 廣州 510760; 4.中國地質科學院, 北京 100037; ?通信作者, E-mail: taohe@pku.edu.cn

針對目前用于建立二維速度結構模型的 RAYINVR 軟件對四分量海底地震儀(OBS)記錄的轉換橫波無法自動反演, 建模效率低的問題, 采用 MATLAB 遺傳反演算法, 對 RAYINVR 軟件進行改進, 實現對橫波速度結構模型各層和層中區塊泊松比的自動同步反演, 可以為天然氣水合物勘探調查提供泊松比和楊氏模量等重要彈性參數。利用采自 Slipstream 海底滑坡的 OBS 數據, 通過同步反演方法, 獲得該地區較為精細的縱波和橫波速度結構模型, 并可與附近 U1326 鉆孔的測井數據進行對比, 驗證了橫波速度結構同步反演建模方法的有效性。模型揭示了兩個泊松比差異大的結構面: 代表水合物穩定底界的似海底反射界面(BSR)(海底之下230±10m)以及淺部異常高速體(有可能是高飽和度水合物富集的砂體)的底界(海底之下 75～100m), 后者與滑脫面大致重合, 指示水合物與Slipstream海底滑坡的形成有關。

同步反演; 橫波速度結構; Slipstream海底滑坡; 天然氣水合物; 射線追蹤

地震勘探是海域天然氣水合物調查的重要途徑[1?2], 不但可以對水合物儲層進行成像刻畫, 還可以建立速度結構模型, 從而有效地評估水合物資源量, 對后續開采意義重大。常規的海面拖纜地震勘探僅能產生和接收單一的縱波信號, 對沉積物孔隙流體中的氣體較為敏感, 少量的游離氣就可能造成大幅度的縱波能量衰減和速度降低[3], 影響水合物地層的成像效果和飽和度估算準確度。這種現象在作為水合物穩定區底界標志的地震似海底反射層(bottom simulating reflector, BSR)上下(表現為雜亂強反射)或攜帶大量游離氣向上運輸的冷泉通道(氣煙囪)內部(表現為條帶狀空白)尤為嚴重[4]。

橫波主要通過沉積物骨架傳播, 幾乎不受孔隙流體的影響[5]。因此, 在海域水合物重點勘探區, 通常使用四分量海底地震儀(Ocean Bottom Seismo-meter, OBS)來記錄由縱波在地層界面上轉換而來的橫波信號[6?7], 從而建立沉積地層的縱橫波速度結構模型, 以便獲取更加豐富、更加準確的水合物儲層信息。

加拿大西海岸溫哥華島外的北卡斯卡底(Casca-dia)陸緣增生楔主要由板塊俯沖過程中刮削下來的洋盆沉積物堆積而成[8?9], 蘊藏著豐富的天然氣水合物資源[10?11], 20 世紀 80 年代以來一直是水合物研究的重點區域, 不但具有豐富的單道地震(single channel seismic, SCS)、多道地震和 OBS 等地球物理勘探數據[12?13], 而且有專門針對水合物的大洋鉆探計劃(Ocean Drilling Program, ODP) 146 航次[14]和綜合大洋鉆探計劃(Integrated Ocean Drilling Program, IODP)X311 航次[15?17]的鉆井數據。特別之處是, 在該增生楔前緣的條狀變形脊上發育一系列海底滑塌構造。Lopez 等[18]和 Yelisetti 等[19]分別運用 OBS 數據和 RAYINVR 軟件[20], 研究其中規模較大的 Orca和 Slipstream 海底滑塌的地震速度結構, 認為其形成與水合物富集層的位置有關。然而, 受限于轉換橫波處理的復雜性和 RAYINVR 軟件的缺陷, 他們只建立滑塌構造的縱波速度模型, 缺乏橫波速度信息, 無法提供滑塌構造力學分析所需的楊氏模量和泊松比等彈性參數。

廣泛使用的 RAYINVR 軟件是用射線追蹤方法對二維速度模型進行走時擬合, 其內置的阻尼最小二乘(damped least squares)反演算法能夠方便地建立純縱(橫)波速度結構模型。但是, 對于 OBS 記錄的轉換橫波反射事件, 該軟件只能通過人工改變模型的泊松比參數, 將縱波速度轉換為相應的橫波速度, 并且, 一旦存在橫向速度結構不均的現象, 還需要將該結構層進一步劃分區塊, 模型會變得更加復雜。因此, 橫波模型需要花費大量的時間進行反復的手動測試, 效率過低, 嚴重妨礙 OBS 橫波信號的處理。

針對上述問題, 本研究在 RAYINVR 軟件的基礎上, 依托 MATLAB 語言建立 OBS 橫波模型的遺傳反演方法, 實現各結構層和層中區塊參數(泊松比)的同步反演。利用該方法建立的較精細的 Slips-tream 海底滑坡縱、橫波速度結構模型, 指示淺部高速異常體(海底之下 75～100m)和深部 BSR (海底之下 230±10m)兩個軟弱面為潛在的滑脫面位置。

1 RAYINVR軟件

對 OBS 等廣角地震數據, 通常采用射線路徑追蹤的算法進行速度結構層析反演, 采用正演試錯擬合的求解方案對模型參數進行迭代[21?22], 直至模型計算的走時與真實記錄的地震信號走時之間誤差足夠小。對于二維速度結構模型, 目前普遍使用 Zelt等[20,23]編寫的 RAYINVR 軟件, 其模型構建采用不規則四邊形網格, 通過靈活地設置速度以及界面節點的個數和位置, 使得模型層界面與地震事件識別情況一致, 能夠在很大程度上避免過參數化引起的反演假象以及欠參數化引起的細結構丟失問題。模型可以分為若干個不規則四邊形, 預先定義每個四邊形的頂點速度, 然后采用雙線性內插法, 由頂點速度得到四邊形內部的速度, 因此在四邊形內部的射線追蹤存在解析解。該軟件可以對折射波、反射波、界面滑行波、轉換波和層間反射波等進行射線追蹤計算, 具有運算速度快和算法穩定等優點。

當沒有轉換波時, RAYINVR 軟件可以使用自帶的阻尼最小二乘反演工具(DMPLSTSQR 程序), 從上往下逐層對模型節點速度和(或)層面深度參數進行自動的快速的迭代反演, 每一層的射線走時擬合流程見圖 1。以建立縱波速度結構模型為例, 首先用 RAYINVR 軟件正演計算某層的射線路徑和走時, 需要輸入 v.in, tx.in 和 r.in 這 3 個文件。其中, 文件 v.in 包含模型每一層的深度和速度信息, 文件tx.in 是從地震成像數據拾取的真實地震事件走時, 文件r.in則是關于射線追蹤的算法控制以及結果畫圖的參數設置。RAYINVR 正演會輸出射線追蹤走時與實際地震事件走時殘差的均方根(root mean squares, RMS)和標準卡方值(normalized2), 后者的計算公式為

其中,為實際拾取的地震走時數量,T為第個實際拾取的地震走時,t為相應的射線追蹤得到的地震走時,為第個地震走時的拾取誤差。在得到標準卡方值后, 通過人機交互判斷2是否滿足要求: 當2≤1 時, 說明走時殘差達到地震事件拾取誤差水平, 可以結束反演; 如果2顯著大于 1, 則將該值輸入 DMPLSTSQR 程序, 同時輸入反演參數控制文件d.in (在該文件中選擇需要反演的目標界面), 通過偏微分方法求出模型參數調整向量, 據此更新模型文件 v.in 中的速度和(或)深度參數。迭代上述正演和反演過程, 直到獲得模型最優解(保存在最后一步的文件 v.in 中)。

可見, DMPLSTSQR 程序是通過改變文件 v.in 中的速度和地層界面深度參數來實現迭代反演的。然而, 當對轉換橫波進行射線追蹤時, 縱波轉換為橫波的反射界面參數(如轉換發生的層號和泊松比等)均存儲在文件 r.in 中, 無法由 DMPLSTSQR 程序自動更新, 只能手動修改。因此, 橫波速度模型的迭代反演過程非常耗時, 并且, 當研究區的地質構造比較復雜時, 需要對目標層進一步細化分塊, 以至很難手動將模型參數調節到令人滿意的程度。

圖1 RAYINVR軟件迭代反演速度模型流程 (修改自文獻[20])

2 轉換橫波同步反演方法

根據震動形式, 可將地震波分為體波和面波, 其中體波又分為縱波和橫波。當縱波到達一個反射界面時, 其反射波中不僅包含縱波, 還包含轉換而來的橫波(圖 2)。由于海面震源激發的地震波在海水中只能傳播縱波, 因此放置在海底的 OBS 記錄的主要橫波事件(能量最強, 可以識別)是由縱波在其反射界面上轉換產生。因此, RAYINVR 對 OBS 橫波的射線追蹤, 不是直接確定橫波速度, 而是在給定的縱波速度結構模型(地層界面深度和縱波速度梯度保持固定)基礎上, 尋找最優的地層泊松比。在通過 RAYINVR 軟件得到模型的泊松比之后, 橫波速度S則由對縱波速度P的變換得到:

如果研究區域的地質構造較為簡單, 可以通過手動調節的方法, 給模型各層設置泊松比參數。但是, 當地質構造比較復雜時, 地層的泊松比不僅在縱向上可能存在差異, 在橫向上也可能存在較大的差異, 需要在層內進一步劃分區塊。此時, 用人工調節參數的方法就會事倍功半, 因此需要研發一套能夠讓程序實現自動同步反演而得到最優解的建模方法。

1為P波在反射界面的入射角,2為P波的折射角,1為轉換S波的反射角,2為轉換S波的折射角

圖2 縱波在彈性界面發生反射和折射示意圖

Fig. 2 Reflection and refraction diagram of P wave incident on an elastic interface

遺傳算法通過模擬自然進化過程, 對值域空間進行搜索, 得到全局最優解, 是一種可用于復雜系統優化的算法。本文將 RAYINVR 軟件的射線追蹤法與 MATLAB 軟件的遺傳算法相結合, 用 MATLAB語言編寫一套能夠實現 OBS 橫波速度結構模型在各層以及層間不同區塊同步反演的方法, 從而節省人工修改參數的時間成本, 實現復雜地質構造環境下反演橫波速度結構的自動化。該同步反演方法的基本流程見圖 3。

1)給定模型各層和層中區塊的初始泊松比組合(保存在文件 r.in 中)。

2)設定遺傳算法的最大代數和每代的種群數量, 啟動自動反演。

3)RAYINVR 軟件讀入射線追蹤控制文件 r.in、拾取的地震走時文件 tx.in 以及速度和層界面深度控制文件 v.in, 通過正演射線路徑擬合得到2。

圖3 OBS橫波速度模型同步反演流程

4)遺傳算法會判斷此時是否達到預設條件: 如果達到預設的迭代次數, 或者2≤1 (意味著此時的泊松比為全局最優解), 反演終止; 如果沒有達到預設的迭代次數, 并且2>1, 那么程序會生成一組新的泊松比組合, 并自動傳給 RAYINVR 軟件, 進行新一輪正演。

5)重復上述迭代過程, 直到滿足預設條件, 得到的最新泊松比組合將保存在文件 r.in 中。

3 Slipstream海底滑坡區地質背景

在卡斯卡底活動陸緣, 胡安德富卡(Juan de Fuca)板塊以約 46mm/a 的速率俯沖到北美板塊之下, 來自洋盆的大量沉積物在變形前緣位置被刮削下來, 堆積形成增生楔[8]。這些增生體沿著一些巨大的深達洋殼的逆沖斷層發生擠壓和縮短, 并相互重疊, 從而在陸坡底部(變形前緣)形成一系列與板塊邊界近平行的長條狀背斜山脊[19], 本文的研究對象 Slip-stream 山脊就是其中之一(圖 4)。卡斯卡底活動陸緣一直是天然氣水合物研究的聚焦區域, 因為大量地震調查都觀測到從變形前緣到陸坡中部的深水沉積物中廣泛存在作為水合物穩定區底界地震標志的BSR, 海底之下的深度約為 200～300m, 顯示大量天然氣水合物的富集[5,24?26,28], 也得到 ODP146 和 IODP X311 大洋鉆探結果的證實。

近年來, 在高分辨率多波束海底地貌的調查中發現, 變形前緣的條狀背斜山脊上經常發育滑塌構造。例如, IODP X311 U1326 站位所在的山脊發育一個規模較大的 Orca 海底滑塌(圖 4), 其滑脫面深度恰好抵達 BSR, 因此推測該滑塌構造的形成與水合物層有關[18]。作為對比, 其相鄰的 Slipstream 滑塌(U1326 站位東南方向約 15km)規模更大, 推測其形成也與水合物有關, 但滑脫面的深度卻只有 100m左右[19]。

根據美國華盛頓大學 TN175 航次采集的 EM300多波束測深數據[19], Slipstream 海底滑坡體的長度約為 3.3km, 寬度約為 2.5km, 滑坡角度較大(約為25°), 出露于海底的陡壁高度可達近 100m。該滑坡體上發育一系列正斷層, 且部分斷層面出露于海底, 因此滑坡體在地貌上類似多個排列不整齊的梯形箱狀高地。較大的滑坡塊體主要分布在距離山腳 1～2km 的海底平原沉積層之上, 一些粒度較小的滑坡體沉積物則外流到更遠的地區(圖 4)。本研究采用的 SCS 和 OBS 數據來自 SeaJade (Seafloor Earthquake Array – Japan Canada Cascadia Experiment)項目。

修改自文獻[27]。紅色虛線框為 Slipstream 海底滑坡區域; 紅色三角形為 OBS 測點 30, 31, 33 和 34 的位置; 4 號地震測線平行地穿過Slipstream 山脊, 其中綠線為 OBS 地震航次, 黃線為 SCS 航次; 紅色五角星為 IODP X311 航次的U1326 站位; 多波束海底地貌圖上疊加的等深線間隔為 25 m

4 縱波速度結構模型

由于 OBS 記錄的橫波是由縱波在地層界面上轉換而來, 因此需要建立 Slipstream 海底滑坡區域的縱波速度結構模型。

傳統的海面拖纜 SCS 數據能夠很好地針對海底沉積層的結構特點進行成像, 但受限于炮點與接收器之間的偏移距(source-receiver offset)在深海中太小, 無法得到準確的速度和深度信息。布設于海底的 OBS 采集的是廣角(大偏移距)數據, 能夠很好地對海底沉積層的速度信息進行分析, 但是對其結構信息的表征則有所欠缺。因此, 在構建海底沉積層縱波速度結構模型時, 我們同時采用 OBS 廣角數據和 SCS 數據來實現對地層速度和深度的同步約束。需要注意的是, RAYINVR 軟件最初是用于大陸勘探的, 設計的是單震源對應多個接收器, 而海洋的地震勘探則是單個 OBS 記錄多個炮點數據, 因此在建立模型時, 需要將 OBS 作為射線發射點, 將炮點作為接收位置。

前人在卡斯卡底增生楔區域創建的縱波模型一般為 5～7 層[5,19,25], 用于滑塌區力學機制分析時精細度不夠。我們通過仔細地辨識 OBS 和 SCS 的反射和界面滑行波事件(圖 5), 在 SCS 數據中拾取 10 個縱波反射事件(line4A_r1～line4A_r10), 在 OBS 數據中拾取 7 個縱波反射事件(RL0～RL6)和 5 個縱波臨界折射(界面滑行波)事件(RF1～RF5)。

(a)在 4 號地震測線的 SCS 數據中拾取的共 10 個縱波反射事件(line4A_r1～line4A_r10), 其中 r1 對應海底, r2 和 r3 分別對應淺部高速層的頂、底界面, r5 對應BSR; (b)在 OBS 測點 33 的數據中拾取的廣角地震反射(RL)和臨界折射事件(RF), 其中 RL0, RL1 和RL2 分別對應海底、高速層頂部和 BSR, RF1 對應高速層頂界面, RF2～RF5 為 BSR 之下的深部臨界折射事件, 下同。折合走時=走時?偏移距/折合速度, 本文折合速度為 4000 m/s

綜合 SCS 和 OBS 數據的拾取事件, 最終將 Slip-stream 海底滑坡沿 4 號測線方向的山脊剖面分為 14層(不包含海水)。反演得到的最終射線路徑追蹤結果見圖 6(a), 擬合得到的走時與實際拾取的走時相符(圖 6(b)), 模型整體的 RMS 走時殘差為 15ms, 小于拾取誤差(20ms),2=0.901, 達到預期要求。

精細化的縱波速度結構模型實現對最小厚度為25m 左右薄層的識別(圖 7(a)), 其中 L1 與海水直接接觸, 平均速度約為 1.51km/s, 代表非常松散的海底沉積物; L2 的頂部界面在海底之下約 75m, 厚度約為 25m, 縱波速度突然增加至約 2.0km/s, 推測可能是有高飽和度天然氣水合物富集的砂體; L3 ～L5, 縱波速度從 1.7km/s 緩慢增加至 2.1km/s, 代表無(或含少量)水合物沉積物的壓實趨勢; L5 的底部界面對應 BSR, 深度約在海底之下 230m; BSR 之下的L6 縱波速度略有削減, 代表游離氣體的存在; 之后, 縱波速度的增加再次向壓實趨勢線回歸。

(a)縱波速度結構模型擬合的射線路徑, 4個圓圈為OBS測點位置及編號(下同); (b)縱波速度結構模型的走時擬合結果

5 橫波速度結構模型

當沉積物孔隙中含游離氣體時, 不但會影響SCS 反射數據和 OBS 垂直分量的成像數據, 衰減嚴重時還會出現地震反射空白, 大幅度地降低縱波速度, 影響水合物飽和度的計算。然而, 橫波傳播不受孔隙流體影響, 可以反映含氣區域(如 BSR 之下的游離氣層)的內部結構。并且, 橫波速度在 BSR之上的大幅度增加可能顯示固體天然氣水合物的存在。此外, 由于 OBS 記錄的橫波是由縱波轉換而來, 其頻率與縱波相同, 但速度更慢(橫波在純天然氣水合物中的速度約為 1890m/s, 在不含水合物的淺層海底沉積物中的速度為 100～600m/s)[29], 因此波長更短, 從而可以解析比縱波更精細的沉積結構信息。

在上述縱波速度結構模型的基礎上, 保持各層頂、底界面的深度節點和縱波速度節點不變, 通過MATLAB 的遺傳算法對橫波模型的泊松比進行同步反演。模型初始的參考泊松比數值(圖 8(c))通過15 km 之外的 IODP U1326D 鉆孔的縱、橫波測井結果計算得到:

需要指出的是, 通過非線性的遺傳算法得到的滿足條件的泊松比組合是全局最優解。作為對比, RAYINVR 的 DMPLSTSQR 程序是線性優化, 在初始參數與模型最終值偏差較大或地質構造復雜時, 容易掉進局部解的陷阱, 造成反演結果不佳, 需要重新手動調節初始參數再嘗試。圖 9 展示轉換橫波的射線追蹤路徑和走時擬合結果, 得到的 RMS 走時殘差為 27ms, 小于橫波事件的拾取誤差(30ms),2=0.903, 達到預期要求。

將橫波速度結構模型(圖 7(b))與縱波速度結構模型(圖 7(a))進行對比, 可見淺部的高速異常體 L2的橫波速度增加幅度(14.3%)明顯小于縱波(24.5%), 推測該層的高飽和度天然氣水合物很可能以固體填隙物的形式富集, 對沉積物骨架的力學性能增加不明顯(泊松比與其上接觸海水的 L1 相近)。BSR 之上的 L5 橫波速度相對于其下游離氣層的增加幅度(14.0%)明顯大于縱波(5.3%), 說明縱波速度受到游離氣的干擾, 沒有顯示出含水合物層的高速特征。圖 8(a)和(b)以 OBS 測點 34 所在位置的一維深度?速度剖面形式對比上述縱、橫波速度的變化情況, 圖8(c)的泊松比則顯示淺部高速層(L1)與 BSR 之下游離氣層相比, 泊松比明顯偏高, 指示這兩個部位為力學軟弱層。L2 的泊松比大于游離氣層, 可能是Slipstream滑塌體沿著相對更軟的高速層底界發育的原因。作為對比, Orca 滑塌體旁邊的 U1326D 測井數據顯示淺部的高速異常體出現在海底之下 60～ 90m 處, 而 Slipstream 海底滑坡區的高速層是在海底之下 75～100m, 這可能是 U1326D 與 Slipstream滑坡區之間約 15km 的空間距離引起的差異。此外, U1326D 的淺部高速異常體泊松比為低值, 指示力學強度相對較高, 并不是一個軟弱面, 因此其滑脫面出現在泊松比值對比更強烈的 BSR 位置(應力容易在力學性質相差較大的分界部位集中), 顯示與Slipstream 滑塌體不同的水合物控制作用。

圖8 Slipstream海底滑坡區縱橫波速度結構模型與U1326D測井數據[30]對比

6 結論

1)將 RAYINVR 軟件射線追蹤算法與 MATLAB遺傳反演相結合, 實現 OBS 橫波速度結構模型各層和層中區塊泊松比的同步反演, 極大地提高了地質構造復雜區域的橫波速度結構建模效率。

2)建立 Slipstream 海底滑坡體的精細縱、橫波速度結構模型, 顯示作為水合物穩定區底界地震標志的似海底反射層(BSR)位于海底之下 230～250m, 而淺部約 75～100m 的高速異常體很可能是高飽和度天然氣水合物富集的砂體。模型結果可與 IODP X311 航次 U1326D 測井數據進行對比, 驗證了同步反演速度結構建模的有效性。

3)淺部高速異常體的深度與滑坡面大致重合, 顯示水合物造成的力學軟弱面對滑坡構造的控制作用, 為后續建立三維力學模型提供了關鍵證據。

(a) OBS 橫波速度結構模型擬合的射線路徑, S1～S8代表轉換橫波反射事件; (b) OBS 橫波速度結構模型的走時擬合結果

致謝 衷心感謝加拿大維多利亞大學 George Spence 教授提供 OBS, SCS 和多波束等數據, 特別感謝 IHS 全球公司提供地球科學綜合解釋軟件包“The Kingdom Suite”。

[1] Makogon Y F. Natural gas hydrates — a promising source of energy. Journal of Natural Gas Science and Engineering, 2010, 2(1): 49?59

[2] Helgerud M B, Dvorkin J, Nur A, et al. Elastic-wave velocity in marine sediments with gas hydrates: effective medium modeling. Geophysical Research Letters, 1999, 26(13): 2021?2024

[3] Yang J, Davies R J. Gravity-driven faults: migration pathways for recycling gas after the dissociation of marine gas hydrates. Marine Geology, 2013, 336: 1?9

[4] Davies R J, Clarke A L. Methane recycling between hydrate and critically pressured stratigraphic traps, offshore Mauritania. Geology, 2010, 38(11): 963?966

[5] Dash R, Spence G. P-wave and S-wave velocity struc-ture of northern Cascadia margin gas hydrates. Geo-physical Journal International, 2011, 187(3): 1363? 1377

[6] 張光學, 徐華寧, 劉學偉, 等. 三維地震與OBS聯合勘探揭示的神狐海域含水合物地層聲波速度特征. 地球物理學報, 2014, 57(4): 1169?1176

[7] 劉伊克, 朱偉林, 米立軍, 等. 南海深水多次波成像. 中國科學: 地球科學, 2015, 45(2): 152?160

[8] Davis E E, Hyndman R D. Accretion and recent defor-mation of sediments along the northern Cascadia sub-duction zone. Geological Society of America Bulle-tin, 1989, 101(11): 1465?1480

[9] Hyndman R D. The Lithoprobe corridor across the Vancouver Island continental margin: the structural and tectonic consequences of subduction. Canadian Journal of Earth Sciences, 1995, 32(10): 1777?1802

[10] Mackay M E, Jarrard R D, Westbrook G K, et al. Origin of bottom-simulating reflectors: geophysical evidence from the Cascadia accretionary prism. Geo-logy, 1994, 22(5): 459

[11] Pohlman J W, Kaneko M, Heuer V B, et al. Methane sources and production in the northern Cascadia mar-gin gas hydrate system. Earth and Planetary Science Letters, 2009, 287(3/4): 504?512

[12] Riedel M, Tréhu A, Spence G. Characterizing the thermal regime of cold vents at the northern Cascadia margin from bottom-simulating reflector distributions, heat-probe measurements and borehole temperature data. Marine Geophysical Researches, 2010, 31(1/2): 1?16

[13] Hyndman R D, Riedel M, Spence G D. Hydrate studies of northern Cascadia margin off Vancouver Island: a reference source // 6th International Con-ference on Gas Hydrates (ICGH 2008). Vancouver, 2008, doi: 10.14288/1.0040999

[14] Westbrook G K, Carson B, Musgrave R J. Proceeding of Ocean Drilling Program, volumn 146 part 1: initial reports [R]. Texas: Ocean Drilling Program, Texas A&M University, 1994

[15] Torres M E, Tréhu A M, Cespedes N, et al. Methane hydrate formation in turbidite sediments of northern Cascadia, IODP Expedition 311. Earth and Planetary Science Letters, 2008, 271: 170?180

[16] Riedel M, Collett T S, Malone M J. Gas hydrate drilling transect across northern Cascadia margin — IODP Expedition 311. Geological Society London Special Publications, 2009, 319(1): 11?19

[17] Expedition 311 Scientists. Cascadia margin gas hyd-rates [R]. Washington, DC: Integrated Ocean Drilling Program Management International, 2005

[18] Lopez C, Spence G, Hyndman R, et al. Frontal ridge slope failure at the northern Cascadia margin: margin-normal fault and gas hydrate control. Geology, 2010, 38(11): 967?970

[19] Yelisetti S, Spence G D, Riedel M. Role of gas hydrates in slope failure on frontal ridge of northern Cascadia margin. Geophysical Journal International, 2014, 199(1): 441?458

[20] Zelt C A, Smith R B. Seismic traveltime inversion for 2-D crustal velocity structure. Geophysical Journal International, 1992, 108(1): 16?34

[21] Catchings R D, Mooney W D. Crustal structure of the Columbia Plateau: evidence for continental rifting. Journal of Geophysical Research: Solid Earth, 1988, 93(B1): 459?474

[22] Boland A V, Ellis R M. Velocity structure of the Kapuskasing Uplift, northern Ontario, from seismic refraction studies. Journal of Geophysical Research: Solid Earth, 1989, 94(B6): 7189?7204

[23] Zelt C A, Ellis R M. Practical and efficient ray tracing in two-dimensional media for rapid traveltime and amplitude forward modeling. Canadian Journal of Exploration Geophysics, 1988, 24(1): 16?31

[24] Hyndman R D, Spence G D. A seismic study of methane hydrate marine bottom simulating reflectors. Journal of Geophysical Research: Solid Earth, 1992, 97(B5): 6683?6698

[25] Chen M. Northern Cascadia marine gas hydrate: con-straints from resistivity, velocity, and AVO [D]. Vic-toria: University of Victoria, 2006

[26] Lopez C. Seismic velocity structure associated with gas hydrate at the frontal ridge of Northern Cascadia Margin [D]. Victoria: University of Victoria, 2007

[27] 隆松伯, 何濤, 梁前勇, 等. 基于ANSYS Work-bench的高精度自動化三維地質建模方法——以天然氣水合物相關的Slipstream海底滑坡為例. 北京大學學報(自然科學版), 2018, 54(5): 994?1002

[28] Hyndman R D, Davis E E. A mechanism for the formation of methane hydrate and seafloor bottom-simulating reflectors by vertical fluid expulsion. Journal of Geophysical Research: Solid Earth, 1992, 97(B5): 7025?7041

[29] Waite W F, Helgerud M B, Nur A, et al. Labora- tory measurements of compressional and shear wave speeds through methane hydrate. Annals of the New York Academy of Sciences, 1999, 912(1): 1003?1010

[30] Riedel M, Collett T S, Malone M J. Expedition 311 summary [R]. Washington, DC: Integrated Ocean Drilling Program Management International, 2006

Model Inversion of Velocity Structure for Slipstream Submarine Slide Related to Gas Hydrate

LAN Kun1,2, ZHU He1,2, HE Tao1,2,?, LIANG Qianyong3, WU Xuemin3, DONG Yifei3, ZHANG Yi4

1. Key Laboratory of Orogenic Belts and Crustal Evolution (MOE), School of Earth and Space Sciences, Peking University, Beijing 100871; 2. Beijing International Center for Gas Hydrate, Peking University, Beijing 100871; 3. Guangzhou Marine Geological Survey, China Geological Survey, Guangzhou 510760; 4. Chinese Academy of Geological Sciences, Beijing 100037; ? Corresponding author, E-mail: taohe@pku.edu.cn

The RAYINVR software, which is widely used in academia for 2-D velocity structure model from four-component OBS (ocean bottom seismometer) data, is unable to invert for the converted S-wave automatically, resulting in the low-efficiency of modeling process. Using MATLAB’s genetic algorithm, the RAYINVR software is improved and able to automatically and synchronously invert for Poisson’s ratios of each layer with all sub-blocks for the S-wave velocity structure model, and thus can provide Young’s modulus, Poisson’s ratio and other important mechanical information for gas hydrate survey. This method is applied to process the OBS data collected at the Slipstream submarine slide, and a fine P- and S-wave velocity structure model is obtained, which is comparable to the logging data of nearby borehole U1326. Therefore, the validity of the auto-synchronous inversion method is verified for the S-wave velocity structure modeling. The optimal velocity model reveals two structural interfaces with large Poisson’s ratio contrast. One is BSR (bottom simulating reflector) at 230±10 mbsf (meter beneath sea floor), which represents the bottom boundary of the gas hydrate stability zone, and the other is the basal boundary of a shallow abnormal high-speed body (possibly a sand body enriched with high saturation gas hydrate) at 75?100 mbsf. The latter agrees roughly with the glide plane of Slipstream submarine slide, indicating that the hydrate is related to the formation of submarine landslide.

synchronous inversion; S-wave velocity structure; Slipstream submarine slide; gas hydrate; ray tracing

10.13209/j.0479-8023.2021.005

國家自然科學基金(41676032)和中國地質調查局國家天然氣水合物專項基金(DD20190234, DD20190218, DD20189320, DD20160217, HD-JJHT-20)資助

2020–03–10;

2020–05–15