(57)【特許請求の範囲】
【請求項１】
地下貯留層構造の貯留層特性を、前記貯留層から利用可能な入力データのコンピュータにおけるサポートベクターマシン処理によってモデリングし、前記利用可能な入力データの回帰分析によって、前記地下貯留層中の対象となる領域での前記貯留層特性の尺度を形成するための、コンピュータ実施方法であって、
（ａ）前記貯留層の地震調査から得られる地震調査データから、地下属性に関するトレーニング入力データを受信するステップと、
（ｂ）前記貯留層中の井戸から得られるデータから、形成岩特徴に関するトレーニングターゲットデータを受信するステップと、
（ｃ）地下属性トレーニングデータ及び形成岩特徴トレーニングターゲットデータを複数のサブセットに区分するステップと、
（ｄ）サポートベクターマシンモデリングに対する形成属性パラメータを選択するステップであり、
（１）カーネルパラメーター値とペナルティーパラメータ値の対を有する動径基底カーネル関数対に対して、地下属性トレーニングデータの各サブセットを、前記複数のサブセットの他のサブセットと相互検証するステップと、
（２）相互検証した各サブセットの誤差関数を形成するステップと、
（３）前記地下属性トレーニングデータのサブセットを前記相互検証するステップを繰り返し、複数の異なる動径基底カーネル関数対に対して誤差関数を形成するステップと、
を実行することによって、選択するステップと、
（ｅ）前記複数の異なる動径基底カーネル関数対に対して前記形成された誤差関数の最小誤差関数を決定することによって、前記選択された形成属性パラメータを最適化するステップと、
（ｆ）前記トレーニングターゲットデータ、前記選択された形成属性パラメータ、地下属性トレーニングデータの前記相互検証されたサブセット、および前記複数の動径基底カーネル関数対に対する前記誤差関数を、サポートベクターマシンモデリングに対するトレーニング入力として提供するステップと、
（ｇ）前記提供されたトレーニング入力の前記誤差関数の最小誤差関数を決定するための回帰分析により、サポートベクターマシンモデリングを実行するステップと、
（ｈ）前記トレーニング入力の前記サポートベクターモデリングに基づいて前記貯留層特性を予測するステップと、
（ｉ）前記予測された貯留層特性の出力表示を形成するステップと、
の前記コンピュータ処理ステップを含む、方法。
【請求項２】
前記コンピュータは、グラフ表示デバイスを含み、さらに、
前記予測された貯留層特性の記録の記憶
を含む、請求項１に記載のコンピュータ実施方法。
【請求項３】
前記貯留層特性は、前記貯留層中の対象となる領域の空隙率を含む、請求項１に記載のコンピュータ実施方法。
【請求項４】
前記貯留層特性は、前記貯留層中の対象となる領域の音響インピーダンスを含む、請求項１に記載のコンピュータ実施方法。
【請求項５】
前記トレーニング入力データの前記地下属性は、振幅、周波数、第１のエンベロープデータ、第２のエンベロープデータ、位相、コヒーレンス、帯域幅、および音響インピーダンスを含む群から選択される、請求項１に記載のコンピュータ実施方法。
【請求項６】
前記トレーニングターゲットデータの前記形成岩特徴は、空隙率および音響インピーダンスを含む群から選択される、請求項１に記載のコンピュータ実施方法。
【請求項７】
前記予測された貯留層特性は、空隙率を含む、請求項１に記載のコンピュータ実施方法。
【請求項８】
地下貯留層構造の貯留層特性を、前記貯留層から利用可能な入力データのコンピュータにおけるサポートベクターマシン処理によってモデリングし、前記利用可能な入力データの回帰分析によって、前記地下貯留層中の対象となる領域での前記貯留層特性の尺度を形成するための、データ処理システムであって、
（ａ）
（１）前記貯留層の地震調査から得られる地震調査データから、地下属性に関するトレーニング入力データを受信するステップと、
（２）前記貯留層中の井戸から得られるデータから、形成岩特徴に関するトレーニングターゲットデータを受信するステップと、
（３）地下属性トレーニングデータ及び形成岩特徴トレーニングターゲットデータを複数のサブセットに区分するステップと、
（４）前記トレーニング入力データの前記トレーニングターゲットデータとの相互相関によってサポートベクターマシンモデリングに対する形成属性パラメータを選択するステップであり、
（ｉ）カーネルパラメーター値とペナルティーパラメータ値の対を有する動径基底カーネル関数対に対して、前記地下属性トレーニングデータの各サブセットを、前記複数のサブセットの他のサブセットと相互検証することと、
（ｉｉ）前記相互検証した各サブセットの誤差関数を形成することと、
（ｉｉｉ）前記地下属性トレーニングデータのサブセットを前記相互検証するステップを繰り返し、複数の異なる動径基底カーネル関数対に対する誤差関数を形成することと、
によって、選択するステップと、
（５）前記複数の異なる動径基底カーネル関数対に対する前記形成された誤差関数の最小誤差関数を決定することにより、前記選択された形成属性パラメータを最適化するステップと、
（６）前記トレーニングターゲットデータ、前記選択された形成属性パラメータ、地下属性トレーニングデータの前記相互検証されたサブセット、および複数の動径基底カーネル関数対に対する前記誤差関数を、サポートベクターマシンモデリングに対するトレーニング入力として提供するステップと、
（７）前記提供されたトレーニング入力の前記誤差関数の最小誤差関数を決定するための回帰分析により、サポートベクターマシンモデリングを実行するステップと、
（８）前記トレーニング入力の前記サポートベクターマシンモデリングに基づき、前記貯留層特性を予測するステップと、
を実行するためのプロセッサと、
（ｂ）
前記予測された貯留層特性の出力表示を形成するステップ
を実行するための表示部と、
を備える、データ処理システム。
【請求項９】
さらに、
前記予測された貯留層特性の記録を記憶する、データメモリ
を含む、請求項８に記載のデータ処理システム。
【請求項１０】
前記貯留層特性は、前記貯留層中の対象となる領域の空隙率を含む、請求項８に記載のデータ処理システム。
【請求項１１】
前記貯留層特性は、前記貯留層中の対象となる領域の音響インピーダンスを含む、請求項８に記載のデータ処理システム。
【請求項１２】
前記トレーニング入力データの前記地下属性は、振幅、周波数、第１のエンベロープデータ、第２のエンベロープデータ、位相、コヒーレンス、帯域幅、および音響インピーダンスを含む群から選択される、請求項８に記載のデータ処理システム。
【請求項１３】
前記トレーニングターゲットデータの前記形成岩特徴は、空隙率および音響インピーダンスを含む群から選択される、請求項８に記載のデータ処理システム。
【請求項１４】
前記予測された貯留層特性は、空隙率を含む、請求項８に記載のデータ処理システム。
【請求項１５】
データ処理システムに、地下貯留層構造の貯留層特性を、前記貯留層から利用可能な入力データのコンピュータにおけるサポートベクターマシン処理によってモデリングさせ、前記利用可能な入力データの回帰分析によって、前記地下貯留層中の対象となる領域での前記貯留層特性の尺度を形成させるための、動作可能な命令を、コンピュータ可読媒体コンピュータ中に記憶している、データ記憶デバイスであって、前記データ記憶デバイス中に記憶された前記命令は、前記データ処理システムに以下の：
（ａ）前記貯留層の地震調査から得られる地震調査データから、地下属性に関するトレーニング入力データを受信するステップと、
（ｂ）前記貯留層中の井戸から得られるデータから、形成岩特徴に関するトレーニングターゲットデータを受信するステップと、
（ｃ）地下属性トレーニングデータ及び形成岩特徴トレーニングターゲットデータを複数のサブセットに区分するステップと、
（ｄ）サポートベクターマシンモデリングに対する形成属性パラメータを選択するステップであり、
（１）カーネルパラメーター値とペナルティーパラメーター値の対を有する動径基底カーネル関数対に対して、前記地下属性トレーニングデータの各サブセットを、前記複数のサブセットの他のサブセットと相互検証するステップと、
（２）前記相互検証した各サブセットの誤差関数を形成するステップと、
（３）前記地下属性トレーニングデータのサブセットを前記相互検証するステップを繰り返し、複数の異なる動径基底カーネル関数対に対する誤差関数を形成するステップと、
を実行することによって、選択するステップと、
（ｅ）前記複数の異なる動径基底カーネル関数対に対する前記形成された誤差関数の最小誤差関数を決定することによって、前記選択された形成属性パラメータを最適化するステップと、
（ｆ）前記トレーニングターゲットデータ、前記選択された形成属性パラメータ、地下属性トレーニングデータの前記相互検証されたサブセット、および前記複数の動径基底カーネル関数対に対する前記誤差関数を、サポートベクターマシンモデリングに対するトレーニング入力として提供するステップと、
（ｇ）前記提供されたトレーニング入力の前記誤差関数の最小誤差関数を決定するための回帰分析により、サポートベクターマシンモデリングを実行するステップと、
（ｈ）前記トレーニング入力の前記サポートベクターマシンモデリングに基づき、前記貯留層特性を予測するステップと、
（ｉ）前記予測された貯留層特性の出力表示を形成するステップと、
を実行させる、データ記憶デバイス。
【請求項１６】
前記データ処理システムは、データメモリを含み、前記命令はさらに、前記データ処理システムに、
前記予測された貯留層特性の記録を記憶するステップ
を実行させる命令を含む、請求項１５に記載のデータ記憶デバイス。
【請求項１７】
前記貯留層特性は、前記貯留層中の対象となる領域の空隙率を含む、請求項１５に記載のデータ記憶デバイス。
【請求項１８】
前記貯留層特性は、前記貯留層中の対象となる領域の音響インピーダンスを含む、請求項１５に記載のデータ記憶デバイス。
【請求項１９】
前記トレーニング入力データの前記地下属性は、振幅、周波数、第１のエンベロープデータ、第２のエンベロープデータ、位相、コヒーレンス、帯域幅、および音響インピーダンスを含む群から選択される、請求項１５に記載のデータ記憶デバイス。
【請求項２０】
前記トレーニングターゲットデータの前記形成岩特徴は、空隙率および音響インピーダンスを含む群から選択される、請求項１５に記載のデータ記憶デバイス。
【請求項２１】
前記予測された貯留層特性は、空隙率を含む、請求項１５に記載のデータ記憶デバイス。
【請求項２２】
前記カーネル対を、前記地下属性トレーニングデータのサブセットに適用するステップを更に含む請求項１に記載のコンピュータ実施方法。
【請求項２３】
前記選択された形成属性パラメータを最適化するステップが、
各属性に対する前記トレーニングデータセット及び試験データセットを、同一のスケールで定義された範囲の値にスケーリングするステップを含む請求項１に記載のコンピュータ実施方法。
【請求項２４】
前記プロセッサは、前記カーネル対を、前記地下属性トレーニングデータのサブセットに適用するステップを更に実行する請求項８に記載のデータ処理システム。
【請求項２５】
前記選択された形成属性パラメータを最適化する前記プロセッサは、
各属性に対する前記トレーニングデータセット及び試験データセットを、同一のスケールで定義された範囲の値にスケーリングするステップを実行する請求項８に記載のデータ処理システム。
【請求項２６】
前記命令は、前記カーネル対を、前記地下属性トレーニングデータのサブセットに適用するステップを実行するための命令を更に含む請求項１５に記載のデータ記憶デバイス。
【請求項２７】
前記選択された形成属性パラメータを最適化する命令が、
各属性に対する前記トレーニングデータセット及び試験データセットを、同一のスケールで定義された範囲の値にスケーリングするステップを実行するための命令を更に含む請求項１５に記載のデータ記憶デバイス。
【発明の詳細な説明】
【技術分野】
【０００１】
本願は、２０１１年９月２８日出願の米国仮出願第６１／５４０，２６３号からの優先権を主張する。米国特許実務の目的で、本願は、参照により全体にこの仮出願の内容を組み込んでいる。
【０００２】
本発明は、地中の炭化水素貯留層のコンピュータ化シミュレーションに関し、具体的には、利用可能な限られた量の検層データがあるときの、貯留層モデルとしての貯留層属性または特性の決定に関する。
【背景技術】
【０００３】
ほんの数個の井戸しか貯留層中に存在しない場合、地震属性から貯留層特性を予測することは、貯留層モデリングおよび評価のための初期の貯留層特性化を提供するための石油産業における共通の課題である。従来、地球統計学およびニューラルネットワークは、種々の地震属性の組み合わせから検層特性を予測するために使用されてきた。既知の井戸と地震との結び付きは、地震データと井戸の値との間の関係を学習するために使用されてきた。近年、多重地震属性は、井戸と地震属性との間の関係の学習に基づき、モデリング技術によって、検層特性を予測するために使用されてきた。しかしながら、小さい母集団（すなわち、ほんの数個の井戸と地震属性との対）に対して、統計的有意性は、場合によっては達成することが不可能であった。しかしながら、この技術の使用は、貯留層中に実際に存在する井戸の数によって決まったか、またはそれによって制限された。
【０００４】
掘削された井戸の限られた利用可能性のため、貯留層石油物理学的特徴モデリングは典型的に、不確実性に悩まされてきた。ニューラルネットワーク方法は、地震属性を使用した貯留層予測のために開発されてきた。かかる方法は、貯留層炭化水素を予測するために、例えば、誤差逆伝搬またはＢＰニューラルネットワークおよび自己組織化マップまたはＳＯＭニューラルネットワークに基づいてきた。しかしながら、ニューラルネットワークモデリング方法に基づき、貯留層特性を予測するという試みは、容易に「過学習」として知られているネットワークをもたらし、それが次に「オーバーフィッティング」をもたらし、このように検証試験において不十分な予測を提供してきた。
【０００５】
地震属性は、地中の特定の間隔または形成層にわたる地震派トレースの特徴の定量的尺度である。地震属性は、地下構造の統合ならびに炭化水素の存在および位置の予測に対して、可能な限りの情報を提供することができる。地震属性は、炭化水素探査識別およびリスク、炭化水素挙動評価、貯留層特性化等に一般的に使用されてきた。地震属性の利点は、それらが井戸においておよび井戸から離れて予測することができ、その上なお井戸データを評価できることである。しばしば、予測は、単純に井戸データを補間するよりも詳細である。地震データおよびそれらの変換の両方から計算されている、多くの異なる地震属性がある。しかしながら、地下構造を統合し、貯留層を予測するためのこれらの属性の使用は、今までに知られている限り、適切に対処されていない問題である。
【発明の概要】
【０００６】
簡潔に、本発明は、地下貯留層構造の貯留層特性を、貯留層から利用可能な入力データのコンピュータにおけるサポートベクターマシン処理によってモデリングし、利用可能な入力データの回帰分析によって、地下貯留層中の対象となる領域での貯留層特性の尺度を形成するための、新しいかつ改良されたコンピュータ実施方法を提供し、方法は、（ａ）貯留層の地震調査から得られる地震調査データから、地下属性に関するトレーニング入力データを受信するステップ、（ｂ）貯留層中の井戸から得られるデータから、形成岩特徴に関するトレーニングターゲットデータを受信するステップ、トレーニング入力データのトレーニングターゲットデータとの相互相関によって、サポートベクターマシンモデリングに対する形成属性パラメータを選択するステップ、選択された形成属性パラメータを最適化するステップ、トレーニング入力データおよび最適化された選択された形成属性パラメータを、サポートベクターマシンモデリングに対する入力として提供するステップ、最適化された選択された形成属性パラメータに基づき、貯留層特性を予測するために、サポートベクターマシンモデリングを実行するステップ、ならびに予測された貯留層特性の出力表示を形成するステップのコンピュータ処理ステップを含む。
【０００７】
本発明はまた、地下貯留層構造の貯留層特性を、貯留層から利用可能な入力データのコンピュータにおけるサポートベクターマシン処理によってモデリングし、利用可能な入力データの回帰分析によって、地下貯留層中の対象となる領域での貯留層特性の尺度を形成する、新しいかつ改良されたデータ処理システムを提供し、データ処理システムは、貯留層の地震調査から得られる地震調査データから、地下属性に関するトレーニング入力データを受信するステップ、貯留層中の井戸から得られるデータから、形成岩特徴に関するトレーニングターゲットデータを受信するステップ、トレーニング入力データのトレーニングターゲットデータとの相互相関によって、サポートベクターマシンモデリングに対する形成属性パラメータを選択するステップ、選択された形成属性パラメータを最適化するステップ、トレーニング入力データおよび最適化された選択された形成属性パラメータを、サポートベクターマシンモデリングに対する入力として提供するステップ、最適化された選択された形成属性パラメータに基づき、貯留層特性を予測するために、サポートベクターマシンモデリングを実行するステップを実行するためのプロセッサと、予測された貯留層特性の出力表示を形成するステップを実行するための表示部とを備える。
【０００８】
本発明はさらに、データ処理システムに、地下貯留層構造の貯留層特性を、貯留層から利用可能な入力データのコンピュータにおけるサポートベクターマシン処理によってモデリングさせ、利用可能な入力データの回帰分析によって、地下貯留層中の対象となる領域での貯留層特性の尺度を形成させるための、動作可能な命令を、コンピュータ可読媒体コンピュータ中に記憶している、新しいかつ改良されたデータ記憶デバイスを提供し、データ記憶デバイス中に記憶された命令は、データ処理システムに以下のステップを実行させる：貯留層の地震調査から得られる地震調査データから、地下属性に関するトレーニング入力データを受信するステップ、貯留層中の井戸から得られるデータから、形成岩特徴に関するトレーニングターゲットデータを受信するステップ、トレーニング入力データのトレーニングターゲットデータとの相互相関によって、サポートベクターマシンモデリングに対する形成属性パラメータを選択するステップ、選択された形成属性パラメータを最適化するステップ、トレーニング入力データおよび最適化された選択された形成属性パラメータを、サポートベクターマシンモデリングに対する入力として提供するステップ、最適化された選択された形成属性パラメータに基づき、貯留層特性を予測するために、サポートベクターマシンモデリングを実行するステップ、ならびに予測された貯留層特性の出力表示を形成するステップ。
【図面の簡単な説明】
【０００９】
【図１】本発明に係る、地下地中形成の貯留層特性予測のためのデータ処理システムで実行される一連のデータ処理ステップの機能ブロック図である。
【００１０】
【図２】本発明に係る、地下地中形成の貯留層特性予測のためのデータ処理システムの概略ブロック図である。
【００１１】
【図３】本発明に係る、地下地中形成の貯留層特性予測中の合成データから得られる試験結果のプロットである。
【００１２】
【図４】本発明に係る、地下地中形成の貯留層特性予測中の合成データから得られる試験結果のプロットである。
【００１３】
【図５】本発明に係る、地下地中形成の貯留層特性予測中に得られる井戸に対する空隙率予測のプロットである。
【００１４】
【図６】本発明に係る、地下地中形成の貯留層特性予測中に得られる別の井戸に対する空隙率予測のプロットである。
【００１５】
【図７】本発明に係る、地下地中形成の貯留層特性予測中に得られる別の井戸に対する空隙率予測のプロットである。
【００１６】
【図８】本発明に係る、地下地中形成の貯留層特性予測中に得られる地下地中形成の縦断面上に重ね合わせた音響インピーダンスのプロットである。
【００１７】
【図９】本発明に係る、地下地中形成の貯留層特性予測中に得られ、かつ図８と同一の地下地中形成の縦断面上に重ね合わせた空隙率のプロットである。
【００１８】
【図１０】本発明に係る、地下地中形成の貯留層特性予測に対し、かつ地下地中形成の縦断面上に重ね合わせた入力周波数属性地震データのプロットである。
【００１９】
【図１１】本発明に係る、地下地中形成の貯留層特性予測中に得られ、かつ図１０と同一の地下地中形成の縦断面上に重ね合わせた空隙率のプロットである。
【発明を実施するための形態】
【００２０】
本発明では、本発明に係るある特定の処理方法もまた実行されるという条件で、サポートベクターマシン（ＳＶＭ）として知られるコンピュータ化モデリング方法が貯留層属性または特性予測に使用され得ることがわかっている。サポートベクターマシンは、教科書「Ｓｔａｔｉｓｔｉｃｓ  Ｌｅａｒｎｉｎｇ  Ｔｈｅｏｒｙ」（１９９８年）に記載されるように、Ｖ．Ｖａｐｎｉｋによって考案された統計的学習理論に基づく。サポートベクターマシンは、機械学習モデルの間で既知の技術である。本発明では、一形態のＳＶＭ回帰が実行される。本発明によるＳＶＭ回帰はまた、処理のための種々のカーネルの適切な形態のうちの選択された一形態を用いた処理、余剰属性またはターゲットに対してほとんど寄与しない可能性がある属性の除去、各属性が均等に寄与することを定めるための属性正規化、ならびにカーネルパラメータおよびペナルティパラメータの自動決定を組み込む。
【００２１】
サポートベクターマシンは、スパースな決定関数展開を作成し、サポートベクターとして知られる選択された数のトレーニング点のみを選択することによって、コンピュータ中に形成される。カーネルの使用により、入力空間中のデータ点間の明示的な内積を伴う線形関数近似アルゴリズムは、それらの非線形一般化に便利にかつ効率的に変換され得る。ＳＶＭは、経験的エラー（リスク）とモデル複雑度との間のバランスのとれたトレードオフを通して、Ｖａｐｎｉｋの構造的リスク最小化原理をほぼ実施する（ＶＣ次元として既知であるものを用いた、可能性の高いモデル試験エラーバウンドの予測を通して測定される）。
【００２２】
まず、最大マージン分類子が、以下のように２クラス分類線形問題を示すことによって確立される：
【数１】
0006006799.tif 000002
式中、φ（χ）は、固定した特徴空間変換を示し、バイアスパラメータｂは、明示的にされている。Ｗは、カーネル関数を示し、ここで、ガウス型動径基底関数カーネルが使用される。
【００２３】
トレーニングデータセットは、Ｎ入力ベクトルｘ１，．．．ｘｎを含み、対応するターゲット値は、ｔ１，．．．ｔｎであり、
0006006799.tif 000003
であり、新しいデータ点ｘは、ｙ（ｘ）の記号に従って分類される。
【００２４】
線形分離可能な特徴空間に対して、パラメータＷおよびｂの少なくとも１つの選択があり、以下を満たす。ｔｎ＝＋１を有する点に対して、
0006006799.tif 000004
であり、ｔｎ＝−１を有する点に対して、ｙ（ｘｎ）＜０であり、全てのトレーニングデータ点に対してｙ（ｘｎ）ｔｎ＞０になるようにする。ｙ（ｘ）の新しいデータ点ｘを予測するためのベクターマシンを表すために使用されるように、最適なＷおよびｂを探す。
【００２５】
回帰モデリング
本発明では、サポートベクターマシン方法は、回帰問題まで及ぶと同時に、スパースネスの特性を維持する。方程式（１）と同一の形態のＳＶＭ回帰モデリングの問題を考慮すると、トレーニングデータセットはさらに、Ｎ入力ベクトルｘ１，．．．ｘｎを含み、対応するターゲット値は、ｔ１，．．．ｔｎ，であり、ここで、ｔｎは、任意のターゲット数列であり得る。サポートベクターマシンは、
【数２】
0006006799.tif 000005
によって与えられる正規化誤差関数を最小化することによって、入力と出力との間の関係をモデリングしようとする。
【００２６】
スパースな解を得るために、上記の二次誤差関数は、以下のように許容誤差関数Ｅに置き換えられる。
【数３】
0006006799.tif 000006
【００２７】
次いで、
【数４】
0006006799.tif 000007
によって与えられる誤差関数等、正規化誤差関数が最小化され得、ここで、ｙ（ｘ）は、方程式（１）によって与えられ、Ｃは、ペナルティパラメータとしても知られる逆正規化パラメータを表す。
【００２８】
次いで、上記で考察された最適化問題は、スラック変数を導入することによって言い換えられ得る。各データ点ｘｎに対して、２つのスラック変数が導入される。
0006006799.tif 000008
【００２９】
ここで、回帰に対する解は、
0006006799.tif 000009
になるように、ほとんどのターゲット点が∈チューブの内にあることを確実にするために、最適なパラメータを見つけることであり、一方で、スラック変数の導入は、スラック変数が非ゼロであるという条件で、点のうちのいくつかがフィッティングの柔軟性のためにチューブの外にあることを可能にする。
【００３０】
次いで、サポートベクター回帰に対する誤差関数は、最終的に、
0006006799.tif 000010
として表され得る。
【００３１】
そのため、
0006006799.tif 000011
および
0006006799.tif 000012
である。ラグランジュ乗数およびカルーシュ・キューン・タッカー（ＫＫＴ）条件として知られるものを使用して、方程式（５）は、バイアス項ｂを計算するためのトレーニングデータ点に使用され得る。非スパースネス最小二乗法に対して、方程式（５）は、
【数６】
0006006799.tif 000013
として表され得る。
【００３２】
したがって、動径基底関数（ＲＢＦ）カーネルを有する最小二乗法は、線形および非線形回帰の両方に対して、方程式（６）の誤差関数を最小化するために組み込まれる。本発明に係る最小二乗ＳＶＭモデリングは、大域的および一意的な解を提供するが、それは、スパースネスおよびサポートベクターの解釈に欠ける。
【００３３】
管理されたトレーニングの前に、自動属性選択が、線形係数重み付け理論の代わりに相互相関によって適用される。相互検証およびグリッド検索技術の２つのパスは、オーバーフィッティングを回避するために、最適なＲＢＦパラメータを達成するように開発される。最小二乗法は、回帰誤差を最小化するために適用される。本発明に係る貯留層特性の予測は、ロバストかつ確実である。
【００３４】
貯留層特性予測のためのコンピュータ処理
図面中、図１に示されるフローチャートＦは、本発明に係る、最小二乗回帰サポートベクターマシンモデリング空隙率を用いて、対象となる地下貯留層の予測された貯留層特性のモデルを形成するために、本発明の基本的なコンピュータ処理シーケンスを示す。フローチャートＦの処理シーケンスは、貯留層中の井戸から得られるデータ、および貯留層を覆う領域中の地震調査の処理から得られる地震属性データに基づく。予測された貯留層特性のモデルを取得するための図１に係るデータの処理は、同様に説明されるデータ処理システムＤ（図２）で実施される。
【００３５】
図１を参照すると、データ処理システムＤにおける処理は、本発明に係る処理のために、入力パラメータまたはデータのデータメモリからの転送によって、ステップ１０（図１）中に開始する。ステップ１０に示されるように、対象となる貯留層を覆う領域の地震調査データの従来の処理から得られる地震属性の形態の入力パラメータまたはデータは、トレーニング入力データとして提供される。トレーニング入力として提供される地震属性データの実施例としては、例えば、振幅、周波数、第１のエンベロープデータ、第２のエンベロープデータ、位相、コヒーレンス、帯域幅、および音響インピーダンスである。他の種類の地震属性データもまた使用され得ることを理解されたい。
【００３６】
ステップ１２に示されるように、サポートベクターマシンモデリングに対するトレーニングターゲットデータの形態の入力データは、データメモリから転送される。ステップ１２中に提供されるトレーニングターゲットデータは、空隙率、音響インピーダンス等、属性の石油物理学的尺度またはパラメータおよび値である。トレーニングターゲットデータは、貯留層中に実際に存在している井戸の検層から得られるデータの従来の処理から得られる。場合によっては、利用可能な場合、コアサンプルデータの処理から得られるトレーニングターゲットデータもまた、利用され得る。
【００３７】
データ処理システムＤにおける処理のステップ１４中、記憶された入力データは、本発明に係る、属性選択を受ける。ステップ１４中の属性選択は、さらなる処理の前に、ステップ１０からのトレーニング入力データの利用可能な地震属性間の選択を提供する。属性選択は、計算作業負荷を最小限に抑えるために、かつ処理効率を高めるために、関連した属性を取得するように実行される。ステップ１４中の属性選択はまた、不適切なデータの悪影響のリスクを低減する。
【００３８】
本発明に係るＳＶＭモデリングの形態を用いて、属性選択に対する相互相関アプローチがステップ１４中に提供される。相互相関アプローチは、費用効率が高く、貯留層予測のためのそれらの関連した属性を選択することがわかっている。ステップ１４中、データ処理は、トレーニング入力データ中の各利用可能な属性と、貯留層中の利用可能な井戸の各トレーニングターゲット入力との間の相関係数を決定する。次いで、指定された閾値よりも低い絶対相関係数を有する属性を除外するために、ユーザ指定の閾値が使用される。
【００３９】
処理は次に、パラメータ最適化のためのステップ１６に進む。パラメータ最適化ステップ１６中、閾値要件を満たすステップからのそれらの選択された属性は、データまたは振幅スケーリングを受ける。データスケーリングは、より小さい数値範囲内の属性より優位である、より大きい数値範囲を有する属性を回避するために実行される。データスケーリングの別の利点は、処理計算中の数値的な問題を回避することである。カーネル値が通常、例えば、線形カーネルおよび多項式カーネル等、特徴ベクトルの内積によるため、大きい属性値は、数値的な問題を引き起こす可能性がある。本発明では、各属性に対するトレーニングおよび試験データセットの両方の線形スケーリングが、同一のスケールで、かつ［−１，１］または［０，１］のいずれかである定義された範囲内の値に提供するように実行される。
【００４０】
ＳＶＭモデリングステップ２４中、データ処理システムＤは、サポートベクターマシン方法に従って動作し、入力データとターゲットデータとの間の関係のモデルを形成する。モデリングは、上記で考察されたように、回帰によって行われ、誤差関数を最小化する。回帰処理中のデータ処理システムＤは、カーネルパラメータγおよびペナルティパラメータＣを有する動径基底関数（ＲＢＦ）カーネルを使用する。
【００４１】
サポートベクターマシンモデリングステップ２４中、かかるＳＶＭ回帰を実行するために、ある特定の生成要件が満たされる必要がある。これらは、種々のカーネルからの適切な処理カーネルの選択、余剰かつターゲットにほとんど寄与しない可能性がある属性の除去、ならびにステップ２４中の使用に選択された各属性が均等に寄与することの確認、ならびにカーネルパラメータγおよびペナルティパラメータＣの自動決定を含む。
【００４２】
以下は、本発明に係る、コンピュータ化貯留層特性予測のためのシーケンスの例示的なリストである。
0006006799.tif 000014
0006006799.tif 000015
【００４３】
方程式（７）に係る動径基底関数（ＲＢＦ）パラメータは、本発明に係るＳＶＭモデリングステップ２４に好ましいことがわかった。
【数７】
0006006799.tif 000016
【００４４】
動径基底関数（ＲＢＦ）カーネルが使用される１つの理由は、それが、サンプルをより高い次元空間に非線形的にマッピングし、そのため、線形カーネルとは違い、クラスラベルと属性との間の関係が非線形である場合に対処することができるからである。さらに、ＲＢＦカーネルは、より単純な線形またはシグモイドカーネルとしての結果を得ることができ、それらは、非線形関係に適合されない。別の理由は、モデル選択の複雑性に影響を与えるハイパーパラメータの数である。多項式カーネルは、ＲＢＦカーネルよりも多くのハイパーパラメータを有する。最終的に、ＲＢＦカーネルは、数値的問題が少ない。
【００４５】
２つのパラメータ、カーネルパラメータγおよびペナルティパラメータＣは、ＳＶＭ予測精度において重要な役割を果たす。実際には、ユーザは、特定のデータセットに対して、最適なパラメータ対（Ｃ，γ）が何なのか前もって知らない。不確実性およびユーザの試行の負担を低減するために、ＳＶＭステップ２４の前に、パラメータグリッド検索ステップ２０および相互検証ステップ２２が、（Ｃ，γ）の最適な対を得るために実行され、その一方で、ＳＶＭが過学習およびオーバーフィッティングされるのを回避する。
【００４６】
ステップ２２中の相互検証は、トレーニングセットがステップ２０中にｎ個のサブセットに分割または区分されるため、ｎ重である。ｎ個のサブセットは必ずしも同一サイズである必要はなく、ｎ重は、ｎ個の検層であり得る。連続的に、トレーニングセット全体の各インスタンスが一度予測されるまで、１個のサブセットが、残りのｎ−１個のサブセットに対してトレーニングされる（Ｃ，γ）の所与のパラメータ対に対して試験される。次いで、最小平均二乗誤差Ｅが、（Ｃ，γ）のこの試験パラメータ対に対して計算および記録される。
【００４７】
ステップ２４中、グリッド検索フレームワークは、（Ｃ，γ）のパラメータ対にわたってループするように提供される。最小誤差Ｅを有する最適なフィッティングパラメータ対（Ｃ，γ）は、次いで、ステップ２６中の生成試験予測として使用される。
【００４８】
計算時間を節約し、なおかなり正確な（Ｃ，γ）パラメータ対を得るために、Ｃの指数関数的に増加しているシーケンスを試すことは、良好なパラメータを識別するための実用的な方法であることがわかっている（例えば、
0006006799.tif 000017
、および
0006006799.tif 000018
）。
【００４９】
また、２回の検索がステップ２４中に実行されることが推奨される。１回目は、ｅ１の指数因子による粗検索であり、２回目のより精細な検索は、ｅ０，１の検索間隔で、１回目の（Ｃ，γ）の周囲で実行される。粗いおよび精細な検索間隔の両方が特定の計算要求を満たすように調節され得ると考えられることを理解されたい。この技術は、最適な（Ｃ，γ）パラメータ対を提供することがわかっている。次いで、ステップ２４中に決定された最適なパラメータ対は、試験データセット全体で対象となる貯留層特性を予測するために、ステップ２６中に、トレーニングデータセット全体に適用される。次いで、予測された貯留層特性は、ステップ２８中に、データ処理Ｄのメモリへの出力として提供され、出力表示を形成するために利用可能である。
【００５０】
データ処理システム
図２に示されるように、本発明に係るデータ処理システムＤは、プロセッサ４０と、その中に動作命令、制御情報、およびデータベース記録を記憶するための、プロセッサ４０に連結されたメモリ４２とを有する、コンピュータ３０を含む。コンピュータ３０は、望みであれば、ラップトップ型コンピュータ、ノート型コンピュータの形態のパーソナルコンピュータ等の携帯デジタルプロセッサ、またはデスクトップコンピュータ等の他の好適なプログラム化もしくはプログラム可能デジタルデータ処理装置であってもよい。また、コンピュータ３０は、Ｉｎｔｅｌ  ＣｏｒｐｏｒａｔｉｏｎもしくはＡｄｖａｎｃｅｄ  Ｍｉｃｒｏ  Ｄｅｖｉｃｅｓ（ＡＭＤ）等からのノードを有するマルチコアプロセッサ、ＨＰＣ  Ｌｉｎｕｘクラスタコンピュータ、もしくはＡｒｍｏｎｋ，Ｎ．Ｙ．のＩｎｔｅｒｎａｔｉｏｎａｌ  Ｂｕｓｉｎｅｓｓ  Ｍａｃｈｉｎｅｓ（ＩＢＭ）から入手可能なもの等、任意の従来の種類の好適な処理能力のメインフレームコンピュータ、または他のソースであってもよいことを理解されたい。
【００５１】
コンピュータ３０は、本発明に係る、岩相および貯留層属性の出力データまたは記録を表示するためのユーザインターフェース４６および出力データ表示部４８を有する。出力表示部４８は、印刷された出力情報を提供することが可能なプリンタおよび出力表示スクリーン、または出力記録もしくは画像としてグラフ、データシート、グラフィックイメージ、データプロット等の形態の可視表示部等のコンポーネントを含む。
【００５２】
コンピュータ３０のユーザインターフェース４６はまた、情報およびデータベース記録を制御し、またはそれらにアクセスし、コンピュータ３０を操作するためのアクセスをユーザに提供するために、好適なユーザ入力デバイスまたは入力／出力制御ユニット５０を含む。データ処理システムＤはさらに、コンピュータメモリ中に記憶されるデータベース５２を含み、それは、内部メモリ４２、または関連データベースサーバ５８中の５６で示されるような外部、ネットワーク化、もしくは非ネットワーク化メモリであってもよい。
【００５３】
データ処理システムＤは、コンピュータ３０のメモリ５４中に記憶されるプログラムコード６０を含む。本発明に係るプログラムコード６０は、データプロセッサ４０に、上記のおよび図１に示される方法で、本発明のコンピュータ実施方法を実施させる、コンピュータで動作可能な命令の形態である。
【００５４】
プログラムコード６０は、データ処理システムＤの機能を制御し、その動作を指示する、順序付けられた動作の特定のセットを提供する、マイクロコード、プログラム、ルーチン、またはコンピュータで動作可能な記号言語の形態であってもよいことに留意されたい。プログラムコード６０の命令は、コンピュータ３０のメモリ５４中、またはコンピュータディスケット、磁気テープ、従来のハードディスクドライブ、電子読み取り専用メモリ、光記憶デバイス、もしくはその上に記憶されるコンピュータ可用媒体を有する他の適切なデータ記憶デバイス上に記憶されてもよい。プログラムコード６０はまた、示されるように、コンピュータ可読媒体として、サーバ５８等のデータ記憶デバイス上に含まれてもよい。
【００５５】
コンピュータ３０で実行される本発明の方法は、メモリ５４中に記憶され、コンピュータ３０のシステムプロセッサ４０によって実行可能である、図１のコンピュータプログラムステップを利用して実施され得る。処理システムＤへの入力データは、検層データおよび上記の貯留層に関する他のデータである。
【００５６】
処理試験データの結果
Ｘが０．１の間隔でサンプルされる理論関数
【数８】
0006006799.tif 000019
を所与として、対応するｙ数列が、図３に対する凡例によって識別されるように、理論値に対して図３にプロットされる。
【００５７】
数学的に定義されるが非線形である、方程式（８）で表される関数に係る合成データを使用した試験モデルにおいて、関数中４回に１回のｘ点を使用して、ＳＶＭをトレーニングした。次いで、ＳＶＭは、残りのｘ点を予測しようと試みた。次いで、その一致がどれほど正確であるかを決定するために、予測された点と理論点との間の比較を行った。
【００５８】
図３は、０〜４πの４ｘ点ごとのトレーニングデータを用いた、第１のかかるＳＶＭ合成試験のプロットである。図３中、理論データ点は、上記のように凡例が示され、トレーニングサブセットは、凡例において識別されるように示され、予測されたＳＶＭ回帰点もまた、凡例において識別されるように示される。
【００５９】
明らかなように、相当数の場所におけるＳＶＭ回帰点が、理論データ点との実質的に正確な重なりのため見えない。ＳＶＭ回帰点のほとんどが理論データ点と実質的に重なり合うため、ＳＶＭ予測の有効性は推定され得ることに留意されたい。
【００６０】
第１の試験と同一の方法で、同一のデータセットに対して、第２の試験を行った。しかしながら、ＳＶＭ回帰のロバスト性をチェックするために、その中のトレーニングデータ点を０〜２πと範囲を半分に減少させた。図４は、０〜２πの４ｘ点ごとのトレーニングデータを用いた、第２のＳＶＭ合成試験のプロットである。図４において、理論データ点は、図４に対する凡例において識別されるように示され、トレーニングサブセットおよび予測されたＳＶＭ回帰点もまた、識別されるように示される。この場合もやはり、相当数の場合における予測されたＳＶＭ回帰点が、理論データ点との実質的に正確な重なりのため見えない。
【００６１】
試験範囲内の予測は、理論データに十分良好に適合したことが観察される。しかしながら、トレーニング範囲外の点に対して、その予測は、十分に近似してはいなかった。これは、方程式（８）の式が線形関係ではないため、予想外ではない。かかる非線形関係は、より完全なトレーニングデータセットを必要とすることが予想され得る。さらに、両方の試験において入力される特徴はたった１つ、すなわちｘであった。貯留層特性予測において、ＳＶＭモデリングに対して複数の属性（特徴）空間が存在した。第１および第２の合成データ試験の結果は、多重地震属性を用いた実用的な貯留層特性予測試験を実行するための実証および十分な基礎を提供した。
【００６２】
実際の貯留層データの処理結果
既知の油田における８つの地震属性の３次元ボリュームをトレーニングデータ入力として提供した。トレーニング井戸ターゲット入力のソースとして使用した、掘削された井戸は３つであった。各井戸は、異なる深度で貯留層空隙率特性を記録した。９つの属性から、適切に十分な数の地震属性を選択し、次いで、ＳＶＭモデルを本発明に従ってトレーニングし、次いで、井戸が存在しない場所でさえ、貯留層の全体積にわたって、本発明に従って、予測された貯留層特性（空隙率）を決定した。
【００６３】
属性は、振幅、周波数、第１のエンベロープおよび第２のエンベロープ、位相、コヒーレンス、帯域幅、ならびに音響インピーダンスであった。最初に、３つの井戸の位置で属性を選択した。属性と空隙率との間の相互相関係数を計算した。具体的には、すなわち、相互相関係数は、０．０５１、０．００６７、０．０９９４、０．０１２４、０．０４４２、０．１８１９、０．０８２、および０．７０９であった。０．０５の拒否閾値を選択し、その結果、５つの属性が、トレーニングおよび試験に対する相互相関閾値に達した。５つの属性は、振幅、周波数、コヒーレンス、帯域幅、および音響インピーダンスであった。
【００６４】
相互検証に対して、各井戸をサブセットとして扱った。２回の検索および検証の後、最終γは４２．２であり、最終Ｃは４１９４３０４であった。０．２５〜０．３５の空隙率範囲と比較して、選択されたγおよびＣを用いた最小相互検証フィッティング誤差は、０．０００４７２であった。予測の検証をチェックするために、３つの井戸の位置で予測された空隙率を比較し、図５、図６、および図７に示した。
【００６５】
図５は、貯留層中の既知の貯留層位置における井戸に対する予測された空隙率予測のプロットである。図５のプロットにおける、井戸データから得られる測定された空隙率に対するデータ点は、凡例によって示されるように識別され、本発明に係るＳＶＭモデリングから予測された空隙率を表すデータ点もまた、その図に対する凡例に従って示される。
【００６６】
図６は、図５の井戸とは異なるが、同一の貯留層中の井戸に対する、予測された空隙率予測のプロットである。この場合もやはり、図６の井戸データからの測定された空隙率に対するデータ点は、凡例によって識別されるように示され、本発明に係るＳＶＭモデリングから予測された空隙率を表すデータ点も、凡例に従って示される。
【００６７】
図７は、図５および６の井戸とは異なるが、同一の貯留層中の井戸に対する、予測された空隙率予測のプロットである。この場合もやはり、井戸データから得られる測定された空隙率に対するデータ点は、凡例に従って示され、本発明に係るＳＶＭモデリングから予測された空隙率に対するデータ点もまた、凡例に従って示される。
【００６８】
図５、６、および７の予測された空隙率点が、実際の空隙率と正確に一致はしないが、測定された空隙率と自然に適合することに留意されたい。本発明では、正確な一致はしないが、かかる自然な適合においてこれらの例と比較できる予測結果を得ることが予想される。このように、井戸が存在しない位置での貯留層特性値を結果として予測し損なうことなく、オーバーフィッティングが回避され得る。
【００６９】
図８は、図５、６、および７に示される結果を得るために、地震調査中に得られ、上記の処理中に５つの属性のうちの１つとして使用される、音響インピーダンスの地震属性入力データの値の表示部である。図８の音響インピーダンスデータは、実際の貯留層にわたる断面の地震測線に沿ってプロットされる。音響インピーダンスデータはまた、地震測線の約一千のトレースにわたる時間／深度の関数としてプロットされ、実際の表示部における属性値は、色の変化によって示されるように、色で示される。図８の最上部分８０は、より高い音響インピーダンスの下部分８２に向かって下方に延在する、より低い音響インピーダンスの領域を表す。
【００７０】
図９は、貯留層にわたる図８と同一の断面の測線に沿って、図５、６、および７に示される結果を得るために使用される他の４つの属性と共に、図８の音響インピーダンス属性データを使用して、本発明に係るＳＶＭモデリングの結果として得られる貯留層特性としての、空隙率の予測された値の表示部である。図９の予測された空隙率はまた、図８と同一のトレースにわたる時間／深度の関数としてプロットされ、実際の実行における予測された空隙率値は、色の変化によって示されるように、色で示される。図９の上部分９０は、高い空隙率の領域を表し、一方で、図９の下部分９２は、より低い空隙率の領域を示す。
【００７１】
本発明に係るＳＶＭモデリングにおける相関処理の間、音響インピーダンスが空隙率との非常に強い逆相関を有し、それが、より低い音響インピーダンスが、より高い空隙率に対応し、逆もまた同様であることを意味することに留意された。図９の予測された空隙率が、予想されたように、図８の音響インピーダンスと同様の逆相関関係を示すことに留意されたい。本発明に従って得られる、貯留層の一部分に対する図９に示される空隙率等、貯留層に対する最終的な完全な予測された貯留層特性データは、データメモリに記憶され、図９のようなプロットは、貯留層特徴のさらなる研究および評価のために、貯留層分析者に利用可能である。
【００７２】
第２の実施例として、１１個の地震属性が既知の油田から得られた別の実際の貯留層３次元ボリュームが、本発明に係るＳＶＭモデリングの対象であった。１１個の属性はすなわち、瞬間振幅、瞬間帯域幅、第１のエンベロープトレース導関数、瞬間周波数、周波数加重エンベロープ、正規化振幅、瞬間位相、瞬間Ｑファクタ、第２のエンベロープトレース導関数、地震振幅、および薄層であった。１１個の属性からの地震データをトレーニング入力として提供した。この貯留層に対する空隙率値を有する井戸はまた、１０個あった。音響インピーダンス属性データは利用可能ではなく、したがって、３つの他の属性が相互相関選択後に選択された。選択された属性は、周波数、正規化振幅、および位相であった。０．１の係数閾値を拒否閾値として使用した。１０個の井戸を相互検証のための４つのサブセットに割り当てた。ＳＶＭモデリングに対して、最終γカーネルパラメータは１０２４であり、最終ペナルティパラメータＣは１８２５６７７であった。０．０９〜０．１１の正常な空隙率値に対して、選択されたγおよびＣを用いた最小相互検証フィッティング誤差は０．００４９であった。空隙率の予測された空隙率貯留層特性が、実際の貯留層にわたる断面の地震測線に沿って図１１に示される一方で、同一線に沿った周波数の地震属性は、図１０に示される。図１１の空隙率表示部が、図９の周波数の入力属性と同様の構造詳細を示すことに留意されたい。予測のさらなる検証は、貯留層シミュレーションまたは追加の検層データによって確認され得る。
本発明では、ＳＶＭモデリングは、貯留層特性予測の目的で回帰を含むように修正されている。記載されたように、本発明に係るコンピュータ処理は、相互相関による属性選択、データ再スケーリングおよび単一化、ＲＢＦカーネルに対するグリッド検索および相互検証の方法によるパラメータ最適化、ならびにＳＶＭモデリングおよび予測を含む。
【００７３】
炭化水素貯留層予測の管理された回帰の問題に対して、ニューラルネットワーク等の従来の学習方法と比較して、ＳＶＭには３つの利点がある。
【００７４】
ＳＶＭは、それにオーバーフィッティングを回避させる構造リスク最小化の原理に基づき、またＳＶＭは、それに優れた一般化能力を有させる経験的リスクおよび信頼限界を考慮に入れる。ＳＶＭは、局所的最適解の代わりに、大域的最適解を与えることができる。
【００７５】
ＳＶＭは、非線形問題を解決するためにカーネルマシンを使用する。低次元特徴空間における非線形問題は常に、何らかの高次元特徴空間における線形問題に変換され得る。また、その一般化能力は、入力空間の次元性ではなく、トレーニングデータの数による。
【００７６】
本発明は、当該分野の平均的な知識を有する者が、本発明において言及された結果を再現および獲得できるように十分に説明されている。それにもかかわらず、本発明の当業者は、本出願において記載されていない修正を行ってもよく、決定された構造にこれらの修正を適用するために、またはその製造プロセスにおいて、以下の特許請求の範囲の請求された事項を必要とし、かかる構造が、本発明の範囲内に含まれるものとする。
【００７７】
添付の特許請求の範囲に記載される発明の精神または範囲から逸脱することなく、改良および修正が上記に詳述される本発明に行われ得ることが、留意および理解されるべきである。
【図１】
0006006799.tif 000020
【図２】
0006006799.tif 000021
【図３】
0006006799.tif 000022
【図４】
0006006799.tif 000023
【図５】
0006006799.tif 000024
【図６】
0006006799.tif 000025
【図７】
0006006799.tif 000026
【図８】
0006006799.tif 000027
【図９】
0006006799.tif 000028
【図１０】
0006006799.tif 000029
【図１１】
0006006799.tif 000030