JP2008014830A - ハイドレートの存在領域探査方法及び探査システム - Google Patents
ハイドレートの存在領域探査方法及び探査システム Download PDFInfo
- Publication number
- JP2008014830A JP2008014830A JP2006187073A JP2006187073A JP2008014830A JP 2008014830 A JP2008014830 A JP 2008014830A JP 2006187073 A JP2006187073 A JP 2006187073A JP 2006187073 A JP2006187073 A JP 2006187073A JP 2008014830 A JP2008014830 A JP 2008014830A
- Authority
- JP
- Japan
- Prior art keywords
- hydrate
- spectrum
- exploration
- frequency band
- seismic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000001228 spectrum Methods 0.000 claims abstract description 55
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims description 29
- 230000015572 biosynthetic process Effects 0.000 claims description 23
- 238000000605 extraction Methods 0.000 claims description 14
- 238000006116 polymerization reaction Methods 0.000 claims description 13
- 150000004677 hydrates Chemical class 0.000 claims description 2
- NMJORVOYSJLJGU-UHFFFAOYSA-N methane clathrate Chemical compound C.C.C.C.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O.O NMJORVOYSJLJGU-UHFFFAOYSA-N 0.000 abstract description 69
- 238000005755 formation reaction Methods 0.000 description 22
- 230000005540 biological transmission Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 150000004678 hydrides Chemical class 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 3
- 239000013535 sea water Substances 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 239000000956 alloy Substances 0.000 description 2
- 229910045601 alloy Inorganic materials 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- VUZPPFZMUPKLLV-UHFFFAOYSA-N methane;hydrate Chemical class C.O VUZPPFZMUPKLLV-UHFFFAOYSA-N 0.000 description 1
- 238000001028 reflection method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000003381 stabilizer Substances 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B41/00—Equipment or details not covered by groups E21B15/00 - E21B40/00
- E21B41/0099—Equipment or details not covered by groups E21B15/00 - E21B40/00 specially adapted for drilling for or production of natural hydrate or clathrate gas reservoirs; Drilling through or monitoring of formations containing gas hydrates or clathrates
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
【課題】海底疑似反射面が存在する海底の地層において、メタンハイドレート存在領域を探査できるハイドレートの存在領域探査方法及び探査システムを提供する。
【解決手段】広周波数帯域を有する震源を用いた地震探査の計測結果から低周波数帯域を用いて海底の地層構造図を作成し、この地層構造図から特定された海底疑似反射面の場所において、海底疑似反射面より上方の計測結果の信号のスペクトルを求め、探査の全周波数に渡ってスペクトルが大きく、かつ、高周波数帯域のスペクトル強度が低周波数帯域のスペクトル強度より大きい部位を抽出し、抽出部位が無ければハイドレートは存在しないと判定し、抽出部位があればこの抽出部位にハイドレートが存在すると判定する。
【選択図】図2
【解決手段】広周波数帯域を有する震源を用いた地震探査の計測結果から低周波数帯域を用いて海底の地層構造図を作成し、この地層構造図から特定された海底疑似反射面の場所において、海底疑似反射面より上方の計測結果の信号のスペクトルを求め、探査の全周波数に渡ってスペクトルが大きく、かつ、高周波数帯域のスペクトル強度が低周波数帯域のスペクトル強度より大きい部位を抽出し、抽出部位が無ければハイドレートは存在しないと判定し、抽出部位があればこの抽出部位にハイドレートが存在すると判定する。
【選択図】図2
Description
本発明は、海底の地層内等にあるハイドレートの存在領域探査方法及び探査システムに関し、より詳細には、海底疑似反射面が存在する海底の地層において、ハイドレードの存在領域を検出できるハイドレートの存在領域探査方法及び探査システムに関する。
日本周辺の海底にはエネルギー資源となりうるメタンハイドレート等のハイドレートが相当量存在しており、これを有効利用することが望まれている。そのためには、これらのハイドレートの資源の探査が必要となる。
この探査手段の一つに反射法地震探査がある。この反射法地震探査では、海面付近や海中や海底で、圧縮した空気を瞬時に海中に放出するエアガンや圧電素子や超磁歪合金等の振動源(送波器)で、振動(音波)を発生させ、海底面や海底下の地層の境界面等で反射して戻って来る音波を、ストリーマーとよばれるケーブルに間隔を設けて配置されたハイドロホン(受波器)で受信し、この受信記録、即ち、海底面下の地震探査記録に現れる物性が異なる地層同士の境界面を探り、地質構造を解析している。
この音波は、海中を伝わり、海中を伝わった音波は海底面で反射して、ある程度の音波が受波器まで戻って来る。しかし、海底面の下の地層まで入っていく音波もあり、この地層の中に入った音波は地層同士の境界面でも反射して受波器まで戻って来る。海底面下の地震探査記録では、物性が異なる地層同士の境界面は黒い線で現れ、この境界面は地質構造を示す。この反射してきた音波を使って探査をするので反射法地震探査と呼ぶ。
海底面下、どの程度の深さまで音波が到達し、受波器に戻ってくるかは、使用する機器によって異なるが、石油探鉱のための地震探査では、海底面下10,000m程度までの記録をとることもある。なお、最近の石油探鉱地震探査では、ストリーマーケーブルも6,000m以上となったり、複数本のストリーマーケーブルを同時に曳航して調査する場合もある。このストリーマーケーブルを複数本使用する場合には、立体的な地震探査記録を得ることができるので、3次元地震探査という。
この反射法地震探査に関して、海上探査では、距離が長く音波が減衰するため、100Hz程度の低周波数を使用するので、反射波の分解能が低く高精度の地層探査が困難となるため、海底近傍を航走させる自律型無人航走体から振動波を発生し、この振動波を、この自律型無人航走体が曳航しているストリーマケーブルにより受信する海底地層探査システムが提案されている(例えば、特許文献1参照。)。
そして、特に、メタンハイドレートを対象とした反射法地震探査では、エアガン等を用いた低周波(例えば、100Hz)による探査でハイドレート層の存在を示唆する海底疑似反射面を確認できるが、深海曳航体などによる高周波(例えば、600Hz)による探査では明確に確認できない場合があり、また、海底疑似反射面(BSR:Bottom Simulated Reflector)の直上に、振幅が小さくなるブランキング現象が確認される等の特徴が知られている。
この海底下のハイドレート等の存在量を推定するために、正確な資源量評価手法が要求されているが、従来の評価手法では、地震探査記録に現れる海底疑似反射面の直上にハイドレートが存在すると考えられており、地震探査記録に現れる海底疑似反射面を検出することにより、ハイドレート存在領域を特定している。つまり、ハイドレートの存在指標とされている海底疑似反射面の分布に基づいて、地震探査記録上に現れる海底疑似反射面をハイドレートの下限とし、また、上限は、検層データ、地震探査結果の解釈(ブランキング、砂岩層分布)、精密重力探査等から推定し、ハイドレート層の層厚推定を行っている。
また、本発明者らは、特願2004−327728号において、300Hz〜1KHzの高周波数で探査して得た地震探査記録上に現れる短い間隔での反射波群の内、その一番早い時間の反射波をハイドレート層の上限とすることにより、従来技術では探査できなかったハイドレート層の上面を決めて、ハイドレート層の層厚を計測できるハイドレート層の層厚推定方法を提案した。
このハイドレート層の層厚推定方法では、薄いハイドレート層と薄い地層の互層構造を対象としているため、薄層互層構造の場合には高周波数帯域において境界面からの反射波が明瞭に観測されるので極めて有効であるが、適用できる範囲が薄層互層構造に限定されるという問題がある。
しかしながら、地震探査記録にはハイドレート層の上面からの明瞭な反射波が観測されない場合が多く、また、海底疑似反射面が存在する複数の地点において試掘したところ、海底疑似反射面が存在してもハイドレート層が存在しない場合もあることが明らかになってきた。このため、このような場合にも適用できる、より汎用的なハイドレートの存在領域探査方法の開発が望まれている。
現状では、ハイドレート層の構造が解明されていないこともあり、有効な探査法は確立されていない。しかし、ハイドレート層を対象とした地震探査の結果から、低周波数の地震探査では海底疑似反射面が観測されるが、高周波数の地震探査では、海底疑似反射面が観測されず、また、どの周波数においてもハイドレート層の上面を特定できないことが分かってきている。
これらの結果から、本発明者は、次のような考察を行った。ハイドレート層が均質であれば、地層とハイドレート層内の音速が異なることからハイドレート層の上面から反射波が生じるが、反射波からハイドレート層の上面の検出ができないので、ハイドレート層は均質ではないと考えられる。
また、低周波数の地震波により、ハイドレート層の下面に存在するフリ−ガス又は水からの反射波である海底疑似反射面を検出できるので、低周波数の地震波はハイドレート層からは大きな影響を受けないもと考えられる。また高周波数の地震波は、フリーガス層まで到達しないか、フリーガス層からの反射波が受信点まで到達していないので、ハイドレート層から影響を大きく受けるものと考えられる。
そして、ハイドレート層が多数の散在する小ハイドレート塊で構成されていると仮定すると、低周波数の地震波はこれらのハイドレート塊では反射され難く、高周波数の地震波はこれらのハイドレート塊で反射され易いので、低周波数の地震波ではフリーガス層まで容易に到達できて、海底疑似反射面を観測することができるが、高周波数の地震波ではフリーガス層まで到達できず海底疑似反射面が観測できないと言うことを説明できる。
また、このハイドレート層の構成の場合には、高周波数の地震波を用いても、個々のハイドレート塊からの反射は微小であり、明瞭な反射波としては観測できないので、ハイドレート層の上面を特定するのは非常に難しいと考えられる。
特開2003−19999号公報
本発明は、上記の考察を基に、上記の問題を解決するためになされたものであり、その目的は、海底疑似反射面が存在する海底の地層において、ハイドレート存在領域を検出できるハイドレートの存在領域探査方法及び探査システムを提供することにある。
上記の目的を達成するための本発明のハイドレートの存在領域探査方法は、反射法地震探査方法によって海底下のハイドレードを探査するためのハイドレートの存在領域探査方法であり、広周波数帯域を有する震源を用いて地震探査を行うステップと、該地震探査の計測結果から低周波数帯域を用いて海底の地層構造図を作成し、該地層構造図から海底疑似反射面の存在領域を特定するステップと、該海底疑似反射面が存在する場所において、広周波数帯域におけるスペクトルを求めるステップと、前記スペクトルにおいて、探査の全周波数に渡ってスペクトルが大きく、かつ、高周波数帯域のスペクトル強度が低周波数帯域のスペクトル強度より大きい部位を抽出するステップと、前記抽出部位が無ければ、メタンハイドレートは存在しないとするステップと、前記抽出部位があれば該抽出部位にメタンハイドレートが存在すると判定するステップとを有してなることを特徴とする。
この広周波数帯域とは、100Hz〜800Hz、好ましくは100Hz〜1000Hz程度の範囲内の周波数幅が700Hz〜900Hzの周波数帯域のことを言い、低周波帯域とは、100Hz〜200Hz程度の範囲に入る周波数帯域のことを言う。また、高周波帯域とは、低周波帯域よりも高い周波数帯域のことを言う。
また、海底疑似反射面とは、BSR(Bottom Simulated Reflector)とも呼ばれる、周辺の順序良く並んだ地層境界面を示す記録を無視して伸びる音波境界面のことで、海底面を示す記録にほぼ平行で、反射の性質が異なる(位相が逆転している)、反射の強度が強い(記録としては太い線として見える)もののことを言う。
このスペクトルは、高速フーリエ変換(FFT)法や最大エントロピー(MEM)法、ウェーブレット(WT)など任意の周波数分析法を適用して求めることができる。このスペクトルは、海底疑似反射面が存在する場所で、できるだけ広い範囲の周波数帯域のデータ、好ましくは計測した全ての周波数帯域のデータを用いて計算する。
また、探査の全周波数に渡ってスペクトルが大きいとは、周囲のハイドレートが存在しない領域におけるスペクトルと比較によって判定する。このハイドレートが存在しない領域は、ハイドレートが存在するための温度範囲、圧力範囲が分かっているため、ハイドレートが存在しないことは高い精度で判定できる。
また、上記のハイドレートの存在領域探査方法において、前記広周波数帯域におけるスペクトルを求める際に、共通反射点重合法を用いると、計算精度を高めるができる。この共通反射点重合法はCDP(Common Depth Point)重合法とも呼ばれ、経路の異なった反射波をその共通反射点(CDP)位置での垂直走時に変換(NMO:Normal Move out 補正)して重合(加算)し、反射点が共通となる距離の異なる反射記録の到達時間を補正し足し合わせて、反射波を強調する方法である。
さらには、上記のハイドレートの存在領域探査方法において、前記スペクトルを深さ方向に移動させる解析窓毎に求め、スペクトルに変化が生じ始める解析窓がある深さをハイドレート層の上面とするステップと、前記海底疑似反射面をハイドレート層の下面とするステップとを有して構成すると、ハイドレート層の上下面を特定できるようになる。この切り出し部位である解析窓は、深さ、言い換えれば、反射波が受波器に到達する時間に関する窓であり、オーバーラップしても、連続しても、また、間隔があいていてもよい。
そして、上記の目的を達成するためのハイドレートの存在領域探査システムは、反射法地震探査方法によって海底下のハイドレードを探査するためのハイドレートの存在領域探査システムであり、広周波数帯域を有する音波を送信する震源となる送波器と、該送波器が発信した音波の地層からの反射波を受信する受信器と、該受信器で受信した反射波の信号を解析する解析手段とを有し、該解析手段が、地震探査の計測結果の低周波数帯域を用いて作成した海底の地層構造図から海底疑似反射面の存在領域を特定すると共に、該海底疑似反射面が存在する場所において、地震探査の計測結果の広周波数帯域において求めたスペクトルが、探査の全周波数に渡って大きく、かつ、高周波数帯域のスペクトル強度が低周波数帯域のスペクトル強度より大きい部位を抽出し、前記抽出部位が無ければ、ハイドレートは存在しないと判定し、前記抽出部位があれば該抽出部位にハイドレートが存在すると判定するように構成される。
また、上記のハイドレート層の存在領域探査システムにおいて、前記解析手段が、前記スペクトルを求める際に、共通反射点重合法を用いるように構成される。
さらには、上記のハイドレート層の存在領域探査システムにおいて、前記解析手段が、前記スペクトルを深さ方向に移動させる解析窓毎に求め、スペクトルにおいて前記抽出部位が生じ始める解析窓がある深さをハイドレート層の上面とし、前記海底疑似反射面をハイドレート層の下面とするように構成される。
これらのハイドレートの存在領域探査システムによれば、上記のハイドレートの存在領域探査方法を実施でき、同様の効果を奏することができる。
本発明のハイドレートの存在領域探査方法及び探査システムによれば、ハイドレートの存在領域を海底地震探査により推定できるようになるので、コストの高い試掘を行う必要がなくなり、経済的にハイドレートを利用することが可能となる。
以下図面を参照して本発明に係るハイドレートの存在領域探査方法及び探査システムの実施の形態について、自律型無人航走体を使用してメタンハイドレートを探索する場合を例にして説明する。
最初に、このハイドレードの存在領域探査方法を実施するためのハイドレードの存在領域探査システムについて説明する。図1に示すように、この探査システムは、自律型無人航走体1と母船6と海底5に設置されたトランスポンダ4等から構成される。
自律型無人航走体1は、その底部に可変周波数振動源(送波器)2を、船尾に送受信装置7をそれぞれ備え、また、船尾からストリーマケーブル3を曳航する。この可変周波数振動源2は、自律型無人航走体1の規模と動力源などの関係からエアガンを使用できないため、小型軽量で省電力型の圧電素子や超磁歪合金などで構成される。これらは、連続的に発信でき、また、波形などを電気的に処理できる等のメリットがある。また、ストリーマケーブル3は、反射波を捉えるハイドロホン(受波器)、深度、方位、温度などを検出するセンサモジュール、データの伝送を行う伝送モジュール、自律型無人航走体1の振動をストリーマケーブル3に伝えないようにする振動制御モジュール、ストリーマケーブル3の姿勢制御を行うスタビライザモジュール等を備えて構成される。
この自律型無人航走体1は、母船6によって探査海域まで運搬された後、海中に降ろされて、予め海底5に設置された複数のトランスポンダ4から発信される信号を受信しながら、また、船尾の送受信装置7と母船6の送受信器8により、交信したり相対位置を検出しながら、所定の海域を所定の水深で航走する。例えば、海底より数m〜数百mの高さを保って航走する。
この自律航走中に、自律型無人航走体1の底部の可変周波数振動源2から海底5に向けて振動波P1を発する。この振動波P1の反射波は、自律型無人航走体1によって曳航されているストリーマケーブル3のハイドロホンによって受信され、記録装置に記録される。この場合において、振動波P1が比較的高い周波数であっても、海底5の近傍から発射され、その反射波を海底近傍で曳航されているストリーマケーブル3のハイドロホンで受信するため高感度で反射波を受信できる。そのため、高い分解能を得ることができる。
この地震探査後に、自律型無人航走体1を母船6に回収して、記録装置に記録されたデータを解析して地層構造等の解析を行う。これによって、振動波P1による地震探査が行われる。
なお、ここでは、このハイドレードの存在領域探査方法を実施するためのハイドレードの存在領域探査システムとして、自律型無人航走体を使用する反射法地震探査システムを例にしているが、その他の反射法地震探査システムとすることができ、所定の低周波と高周波による探査ができるシステムであればよい。
このハイドレートの存在領域探査システムの解析装置(解析手段)は、データ収集部とデータ処理部で構成される。データ収集部はハイドロホンで受信したデータを記録装置に記録する部分であり、データ処理部は、地震探査処理部とハイドレート存在領域探査部で構成される。この地震探査処理部で地層の構造を描画し、メタンハイドレート層の下面と推定される海底疑似反射面の検出を行い、ハイドレート存在領域探査部では、メタンハイドレートの存否の判定と、メタンハイドレート層の上面を推定して、ハイドレート存在領域を特定するための処理を行う。
このハイドレートの存在領域探査方法は、メタンハイドレートの探査においては、図2のフローに示すように、最初のS1ステップで、100〜800Hz好ましくは100〜1000Hz程度の広周波数帯域の地震波を用いて、広周波数地震探査を行う。
次のS2ステップで、この計測結果のデータにおいて、100〜200Hzの範囲に入るような低周波数に対する応答データから海底疑似反射面(BSR)を検出し、メタンハイドレート層の下面を得る。つまり、応答波形の100Hz付近の低周波成分を用いて、通常の地震探査処理により海底の地層構造図を描画し、この地層構造図から海底疑似反射面を探索し、検出した海底疑似反射面をメタンハイドレート層の下面とする。
次のS3ステップで、海底疑似反射面の存在領域を対象として、海底疑似反射面の直上の範囲で、共通反射点(CDP)重合法等の重合法を用いて、経路の異なった反射波をその共通反射点位置での垂直走時に変換して重合(加算)し、この重合波形について周波数分析を行う。つまり、海底疑似反射面が検出された部位の上部において、広い範囲の、好ましくは、計測した全ての範囲の周波数帯域を用いて重合波形を計算し、この計算された重合波形の周波数分析を行う。
この周波数分析は、図3に示すように、解析窓(時間窓)Wを設けて、この解析窓の中にある信号の周波数分析を行い、一種のソナグラムを得る。この解析窓Wを縦軸の深さ方向(時間軸の時間方向)に移動させて、各深さにおけるスペクトルを得る。この解析窓Wはオーバーラップしても、連続しても、また、間隔があいていてもよく、要は深さ方向のスペクトルの変化の情報を得られればよい。また、この周波数分析においては、高速フーリエ変換法(FFT)、最大エントロピー法(MEM)、ウェーブレット(WT)などの任意の周波数分析法を適用することができ、これにより、周波数別の強度分布、あるいはエネルギー分布であるスペクトルを得る。
次のS4ステップでは、スペクトルのチェックを行い、解析窓Wの移動に応じて順次得られるスペクトルの態樣を検査し、全体的にスペクトルが大きいか、低周波数成分よりも高周波数成分の方が強いかを判定する。なお、メタンハイドレートの有無のみの判定であれば、海底疑似反射面より上のメタンハイドレートが存在しそうな範囲をカバーする大きな解析窓一つのスペクトラムでチェックしてもよい。
このS4ステップの判定で、通常と比較してスペクトルが大きくなく、又は、図4に模式的に示すように、低周波成分と比較して高周波成分のスペクトル値が大きくない場合には、S5ステップに行き、メタンハイドレートが存在しないものと判定する。これにより、探査を終了し、必要に応じて次の領域における探査を行う。
また、このS4ステップの判定で、通常のスペクトラルと比較してスペクトルが大きく、かつ、図5に模式的に示すように、低周波成分と比較して高周波成分のスペクトル値が大きい場合には、S6ステップに行き、メタンハイドレートが存在するものと判定する。
この通常のスペクトルとは、海底疑似反射面よりもかなり上の領域や周囲のメタンハイドレートが存在しない領域におけるスペクトルのことである。このメタンハイドレートが存在可能かどうかは、メタンハイドレートが存在するための温度範囲、圧力範囲が分かっているため、これらのデータを併用することにより、高い精度で判定できる。また、実用上は、海底疑似反射面よりもかなり上の領域や周囲のメタンハイドレート層が存在しない領域におけるスペクトルとの比較によって判定する。
なお、下記の地層モデルによるシミュレーションの結果から考えると、メタンハイドレートが存在しない領域と、メタンハイドレートが存在する領域で、模式的には、図4、図5のようなスペクトルが得られることが予想されるが、実際にはスペクトルは震源の周波数特性やその他の影響を受けて、曲線になり、その形態も状況に応じて変化すると考えられる。
次のS7ステップに行き、深さ方向に移動させた解析窓W毎のスペクトルを比較して、スペクトルに変化が生じ始める深さをメタンハイドレート層の上面とする。この最初にメタンハイドレートが存在すると判定される深さから、海底疑似反射面までをメタンハイドレート存在領域とする。これにより、探査を終了し、必要に応じて次の領域における探査を行う。
次に、多数の散在する微小なメタンハイドレートの塊で構成されるメタンハイドレート層を有する地層モデルに対する地震応答をシミュレーション計算により求めた結果を示す。図6にこのメタンハイドレート層15を有する地層モデル10を示し、図7にメタンハイドレートの小塊群の詳細を示す。また、比較のために、メタンハイドレート塊群を取り除いた地層モデルに対する地震応答も求めた。図8にこのメタンハイドレート塊群15を取り除いた地層モデル10Aを示す。
これらの地層モデル10,10Aでは、海水11の下に第1地層12とフリーガス層13と第2地層14を設け、地層モデル10では、このフリーガス層13の上に、メタンガスハイドレート塊群15を配置した。この地層モデル10は、実物大とし、幅390mで、水深320mの海水部分を含めた全深さが440mであり、第1地層12は厚さ55mので、フリーガス層13は、厚さ5mで、第2地層14は厚さ60mとした。また、メタンハイドレート塊群15の小塊15aは、実際の寸法では上底の長さ1m、下底の長さ3m、高さ1mの台形形状のものに相当する寸法とした。そして、それぞれの層に対してそれぞれの層の物性値を与えている。また、メタンハイドレート内の音速は3000m/s(縦波)、1175m/s(横波)とした。
この地震応答試験では、水深280mの送波点のトランスデューサ(送波器)16から音波を発信し、水深280mの受波点のハイドロホン(受波器)17で反射波を受信した。この送受波点を移動して、送波点と受波点との間の距離Dを50mから90mに変更しながら反射波を受信した。
この地層モデル10,10Aにおいては、地層の縦波速度を1800m/sと仮定すると、1mは1800Hzの波長に相当することになる。また、100Hzの波長は18m、500Hzの波長は3.6m、1000Hzの波長は1.8mである。
図9〜図12に、それぞれ100Hz、200Hz、400Hz、500Hzにおける地震応答を示す。それぞれ、図の(a)はメタンハイドレートの塊15aがある場合を示し、図の(b)はメタンハイドレートの塊15aがない場合を示す。これらの図において、横軸は送受波点間の距離(オフセット)Dであり、縦軸は時刻tであり、深さに対応する。この時刻tは震源であるトランスデューサ16が音波を発生した時刻をゼロとして表示している。
この図9〜図12によると、主として反射波は2つ現れており、図中のAで示す反射波は、海底面からの反射波であり、Bで示す反射波はフリーガス層からの反射波であり、海底疑似反射面に相当する。図9からは、周波数が100Hzの場合には、メタンハイドレートの塊15aの影響が小さく、メタンハイドレートの塊15aの有無による違いが殆ど無いことが分かる。一方、図10〜図12に示すように周波数が高くなれば、海底疑似反射面に相当する反射が小さくなり、この海底疑似反射面の直前に微小な反射波が発生することが分かり、これらの変化は周波数が高くなるほど大きくなることも分かった。これらの地震応答結果は、従来の知見とよく符合している。
次に、図9〜図12の計測データに対して、共通反射点(CDP:Common Depth Point)重合を行った結果を図13〜図16に示す。図13では、重合波形においてもメタンハイドレートの塊15aの有無による違いが殆ど見られないが、図14では海底疑似反射面の強度が若干弱くなっており、矢印Cで示すように海底疑似反射面の前に微小な反射が見られる。この傾向は図15及び図16では強くなっている。
上記の図9〜図12に示した地震応答波形、及び、図13〜図16に示した重合波形の検討から、多数散在するメタンハイドレートの小塊15aでメタンハイドレート層15が構成されている地層モデル10において、そのメタンハイドレート層15の探査に有効であることが確認され、本発明は、多数散在するメタンハイドレートの小塊15aで構成されたメタンハイドレート層15に対しては、有効なハイドレートの存在領域探査法となることが分かった。
1 自律型無人航走体
2 振動源(送波器)
3 ストリーマケーブル
4 トランスポンダ
5 海底
6 母船
7 送受信装置
8 送受信器
10 メタンハイドレート層を有する地層モデル
10A メタンハイドレート層が無い地層モデル
11 海水
12 第1地層
13 フリーガス層
14 第2地層
15 メタンハイドレート層
15a メタンハイドレート小塊
15b 地層
16 トランスデューサ(送波器)
17 ハイドロホン(受波器)
2 振動源(送波器)
3 ストリーマケーブル
4 トランスポンダ
5 海底
6 母船
7 送受信装置
8 送受信器
10 メタンハイドレート層を有する地層モデル
10A メタンハイドレート層が無い地層モデル
11 海水
12 第1地層
13 フリーガス層
14 第2地層
15 メタンハイドレート層
15a メタンハイドレート小塊
15b 地層
16 トランスデューサ(送波器)
17 ハイドロホン(受波器)
Claims (6)
- 反射法地震探査方法によって海底下のハイドレードを探査するためのハイドレートの存在領域探査方法であり、広周波数帯域を有する震源を用いて地震探査を行うステップと、該地震探査の計測結果から低周波数帯域を用いて海底の地層構造図を作成し、該地層構造図から海底疑似反射面の存在領域を特定するステップと、該海底疑似反射面が存在する場所において、広周波数帯域におけるスペクトルを求めるステップと、前記スペクトルにおいて、探査の全周波数に渡ってスペクトルが大きく、かつ、高周波数帯域のスペクトル強度が低周波数帯域のスペクトル強度より大きい部位を抽出するステップと、前記抽出部位が無ければ、ハイドレートは存在しないとするステップと、前記抽出部位があれば該抽出部位にハイドレートが存在すると判定するステップとを有してなることを特徴とするハイドレートの存在領域探査方法。
- 前記スペクトルを求める際に、共通反射点重合法を用いることを特徴とする請求項1記載のハイドレートの存在領域探査方法。
- 前記スペクトルを深さ方向に移動させる解析窓毎に求め、スペクトルにおいて前記抽出部位が生じ始める解析窓がある深さをハイドレート層の上面とするステップと、前記海底疑似反射面をハイドレート層の下面とするステップとを有してなることを特徴とする請求項1又は2記載のハイドレートの存在領域探査方法。
- 反射法地震探査方法によって海底下のハイドレードを探査するためのハイドレートの存在領域探査システムであり、広周波数帯域を有する音波を送信する震源となる送波器と、該送波器が発信した音波の地層からの反射波を受信する受信器と、該受信器で受信した反射波の信号を解析する解析手段とを有し、該解析手段が、地震探査の計測結果の低周波数帯域を用いて作成した海底の地層構造図から海底疑似反射面の存在領域を特定すると共に、該海底疑似反射面が存在する場所において、地震探査の計測結果の広周波数帯域において求めたスペクトルが、探査の全周波数に渡って大きく、かつ、高周波数帯域のスペクトル強度が低周波数帯域のスペクトル強度より大きい部位を抽出し、前記抽出部位が無ければ、ハイドレートは存在しないと判定し、前記抽出部位があれば該抽出部位にハイドレートが存在すると判定することを特徴とするハイドレートの存在領域探査システム。
- 前記解析手段が、前記スペクトルを求める際に、共通反射点重合法を用いることを特徴とする請求項4記載のハイドレートの存在領域探査システム。
- 前記解析手段が、前記スペクトルを深さ方向に移動させる解析窓毎に求め、スペクトルにおいて前記抽出部位が生じ始める解析窓がある深さをハイドレート層の上面とし、前記海底疑似反射面をハイドレート層の下面とすることを特徴とする請求項4又は5記載のハイドレートの存在領域探査システム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006187073A JP2008014830A (ja) | 2006-07-06 | 2006-07-06 | ハイドレートの存在領域探査方法及び探査システム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006187073A JP2008014830A (ja) | 2006-07-06 | 2006-07-06 | ハイドレートの存在領域探査方法及び探査システム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2008014830A true JP2008014830A (ja) | 2008-01-24 |
Family
ID=39071982
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006187073A Pending JP2008014830A (ja) | 2006-07-06 | 2006-07-06 | ハイドレートの存在領域探査方法及び探査システム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2008014830A (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014112261A1 (ja) | 2013-01-18 | 2014-07-24 | 株式会社Ihi | 水域地中探査システム及び水域地中探査方法 |
CN107817522A (zh) * | 2017-10-31 | 2018-03-20 | 西安恒泰艾普能源发展有限公司 | 一种海上定性识别水合物顶底界面的确定性方法 |
CN113153232A (zh) * | 2021-03-11 | 2021-07-23 | 浙江大学 | 一种基于小型阵列的海底天然气水合物气泡泄露定位装置及方法 |
CN114194364A (zh) * | 2021-09-17 | 2022-03-18 | 中国人民解放军海军工程大学 | 一种基于水下滑翔机的声电联合探测装置及方法 |
CN114720988A (zh) * | 2022-03-10 | 2022-07-08 | 西北工业大学 | 一种基于模态相位的单水听器低频宽带匹配场测距方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000249760A (ja) * | 1999-03-01 | 2000-09-14 | Nec Corp | 埋没物体探知ソーナーシステムおよびその探知方法 |
JP2003019999A (ja) * | 2001-07-09 | 2003-01-21 | Mitsui Eng & Shipbuild Co Ltd | 海底地層探査システム |
JP2003084069A (ja) * | 2001-09-14 | 2003-03-19 | Mitsubishi Heavy Ind Ltd | ガスハイドレート探査システム及び探査方法 |
JP2006138706A (ja) * | 2004-11-11 | 2006-06-01 | Mitsui Eng & Shipbuild Co Ltd | ハイドレート層の層厚推定方法 |
-
2006
- 2006-07-06 JP JP2006187073A patent/JP2008014830A/ja active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000249760A (ja) * | 1999-03-01 | 2000-09-14 | Nec Corp | 埋没物体探知ソーナーシステムおよびその探知方法 |
JP2003019999A (ja) * | 2001-07-09 | 2003-01-21 | Mitsui Eng & Shipbuild Co Ltd | 海底地層探査システム |
JP2003084069A (ja) * | 2001-09-14 | 2003-03-19 | Mitsubishi Heavy Ind Ltd | ガスハイドレート探査システム及び探査方法 |
JP2006138706A (ja) * | 2004-11-11 | 2006-06-01 | Mitsui Eng & Shipbuild Co Ltd | ハイドレート層の層厚推定方法 |
Non-Patent Citations (1)
Title |
---|
門元之郎、他: "メタンハイドレート層厚探査手法の研究", 三井造船技報, vol. 第187号, JPN6010065203, 28 February 2006 (2006-02-28), pages 8 - 14, ISSN: 0001775734 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014112261A1 (ja) | 2013-01-18 | 2014-07-24 | 株式会社Ihi | 水域地中探査システム及び水域地中探査方法 |
US9726773B2 (en) | 2013-01-18 | 2017-08-08 | Ihi Corporation | Subaqueous underground survey system and subaqueous underground survey method |
CN107817522A (zh) * | 2017-10-31 | 2018-03-20 | 西安恒泰艾普能源发展有限公司 | 一种海上定性识别水合物顶底界面的确定性方法 |
CN113153232A (zh) * | 2021-03-11 | 2021-07-23 | 浙江大学 | 一种基于小型阵列的海底天然气水合物气泡泄露定位装置及方法 |
CN113153232B (zh) * | 2021-03-11 | 2024-04-02 | 浙江大学 | 一种基于小型阵列的海底天然气水合物气泡泄露定位装置及方法 |
CN114194364A (zh) * | 2021-09-17 | 2022-03-18 | 中国人民解放军海军工程大学 | 一种基于水下滑翔机的声电联合探测装置及方法 |
CN114720988A (zh) * | 2022-03-10 | 2022-07-08 | 西北工业大学 | 一种基于模态相位的单水听器低频宽带匹配场测距方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7539079B2 (en) | System and method for determining positions of towed marine source-array elements | |
US10690792B2 (en) | Amplitude-versus-angle analysis for quantitative interpretation | |
AU2009293586B2 (en) | Cetacean protection system | |
US9477000B2 (en) | System and method for the removal of shallow water multiples using a hybrid multi-channel prediction method | |
US9052380B2 (en) | System for doppler positioning of seismic sensors and method | |
EP2530491B1 (en) | Methods and apparatus for seismic exploration using pressure changes caused by sea-surface variations | |
US20120224454A1 (en) | Method and device to acquire marine seismic data | |
CN113534260B (zh) | 一种宽频带近海底深海地质结构声学探测系统与方法 | |
EP2607920A1 (en) | Method and device for estimating an inter-node distance between nodes arranged along towed acoustic linear antennas. | |
MX2011010162A (es) | Procesamiento de datos sismicos. | |
MX2011008293A (es) | Sistema de posicinamiento de cable marino basado en un sensor de movimiento de particulas. | |
US9759828B2 (en) | Determining a streamer position | |
JP2008014830A (ja) | ハイドレートの存在領域探査方法及び探査システム | |
CN104181586A (zh) | 一种水陆检波器数据海底反射系数反演的方法 | |
KR101864307B1 (ko) | 그래디언트 분석을 통한 탄성파 너울보정 및 지하지층구조 탐사 자료 처리방법 및 그 기록매체 | |
US20130046472A1 (en) | Method of determining the relative position of two detectors at the bottom of the sea | |
KR101416278B1 (ko) | 첩 방식 천부지층 탐사 장비 및 이를 이용한 천부지층 탐사 방법 | |
CN102939548B (zh) | 确定放置于海底的探测器的位置的方法 | |
US9606251B2 (en) | Hyperbolic shooting method and device | |
JP2006138706A (ja) | ハイドレート層の層厚推定方法 | |
US12050200B2 (en) | Method of monitoring the interior of a subsea pipeline | |
US20170248710A1 (en) | Device and method for mitigating seismic survey interference | |
US20140254308A1 (en) | Identifying reflection acoustic signals | |
US20210048545A1 (en) | Surveying with low frequency impulse sources | |
Nasralla | Distributed Acoustic Sensing (DAS) As a New Tool For Subsurface Imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20090310 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20101027 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101116 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20110315 |