ガス水和物を含む堆積物のＩＮＳＩＴＵ機械特性の推定
【課題】多成分粒子系の弾塑性特性相関をとる方法を提供する。
【解決手段】堆積物の水和物系の地球物理学的データから、ｓｈｙｄおよびｐｃを有するパラメータを取得する工程を含み、ここで、ｓｈｙｄは水和物飽和度で、ｐｃは封圧である。さらに、前記パラメータおよび相関を使って前記堆積物の水和物系のヤング率（Ｅ）を決定する工程を含み、ここで、前記相関は
であり、α＝ｐｃ／ｃで、ｃは圧力の次元を伴ったスケール変数である。当該方法は、さらに、前記堆積物の水和物系のヤング率に基づいて現場作業を調整する工程を含む。
【解決手段】堆積物の水和物系の地球物理学的データから、ｓｈｙｄおよびｐｃを有するパラメータを取得する工程を含み、ここで、ｓｈｙｄは水和物飽和度で、ｐｃは封圧である。さらに、前記パラメータおよび相関を使って前記堆積物の水和物系のヤング率（Ｅ）を決定する工程を含み、ここで、前記相関は
であり、α＝ｐｃ／ｃで、ｃは圧力の次元を伴ったスケール変数である。当該方法は、さらに、前記堆積物の水和物系のヤング率に基づいて現場作業を調整する工程を含む。
【発明の詳細な説明】
【技術分野】
【０００１】
本願は、２００８年６月２日付けで出願された米国仮特許出願第６１／０５８，１５５号「Ｍｅｔｈｏｄ ａｎｄ Ｓｙｓｔｅｍ ｆｏｒ Ｅｓｔｉｍａｔｉｎｇ Ｉｎ Ｓｉｔｕ Ｍｅｃｈａｎｉｃａｌ Ｐｒｏｐｅｒｔｉｅｓ ｏｆ Ｓｅｄｉｍｅｎｔｓ Ｃｏｎｔａｉｎｉｎｇ Ｇａｓ Ｈｙｄｒａｔｅｓ」（堆積物を含むガス水和物のＩｎ Ｓｉｔｕ機械特性を推定する方法およびシステム）に基づく利益を主張するものである（参照により本明細書に組み込むものとする）。
【０００２】
さらに、本願は、２００８年１２月３０日付けで交付された米国特許第７，４７２，０２２号「Ｍｅｔｈｏｄ ａｎｄ Ｓｙｓｔｅｍ ｆｏｒ Ｍａｎａｇｉｎｇ ａ Ｄｒｉｌｌｉｎｇ Ｏｐｅｒａｔｉｏｎ ｉｎ ａ Ｍｕｌｔｉｃｏｍｐｏｎｅｎｔ Ｐａｒｔｉｃｕｌａｔｅ Ｓｙｓｔｅｍ」（多成分粒子系における掘削作業を管理する方法およびシステム）に関連する（参照により本明細書に組み込むものとする）。
【背景技術】
【０００３】
図１は掘削現場作業を示した図であり、掘削リグ（１００）を使って、ボアホール（ボーリング孔）（１４５）内のドリルパイプ（１４０）遠端に連結されたドリルビット（１５０）を回転させている。この現場作業により、石油、天然ガス、水、または掘削で取得可能なまたは取得される他の任意タイプの物質を取得することができる。図１に示した現場作業は地表から地層を直接掘削していくものであるが、当業者であれば、湖底の掘削、深海底の掘削、探鉱作業、開発作業、仕上げ作業、生産作業、処理作業など、他タイプの現場作業が存在することも理解されるであろう。
【０００４】
図１に示すように、ロータリーテーブル（１２５）で生成される回転力は、前記掘削リグ（１００）から前記ドリルパイプ（１４０）を介して前記ドリルビット（１５０）へ伝達される。さらに、掘削流体（「マッド（ｍｕｄ）」とも呼ばれる）が前記ドリルパイプ（１４０）の中空コアを通じて前記ドリルビット（１５０）へ送られる。具体的には、マッドポンプ（１８０）を使い、スタンドパイプ（１６０）とホース（１５５）とケリー（１２０）とを通じて、前記ドリルパイプ（１４０）へマッドを送る。暴噴の危険性を軽減するには、防噴装置（１３０）を使って、前記ボアホール（１４５）内の流体圧力を制御する。さらに、ケーシング（１３５）を使って前記ボアホール（１４５）を強化することにより、暴噴、または当該ボアホール（１４５）に作用する他の力による倒壊を防ぐことができる。当該掘削リグ（１００）には、クラウンブロック（１０５）、トラベリングブロック（１１０）、スイベル（１１５）、および他の構成要素（図示せず）を含めることもできる。
【０００５】
前記ボアホール（１４５）から地表へ返ってくるマッドは、マッド戻り配管（１６５）を通じてマッド処理用の機器へ送られる。例えば、マッドは、掘削された固体を取り除くよう構成されたシェーカー（１７０）へ送られる。除去された固体はリザーブ（貯留）ピット（１７５）へ移動され、マッドはマッドピット（１９０）に堆積する。前記マッドポンプ（１８０）は、ろ過されたマッドを前記マッドピット（１９０）からマッド吸引管（１８５）で吸引し、前記掘削リグ（１００）内へ再注入する。当業者であれば、デガッサー、デサンダー、デシルター、遠心分離機、混合ホッパーなど他のマッド処理装置も使用できることが理解されるであろう。さらに、掘削作業には、流体技術実施、掘削シミュレーション、圧力制御、抗井からの不要物除去、および廃棄物管理などの作業に使用される他のタイプの掘削構成要素が含まれることがある。
【０００６】
所与の現場作業では（図１に示した掘削作業など）、地層の地質力学特性に関する知識を使って、現場に関連した種々の問題を軽減することができる。例えば、一部の地層では、岩石が変形し若しくは破壊されるおそれがある。塑性パラメータは、通常、コアに機械的試験を行って直接測定されるが、数学的に相関を構築する（相関をとる）と脆性−弾性特性を予測することができる。
【０００７】
例えば、静的ヤング率を推定する方法は、Ｅｉｓｓａ，Ｅ．Ａ．& Ｋａｚｉ，Ａ．らの「Ｒｅｌａｔｉｏｎ ｂｅｔｗｅｅｎ ｓｔａｔｉｃ ａｎｄ ｄｙｎａｍｉｃ Ｙｏｕｎｇ’ｓ ｍｏｄｕｌｉ ｆｏｒ ｒｏｃｋｓ」（岩石に関する静的ヤング率および動的ヤング率の関係）（Ｉｎｔ．Ｊ．Ｒｏｃｋ Ｍｅｃｈ．Ｍｉｎ．Ｓｃｉ．&Ｇｅｏｍｅｃｈ．Ａｂｓｔｒ．２５（１９８８年）、ページ４７９−４８２）、Ｍｏｎｔｍａｙｏｕｒ，Ｈ．& Ｇｒａｖｅｓ，Ｒ．Ｍ．らの「Ｐｒｅｄｉｃｔｉｏｎ ｏｆ Ｓｔａｔｉｃ Ｅｌａｓｔｉｃ／Ｍｅｃｈａｎｉｃａｌ Ｐｒｏｐｅｒｔｉｅｓ ｏｆ Ｃｏｎｓｏｌｉｄａｔｅｄ ａｎｄ Ｕｎｃｏｎｓｏｌｉｄａｔｅｄ Ｓａｎｄｓ Ｆｒｏｍ Ａｃｏｕｓｔｉｃ Ｍｅａｓｕｒｅｍｅｎｔｓ： Ｃｏｒｒｅｌａｔｉｏｎｓ」（音響測定による固結砂および未固結砂の静的な弾性特性／機械特性の予測：相関）（６１ｓｔ Ａｎｎｕａｌ Ｔｅｃｈｎｉｃａｌ Ｃｏｎｆｅｒｅｎｃｅ ａｎｄ Ｅｘｈｉｂｉｔｉｏｎ ｏｆ ｔｈｅ Ｓｏｃｉｅｔｙ ｏｆ Ｐｅｔｒｏｌｅｕｍ Ｅｎｇｉｎｅｅｒｓ、Ｎｅｗ Ｏｒｌｅａｎｓ、ＬＡ．（１９８６年）ＳＰＥ １５６４４）、Ｍｏｒａｌｅｓ，Ｒ．Ｈ．& Ｍａｒｃｉｎｅｗ，Ｒ．Ｐ．らの「Ｆｒａｃｔｕｒｉｎｇ ｏｆ ｈｉｇｈ−ｐｅｒｍｅａｂｉｌｉｔｙ ｆｏｒｍａｔｉｏｎｓ： Ｍｅｃｈａｎｉｃａｌ ｐｒｏｐｅｒｔｉｅｓ ｃｏｒｒｅｌａｔｉｏｎｓ」（浸透率の高い地層の破砕：機械特性の相関）（ＳＰＥ ２６５６１（１９９３年））、Ｙａｌｅ，Ｄ．Ｐ．& Ｊａｍｉｅｓｏｎ，Ｗ．Ｈ．らの「Ｓｔａｔｉｃ ａｎｄ Ｄｙｎａｍｉｃ Ｒｏｃｋ Ｍｅｃｈａｎｉｃａｌ Ｐｒｏｐｅｒｔｉｅｓ ｉｎ ｔｈｅ Ｈｕｇｏｔｏｎ ａｎｄ Ｐａｎｏｍａ Ｆｉｅｌｄｓ， Ｋａｎｓａｓ」（米国カンザス州Ｈｕｇｏｔｏｎ ＦｉｅｌｄおよびＰａｎｏｍａ Ｆｉｅｌｄにおける静的および動的な岩石機械特性）（ＳＰＥ Ｍｉｄ−Ｃｏｎｔｉｎｅｎｔ Ｇａｓ Ｓｙｍｐｏｓｉｕｍ，Ａｍａｒｉｌｌｏ，ＴＸ．、ＳＰＥ ２７９３９（１９９４年））、およびＴｕｔｕｎｃｕ，Ａ．Ｎ．& Ｓｈａｒｍａ，Ｍ．Ｍ．らの「Ｒｅｌａｔｉｎｇ Ｓｔａｔｉｃ ａｎｄ Ｕｌｔｒａｓｏｎｉｃ Ｌａｂｏｒａｔｏｒｙ Ｍｅａｓｕｒｅｍｅｎｔｓ ｔｏ Ａｃｏｕｓｔｉｃ Ｌｏｇ Ｍｅａｓｕｒｅｍｅｎｔｓ ｉｎ Ｔｉｇｈｔ Ｇａｓ Ｓａｎｄｓ」（タイトな（硬質）ガス砂における静的実験および超音波実験の測定値の音響検層測定値との関係付け）（６７ｔｈ Ａｎｎｕａｌ Ｔｅｃｈｎｉｃａｌ Ｃｏｎｆｅｒｅｎｃｅ ａｎｄ Ｅｘｈｉｂｉｔｉｏｎ ｏｆ ｔｈｅ Ｓｏｃｉｅｔｙ ｏｆ Ｐｅｔｒｏｌｅｕｍ Ｅｎｇｉｎｅｅｒｓ，Ｗａｓｈｉｎｇｔｏｎ， ＤＣ、ＳＰＥ ２４６８９（１９９２年））の諸論文に説明されている。
【０００８】
さらに、動的ポアソン比から静的ポアソン比を推定する方法は、Ｔｕｔｕｎｃｕ&Ｓｈａｒｍａ（１９９２）（上記で参照）およびＹａｌｅ&Ｊａｍｉｅｓｏｎ（１９９４）（上記で参照）で説明されている。ビオ数を評価するための経験的相関については、Ｋｒｉｅｆ，Ｍ．、Ｇａｒａｔ，Ｊ．、Ｓｔｅｌｌｉｎｇｗｅｒｆｆ，Ｊ．、Ｖｅｎｔｒｅ，Ｊ．らの論文「Ａ ｐｅｔｒｏｐｈｙｓｉｃａｌ ｉｎｔｅｒｐｒｅｔａｔｉｏｎ ｕｓｉｎｇ ｔｈｅ ｖｅｌｏｃｉｔｉｅｓ ｏｆ Ｐ ａｎｄ Ｓ ｗａｖｅｓ （ｆｕｌｌ−ｗａｖｅｆｏｒｍ ｓｏｎｉｃ）」（Ｐ波およびＳ波（完全波形の音波）（Ｔｈｅ Ｌｏｇ Ａｎａｌｙｓｔ ３１、１９９０年１１月、ページ３５５−３６９）の速度を使った解釈）に説明されている。
【０００９】
さらに、岩石の一軸圧縮強度を推定するための相関をとる方法が複数の著者により考案されており、Ｃｈａｎｇ，Ｃ．らの論文「Ｅｍｐｉｒｉｃａｌ Ｒｏｃｋ Ｓｔｒｅｎｇｔｈ Ｌｏｇｇｉｎｇ ｉｎ Ｂｏｒｅｈｏｌｅｓ Ｐｅｎｅｔｒａｔｉｎｇ Ｓｅｄｉｍｅｎｔａｒｙ Ｆｏｒｍａｔｉｏｎｓ」（堆積岩層を貫通したボアホールに関する経験的岩石強度の孔内検層）（ＭＵＬＬＩ−ＴＡＭＳＡ（Ｇｅｏｐｈｙｓｉｃａｌ Ｅｘｐｌｏｒａｔｉｏｎ）７、２００４年、ページ１７４−１８３）でレビューされている。この目的での付加的な相関については、Ｐｌｕｍｂ，Ｒ．Ａ．、Ｈｅｒｒｏｎ Ｓ．Ｌ．&Ｏｌｓｅｎ，Ｍ．Ｐ．「Ｉｎｆｌｕｅｎｃｅ ｏｆ Ｃｏｍｐｏｓｉｔｉｏｎ ａｎｄ Ｔｅｘｔｕｒｅ ｏｎ Ｃｏｍｐｒｅｓｓｉｖｅ Ｓｔｒｅｎｇｔｈ Ｖａｒｉａｔｉｏｎｓ ｉｎ ｔｈｅ Ｔｒａｖｉｓ Ｐｅａｋ Ｆｏｒｍａｔｉｏｎ」（組成およびテクスチャーがＴｒａｖｉｓ Ｐｅａｋ地層の圧縮強度変動に及ぼす影響）（６７ｔｈ Ａｎｎｕａｌ Ｔｅｃｈｎｉｃａｌ Ｃｏｎｆｅｒｅｎｃｅ ａｎｄ Ｅｘｈｉｂｉｔｉｏｎ ｏｆ ｔｈｅ Ｓｏｃｉｅｔｙ ｏｆ Ｐｅｔｒｏｌｅｕｍ Ｅｎｇｉｎｅｅｒｓ，Ｗａｓｈｉｎｇｔｏｎ，Ｄ．Ｃ．（１９９２年）ＳＰＥ ２４７５８）およびＱｉｕ，Ｋ．、Ｍａｒｓｄｅｎ，Ｊ．Ｒ．、Ｓｏｌｏｖｙｏｖ，Ｙ．、Ｓａｆｄａｒ，Ｍ．& Ｃｈａｒｄａｃ，Ｏ．「Ｄｏｗｎｓｃａｌｉｎｇ Ｇｅｏｍｅｃｈａｎｉｃｓ Ｄａｔａ ｆｏｒ Ｔｈｉｎ Ｂｅｄ Ｕｓｉｎｇ Ｐｅｔｒｏｐｈｙｓｉｃａｌ Ｔｅｃｈｎｉｑｕｅｓ」（石油物理学的技術を使った薄層理に関する地質力学データのダウンスケーリング）（１４ｔｈ ＳＰＥ Ｍｉｄｄｌｅ Ｅａｓｔ Ｏｉｌ ａｎｄ Ｇａｓ Ｓｈｏｗ ａｎｄ Ｃｏｎｆｅｒｅｎｃｅ， Ｂａｈｒａｉｎ（２００５年）ＳＰＥ ９３６０５）の論文で説明されている。
【発明の概要】
【課題を解決するための手段】
【００１０】
多成分粒子系の弾塑性特性相関をとる方法。当該方法は、堆積物の水和物系の地球物理学的データから、ｓｈｙｄおよびｐｃを有するパラメータを取得する工程を含み、ここで、ｓｈｙｄは水和物飽和度で、ｐｃは封圧である。当該方法は、さらに、前記パラメータおよび相関を使って前記堆積物の水和物系のヤング率（Ｅ）を決定する工程を含み、ここで、前記相関は
【数１】
【００１１】
であり、α＝ｐｃ／ｃで、ｃは圧力の次元を伴ったスケール変数である。当該方法は、さらに、前記堆積物の水和物系のヤング率に基づいて現場作業を調整する工程を含む。
【００１２】
ガス水和物を含む堆積物のｉｎ ｓｉｔｕ（測定対象を動かさず、原位置で）機械特性を推定する工程の他の態様は、以下の説明および添付の請求項から明確に理解されるであろう。
【図面の簡単な説明】
【００１３】
【図１】図１は、現場作業の図を示したものである。
【図２】図２は、１若しくはそれ以上の実施形態に係るシステムを示した図である。
【図３】図３は、１若しくはそれ以上の実施形態に従って現場作業を管理する方法のフローチャートを示したものである。
【図４】図４は、１若しくはそれ以上の実施形態に係る相関図を示したものである。
【図５】図５は、１若しくはそれ以上の実施形態に係るコンピュータシステムを示した図である。
【発明を実施するための形態】
【００１４】
以下、添付の図面を参照して、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ（測定対象を動かさず、原位置で）機械特性を推定する具体的な実施形態を詳しく説明する。種々の図にわたる同様な要素は、一貫性のため同様な参照番号で示している。
【００１５】
以下の詳細な説明では、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態について、それらの実施形態がより完全に理解されるよう具体的な詳細事項を多数記載する。ただし当業者であれば、それらの具体的な詳細事項がなくとも前記実施形態が実施可能であることは明確に理解されるであろう。他の場合、周知の特徴（機能）については、説明を不要に複雑化しないよう詳しい説明を省いている。
【００１６】
一般に、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態は、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を地震観測データまたは検層データから推定する方法およびシステムに関し、掘削作業中、仕上げ作業中、および生産作業中の機械的破壊の危険性を評価する上で極めて重要である。より具体的には、１若しくはそれ以上の実施形態は、ガス水和物含有堆積物の機械特性と地球物理学的データとの相関をとる方法に関する。
【００１７】
一般に、前記方法の実施形態では、マイクロメカニクスモデルを使った技術を適用して堆積物の物理的挙動を決定するパラメータの選択を誘導する。上述の選択に基づき、弾塑性特性（従属変数）と、検層データまたは地震観測データから推定できる特性（独立変数）との無次元化した関係のセットが導出される。前記方法では、これらの関係を使って、従属変数に関する弾塑性特性の代数展開を提案し、多種多様な情報源のデータを使って、メタンを含む砂質土のヤング率（Ｅ）と、ＴＨＦ（テトラヒドロフラン）水和物との相関を導出する。当業者であれば、上記の相関を導出する技術は、メタンを伴った砂質土およびＴＨＦ水和物だけに限定されず、粒子状物質全般に拡張できることが理解されるであろう。当業者であれば、前記技術を使用すると、機械強度、摩擦角、ダイレーション角（ダイレイタンシー角）に関する相関など、粒子状物質に関する他種の弾塑性特性相関を導出できることも理解されるであろう。
【００１８】
本発明の実施形態は、全般的に、多成分粒子系（２種若しくはそれ以上の粒子を含んだ地層）における現場作業（掘削作業、測量作業、仕上げ作業、生産作業など）を管理する方法およびシステムを提供する。ただし、本願で説明する方法は、単一成分系に適用することもできる。粒子系の地球物理学的特性、応力特性、および石油物理学的特性は、ｉｎ ｓｉｔｕ（測定対象を動かさず、原位置で）または実験室での測定値から得られる。適切な特性の選択は、マイクロメカニクスにより推定される。本明細書において、地球物理学的特性および石油物理学的特性とは、一般に岩石フレーム（構造）に関するものであるが、個々の粒子タイプの特性を含めていう場合もある。これらの特性は、適切なスケーリングにより無次元数で表せる。また、粒子系の排水条件での静的弾塑性変形を特徴付ける特性も同様に得られ、無次元形態に変換できる。弾塑性特性は、１若しくはそれ以上の経験的相関を介して地球物理学的特性、応力特性、および石油物理学的特性に関係している。前記経験的相関は、堆積系のｉｎ ｓｉｔｕ弾塑性特性を推定するため使用される。これらの特性は、岩石の変形および破壊モデルにおいて、堆積系で掘削されたボアホール周辺の岩石の機械的完全性を評価するため使用され、この評価に基づき現場作業が調整される。
【００１９】
マイクロメカニクスは複合系に関する学問であり、そのような系の機械的挙動を、構造部分の挙動により説明することを狙いとしている。複合系とは、２若しくはそれ以上の部分から成るものをいう。
【００２０】
フレーム特性は、多孔質集合体に関連するもので、その多孔質集合体内に常駐する流体（間隙流体）は除く。フレーム特性を測定する際は、試験前に間隙流体が除去され、または間隙流体の影響を最小限に抑えるような態様で試験が行われる。
【００２１】
静的特性は、一定の荷重または緩慢に変化する荷重を物質に与えた場合の変形を特徴付けるため使用される。これと対照的に、動的特性は、急速に変化する荷重を物質が受けた場合の変形を特徴付けるものである。
【００２２】
排水特性（排水条件下の特性）は、物質の変形を特徴付けるもので、その場合、当該変形は十分緩慢であるため、当該変形に影響を及ぼすことなく間隙流体が試料から流出できる。
【００２３】
弾塑性特性は、弾塑性物質の応力−ひずみ挙動をモデル化するため使用される特性である。弾塑性物質は、ひずみが小さいと弾性応答を呈し、ひずみが大きいと塑性応答を呈する。弾性変形中、当該物質に永続的な損傷は生じない。これは、当該物質への荷重が取り除かれると、荷重をかける元の状態に戻ることを意味している。ただし、荷重が十分大きくなると当該物質は塑性を呈し、永続的な損傷が生じて、荷重を取り除いても元の状態には戻らなくなる。弾塑性特性は、弾性変形および塑性変形の双方の体系において物質の応力−ひずみ応答をモデル化する構成則（構成式）に使用される。
【００２４】
ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態は、種々の多成分粒子系について実施できる。例えば、現場作業では、クラスレート（包接化合物）水和物（ハイドレート）（当該技術分野では、ガスクラスレート水和物、ガス水和物、またはクラスレートとも呼ばれる）を含んだ地層を標的にすることがある。クラスレート水和物は、水分子の格子から構成された氷状の結晶性固体である。この格子にガスまたは液体の分子が取り込まれて、この水和物構造を安定なものにしている。メタンガス水和物（メタンハイドレート）は、クラスレート水和物として自然発生する最も一般的な種である。ただし、プロパンやエタンなど他の炭化水素を含んだクラスレート水和物も存在する。さらに、テトラヒドロフラン（ｔｅｔｒａｈｙｄｒｏｆｕｒａｎ、略称ＴＨＦ）、二酸化炭素（ＣＯ２）、および硫化水素（Ｈ２Ｓ）などの非炭化水素物質が水和物構造に組み込まれることもある。クラスレート水和物は、永久凍土地域および大陸縁辺の深海環境に広く存在してる。
【００２５】
クラスレート水和物ゾーンの掘削に対する関心は増している。例えば、クラスレート水和物は、エネルギー源として使用することができる。しかし、クラスレート水和物ゾーンでの掘削および生産には種々の課題が存在する。掘削工程に起因した温度および圧力の外乱によりクラスレート水和物中のガスが解離し、抗井中への放出が制御できず、火災または暴噴が起こるおそれがある。さらに、解放されたガスが掘削用マッドをガス化させるおそれもある。場合により、解離中のクラスレート水和物を含んだ堆積層部分の崩壊で抗井が不安定になり、結果的に穴または枝掘り部が失われてしまうおそれもある。例えば、Ｃｏｌｌｅｔｔ，Ｔ．Ｓ．& Ｄａｌｌｉｍｏｒｅ，Ｓ．Ｒ．（２００２年）「Ｄｅｔａｉｌｅｄ ａｎａｌｙｓｉｓ ｏｆ ｇａｓ ｈｙｄｒａｔｅ ｉｎｄｕｃｅｄ ｄｒｉｌｌｉｎｇ ａｎｄ ｐｒｏｄｕｃｔｉｏｎ ｈａｚａｒｄｓ」（ガス水和物により誘発される掘削および生産時の危険に関する詳細解析）（４ｔｈ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｇａｓ Ｈｙｄｒａｔｅｓ（ＩＣＧＨ−４）、Ｙｏｋｏｈａｍａ、２００２年５月）を参照。その影響が壊滅的でない場合でさえ、穴の状態が良好でなければクラスレート水和物ゾーンの孔内検層結果の質が落ちるおそれがある。
【００２６】
クラスレート水和物ゾーンにおける坑井の完成およびそこでの生産にも課題が存在する。穴の状態が良好でないとセメントが効果的に結合（ボンディング）せず、ケーシング外部にガスが漏れてしまう。さらに、クラスレート水和物を除去したことにより地層の維持能力が失われると、砂質土が生じ、地層が地盤沈下を起こし、またはケーシングが破壊されるおそれがある。また、クラスレート水和物ゾーンで完成した坑井内では、高温の炭化水素により抗井ケーシングの周囲で水和物の解離が生じ、地層への圧力が高まるおそれがある。これについては、Ｆｒａｎｋｌｉｎ，Ｌ．Ｊ．（１９８１年）「Ｈｙｄｒａｔｅｓ ｉｎ Ａｒｔｉｃ Ｉｓｌａｎｄｓ」（北極諸島における水和物）（Ａ．Ｌ． Ｂｒｏｗｓｅｒ（ｅｄ．）Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ａ Ｗｏｒｋｓｈｏｐ ｏｎ Ｃｌａｔｈｒａｔｅｓ（ｇａｓ ｈｙｄｒａｔｅｓ）ｉｎ ｔｈｅ Ｎａｔｉｏｎａｌ Ｐｅｔｒｏｌｅｕｍ Ｒｅｓｅｒｖｅ ｉｎ Ａｌａｓｋａ．ＵＳＧＳ Ｏｐｅｎ−Ｆｉｌｅ Ｒｅｐｏｒｔ ８１−１２５９、１８−２１）を参照。
【００２７】
当業者であれば、上記で具体的に説明した課題が単なる例であり、クラスレート水和物ゾーンでの掘削には他の課題が多種存在することが理解されるであろう。さらに、本明細書の例ではクラスレート水和物ゾーンでの現場作業について強調したが、他種の多成分粒子系でも同様な課題が存在する。また、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態は、単一成分系にも適用でき、２より多くの成分を含む多成分粒子系に拡張することもできる。
【００２８】
図２は、１若しくはそれ以上の実施形態に係るシステム（２００）の図を示したものである。このシステム（２００）には、掘削対象である多成分粒子系（２２０）の特性を測定するよう構成された１若しくはそれ以上の測定機構（測定機構Ａ（２１０）、測定機構Ｎ（２１５）など）が含まれる。測定は、実験室試料に対し、またはｉｎ ｓｉｔｕで（測定対象を動かさず、原位置で）行うことができる。あるいは、前記多成分粒子系（２２０）のシミュレーションから、または類似組成を有する実際の地層若しくはシミュレートした地層から測定値を得ることもできる。前記多成分粒子系（２２０）は２若しくはそれ以上のタイプの粒子を含むため、１若しくはそれ以上の異種の粒子について別個の測定を行う場合もある。
【００２９】
当業者であれば、多成分粒子系には異なる多数の測定可能な特性と、それらの測定値を得る上で種々の物理機構および論理機序とが存在することが理解されるであろう。その一例として、前記多成分粒子系（２２０）について、三軸荷重試験または他の同様な機械的試験で、排水条件下の弾塑性特性を測定するよう測定機構を構成することができる。別の例としては、表面形状を取得する方法または孔内検層ツールを使って、前記多成分粒子系（２２０）の地球物理学的特性（音波速度またはかさ密度など）を測定するよう測定機構を構成することができる。
【００３０】
１若しくはそれ以上の実施形態では、測定機構がデータ解析システム（２０５）に通信可能に連結される。このデータ解析システム（２０５）は、測定された特性を処理することにより、前記多成分粒子系（２２０）の弾塑性変形に伴う巨視的特性を計算するよう構成されている。これらの特性としては、弾塑性物質パラメータ自体、またはそれに関する特性、例えば地球物理学的属性、応力属性、または石油物理学的属性などがある。巨視的特性とは、多数の粒子の総体的な挙動を特徴付けるものをいう。粒子の数が十分大きいため、堆積系は（離散粒子（不連続粒子）のセットではなく）連続的な物質として扱うことができ、連続的な系と同様に総体的な特性を伴ったものと仮定される。これらの特性は、当該系の構成要素全体の「平均」を表す。さらに、前記データ解析システム（２０５）は、当該系について、弾塑性特性と、地球物理学的属性、応力属性、または石油物理学的属性との相関をとる（または数学的関係を構築する）よう構成されている。巨視的特性の計算および弾塑性特性相関の構築については、以下で詳しく説明する。
【００３１】
１若しくはそれ以上の実施形態では、前記データ解析システム（２０５）により構築された弾塑性特性相関を使って、ボアホール周囲の堆積物の機械的完全性を維持するため、現場作業における１若しくはそれ以上の現場作業構成要素（２０２）を調整をする。１若しくはそれ以上の実施形態では、ドリルビット、マッド処理構成要素、ケーシングなど上述した構成要素が、前記現場作業構成要素（２０２）に１若しくはそれ以上含まれる。まず、弾塑性特性相関を使って、掘削を計画している地層の弾塑性特性を推定する。これらの特性は、ボアホール周囲の岩石の機械的完全性を評価するよう設計された抗井安定度モデルへの入力として使う。そのようなモデルの結果に基づき、１若しくはそれ以上の現場作業が調整される。そのような調整の例としては、ボアホール掘削に使用するマッドの密度または粘度の変更、ボアホールへのマッド循環率の変更、ボアホールの掘削率の変更、坑井軌跡の変更などがある。
【００３２】
引き続き図２を説明すると、前記データ解析システム（２０５）は、前記現場作業構成要素（２０２）に通信可能に連結できる。そのような場合、当該データ解析システム（２０５）は、電気的切り替え信号または遠隔手続き呼び出しなどの自動工程で前記現場作業構成要素（２０２）を調整するよう構成することができる。あるいは、前記データ解析システム（２０５）は、プリントアウトやコンピュータ表示など人間が読んで理解できる形態で計算結果を提示するよう構成でき、人間である作業者は、それらの結果に基づいて前記現場作業構成要素（２０２）を手動で調整することができる。
【００３３】
図３は、１若しくはそれ以上の実施形態に従って現場作業を管理する方法のフローチャートを示したものである。１若しくはそれ以上の実施形態では、図３に示す１若しくはそれ以上のブロックを省略し、反復し、または異なる順序で実施することができる。そのため、図３に示した方法の具体的な構成は、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程の範囲を限定するものと解釈すべきではない。
【００３４】
ブロック３００では、多成分系の地球物理学的特性、応力特性、石油物理学的特性（すなわちパラメータ）を取得する。そのような特性は、フレームまたは個々の粒子タイプの各特性に対応し、上記の測定機構により測定できる。特に、地球物理学的または石油物理学的なデータは、多成分粒子系における２若しくはそれ以上のタイプの粒子の特性を含むことがあり、多成分粒子系で測定される剪断（せん断）音波をさらに含むこともある。ブロック３０５では、多成分系の排水条件下の静的弾塑性特性を取得する。そのような特性は、上記の測定機構で測定できる。
【００３５】
上記のとおり、地球物理学的特性、応力特性、および石油物理学的特性の選択は、マイクロメカニカルな（マイクロメカニクスによる）原則により決定される。ブロック３１０では、それらの特性を無次元形態に変換するか、元の次元付き形態のまま残すことができる。同様に、静的弾塑性特性も、無次元形態に変換するか、元の次元付き形態のまま残すことができる。ブロック３１５では、この弾塑性特性を、１若しくはそれ以上の数学的相関を介して地球物理学的特性、応力特性、および石油物理学的特性に関係付ける。
【００３６】
１若しくはそれ以上の実施形態では、特定形態の相関が多成分粒子系の水和物飽和レベルに依存することがある。例えば、堆積物の水和物系では、水和物飽和度が低い場合および水和物飽和度が高い場合に、異なる形態を使用できる。そのため１若しくはそれ以上の実施形態では、ブロック３１５で、多成分粒子系での水和物飽和レベルに基づき、弾塑性相関の形態を選択する。異なる水和物飽和レベルに適用できる各種形態例については、以下で詳しく説明する。
【００３７】
ブロック３２０では、ブロック３１５で導出された相関の１つを使って、１若しくはそれ以上の弾塑性特性を推定する。ブロック３１５で飽和度固有の相関が選択された場合はその相関を使い、ブロック３２０で、１若しくはそれ以上の弾塑性特性を多成分系について推定することができる。ブロック３２５では、ブロック３２０で推定された１若しくはそれ以上の弾塑性特性を使って、抗井の機械的完全性をモデル化する。ブロック３３０では、そのようなモデル化の結果に基づき、現場作業が調整される。
【００３８】
以下の説明では、例えばブロック３１５で行われる弾塑性特性相関の構築例を挙げる。当業者であれば、本開示に基づき、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程を全体として逸脱しない範囲で以下の説明と異なる本発明の実施形態を想定できることが理解されるであろう。
【００３９】
１若しくはそれ以上の実施形態において、弾塑性特性相関は、粒状物質の変形を決定するマイクロメカニカルなモデルに基づいている。以下の例では、乾燥した同位体の堆積物−水和物集合体を仮定しており、クラスレート水和物は自然状態で粒子状であると見なされる。この仮定は、ＴＨＦ水和物を含んだガラスビーズについて行われた実験によると、低〜中程度の水和物飽和度について妥当といえることが示唆されている。これについては、Ｙａｎｇ，Ｊ．、Ｔｏｈｉｄｉ，Ｂ．、Ｃｌｅｎｎｅｌｌ，Ｂ．Ｍ．（２００４年）「Ｍｉｃｒｏ ａｎｄ ｍａｃｒｏ−ｓｃａｌｅ ｉｎｖｅｓｔｉｇａｔｉｏｎ ｏｆ ｃｅｍｅｎｔｉｎｇ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｇａｓ ｈｙｄｒａｔｅｓ」（ガス水和物のセメンチング特性に関するマイクロスケールおよびマクロスケールの研究）（ＡＡＰＧ Ｈｅｄｂｕｒｇ Ｒｅｓｅａｒｃｈ Ｃｏｎｆｅｒｅｎｃｅ、Ｖａｎｃｏｕｖｅｒ、２００４年９月１２〜１６日）を参照。
【００４０】
水和物の飽和度が高い場合、クラスレート水和物は連続的なマトリックス（基質）を形成する。したがって、そのような場合、後述する理論分析の態様は厳密に該当しない場合がある。ただし、水和物飽和度が高い場合であっても、本明細書で説明する理論に基づき考案される相関が有効な場合はある。
【００４１】
以下の説明において、別段の断りがない限り、用語「粒」（ｇｒａｉｎ）は鉱物粒をいい、用語「粒子」（ｐａｒｔｉｃｌｅ）は多成分粒子系の任意の成分（水和物や粒など）をいう。さらに、以下の説明では単純化した球形粒子の２成分（二元）モデルを仮定し、各タイプの粒子（水和物や粒など）は均一な直径および同一の物理特性を有する。当業者であれば、粒子が非球形またはマルチサイズで、あるいは他の物理特性が可変である、より複雑な粒子系モデルを使用できることが理解されるであろう。さらに上記のとおり、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態は、２種より多くのタイプの粒子を含む多成分粒子系へと拡張できる。
【００４２】
マイクロメカニカルモデルの態様については、Ｄａｒｖｅ，Ｆ．& Ｎｉｃｏｔ，Ｆ．（２００５年）「Ｏｎ ｉｎｃｒｅｍｅｎｔａｌ ｎｏｎ−ｌｉｎｅａｒｉｔｙ ｉｎ ｇｒａｎｕｌａｒ ｍｅｄｉａ： ｐｈｅｎｏｍｅｎｏｌｏｇｉｃａｌ ａｎｄ ｍｕｌｔｉ−ｓｃａｌｅ ｖｉｅｗｓ （Ｐａｒｔ Ｉ）」（粒状媒質における増分非線形について：現象論的およびマルチスケールの見解（パートＩ））（Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｆｏｒ Ｎｕｍｅｒｉｃａｌ ａｎｄ Ａｎａｌｙｔｉｃａｌ Ｍｅｔｈｏｄｓ ｉｎ Ｇｅｏｍｅｃｈａｎｉｃｓ、２９、１３８７−１４０９）、Ｅｍｅｒｉａｕｌｔ，Ｆ．& Ｃｌａｑｕｉｎ，Ｃ．（２００４年）「Ｓｔａｔｉｓｔｉｃａｌ ｈｏｍｏｇｅｎｉｚａｔｉｏｎ ｆｏｒ ａｓｓｅｍｂｌｉｅｓ ｏｆ ｅｌｌｉｐｔｉｃａｌ ｇｒａｉｎｓ： ｅｆｆｅｃｔ ｏｆ ｔｈｅ ａｓｐｅｃｔ ｒａｔｉｏ ａｎｄ ｐａｒｔｉｃｌｅ ｏｒｉｅｎｔａｔｉｏｎ」（楕円形粒の集合体の統計的均質化：アスペクト比および粒子配向の効果）（Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｓｏｌｉｄｓ ａｎｄ Ｓｔｒｕｃｔｕｒｅｓ、４１、５８３７−５８４９）、Ｇａｒｄｉｎｅｒ，Ｂ．Ｓ．& Ｔｏｒｄｅｓｉｌｌａｓ，Ａ．（２００４年）「Ｍｉｃｒｏｍｅｃｈａｎｉｃｓ ｏｆ ｓｈｅａｒ ｂａｎｄｓ」（剪断帯のマイクロメカニクス）（Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｓｏｌｉｄｓ ａｎｄ Ｓｔｒｕｃｔｕｒｅｓ、４１、５８８５−５９０１）、Ｓａｙｅｒｓ，Ｃ．Ｍ．（２００２年）「Ｓｔｒｅｓｓ−ｄｅｐｅｎｄｅｎｔ ｅｌａｓｔｉｃ ａｎｉｓｏｔｒｏｐｙ ｏｆ ｓａｎｄｓｔｏｎｅｓ」（応力に依存する砂岩の弾性異方性）（Ｇｅｏｐｈｙｓｉｃａｌ Ｐｒｏｓｐｅｃｔｉｎｇ，５０，８５−９５）、およびＴｏｒｄｅｓｉｌｌａｓ，Ａ．& Ｗａｌｓｈ，Ｓ．（２００２年）「Ｉｎｃｏｒｐｏｒａｔｉｎｇ ｒｏｌｌｉｎｇ ｒｅｓｉｓｔａｎｃｅ ａｎｄ ｃｏｎｔａｃｔ ａｎｉｓｏｔｒｏｐｙ ｉｎ ｍｉｃｒｏｍｅｃｈａｎｉｃａｌ ｍｏｄｅｌｓ ｏｆ ｇｒａｎｕｌａｒ ｍｅｄｉａ」（粒状媒質のマイクロメカニカルモデルへの転がり抵抗および接触異方性の導入）（Ｐｏｗｄｅｒ Ｔｅｃｈｎｏｌｏｇｙ、１２４、１０６−１１１）の論文に説明されている。
【００４３】
マイクロメカニカルモデルに関する上記文献の調査および論理的分析によれば、排水条件下の動的および静的な変形を考慮することができ、その場合以下の物理属性が必要になることが示唆されている。以下に記載したこれらの物理属性の番号は、単に読者の便宜を考慮したもので、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程の範囲を限定するものと解釈すべきではない。
【００４４】
１）粒の直径（ｄｇ）および水和物の粒径（ｄｈ）。
【００４５】
２）粒の密度（ρｇ）、水和物の密度（ρｈ）、および乾燥した集合体の密度（ρｂ）。乾燥した集合体は、クラスレート水和物および鉱物粒から成るが、間隙流体は含まない。
【００４６】
３）粒−粒接触で摩擦を決定する粒界滑り摩擦係数（μｇｇ）、水和物−水和物接触での粒界滑り摩擦係数（μｈｈ）、および粒−水和物接触での粒界滑り摩擦係数（μｇｈ）。１若しくはそれ以上の実施形態において、これらの摩擦係数は、堆積物および水和物の所与の組み合わせについて各々一定と見なされる。あるいは、より複雑な摩擦モデルも使用できる。
【００４７】
４）水和物−水和物接触における凝集力（ｃｈｈ）および粒−水和物接触における凝集力（ｃｇｈ）。１若しくはそれ以上の実施形態において、これら凝集力値は、堆積物および水和物の所与の組み合わせについて各々一定と見なされる。さらに、未固結堆積物について一般に仮定されるように、粒−粒接触における凝集力は無視できるか存在しないものとして扱うことができる。あるいは、より複雑な凝集モデルも使用できる。
【００４８】
５）所与の方向で接点が見つかる期待値を記述する接触密度分布関数。この例では、このような分布が粒−粒、粒−水和物、水和物−粒、および水和物−水和物の接触について４つ使用される。等方性物質についての三軸試験では、接触密度分布関数は、軸方向に圧縮がかかる前の、初期等方性応力状態の配向とは無関係である。軸圧縮が起こると、接触密度分布関数は一般に変化（進化）し、角度方向に依存するようになる。これについては、Ｄａｒｖｅ，Ｆ．& Ｎｉｃｏｔ，Ｆ．（２００５年）「Ｏｎ ｉｎｃｒｅｍｅｎｔａｌ ｎｏｎ−ｌｉｎｅａｒｉｔｙ ｉｎ ｇｒａｎｕｌａｒ ｍｅｄｉａ： ｐｈｅｎｏｍｅｎｏｌｏｇｉｃａｌ ａｎｄ ｍｕｌｔｉ−ｓｃａｌｅ ｖｉｅｗｓ （Ｐａｒｔ Ｉ）」（粒状媒質における増分非線形について：現象論的およびマルチスケールの見解（パートＩ））（Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｆｏｒ Ｎｕｍｅｒｉｃａｌ ａｎｄ Ａｎａｌｙｔｉｃａｌ Ｍｅｔｈｏｄｓ ｉｎ Ｇｅｏｍｅｃｈａｎｉｃｓ、２９、１３８７−１４０９）を参照。ただし上記の変化（進化）は、初期接触密度分布関数と、与えられる応力と、本明細書で定義する特性とにより完全に決定される。このため、初期等方性応力状態の接触密度分布関数を指定するだけでよい。
【００４９】
そのような状態では、平均部分配位数および種々の粒子の数量が接触密度分布関数について置換可能になる。そのため、変形過程は、単位体積あたりの粒の数（ｎｇ）と、単位体積あたりの水和物粒子の数（ｎｈ）と、平均部分配位数（Ｎｇｇ、Ｎｇｈ、Ｎｈｇ、およびＮｈｈ）とに依存すると仮定することができ、ここで、上記の下付き文字は接触の性質を表す。ｎｇＮｇｈ＝ｎｈＮｈｇであるため、１若しくはそれ以上の一実施形態では、Ｎｈｇの指定が不要になる。あるいは、別の変数を指定から省略することができる。
【００５０】
６）各種の粒子接触（粒−粒、粒−水和物、水和物−粒、および水和物−水和物）の配位数に関する分布関数。本明細書を執筆している時点で、多成分粒子系の配位数分布について構築された理論はなく、この分布を測定するための実験もほとんど行われていない。これについては、Ｐｉｎｓｏｎ，Ｄ．、Ｚｏｕ，Ｒ．Ｐ．、Ｙｕ，Ａ．Ｂ．、Ｚｕｌｌｉ，Ｐ．& ＭｃＣａｒｔｈｙ，Ｍ．Ｊ．（１９９８年）「Ｃｏｏｒｄｉｎａｔｉｏｎ ｎｕｍｂｅｒ ｏｆ ｂｉｎａｒｙ ｍｉｘｔｕｒｅｓ ｏｆ ｓｐｈｅｒｅｓ」（球体の２成分混合物（二元混合物）の配位数）（Ｊ．Ｐｈｙｓ．Ｄ、３１、４５７−４６２）などを参照。このため、各粒子接触の分布の平均値（Ｎｇｇ、Ｎｇｈ、およびＮｈｈ）だけが指定可能である。言い換えると、より高次のモーメントの役割は無視できる。ただし、多成分粒子系の配位数分布について理論が考案された場合は、その理論を使って、配位数分布を得ることができる。
【００５１】
７）粒子接触の法線方向および剪断方向の弾性剛性（弾性スティフネス）。一般に、これらの弾性スティフネスは接触ごとに異なり、与えられる応力と、接触しあう粒子の物質特性とに依存する。ただし、これらの弾性スティフネスは、粒の弾性特性がわかっており、各接触における応力を計算する上で十分な情報がある場合、既知のものと仮定できる。例えば、この仮定が可能になるのは、多成分粒子系の幾何学的構造および封圧が定義済みの場合である。１若しくはそれ以上の実施形態では、多成分粒子系の幾何学的構造を定義するには、上記の変数（Ｎｇｇ、Ｎｇｈ、Ｎｈｈ、ｎｈ、ｎｇ、ｄｇ、およびｄｈ）だけで十分と見なせる。これにより、粒子接触の法線方向および剪断方向の弾性スティフネスを指定する上で必要な追加変数は、粒および水和物のヤング率（Ｅｇ、Ｅｈ）と、粒および水和物のポアソン比（ｖｇ、ｖｈ）と、封圧（ｐｃ）とになる。
【００５２】
この議論では、考慮する物理属性のリストからクラスレート水和物の含有量を除外する。クラスレート水和物の含有量は、すでに前記リストに含まれている物理特性に作用する範囲で変形に影響を及ぼす。
【００５３】
１若しくはそれ以上の実施形態では、多成分粒子系が、粒接触によるスライドで可塑的に変形すると仮定できる。そのため、変形に対する粒子回転の寄与は無視できる。剪断帯の最終段階は、滑り接触ではなく転がり接触により決定されると期待される。これについては、Ｂａｒｄｅｔ，Ｊ．Ｐ．& Ｐｒｏｕｂｅｔ，Ｊ．（１９９１年）「Ａ ｎｕｍｅｒｉｃａｌ ｉｎｖｅｓｔｉｇａｔｉｏｎ ｏｆ ｔｈｅ ｓｔｒｕｃｔｕｒｅ ｏｆ ｐｅｒｓｉｓｔｅｎｔ ｓｈｅａｒ ｂａｎｄｓ ｉｎ ｇｒａｎｕｌａｒ ｍｅｄｉａ」（粒状媒質において持続する剪断帯構造の数値解析）（Ｇｅｏｔｅｃｈｎｉｑｕｅ、４１、５９９−６１３）およびＴｏｒｄｅｓｉｌｌａｓ，Ａ．& Ｗａｌｓｈ，Ｓ．（２００２年）「Ｉｎｃｏｒｐｏｒａｔｉｎｇ ｒｏｌｌｉｎｇ ｒｅｓｉｓｔａｎｃｅ ａｎｄ ｃｏｎｔａｃｔ ａｎｉｓｏｔｒｏｐｙ ｉｎ ｍｉｃｒｏｍｅｃｈａｎｉｃａｌ ｍｏｄｅｌｓ ｏｆ ｇｒａｎｕｌａｒ ｍｅｄｉａ」（粒状媒質のマイクロメカニカルモデルにおける転がり抵抗および接触異方性の導入）（Ｐｏｗｄｅｒ Ｔｅｃｈｎｏｌｏｇｙ、１２４、１０６−１１１）などを参照。あるいは、接点での転がり抵抗および回転剛性を記述する付加的なパラメータを指定することにより、粒子回転を解析に組み込むこともできる。
【００５４】
上記の議論に基づくと、弾塑性変形を決定する巨視的特性ηは、以下のパラメータに依存すると見なすことができる。
【数２】
【００５５】
さらに、式（１）は以下の関係に基づいて修正でき、ここで、φは多孔率、ｓｈｙｄは水和物の飽和度である。
【数３】
【００５６】
これにより、式（２）、（３）、および（４）を使って式（１）からｎｇ、ｄｈ、およびρｈを排除すると、次式が得られる。
【数４】
【００５７】
式（５）でどの変数を排除するかは、上記だけに限定されるものではない。式（５）は、次のように表現できる音速も含め、多成分粒子系の弾塑性特性を決定する変数を定義する。
【数５】
【００５８】
ここで、Ｖｐは疎密波（圧縮波）の速度、Ｖｓは剪断波の速度である。
【００５９】
この例では、どちらの音速もフレーム速度である。さらに、必要な場合は、流体が飽和した堆積物中で測定された音速をガスマンの関係または他の任意の同様な方法を使って補正できる。これについては、Ｍａｖｋｏ，Ｇ．、Ｍｕｋｅｒｊｉ，Ｔ．，& Ｄｖｏｒｋｉｎ，Ｊ．（１９９８年）「Ｔｈｅ Ｒｏｃｋ Ｐｈｙｓｉｃｓ Ｈａｎｄｂｏｏｋ」（岩石物理学ハンドブック）（Ｃａｍｂｒｉｄｇｅ Ｕｎｉｖｅｒｓｉｔｙ Ｐｒｅｓｓ、ｐ．３２９）などを参照。式（６）および（７）では、音波が誘発するひずみにより接点で生じる滑りは、ひずみが十分小さいため無視できると仮定している。したがって、ＶｐおよびＶｓは、接点における凝集特性および摩擦特性とは無関係であると見なせる。
【００６０】
引き続き上記の例を議論すると、式（６）はＮｇｇについての式に逆変換できると仮定できる。具体的には、式（６）の他のすべての変数を固定すると、ＶｐとＮｇｇとに１対１の関係を仮定することができる。この仮定は、Ｎｇｇの増加ともにＶｐが単調に増加すると期待されるため、物理的に妥当である。したがって、式（５）のＮｇｇをＶｐで置き換えて、次式のようにできる。
【数６】
【００６１】
１若しくはそれ以上の実施形態では、ＶｓをＶｐで置き換えることができる。さらに、水和物飽和度が低い場合は、機械特性が水和物の特性と無関係であると見なせる。この仮定は、核形成後、水和物が粒同士の接触領域ではなく微細孔の空間に成長する場合に適用できると期待される。したがって、式（８）は、次のように単純化できる。
【数７】
【００６２】
ただし、中程度〜高い水和物飽和度の場合は、水和物粒子に荷重がかかり、水和物粒子の幾何学的特性および物理特性が重要になる。そのような場合、Ｖｓの感度は、これらの特性に対し高まりうる。
【００６３】
さらに、式（７）はＮｈｈについての式に逆変換できると仮定できる。すなわち、式（７）の他のすべての変数を固定すると、ＶSとＮｈｈとに１対１の関係を仮定することができる。この１対１の関係は、隣接しあう水和物粒子の接点数が増加するともにＶSが単調に増加すると期待されるため、ここでも物理的に妥当である。したがって、式（８）のＮｈｈをＶＳで置き換えて、次式のようにできる。
【数８】
【００６４】
式（９）および（１０）は、弾塑性特性相関をとる際、使用できる一般的な関係である。ただし、これらの式は、一部の場合に適用するには扱いにくいことがある。その場合は、上記の式を単純化して、適用を容易にできる。以下では、そのような単純化セットの１つについて説明する。以下の単純化は、単なる例示目的で提供するものであり、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程の範囲を限定するものと解釈すべきではない。
【００６５】
一例として、式（９）および（１０）は、所与の堆積物−水和物２成分（二元）系に適用して単純化できる（あるいは、式（９）および（１０）を使用すると、複数の異なる２成分（二元）系のデータを重ね合わせたものの相関をとることができる）。上記の仮定により、堆積物−水和物の所与の組み合わせについて、ｄｇ、ｒｇ、Ｅｇ、Ｅｈ、ｖｇ、ｖｈ、μｇｇ、μｇｈ、μｈｈ、ｃｇｈ、およびｃｈｈを、固定値と仮定することができる。
【００６６】
式（１０）からＮｇｈおよびｎｈを排除できれば望ましい。本明細書を執筆している時点で、３つの平均部分配位数（Ｎｇｇ、Ｎｇｈ、およびＮｈｈ）に関する厳密な理論的表現は存在しない。ただし、多成分（多元）粒子系の平均部分配位数に関するおおよその理論的関係は、Ｄｏｄｄｓ，Ｊ．Ａ．（１９８０年）「Ｔｈｅ ｐｏｒｏｓｉｔｙ ａｎｄ ｃｏｎｔａｃｔ ｐｏｉｎｔｓ ｉｎ ｍｕｌｔｉｃｏｍｐｏｎｅｎｔ ｒａｎｄｏｍ ｓｐｈｅｒｅ ｐａｃｋｉｎｇｓ ｃａｌｃｕｌａｔｅｄ ｂｙ ａ ｓｉｍｐｌｅ ｓｔａｔｉｓｔｉｃａｌ ｇｅｏｍｅｔｒｉｃ ｍｏｄｅｌ」（多成分ランダム球体充填における多孔率および接点の、単純な統計幾何学的モデルによる計算）（Ｊ． Ｃｏｌｌｏｉｄ Ｉｎｔｅｒｆａｃｅ Ｓｃｉ ７７、３１７-３２７）、Ｏｕｃｈｉｙａｍａ，Ｎ．& Ｔａｎａｋａ，Ｔ．（１９８０年）「Ｅｓｔｉｍａｔｉｏｎ ｏｆ ｔｈｅ ａｖｅｒａｇｅ ｎｕｍｂｅｒ ｏｆ ｃｏｎｔａｃｔｓ ｂｅｔｗｅｅｎ ｒａｎｄｏｍｌｙ ｍｉｘｅｄ ｓｏｌｉｄ ｐａｒｔｉｃｌｅｓ」（無作為に混合した固体粒子間の平均接点数の推定）（Ｉｎｄ．Ｅｎｇ．Ｃｈｅｍ．Ｆｕｎｄａｍ．１９、３３８−３４０）、およびＳｕｚｕｋｉ，Ｍ．& Ｏｓｈｉｍａ，Ｔ．（１９８５年）「Ｃｏ−ｏｒｄｉｎａｔｉｏｎ ｎｕｍｂｅｒ ｏｆ ａ ｍｕｌｔｉｃｏｍｐｏｎｅｎｔ ｒａｎｄｏｍｌｙ ｐａｃｋｅｄ ｂｅｄ ｏｆ ｓｐｈｅｒｅｓ ｗｉｔｈ ｓｉｚｅ ｄｉｓｔｒｉｂｕｔｉｏｎ」（サイズ分布を伴った多成分球体を無作為に充填した層の配位数）（Ｐｏｗｄｅｒ Ｔｅｃｈｎｏｌ．４４、２１３）の論文に説明されている。
【００６７】
Ｏｕｃｈｉｙａｍａ & Ｔａｎａｋａ（１９８０年）（上記参照）により開発された式の検査法を使うと、Ｎｇｇ、Ｎｇｈ、およびＮｈｈは、次のように近似的にモデル化できることが明らかになる。
【数９】
【００６８】
さらに、式（２）および（３）を使うと、上記のようにｎｇおよびｄｈを排除することができる。したがって、式（１１）は、次のように単純化できる。
【数１０】
【００６９】
さらに、Ｏｕｃｈｉｙａｍａ & Ｔａｎａｋａ（１９８０年）（上記参照）関係の簡潔性により、式（１２）を逆変換してＮｈｈからｎｈを得ることができる。
【数１１】
【００７０】
式（１３）を式（１２）に代入するとＮｇｈについて、次式が得られる。
【数１２】
【００７１】
この例では、最後に式（１３）および（１４）を式（１０）に代入し、式（７）の逆変換によりＮｈｈをＶｓで表せることを利用して、次式が得られる。
【数１３】
【００７２】
次元解析を行うと、式（９）および（１５）における独立変数の数を減らすことができる。一般に、前記巨視的特性ηは、無次元のもの（ポアソン比やダイレーション角（ダイレイタンシー角）など）としても、圧力次元を有するもの（ヤング率、一軸圧縮強度（ｕｎｃｏｎｆｉｎｅｄ ｃｏｍｐｒｅｓｓｉｏｎ ｓｔｒｅｎｇｔｈ、略称ＵＣＳ）、凝集硬化パラメータなど）としても考慮できる。無次元化を目的とした場合、ｄｇ、ρｇ、およびＶｐをスケール変数として選択できる。ただし、他のスケール変数の組み合わせも使用できる（｛ｐｃ、ｄｇ、ρｇ｝、｛ｐｃ、ｄｇ、Ｖｐ｝など）。さらに、結果的に得られる関係は、排他的にも他の関係と組み合わせて使用することもできる。
【００７３】
上記のスケール変数を選択すると、Ｅｇ、Ｅｈ、ｃｇｈ、およびｃｈｈの項（圧力次元を有する）を
【数１４】
【００７４】
の値でスケーリングして４つの無次元変数を
【数１５】
【００７５】
の形で得ることができる。ここで、ｃは固定された堆積物−水和物組み合わせに関する定数である。ただし、４つの変数は、すべて一般性を失うことなく単一の変数
【数１６】
【００７６】
で置換することができる。ここで、ｃは圧力次元を有した任意のスケール変数である。
【００７７】
したがって、無次元化を行うと、無次元量間に以下の関係が得られる。
【数１７】
【００７８】
ここで、ηが圧力次元を有する場合は
【数１８】
【００７９】
ηが無次元の場合はη*＝η、また
【数１９】
【００８０】
である。
【００８１】
式（１６）および（１７）により、２成分（二元）堆積物の水和物系における機械特性および地球物理学的データの相関を考慮する基礎が得られる。一定の堆積物−水和物組み合わせについては、νｇ、νｈ、μｇｇ、μｇｈ、およびμｈｈも固定されるため、残りの変数についてのみ相関を求めればよい。独立変数、
【数２０】
【００８２】
は、すべて地球物理学的データから測定または推定することができる。
【００８３】
図４は、１若しくはそれ以上の実施形態に係る相関図を示したものである。具体的にいうと、図４は、ダイレーション角（ダイレイタンシー角）と無次元かさ密度ρ＊との例示的な相関（一点鎖線）の図であり、これは上記の例示的な方程式に基づく仮想データを使って構成したものである。ダイレーション角は、剪断力を受けた固体が体積を増す傾向を示す尺度である。図４において、水和物飽和度は、凡例に示したデータ点に対応している。実際には、以下の手順で図４のデータを得ることができる。水和物を含んだ堆積物の試料は、対象となる地層をコアリングすることにより、または実験室で合成代用物質を製造することにより得られる。試料の特性ρｂおよびρｇは実験室で測定され、ρ＊が計算される。次に、試料を、種々の荷重体系下で三軸圧縮試験装置にかける。そして、その試料について応力−ひずみ曲線をプロットする。ダイレーション角（ダイレイタンシー角）は、これらの応力−ひずみ曲線から推定する。これによりダイレーション角とρ＊のプロットを作成し、そのプロットの点に対し曲線フィッティングを行う。この曲線を記述する数学関数が、ダイレーション角およびρ＊の相関関数である。この例ではダイレーション角およびρ＊に良好な相関があると仮定したが、実際には一般に調査を行って、最良の相関が得られる式（１６）および（１７）の変数の組み合わせを見つける。例えば、ダイレーション角は、最終的に積ρ＊γ、または３つの変数ｓｈｙｄ、ρ＊、およびφに最もよく相関する可能性がある。当業者であれば、多種の相関が存在することが理解されるであろう。図４に示した相関は、単に例示目的で提供するものであり、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程の範囲を限定するものと解釈すべきではない。
【００８４】
ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態では、多成分粒子系（例えば、クラスレート水和物を含んだ堆積物など）の機械特性を正確に予測することができる。１若しくはそれ以上の実施形態では、そのような系における機械的な作用を説明することにより、過度に保守的な機械的破壊モデルおよび弾性−脆性モデルに伴う非効率な作業手順を避けることができる。例えば、それらの実施形態を使用すると、岩石の変形または破壊を予測でき、その予測に基づいて現場作業を調整できる。より一般には、前記実施形態により現場作業を調整して、多成分粒子系の掘削に伴う問題を軽減することができる。
【００８５】
１若しくはそれ以上の実施形態では、上記の式（１）〜（１７）を参照して説明し米国特許第７，４７２，０２２号で説明された方法を使って、次式が導出される。
【数２１】
【００８６】
ここで、式（１８）および（１９）でアスタリスクが付いている項は、すべて無次元である。具体的にいうと、
【数２２】
【００８７】
は、ガス水和物飽和度がそれぞれ低い場合および高い場合の無次元の機械特性であり、
【数２３】
【００８８】
は、それぞれ無次元の封圧、疎密波（圧縮波）速度、および剪断波速度であり、ｓｈｙｄは水和物飽和度であり、φは多孔率であり、左辺の機械特性は、圧力単位の値を封圧ｐｃでスケーリングしたものである。無次元化した上記の各項は、
【数２４】
【００８９】
で定義され、ここで、ｃは圧力の次元を伴った任意のスケール変数、ＶｐおよびＶｓは、それぞれ疎密波および剪断波の速度、そしてρｇは密度である。１若しくはそれ以上の実施形態では、これら独立変数の一部が利用できない。この場合、相関は、変数の全補集合サブセットにわたり探索を行うことにより導出できる。例えば、音速ＶｐおよびＶｓがわからない場合は、式（１８）および（１９）を使って、無次元の静的排水ヤング率（Ｅ／ｐｃ）と独立変数ｓｈｙｄ、α＊、およびφとの相関を導出することができる。一実施形態では、これらの独立変数ｓｈｙｄ、α＊、およびφを、地球物理学的データおよびスケーリング定数ｐｃから測定または推定することができる。便宜上、無次元変数の上付き位置に付けるアスタリスクは、次式には含めていない。
【００９０】
本発明の１若しくはそれ以上の一実施形態において、式（１９）は次のように展開される。
【数２５】
【００９１】
この例では、音速Ｖｐ、Ｖｓが利用できない場合に式（２０）が使用される。式（１９）の展開形態を使うと、独立変数（ｓｈｙｄ、α、φ）間の２次結合が可能になる。当業者であれば、それに代わる展開が可能であり、例えばより高次の代数項、または非代数項を使用できることが理解されるであろう。さらに、展開には、音速に関する有次元の変数や項、または式（８）、（９）、および／または（１６）の右辺に見られる項の種々の組み合わせを含めることもできる。次いで、現場から得られるデータ（または試料データ）とともに数値探索スキームを使用すると、展開式から冗長項を排除し、残りの係数について最適な値を決定することができる。この場合、冗長項は、検査値に感度のない項（相関能力の改善による検査値へのフィッティングに効果を示さない項）に対応したものであってよい。１若しくはそれ以上の実施形態では、ベイズ法に基づく逆解析（インバージョン解析）が反復されて、検査値との最適マッチをもたらす係数の値が得られる。不確定性が最も大きい（データへの感度が最も小さい）係数は展開から排除され、逆解析の別の反復が行われる。この工程は、さらに排除できる項がなくなるまで反復できる。
【００９２】
１若しくはそれ以上の実施形態では、上述のベイズ法に基づいた探索スキームを使って以下の相関が決定される。
【数２６】
【００９３】
ここで、Ａ＝９０．５８、Ｂ＝７８．９０、Ｃ＝０．５８３１、Ｄ＝８００．４、Ｅ＝１．３７１、およびＦ＝１．０２２である。この例では、上記の相関が、αおよびｓｈｙｄの２次結合を説明する。一般化した展開手法を使うと、そのような結合が発見可能になる。さらに、一般化した展開を使うと、相関の所定形態を強制することに伴うバイアスなしで、相関を生成することができる。
【００９４】
前記一般化した展開手法とその結果（式（２１）など）には、水和物を含有する堆積物の機械特性に関する比較的最近の実験測定値が利用されることが、当業者であれば理解されるであろう。さらに、ベイズ法に基づく逆解析による展開式からの項の逐次排除には、著しい計算資源が必要とされる。
【００９５】
１若しくはそれ以上の実施形態では、式（２１）の相関を適用して得られた結果により、クラスレート水和物の飽和度が高い場合のヤング率は、実際の値より低く予測される可能性があると示されている。この場合、以下の条件が満たされれば上述の相関（式（２１））を使用することができる。（ｉ）クラスレート水和物の飽和度が５０％を実質的に超えず、（ｉｉ）封圧が
【数２７】
【００９６】
の範囲内である。当業者であれば、これらの条件が満たされなくとも、前記相関を使用できることが理解されるであろう。さらに当業者であれば、砂質土に存在するガス水和物および／または数値探索スキームのデータ源である砂質土若しくは環境の他の物理特性に応じて、上記の相関における定数または当該相関の形態が異なる可能性のあることが理解されるであろう。
【００９７】
実施形態は、使用するプラットフォームにかかわらず、ほぼ全タイプのコンピュータで実施できる。例えば図５に示すように、コンピュータシステム（５００）には、プロセッサ（５０２）と、それに伴うメモリ（５０４）と、記憶装置（５０６）と、今日のコンピュータに一般に備えられた他の多くの要素および機能（図示せず）とが含まれる。前記コンピュータ（５００）には、キーボード（５０８）およびマウス（５１０）などの入力手段と、モニター装置（５１２）などの出力手段とが含まれることもある。前記コンピュータシステム（５００）は、ネットワーク（５１４）（例えば、ローカルエリアネットワーク（ＬＡＮ）、インターネットなどの広域ネットワーク（ＷＡＮ）、または他の任意の類似タイプのネットワーク）に、ネットワークインターフェース接続（図示せず）を通じて接続できる。当業者であれば、これらの入出力手段が他の形態を取ることも可能なことが理解されるであろう。
【００９８】
さらに、当業者であれば、上述のコンピュータシステム（５００）の１若しくはそれ以上の要素を遠隔位置に配置することも、またネットワーク経由で他の要素に接続することも可能なことが理解されるであろう。さらに、１若しくはそれ以上の実施形態は、複数のノードを有した分散システムで実施でき、その場合、本システムの各部分（現場作業構成要素、データ解析システム、測定機構など）を前記分散システム内の異なるノードに配置することができる。１若しくはそれ以上の実施形態において、前記ノードはコンピュータシステムに対応する。あるいは、前記ノードは、物理メモリを伴ったプロセッサに対応する。あるいは、前記ノードは、共有メモリまたは共有資源（リソース）を伴ったプロセッサに対応する。さらに、前記方法の実施形態を実行するソフトウェア命令は、コンパクトディスク（ＣＤ）、ディスケット、テープ、または他の任意のコンピュータ可読記憶装置といったコンピュータ可読媒体に格納できる。
【００９９】
以上、限定された数の実施形態を参照して、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程を説明したが、当業者であれば、本明細書に開示した、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程の範囲を逸脱しない範囲で、本開示に基づき他の実施形態も考案できることが理解されるであろう。そのため、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する当該工程の範囲は、添付の特許請求の範囲によってのみ限定されるべきである。
【技術分野】
【０００１】
本願は、２００８年６月２日付けで出願された米国仮特許出願第６１／０５８，１５５号「Ｍｅｔｈｏｄ ａｎｄ Ｓｙｓｔｅｍ ｆｏｒ Ｅｓｔｉｍａｔｉｎｇ Ｉｎ Ｓｉｔｕ Ｍｅｃｈａｎｉｃａｌ Ｐｒｏｐｅｒｔｉｅｓ ｏｆ Ｓｅｄｉｍｅｎｔｓ Ｃｏｎｔａｉｎｉｎｇ Ｇａｓ Ｈｙｄｒａｔｅｓ」（堆積物を含むガス水和物のＩｎ Ｓｉｔｕ機械特性を推定する方法およびシステム）に基づく利益を主張するものである（参照により本明細書に組み込むものとする）。
【０００２】
さらに、本願は、２００８年１２月３０日付けで交付された米国特許第７，４７２，０２２号「Ｍｅｔｈｏｄ ａｎｄ Ｓｙｓｔｅｍ ｆｏｒ Ｍａｎａｇｉｎｇ ａ Ｄｒｉｌｌｉｎｇ Ｏｐｅｒａｔｉｏｎ ｉｎ ａ Ｍｕｌｔｉｃｏｍｐｏｎｅｎｔ Ｐａｒｔｉｃｕｌａｔｅ Ｓｙｓｔｅｍ」（多成分粒子系における掘削作業を管理する方法およびシステム）に関連する（参照により本明細書に組み込むものとする）。
【背景技術】
【０００３】
図１は掘削現場作業を示した図であり、掘削リグ（１００）を使って、ボアホール（ボーリング孔）（１４５）内のドリルパイプ（１４０）遠端に連結されたドリルビット（１５０）を回転させている。この現場作業により、石油、天然ガス、水、または掘削で取得可能なまたは取得される他の任意タイプの物質を取得することができる。図１に示した現場作業は地表から地層を直接掘削していくものであるが、当業者であれば、湖底の掘削、深海底の掘削、探鉱作業、開発作業、仕上げ作業、生産作業、処理作業など、他タイプの現場作業が存在することも理解されるであろう。
【０００４】
図１に示すように、ロータリーテーブル（１２５）で生成される回転力は、前記掘削リグ（１００）から前記ドリルパイプ（１４０）を介して前記ドリルビット（１５０）へ伝達される。さらに、掘削流体（「マッド（ｍｕｄ）」とも呼ばれる）が前記ドリルパイプ（１４０）の中空コアを通じて前記ドリルビット（１５０）へ送られる。具体的には、マッドポンプ（１８０）を使い、スタンドパイプ（１６０）とホース（１５５）とケリー（１２０）とを通じて、前記ドリルパイプ（１４０）へマッドを送る。暴噴の危険性を軽減するには、防噴装置（１３０）を使って、前記ボアホール（１４５）内の流体圧力を制御する。さらに、ケーシング（１３５）を使って前記ボアホール（１４５）を強化することにより、暴噴、または当該ボアホール（１４５）に作用する他の力による倒壊を防ぐことができる。当該掘削リグ（１００）には、クラウンブロック（１０５）、トラベリングブロック（１１０）、スイベル（１１５）、および他の構成要素（図示せず）を含めることもできる。
【０００５】
前記ボアホール（１４５）から地表へ返ってくるマッドは、マッド戻り配管（１６５）を通じてマッド処理用の機器へ送られる。例えば、マッドは、掘削された固体を取り除くよう構成されたシェーカー（１７０）へ送られる。除去された固体はリザーブ（貯留）ピット（１７５）へ移動され、マッドはマッドピット（１９０）に堆積する。前記マッドポンプ（１８０）は、ろ過されたマッドを前記マッドピット（１９０）からマッド吸引管（１８５）で吸引し、前記掘削リグ（１００）内へ再注入する。当業者であれば、デガッサー、デサンダー、デシルター、遠心分離機、混合ホッパーなど他のマッド処理装置も使用できることが理解されるであろう。さらに、掘削作業には、流体技術実施、掘削シミュレーション、圧力制御、抗井からの不要物除去、および廃棄物管理などの作業に使用される他のタイプの掘削構成要素が含まれることがある。
【０００６】
所与の現場作業では（図１に示した掘削作業など）、地層の地質力学特性に関する知識を使って、現場に関連した種々の問題を軽減することができる。例えば、一部の地層では、岩石が変形し若しくは破壊されるおそれがある。塑性パラメータは、通常、コアに機械的試験を行って直接測定されるが、数学的に相関を構築する（相関をとる）と脆性−弾性特性を予測することができる。
【０００７】
例えば、静的ヤング率を推定する方法は、Ｅｉｓｓａ，Ｅ．Ａ．& Ｋａｚｉ，Ａ．らの「Ｒｅｌａｔｉｏｎ ｂｅｔｗｅｅｎ ｓｔａｔｉｃ ａｎｄ ｄｙｎａｍｉｃ Ｙｏｕｎｇ’ｓ ｍｏｄｕｌｉ ｆｏｒ ｒｏｃｋｓ」（岩石に関する静的ヤング率および動的ヤング率の関係）（Ｉｎｔ．Ｊ．Ｒｏｃｋ Ｍｅｃｈ．Ｍｉｎ．Ｓｃｉ．&Ｇｅｏｍｅｃｈ．Ａｂｓｔｒ．２５（１９８８年）、ページ４７９−４８２）、Ｍｏｎｔｍａｙｏｕｒ，Ｈ．& Ｇｒａｖｅｓ，Ｒ．Ｍ．らの「Ｐｒｅｄｉｃｔｉｏｎ ｏｆ Ｓｔａｔｉｃ Ｅｌａｓｔｉｃ／Ｍｅｃｈａｎｉｃａｌ Ｐｒｏｐｅｒｔｉｅｓ ｏｆ Ｃｏｎｓｏｌｉｄａｔｅｄ ａｎｄ Ｕｎｃｏｎｓｏｌｉｄａｔｅｄ Ｓａｎｄｓ Ｆｒｏｍ Ａｃｏｕｓｔｉｃ Ｍｅａｓｕｒｅｍｅｎｔｓ： Ｃｏｒｒｅｌａｔｉｏｎｓ」（音響測定による固結砂および未固結砂の静的な弾性特性／機械特性の予測：相関）（６１ｓｔ Ａｎｎｕａｌ Ｔｅｃｈｎｉｃａｌ Ｃｏｎｆｅｒｅｎｃｅ ａｎｄ Ｅｘｈｉｂｉｔｉｏｎ ｏｆ ｔｈｅ Ｓｏｃｉｅｔｙ ｏｆ Ｐｅｔｒｏｌｅｕｍ Ｅｎｇｉｎｅｅｒｓ、Ｎｅｗ Ｏｒｌｅａｎｓ、ＬＡ．（１９８６年）ＳＰＥ １５６４４）、Ｍｏｒａｌｅｓ，Ｒ．Ｈ．& Ｍａｒｃｉｎｅｗ，Ｒ．Ｐ．らの「Ｆｒａｃｔｕｒｉｎｇ ｏｆ ｈｉｇｈ−ｐｅｒｍｅａｂｉｌｉｔｙ ｆｏｒｍａｔｉｏｎｓ： Ｍｅｃｈａｎｉｃａｌ ｐｒｏｐｅｒｔｉｅｓ ｃｏｒｒｅｌａｔｉｏｎｓ」（浸透率の高い地層の破砕：機械特性の相関）（ＳＰＥ ２６５６１（１９９３年））、Ｙａｌｅ，Ｄ．Ｐ．& Ｊａｍｉｅｓｏｎ，Ｗ．Ｈ．らの「Ｓｔａｔｉｃ ａｎｄ Ｄｙｎａｍｉｃ Ｒｏｃｋ Ｍｅｃｈａｎｉｃａｌ Ｐｒｏｐｅｒｔｉｅｓ ｉｎ ｔｈｅ Ｈｕｇｏｔｏｎ ａｎｄ Ｐａｎｏｍａ Ｆｉｅｌｄｓ， Ｋａｎｓａｓ」（米国カンザス州Ｈｕｇｏｔｏｎ ＦｉｅｌｄおよびＰａｎｏｍａ Ｆｉｅｌｄにおける静的および動的な岩石機械特性）（ＳＰＥ Ｍｉｄ−Ｃｏｎｔｉｎｅｎｔ Ｇａｓ Ｓｙｍｐｏｓｉｕｍ，Ａｍａｒｉｌｌｏ，ＴＸ．、ＳＰＥ ２７９３９（１９９４年））、およびＴｕｔｕｎｃｕ，Ａ．Ｎ．& Ｓｈａｒｍａ，Ｍ．Ｍ．らの「Ｒｅｌａｔｉｎｇ Ｓｔａｔｉｃ ａｎｄ Ｕｌｔｒａｓｏｎｉｃ Ｌａｂｏｒａｔｏｒｙ Ｍｅａｓｕｒｅｍｅｎｔｓ ｔｏ Ａｃｏｕｓｔｉｃ Ｌｏｇ Ｍｅａｓｕｒｅｍｅｎｔｓ ｉｎ Ｔｉｇｈｔ Ｇａｓ Ｓａｎｄｓ」（タイトな（硬質）ガス砂における静的実験および超音波実験の測定値の音響検層測定値との関係付け）（６７ｔｈ Ａｎｎｕａｌ Ｔｅｃｈｎｉｃａｌ Ｃｏｎｆｅｒｅｎｃｅ ａｎｄ Ｅｘｈｉｂｉｔｉｏｎ ｏｆ ｔｈｅ Ｓｏｃｉｅｔｙ ｏｆ Ｐｅｔｒｏｌｅｕｍ Ｅｎｇｉｎｅｅｒｓ，Ｗａｓｈｉｎｇｔｏｎ， ＤＣ、ＳＰＥ ２４６８９（１９９２年））の諸論文に説明されている。
【０００８】
さらに、動的ポアソン比から静的ポアソン比を推定する方法は、Ｔｕｔｕｎｃｕ&Ｓｈａｒｍａ（１９９２）（上記で参照）およびＹａｌｅ&Ｊａｍｉｅｓｏｎ（１９９４）（上記で参照）で説明されている。ビオ数を評価するための経験的相関については、Ｋｒｉｅｆ，Ｍ．、Ｇａｒａｔ，Ｊ．、Ｓｔｅｌｌｉｎｇｗｅｒｆｆ，Ｊ．、Ｖｅｎｔｒｅ，Ｊ．らの論文「Ａ ｐｅｔｒｏｐｈｙｓｉｃａｌ ｉｎｔｅｒｐｒｅｔａｔｉｏｎ ｕｓｉｎｇ ｔｈｅ ｖｅｌｏｃｉｔｉｅｓ ｏｆ Ｐ ａｎｄ Ｓ ｗａｖｅｓ （ｆｕｌｌ−ｗａｖｅｆｏｒｍ ｓｏｎｉｃ）」（Ｐ波およびＳ波（完全波形の音波）（Ｔｈｅ Ｌｏｇ Ａｎａｌｙｓｔ ３１、１９９０年１１月、ページ３５５−３６９）の速度を使った解釈）に説明されている。
【０００９】
さらに、岩石の一軸圧縮強度を推定するための相関をとる方法が複数の著者により考案されており、Ｃｈａｎｇ，Ｃ．らの論文「Ｅｍｐｉｒｉｃａｌ Ｒｏｃｋ Ｓｔｒｅｎｇｔｈ Ｌｏｇｇｉｎｇ ｉｎ Ｂｏｒｅｈｏｌｅｓ Ｐｅｎｅｔｒａｔｉｎｇ Ｓｅｄｉｍｅｎｔａｒｙ Ｆｏｒｍａｔｉｏｎｓ」（堆積岩層を貫通したボアホールに関する経験的岩石強度の孔内検層）（ＭＵＬＬＩ−ＴＡＭＳＡ（Ｇｅｏｐｈｙｓｉｃａｌ Ｅｘｐｌｏｒａｔｉｏｎ）７、２００４年、ページ１７４−１８３）でレビューされている。この目的での付加的な相関については、Ｐｌｕｍｂ，Ｒ．Ａ．、Ｈｅｒｒｏｎ Ｓ．Ｌ．&Ｏｌｓｅｎ，Ｍ．Ｐ．「Ｉｎｆｌｕｅｎｃｅ ｏｆ Ｃｏｍｐｏｓｉｔｉｏｎ ａｎｄ Ｔｅｘｔｕｒｅ ｏｎ Ｃｏｍｐｒｅｓｓｉｖｅ Ｓｔｒｅｎｇｔｈ Ｖａｒｉａｔｉｏｎｓ ｉｎ ｔｈｅ Ｔｒａｖｉｓ Ｐｅａｋ Ｆｏｒｍａｔｉｏｎ」（組成およびテクスチャーがＴｒａｖｉｓ Ｐｅａｋ地層の圧縮強度変動に及ぼす影響）（６７ｔｈ Ａｎｎｕａｌ Ｔｅｃｈｎｉｃａｌ Ｃｏｎｆｅｒｅｎｃｅ ａｎｄ Ｅｘｈｉｂｉｔｉｏｎ ｏｆ ｔｈｅ Ｓｏｃｉｅｔｙ ｏｆ Ｐｅｔｒｏｌｅｕｍ Ｅｎｇｉｎｅｅｒｓ，Ｗａｓｈｉｎｇｔｏｎ，Ｄ．Ｃ．（１９９２年）ＳＰＥ ２４７５８）およびＱｉｕ，Ｋ．、Ｍａｒｓｄｅｎ，Ｊ．Ｒ．、Ｓｏｌｏｖｙｏｖ，Ｙ．、Ｓａｆｄａｒ，Ｍ．& Ｃｈａｒｄａｃ，Ｏ．「Ｄｏｗｎｓｃａｌｉｎｇ Ｇｅｏｍｅｃｈａｎｉｃｓ Ｄａｔａ ｆｏｒ Ｔｈｉｎ Ｂｅｄ Ｕｓｉｎｇ Ｐｅｔｒｏｐｈｙｓｉｃａｌ Ｔｅｃｈｎｉｑｕｅｓ」（石油物理学的技術を使った薄層理に関する地質力学データのダウンスケーリング）（１４ｔｈ ＳＰＥ Ｍｉｄｄｌｅ Ｅａｓｔ Ｏｉｌ ａｎｄ Ｇａｓ Ｓｈｏｗ ａｎｄ Ｃｏｎｆｅｒｅｎｃｅ， Ｂａｈｒａｉｎ（２００５年）ＳＰＥ ９３６０５）の論文で説明されている。
【発明の概要】
【課題を解決するための手段】
【００１０】
多成分粒子系の弾塑性特性相関をとる方法。当該方法は、堆積物の水和物系の地球物理学的データから、ｓｈｙｄおよびｐｃを有するパラメータを取得する工程を含み、ここで、ｓｈｙｄは水和物飽和度で、ｐｃは封圧である。当該方法は、さらに、前記パラメータおよび相関を使って前記堆積物の水和物系のヤング率（Ｅ）を決定する工程を含み、ここで、前記相関は
【数１】
【００１１】
であり、α＝ｐｃ／ｃで、ｃは圧力の次元を伴ったスケール変数である。当該方法は、さらに、前記堆積物の水和物系のヤング率に基づいて現場作業を調整する工程を含む。
【００１２】
ガス水和物を含む堆積物のｉｎ ｓｉｔｕ（測定対象を動かさず、原位置で）機械特性を推定する工程の他の態様は、以下の説明および添付の請求項から明確に理解されるであろう。
【図面の簡単な説明】
【００１３】
【図１】図１は、現場作業の図を示したものである。
【図２】図２は、１若しくはそれ以上の実施形態に係るシステムを示した図である。
【図３】図３は、１若しくはそれ以上の実施形態に従って現場作業を管理する方法のフローチャートを示したものである。
【図４】図４は、１若しくはそれ以上の実施形態に係る相関図を示したものである。
【図５】図５は、１若しくはそれ以上の実施形態に係るコンピュータシステムを示した図である。
【発明を実施するための形態】
【００１４】
以下、添付の図面を参照して、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ（測定対象を動かさず、原位置で）機械特性を推定する具体的な実施形態を詳しく説明する。種々の図にわたる同様な要素は、一貫性のため同様な参照番号で示している。
【００１５】
以下の詳細な説明では、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態について、それらの実施形態がより完全に理解されるよう具体的な詳細事項を多数記載する。ただし当業者であれば、それらの具体的な詳細事項がなくとも前記実施形態が実施可能であることは明確に理解されるであろう。他の場合、周知の特徴（機能）については、説明を不要に複雑化しないよう詳しい説明を省いている。
【００１６】
一般に、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態は、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を地震観測データまたは検層データから推定する方法およびシステムに関し、掘削作業中、仕上げ作業中、および生産作業中の機械的破壊の危険性を評価する上で極めて重要である。より具体的には、１若しくはそれ以上の実施形態は、ガス水和物含有堆積物の機械特性と地球物理学的データとの相関をとる方法に関する。
【００１７】
一般に、前記方法の実施形態では、マイクロメカニクスモデルを使った技術を適用して堆積物の物理的挙動を決定するパラメータの選択を誘導する。上述の選択に基づき、弾塑性特性（従属変数）と、検層データまたは地震観測データから推定できる特性（独立変数）との無次元化した関係のセットが導出される。前記方法では、これらの関係を使って、従属変数に関する弾塑性特性の代数展開を提案し、多種多様な情報源のデータを使って、メタンを含む砂質土のヤング率（Ｅ）と、ＴＨＦ（テトラヒドロフラン）水和物との相関を導出する。当業者であれば、上記の相関を導出する技術は、メタンを伴った砂質土およびＴＨＦ水和物だけに限定されず、粒子状物質全般に拡張できることが理解されるであろう。当業者であれば、前記技術を使用すると、機械強度、摩擦角、ダイレーション角（ダイレイタンシー角）に関する相関など、粒子状物質に関する他種の弾塑性特性相関を導出できることも理解されるであろう。
【００１８】
本発明の実施形態は、全般的に、多成分粒子系（２種若しくはそれ以上の粒子を含んだ地層）における現場作業（掘削作業、測量作業、仕上げ作業、生産作業など）を管理する方法およびシステムを提供する。ただし、本願で説明する方法は、単一成分系に適用することもできる。粒子系の地球物理学的特性、応力特性、および石油物理学的特性は、ｉｎ ｓｉｔｕ（測定対象を動かさず、原位置で）または実験室での測定値から得られる。適切な特性の選択は、マイクロメカニクスにより推定される。本明細書において、地球物理学的特性および石油物理学的特性とは、一般に岩石フレーム（構造）に関するものであるが、個々の粒子タイプの特性を含めていう場合もある。これらの特性は、適切なスケーリングにより無次元数で表せる。また、粒子系の排水条件での静的弾塑性変形を特徴付ける特性も同様に得られ、無次元形態に変換できる。弾塑性特性は、１若しくはそれ以上の経験的相関を介して地球物理学的特性、応力特性、および石油物理学的特性に関係している。前記経験的相関は、堆積系のｉｎ ｓｉｔｕ弾塑性特性を推定するため使用される。これらの特性は、岩石の変形および破壊モデルにおいて、堆積系で掘削されたボアホール周辺の岩石の機械的完全性を評価するため使用され、この評価に基づき現場作業が調整される。
【００１９】
マイクロメカニクスは複合系に関する学問であり、そのような系の機械的挙動を、構造部分の挙動により説明することを狙いとしている。複合系とは、２若しくはそれ以上の部分から成るものをいう。
【００２０】
フレーム特性は、多孔質集合体に関連するもので、その多孔質集合体内に常駐する流体（間隙流体）は除く。フレーム特性を測定する際は、試験前に間隙流体が除去され、または間隙流体の影響を最小限に抑えるような態様で試験が行われる。
【００２１】
静的特性は、一定の荷重または緩慢に変化する荷重を物質に与えた場合の変形を特徴付けるため使用される。これと対照的に、動的特性は、急速に変化する荷重を物質が受けた場合の変形を特徴付けるものである。
【００２２】
排水特性（排水条件下の特性）は、物質の変形を特徴付けるもので、その場合、当該変形は十分緩慢であるため、当該変形に影響を及ぼすことなく間隙流体が試料から流出できる。
【００２３】
弾塑性特性は、弾塑性物質の応力−ひずみ挙動をモデル化するため使用される特性である。弾塑性物質は、ひずみが小さいと弾性応答を呈し、ひずみが大きいと塑性応答を呈する。弾性変形中、当該物質に永続的な損傷は生じない。これは、当該物質への荷重が取り除かれると、荷重をかける元の状態に戻ることを意味している。ただし、荷重が十分大きくなると当該物質は塑性を呈し、永続的な損傷が生じて、荷重を取り除いても元の状態には戻らなくなる。弾塑性特性は、弾性変形および塑性変形の双方の体系において物質の応力−ひずみ応答をモデル化する構成則（構成式）に使用される。
【００２４】
ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態は、種々の多成分粒子系について実施できる。例えば、現場作業では、クラスレート（包接化合物）水和物（ハイドレート）（当該技術分野では、ガスクラスレート水和物、ガス水和物、またはクラスレートとも呼ばれる）を含んだ地層を標的にすることがある。クラスレート水和物は、水分子の格子から構成された氷状の結晶性固体である。この格子にガスまたは液体の分子が取り込まれて、この水和物構造を安定なものにしている。メタンガス水和物（メタンハイドレート）は、クラスレート水和物として自然発生する最も一般的な種である。ただし、プロパンやエタンなど他の炭化水素を含んだクラスレート水和物も存在する。さらに、テトラヒドロフラン（ｔｅｔｒａｈｙｄｒｏｆｕｒａｎ、略称ＴＨＦ）、二酸化炭素（ＣＯ２）、および硫化水素（Ｈ２Ｓ）などの非炭化水素物質が水和物構造に組み込まれることもある。クラスレート水和物は、永久凍土地域および大陸縁辺の深海環境に広く存在してる。
【００２５】
クラスレート水和物ゾーンの掘削に対する関心は増している。例えば、クラスレート水和物は、エネルギー源として使用することができる。しかし、クラスレート水和物ゾーンでの掘削および生産には種々の課題が存在する。掘削工程に起因した温度および圧力の外乱によりクラスレート水和物中のガスが解離し、抗井中への放出が制御できず、火災または暴噴が起こるおそれがある。さらに、解放されたガスが掘削用マッドをガス化させるおそれもある。場合により、解離中のクラスレート水和物を含んだ堆積層部分の崩壊で抗井が不安定になり、結果的に穴または枝掘り部が失われてしまうおそれもある。例えば、Ｃｏｌｌｅｔｔ，Ｔ．Ｓ．& Ｄａｌｌｉｍｏｒｅ，Ｓ．Ｒ．（２００２年）「Ｄｅｔａｉｌｅｄ ａｎａｌｙｓｉｓ ｏｆ ｇａｓ ｈｙｄｒａｔｅ ｉｎｄｕｃｅｄ ｄｒｉｌｌｉｎｇ ａｎｄ ｐｒｏｄｕｃｔｉｏｎ ｈａｚａｒｄｓ」（ガス水和物により誘発される掘削および生産時の危険に関する詳細解析）（４ｔｈ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｇａｓ Ｈｙｄｒａｔｅｓ（ＩＣＧＨ−４）、Ｙｏｋｏｈａｍａ、２００２年５月）を参照。その影響が壊滅的でない場合でさえ、穴の状態が良好でなければクラスレート水和物ゾーンの孔内検層結果の質が落ちるおそれがある。
【００２６】
クラスレート水和物ゾーンにおける坑井の完成およびそこでの生産にも課題が存在する。穴の状態が良好でないとセメントが効果的に結合（ボンディング）せず、ケーシング外部にガスが漏れてしまう。さらに、クラスレート水和物を除去したことにより地層の維持能力が失われると、砂質土が生じ、地層が地盤沈下を起こし、またはケーシングが破壊されるおそれがある。また、クラスレート水和物ゾーンで完成した坑井内では、高温の炭化水素により抗井ケーシングの周囲で水和物の解離が生じ、地層への圧力が高まるおそれがある。これについては、Ｆｒａｎｋｌｉｎ，Ｌ．Ｊ．（１９８１年）「Ｈｙｄｒａｔｅｓ ｉｎ Ａｒｔｉｃ Ｉｓｌａｎｄｓ」（北極諸島における水和物）（Ａ．Ｌ． Ｂｒｏｗｓｅｒ（ｅｄ．）Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ａ Ｗｏｒｋｓｈｏｐ ｏｎ Ｃｌａｔｈｒａｔｅｓ（ｇａｓ ｈｙｄｒａｔｅｓ）ｉｎ ｔｈｅ Ｎａｔｉｏｎａｌ Ｐｅｔｒｏｌｅｕｍ Ｒｅｓｅｒｖｅ ｉｎ Ａｌａｓｋａ．ＵＳＧＳ Ｏｐｅｎ−Ｆｉｌｅ Ｒｅｐｏｒｔ ８１−１２５９、１８−２１）を参照。
【００２７】
当業者であれば、上記で具体的に説明した課題が単なる例であり、クラスレート水和物ゾーンでの掘削には他の課題が多種存在することが理解されるであろう。さらに、本明細書の例ではクラスレート水和物ゾーンでの現場作業について強調したが、他種の多成分粒子系でも同様な課題が存在する。また、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態は、単一成分系にも適用でき、２より多くの成分を含む多成分粒子系に拡張することもできる。
【００２８】
図２は、１若しくはそれ以上の実施形態に係るシステム（２００）の図を示したものである。このシステム（２００）には、掘削対象である多成分粒子系（２２０）の特性を測定するよう構成された１若しくはそれ以上の測定機構（測定機構Ａ（２１０）、測定機構Ｎ（２１５）など）が含まれる。測定は、実験室試料に対し、またはｉｎ ｓｉｔｕで（測定対象を動かさず、原位置で）行うことができる。あるいは、前記多成分粒子系（２２０）のシミュレーションから、または類似組成を有する実際の地層若しくはシミュレートした地層から測定値を得ることもできる。前記多成分粒子系（２２０）は２若しくはそれ以上のタイプの粒子を含むため、１若しくはそれ以上の異種の粒子について別個の測定を行う場合もある。
【００２９】
当業者であれば、多成分粒子系には異なる多数の測定可能な特性と、それらの測定値を得る上で種々の物理機構および論理機序とが存在することが理解されるであろう。その一例として、前記多成分粒子系（２２０）について、三軸荷重試験または他の同様な機械的試験で、排水条件下の弾塑性特性を測定するよう測定機構を構成することができる。別の例としては、表面形状を取得する方法または孔内検層ツールを使って、前記多成分粒子系（２２０）の地球物理学的特性（音波速度またはかさ密度など）を測定するよう測定機構を構成することができる。
【００３０】
１若しくはそれ以上の実施形態では、測定機構がデータ解析システム（２０５）に通信可能に連結される。このデータ解析システム（２０５）は、測定された特性を処理することにより、前記多成分粒子系（２２０）の弾塑性変形に伴う巨視的特性を計算するよう構成されている。これらの特性としては、弾塑性物質パラメータ自体、またはそれに関する特性、例えば地球物理学的属性、応力属性、または石油物理学的属性などがある。巨視的特性とは、多数の粒子の総体的な挙動を特徴付けるものをいう。粒子の数が十分大きいため、堆積系は（離散粒子（不連続粒子）のセットではなく）連続的な物質として扱うことができ、連続的な系と同様に総体的な特性を伴ったものと仮定される。これらの特性は、当該系の構成要素全体の「平均」を表す。さらに、前記データ解析システム（２０５）は、当該系について、弾塑性特性と、地球物理学的属性、応力属性、または石油物理学的属性との相関をとる（または数学的関係を構築する）よう構成されている。巨視的特性の計算および弾塑性特性相関の構築については、以下で詳しく説明する。
【００３１】
１若しくはそれ以上の実施形態では、前記データ解析システム（２０５）により構築された弾塑性特性相関を使って、ボアホール周囲の堆積物の機械的完全性を維持するため、現場作業における１若しくはそれ以上の現場作業構成要素（２０２）を調整をする。１若しくはそれ以上の実施形態では、ドリルビット、マッド処理構成要素、ケーシングなど上述した構成要素が、前記現場作業構成要素（２０２）に１若しくはそれ以上含まれる。まず、弾塑性特性相関を使って、掘削を計画している地層の弾塑性特性を推定する。これらの特性は、ボアホール周囲の岩石の機械的完全性を評価するよう設計された抗井安定度モデルへの入力として使う。そのようなモデルの結果に基づき、１若しくはそれ以上の現場作業が調整される。そのような調整の例としては、ボアホール掘削に使用するマッドの密度または粘度の変更、ボアホールへのマッド循環率の変更、ボアホールの掘削率の変更、坑井軌跡の変更などがある。
【００３２】
引き続き図２を説明すると、前記データ解析システム（２０５）は、前記現場作業構成要素（２０２）に通信可能に連結できる。そのような場合、当該データ解析システム（２０５）は、電気的切り替え信号または遠隔手続き呼び出しなどの自動工程で前記現場作業構成要素（２０２）を調整するよう構成することができる。あるいは、前記データ解析システム（２０５）は、プリントアウトやコンピュータ表示など人間が読んで理解できる形態で計算結果を提示するよう構成でき、人間である作業者は、それらの結果に基づいて前記現場作業構成要素（２０２）を手動で調整することができる。
【００３３】
図３は、１若しくはそれ以上の実施形態に従って現場作業を管理する方法のフローチャートを示したものである。１若しくはそれ以上の実施形態では、図３に示す１若しくはそれ以上のブロックを省略し、反復し、または異なる順序で実施することができる。そのため、図３に示した方法の具体的な構成は、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程の範囲を限定するものと解釈すべきではない。
【００３４】
ブロック３００では、多成分系の地球物理学的特性、応力特性、石油物理学的特性（すなわちパラメータ）を取得する。そのような特性は、フレームまたは個々の粒子タイプの各特性に対応し、上記の測定機構により測定できる。特に、地球物理学的または石油物理学的なデータは、多成分粒子系における２若しくはそれ以上のタイプの粒子の特性を含むことがあり、多成分粒子系で測定される剪断（せん断）音波をさらに含むこともある。ブロック３０５では、多成分系の排水条件下の静的弾塑性特性を取得する。そのような特性は、上記の測定機構で測定できる。
【００３５】
上記のとおり、地球物理学的特性、応力特性、および石油物理学的特性の選択は、マイクロメカニカルな（マイクロメカニクスによる）原則により決定される。ブロック３１０では、それらの特性を無次元形態に変換するか、元の次元付き形態のまま残すことができる。同様に、静的弾塑性特性も、無次元形態に変換するか、元の次元付き形態のまま残すことができる。ブロック３１５では、この弾塑性特性を、１若しくはそれ以上の数学的相関を介して地球物理学的特性、応力特性、および石油物理学的特性に関係付ける。
【００３６】
１若しくはそれ以上の実施形態では、特定形態の相関が多成分粒子系の水和物飽和レベルに依存することがある。例えば、堆積物の水和物系では、水和物飽和度が低い場合および水和物飽和度が高い場合に、異なる形態を使用できる。そのため１若しくはそれ以上の実施形態では、ブロック３１５で、多成分粒子系での水和物飽和レベルに基づき、弾塑性相関の形態を選択する。異なる水和物飽和レベルに適用できる各種形態例については、以下で詳しく説明する。
【００３７】
ブロック３２０では、ブロック３１５で導出された相関の１つを使って、１若しくはそれ以上の弾塑性特性を推定する。ブロック３１５で飽和度固有の相関が選択された場合はその相関を使い、ブロック３２０で、１若しくはそれ以上の弾塑性特性を多成分系について推定することができる。ブロック３２５では、ブロック３２０で推定された１若しくはそれ以上の弾塑性特性を使って、抗井の機械的完全性をモデル化する。ブロック３３０では、そのようなモデル化の結果に基づき、現場作業が調整される。
【００３８】
以下の説明では、例えばブロック３１５で行われる弾塑性特性相関の構築例を挙げる。当業者であれば、本開示に基づき、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程を全体として逸脱しない範囲で以下の説明と異なる本発明の実施形態を想定できることが理解されるであろう。
【００３９】
１若しくはそれ以上の実施形態において、弾塑性特性相関は、粒状物質の変形を決定するマイクロメカニカルなモデルに基づいている。以下の例では、乾燥した同位体の堆積物−水和物集合体を仮定しており、クラスレート水和物は自然状態で粒子状であると見なされる。この仮定は、ＴＨＦ水和物を含んだガラスビーズについて行われた実験によると、低〜中程度の水和物飽和度について妥当といえることが示唆されている。これについては、Ｙａｎｇ，Ｊ．、Ｔｏｈｉｄｉ，Ｂ．、Ｃｌｅｎｎｅｌｌ，Ｂ．Ｍ．（２００４年）「Ｍｉｃｒｏ ａｎｄ ｍａｃｒｏ−ｓｃａｌｅ ｉｎｖｅｓｔｉｇａｔｉｏｎ ｏｆ ｃｅｍｅｎｔｉｎｇ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｇａｓ ｈｙｄｒａｔｅｓ」（ガス水和物のセメンチング特性に関するマイクロスケールおよびマクロスケールの研究）（ＡＡＰＧ Ｈｅｄｂｕｒｇ Ｒｅｓｅａｒｃｈ Ｃｏｎｆｅｒｅｎｃｅ、Ｖａｎｃｏｕｖｅｒ、２００４年９月１２〜１６日）を参照。
【００４０】
水和物の飽和度が高い場合、クラスレート水和物は連続的なマトリックス（基質）を形成する。したがって、そのような場合、後述する理論分析の態様は厳密に該当しない場合がある。ただし、水和物飽和度が高い場合であっても、本明細書で説明する理論に基づき考案される相関が有効な場合はある。
【００４１】
以下の説明において、別段の断りがない限り、用語「粒」（ｇｒａｉｎ）は鉱物粒をいい、用語「粒子」（ｐａｒｔｉｃｌｅ）は多成分粒子系の任意の成分（水和物や粒など）をいう。さらに、以下の説明では単純化した球形粒子の２成分（二元）モデルを仮定し、各タイプの粒子（水和物や粒など）は均一な直径および同一の物理特性を有する。当業者であれば、粒子が非球形またはマルチサイズで、あるいは他の物理特性が可変である、より複雑な粒子系モデルを使用できることが理解されるであろう。さらに上記のとおり、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態は、２種より多くのタイプの粒子を含む多成分粒子系へと拡張できる。
【００４２】
マイクロメカニカルモデルの態様については、Ｄａｒｖｅ，Ｆ．& Ｎｉｃｏｔ，Ｆ．（２００５年）「Ｏｎ ｉｎｃｒｅｍｅｎｔａｌ ｎｏｎ−ｌｉｎｅａｒｉｔｙ ｉｎ ｇｒａｎｕｌａｒ ｍｅｄｉａ： ｐｈｅｎｏｍｅｎｏｌｏｇｉｃａｌ ａｎｄ ｍｕｌｔｉ−ｓｃａｌｅ ｖｉｅｗｓ （Ｐａｒｔ Ｉ）」（粒状媒質における増分非線形について：現象論的およびマルチスケールの見解（パートＩ））（Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｆｏｒ Ｎｕｍｅｒｉｃａｌ ａｎｄ Ａｎａｌｙｔｉｃａｌ Ｍｅｔｈｏｄｓ ｉｎ Ｇｅｏｍｅｃｈａｎｉｃｓ、２９、１３８７−１４０９）、Ｅｍｅｒｉａｕｌｔ，Ｆ．& Ｃｌａｑｕｉｎ，Ｃ．（２００４年）「Ｓｔａｔｉｓｔｉｃａｌ ｈｏｍｏｇｅｎｉｚａｔｉｏｎ ｆｏｒ ａｓｓｅｍｂｌｉｅｓ ｏｆ ｅｌｌｉｐｔｉｃａｌ ｇｒａｉｎｓ： ｅｆｆｅｃｔ ｏｆ ｔｈｅ ａｓｐｅｃｔ ｒａｔｉｏ ａｎｄ ｐａｒｔｉｃｌｅ ｏｒｉｅｎｔａｔｉｏｎ」（楕円形粒の集合体の統計的均質化：アスペクト比および粒子配向の効果）（Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｓｏｌｉｄｓ ａｎｄ Ｓｔｒｕｃｔｕｒｅｓ、４１、５８３７−５８４９）、Ｇａｒｄｉｎｅｒ，Ｂ．Ｓ．& Ｔｏｒｄｅｓｉｌｌａｓ，Ａ．（２００４年）「Ｍｉｃｒｏｍｅｃｈａｎｉｃｓ ｏｆ ｓｈｅａｒ ｂａｎｄｓ」（剪断帯のマイクロメカニクス）（Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｓｏｌｉｄｓ ａｎｄ Ｓｔｒｕｃｔｕｒｅｓ、４１、５８８５−５９０１）、Ｓａｙｅｒｓ，Ｃ．Ｍ．（２００２年）「Ｓｔｒｅｓｓ−ｄｅｐｅｎｄｅｎｔ ｅｌａｓｔｉｃ ａｎｉｓｏｔｒｏｐｙ ｏｆ ｓａｎｄｓｔｏｎｅｓ」（応力に依存する砂岩の弾性異方性）（Ｇｅｏｐｈｙｓｉｃａｌ Ｐｒｏｓｐｅｃｔｉｎｇ，５０，８５−９５）、およびＴｏｒｄｅｓｉｌｌａｓ，Ａ．& Ｗａｌｓｈ，Ｓ．（２００２年）「Ｉｎｃｏｒｐｏｒａｔｉｎｇ ｒｏｌｌｉｎｇ ｒｅｓｉｓｔａｎｃｅ ａｎｄ ｃｏｎｔａｃｔ ａｎｉｓｏｔｒｏｐｙ ｉｎ ｍｉｃｒｏｍｅｃｈａｎｉｃａｌ ｍｏｄｅｌｓ ｏｆ ｇｒａｎｕｌａｒ ｍｅｄｉａ」（粒状媒質のマイクロメカニカルモデルへの転がり抵抗および接触異方性の導入）（Ｐｏｗｄｅｒ Ｔｅｃｈｎｏｌｏｇｙ、１２４、１０６−１１１）の論文に説明されている。
【００４３】
マイクロメカニカルモデルに関する上記文献の調査および論理的分析によれば、排水条件下の動的および静的な変形を考慮することができ、その場合以下の物理属性が必要になることが示唆されている。以下に記載したこれらの物理属性の番号は、単に読者の便宜を考慮したもので、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程の範囲を限定するものと解釈すべきではない。
【００４４】
１）粒の直径（ｄｇ）および水和物の粒径（ｄｈ）。
【００４５】
２）粒の密度（ρｇ）、水和物の密度（ρｈ）、および乾燥した集合体の密度（ρｂ）。乾燥した集合体は、クラスレート水和物および鉱物粒から成るが、間隙流体は含まない。
【００４６】
３）粒−粒接触で摩擦を決定する粒界滑り摩擦係数（μｇｇ）、水和物−水和物接触での粒界滑り摩擦係数（μｈｈ）、および粒−水和物接触での粒界滑り摩擦係数（μｇｈ）。１若しくはそれ以上の実施形態において、これらの摩擦係数は、堆積物および水和物の所与の組み合わせについて各々一定と見なされる。あるいは、より複雑な摩擦モデルも使用できる。
【００４７】
４）水和物−水和物接触における凝集力（ｃｈｈ）および粒−水和物接触における凝集力（ｃｇｈ）。１若しくはそれ以上の実施形態において、これら凝集力値は、堆積物および水和物の所与の組み合わせについて各々一定と見なされる。さらに、未固結堆積物について一般に仮定されるように、粒−粒接触における凝集力は無視できるか存在しないものとして扱うことができる。あるいは、より複雑な凝集モデルも使用できる。
【００４８】
５）所与の方向で接点が見つかる期待値を記述する接触密度分布関数。この例では、このような分布が粒−粒、粒−水和物、水和物−粒、および水和物−水和物の接触について４つ使用される。等方性物質についての三軸試験では、接触密度分布関数は、軸方向に圧縮がかかる前の、初期等方性応力状態の配向とは無関係である。軸圧縮が起こると、接触密度分布関数は一般に変化（進化）し、角度方向に依存するようになる。これについては、Ｄａｒｖｅ，Ｆ．& Ｎｉｃｏｔ，Ｆ．（２００５年）「Ｏｎ ｉｎｃｒｅｍｅｎｔａｌ ｎｏｎ−ｌｉｎｅａｒｉｔｙ ｉｎ ｇｒａｎｕｌａｒ ｍｅｄｉａ： ｐｈｅｎｏｍｅｎｏｌｏｇｉｃａｌ ａｎｄ ｍｕｌｔｉ−ｓｃａｌｅ ｖｉｅｗｓ （Ｐａｒｔ Ｉ）」（粒状媒質における増分非線形について：現象論的およびマルチスケールの見解（パートＩ））（Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｆｏｒ Ｎｕｍｅｒｉｃａｌ ａｎｄ Ａｎａｌｙｔｉｃａｌ Ｍｅｔｈｏｄｓ ｉｎ Ｇｅｏｍｅｃｈａｎｉｃｓ、２９、１３８７−１４０９）を参照。ただし上記の変化（進化）は、初期接触密度分布関数と、与えられる応力と、本明細書で定義する特性とにより完全に決定される。このため、初期等方性応力状態の接触密度分布関数を指定するだけでよい。
【００４９】
そのような状態では、平均部分配位数および種々の粒子の数量が接触密度分布関数について置換可能になる。そのため、変形過程は、単位体積あたりの粒の数（ｎｇ）と、単位体積あたりの水和物粒子の数（ｎｈ）と、平均部分配位数（Ｎｇｇ、Ｎｇｈ、Ｎｈｇ、およびＮｈｈ）とに依存すると仮定することができ、ここで、上記の下付き文字は接触の性質を表す。ｎｇＮｇｈ＝ｎｈＮｈｇであるため、１若しくはそれ以上の一実施形態では、Ｎｈｇの指定が不要になる。あるいは、別の変数を指定から省略することができる。
【００５０】
６）各種の粒子接触（粒−粒、粒−水和物、水和物−粒、および水和物−水和物）の配位数に関する分布関数。本明細書を執筆している時点で、多成分粒子系の配位数分布について構築された理論はなく、この分布を測定するための実験もほとんど行われていない。これについては、Ｐｉｎｓｏｎ，Ｄ．、Ｚｏｕ，Ｒ．Ｐ．、Ｙｕ，Ａ．Ｂ．、Ｚｕｌｌｉ，Ｐ．& ＭｃＣａｒｔｈｙ，Ｍ．Ｊ．（１９９８年）「Ｃｏｏｒｄｉｎａｔｉｏｎ ｎｕｍｂｅｒ ｏｆ ｂｉｎａｒｙ ｍｉｘｔｕｒｅｓ ｏｆ ｓｐｈｅｒｅｓ」（球体の２成分混合物（二元混合物）の配位数）（Ｊ．Ｐｈｙｓ．Ｄ、３１、４５７−４６２）などを参照。このため、各粒子接触の分布の平均値（Ｎｇｇ、Ｎｇｈ、およびＮｈｈ）だけが指定可能である。言い換えると、より高次のモーメントの役割は無視できる。ただし、多成分粒子系の配位数分布について理論が考案された場合は、その理論を使って、配位数分布を得ることができる。
【００５１】
７）粒子接触の法線方向および剪断方向の弾性剛性（弾性スティフネス）。一般に、これらの弾性スティフネスは接触ごとに異なり、与えられる応力と、接触しあう粒子の物質特性とに依存する。ただし、これらの弾性スティフネスは、粒の弾性特性がわかっており、各接触における応力を計算する上で十分な情報がある場合、既知のものと仮定できる。例えば、この仮定が可能になるのは、多成分粒子系の幾何学的構造および封圧が定義済みの場合である。１若しくはそれ以上の実施形態では、多成分粒子系の幾何学的構造を定義するには、上記の変数（Ｎｇｇ、Ｎｇｈ、Ｎｈｈ、ｎｈ、ｎｇ、ｄｇ、およびｄｈ）だけで十分と見なせる。これにより、粒子接触の法線方向および剪断方向の弾性スティフネスを指定する上で必要な追加変数は、粒および水和物のヤング率（Ｅｇ、Ｅｈ）と、粒および水和物のポアソン比（ｖｇ、ｖｈ）と、封圧（ｐｃ）とになる。
【００５２】
この議論では、考慮する物理属性のリストからクラスレート水和物の含有量を除外する。クラスレート水和物の含有量は、すでに前記リストに含まれている物理特性に作用する範囲で変形に影響を及ぼす。
【００５３】
１若しくはそれ以上の実施形態では、多成分粒子系が、粒接触によるスライドで可塑的に変形すると仮定できる。そのため、変形に対する粒子回転の寄与は無視できる。剪断帯の最終段階は、滑り接触ではなく転がり接触により決定されると期待される。これについては、Ｂａｒｄｅｔ，Ｊ．Ｐ．& Ｐｒｏｕｂｅｔ，Ｊ．（１９９１年）「Ａ ｎｕｍｅｒｉｃａｌ ｉｎｖｅｓｔｉｇａｔｉｏｎ ｏｆ ｔｈｅ ｓｔｒｕｃｔｕｒｅ ｏｆ ｐｅｒｓｉｓｔｅｎｔ ｓｈｅａｒ ｂａｎｄｓ ｉｎ ｇｒａｎｕｌａｒ ｍｅｄｉａ」（粒状媒質において持続する剪断帯構造の数値解析）（Ｇｅｏｔｅｃｈｎｉｑｕｅ、４１、５９９−６１３）およびＴｏｒｄｅｓｉｌｌａｓ，Ａ．& Ｗａｌｓｈ，Ｓ．（２００２年）「Ｉｎｃｏｒｐｏｒａｔｉｎｇ ｒｏｌｌｉｎｇ ｒｅｓｉｓｔａｎｃｅ ａｎｄ ｃｏｎｔａｃｔ ａｎｉｓｏｔｒｏｐｙ ｉｎ ｍｉｃｒｏｍｅｃｈａｎｉｃａｌ ｍｏｄｅｌｓ ｏｆ ｇｒａｎｕｌａｒ ｍｅｄｉａ」（粒状媒質のマイクロメカニカルモデルにおける転がり抵抗および接触異方性の導入）（Ｐｏｗｄｅｒ Ｔｅｃｈｎｏｌｏｇｙ、１２４、１０６−１１１）などを参照。あるいは、接点での転がり抵抗および回転剛性を記述する付加的なパラメータを指定することにより、粒子回転を解析に組み込むこともできる。
【００５４】
上記の議論に基づくと、弾塑性変形を決定する巨視的特性ηは、以下のパラメータに依存すると見なすことができる。
【数２】
【００５５】
さらに、式（１）は以下の関係に基づいて修正でき、ここで、φは多孔率、ｓｈｙｄは水和物の飽和度である。
【数３】
【００５６】
これにより、式（２）、（３）、および（４）を使って式（１）からｎｇ、ｄｈ、およびρｈを排除すると、次式が得られる。
【数４】
【００５７】
式（５）でどの変数を排除するかは、上記だけに限定されるものではない。式（５）は、次のように表現できる音速も含め、多成分粒子系の弾塑性特性を決定する変数を定義する。
【数５】
【００５８】
ここで、Ｖｐは疎密波（圧縮波）の速度、Ｖｓは剪断波の速度である。
【００５９】
この例では、どちらの音速もフレーム速度である。さらに、必要な場合は、流体が飽和した堆積物中で測定された音速をガスマンの関係または他の任意の同様な方法を使って補正できる。これについては、Ｍａｖｋｏ，Ｇ．、Ｍｕｋｅｒｊｉ，Ｔ．，& Ｄｖｏｒｋｉｎ，Ｊ．（１９９８年）「Ｔｈｅ Ｒｏｃｋ Ｐｈｙｓｉｃｓ Ｈａｎｄｂｏｏｋ」（岩石物理学ハンドブック）（Ｃａｍｂｒｉｄｇｅ Ｕｎｉｖｅｒｓｉｔｙ Ｐｒｅｓｓ、ｐ．３２９）などを参照。式（６）および（７）では、音波が誘発するひずみにより接点で生じる滑りは、ひずみが十分小さいため無視できると仮定している。したがって、ＶｐおよびＶｓは、接点における凝集特性および摩擦特性とは無関係であると見なせる。
【００６０】
引き続き上記の例を議論すると、式（６）はＮｇｇについての式に逆変換できると仮定できる。具体的には、式（６）の他のすべての変数を固定すると、ＶｐとＮｇｇとに１対１の関係を仮定することができる。この仮定は、Ｎｇｇの増加ともにＶｐが単調に増加すると期待されるため、物理的に妥当である。したがって、式（５）のＮｇｇをＶｐで置き換えて、次式のようにできる。
【数６】
【００６１】
１若しくはそれ以上の実施形態では、ＶｓをＶｐで置き換えることができる。さらに、水和物飽和度が低い場合は、機械特性が水和物の特性と無関係であると見なせる。この仮定は、核形成後、水和物が粒同士の接触領域ではなく微細孔の空間に成長する場合に適用できると期待される。したがって、式（８）は、次のように単純化できる。
【数７】
【００６２】
ただし、中程度〜高い水和物飽和度の場合は、水和物粒子に荷重がかかり、水和物粒子の幾何学的特性および物理特性が重要になる。そのような場合、Ｖｓの感度は、これらの特性に対し高まりうる。
【００６３】
さらに、式（７）はＮｈｈについての式に逆変換できると仮定できる。すなわち、式（７）の他のすべての変数を固定すると、ＶSとＮｈｈとに１対１の関係を仮定することができる。この１対１の関係は、隣接しあう水和物粒子の接点数が増加するともにＶSが単調に増加すると期待されるため、ここでも物理的に妥当である。したがって、式（８）のＮｈｈをＶＳで置き換えて、次式のようにできる。
【数８】
【００６４】
式（９）および（１０）は、弾塑性特性相関をとる際、使用できる一般的な関係である。ただし、これらの式は、一部の場合に適用するには扱いにくいことがある。その場合は、上記の式を単純化して、適用を容易にできる。以下では、そのような単純化セットの１つについて説明する。以下の単純化は、単なる例示目的で提供するものであり、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程の範囲を限定するものと解釈すべきではない。
【００６５】
一例として、式（９）および（１０）は、所与の堆積物−水和物２成分（二元）系に適用して単純化できる（あるいは、式（９）および（１０）を使用すると、複数の異なる２成分（二元）系のデータを重ね合わせたものの相関をとることができる）。上記の仮定により、堆積物−水和物の所与の組み合わせについて、ｄｇ、ｒｇ、Ｅｇ、Ｅｈ、ｖｇ、ｖｈ、μｇｇ、μｇｈ、μｈｈ、ｃｇｈ、およびｃｈｈを、固定値と仮定することができる。
【００６６】
式（１０）からＮｇｈおよびｎｈを排除できれば望ましい。本明細書を執筆している時点で、３つの平均部分配位数（Ｎｇｇ、Ｎｇｈ、およびＮｈｈ）に関する厳密な理論的表現は存在しない。ただし、多成分（多元）粒子系の平均部分配位数に関するおおよその理論的関係は、Ｄｏｄｄｓ，Ｊ．Ａ．（１９８０年）「Ｔｈｅ ｐｏｒｏｓｉｔｙ ａｎｄ ｃｏｎｔａｃｔ ｐｏｉｎｔｓ ｉｎ ｍｕｌｔｉｃｏｍｐｏｎｅｎｔ ｒａｎｄｏｍ ｓｐｈｅｒｅ ｐａｃｋｉｎｇｓ ｃａｌｃｕｌａｔｅｄ ｂｙ ａ ｓｉｍｐｌｅ ｓｔａｔｉｓｔｉｃａｌ ｇｅｏｍｅｔｒｉｃ ｍｏｄｅｌ」（多成分ランダム球体充填における多孔率および接点の、単純な統計幾何学的モデルによる計算）（Ｊ． Ｃｏｌｌｏｉｄ Ｉｎｔｅｒｆａｃｅ Ｓｃｉ ７７、３１７-３２７）、Ｏｕｃｈｉｙａｍａ，Ｎ．& Ｔａｎａｋａ，Ｔ．（１９８０年）「Ｅｓｔｉｍａｔｉｏｎ ｏｆ ｔｈｅ ａｖｅｒａｇｅ ｎｕｍｂｅｒ ｏｆ ｃｏｎｔａｃｔｓ ｂｅｔｗｅｅｎ ｒａｎｄｏｍｌｙ ｍｉｘｅｄ ｓｏｌｉｄ ｐａｒｔｉｃｌｅｓ」（無作為に混合した固体粒子間の平均接点数の推定）（Ｉｎｄ．Ｅｎｇ．Ｃｈｅｍ．Ｆｕｎｄａｍ．１９、３３８−３４０）、およびＳｕｚｕｋｉ，Ｍ．& Ｏｓｈｉｍａ，Ｔ．（１９８５年）「Ｃｏ−ｏｒｄｉｎａｔｉｏｎ ｎｕｍｂｅｒ ｏｆ ａ ｍｕｌｔｉｃｏｍｐｏｎｅｎｔ ｒａｎｄｏｍｌｙ ｐａｃｋｅｄ ｂｅｄ ｏｆ ｓｐｈｅｒｅｓ ｗｉｔｈ ｓｉｚｅ ｄｉｓｔｒｉｂｕｔｉｏｎ」（サイズ分布を伴った多成分球体を無作為に充填した層の配位数）（Ｐｏｗｄｅｒ Ｔｅｃｈｎｏｌ．４４、２１３）の論文に説明されている。
【００６７】
Ｏｕｃｈｉｙａｍａ & Ｔａｎａｋａ（１９８０年）（上記参照）により開発された式の検査法を使うと、Ｎｇｇ、Ｎｇｈ、およびＮｈｈは、次のように近似的にモデル化できることが明らかになる。
【数９】
【００６８】
さらに、式（２）および（３）を使うと、上記のようにｎｇおよびｄｈを排除することができる。したがって、式（１１）は、次のように単純化できる。
【数１０】
【００６９】
さらに、Ｏｕｃｈｉｙａｍａ & Ｔａｎａｋａ（１９８０年）（上記参照）関係の簡潔性により、式（１２）を逆変換してＮｈｈからｎｈを得ることができる。
【数１１】
【００７０】
式（１３）を式（１２）に代入するとＮｇｈについて、次式が得られる。
【数１２】
【００７１】
この例では、最後に式（１３）および（１４）を式（１０）に代入し、式（７）の逆変換によりＮｈｈをＶｓで表せることを利用して、次式が得られる。
【数１３】
【００７２】
次元解析を行うと、式（９）および（１５）における独立変数の数を減らすことができる。一般に、前記巨視的特性ηは、無次元のもの（ポアソン比やダイレーション角（ダイレイタンシー角）など）としても、圧力次元を有するもの（ヤング率、一軸圧縮強度（ｕｎｃｏｎｆｉｎｅｄ ｃｏｍｐｒｅｓｓｉｏｎ ｓｔｒｅｎｇｔｈ、略称ＵＣＳ）、凝集硬化パラメータなど）としても考慮できる。無次元化を目的とした場合、ｄｇ、ρｇ、およびＶｐをスケール変数として選択できる。ただし、他のスケール変数の組み合わせも使用できる（｛ｐｃ、ｄｇ、ρｇ｝、｛ｐｃ、ｄｇ、Ｖｐ｝など）。さらに、結果的に得られる関係は、排他的にも他の関係と組み合わせて使用することもできる。
【００７３】
上記のスケール変数を選択すると、Ｅｇ、Ｅｈ、ｃｇｈ、およびｃｈｈの項（圧力次元を有する）を
【数１４】
【００７４】
の値でスケーリングして４つの無次元変数を
【数１５】
【００７５】
の形で得ることができる。ここで、ｃは固定された堆積物−水和物組み合わせに関する定数である。ただし、４つの変数は、すべて一般性を失うことなく単一の変数
【数１６】
【００７６】
で置換することができる。ここで、ｃは圧力次元を有した任意のスケール変数である。
【００７７】
したがって、無次元化を行うと、無次元量間に以下の関係が得られる。
【数１７】
【００７８】
ここで、ηが圧力次元を有する場合は
【数１８】
【００７９】
ηが無次元の場合はη*＝η、また
【数１９】
【００８０】
である。
【００８１】
式（１６）および（１７）により、２成分（二元）堆積物の水和物系における機械特性および地球物理学的データの相関を考慮する基礎が得られる。一定の堆積物−水和物組み合わせについては、νｇ、νｈ、μｇｇ、μｇｈ、およびμｈｈも固定されるため、残りの変数についてのみ相関を求めればよい。独立変数、
【数２０】
【００８２】
は、すべて地球物理学的データから測定または推定することができる。
【００８３】
図４は、１若しくはそれ以上の実施形態に係る相関図を示したものである。具体的にいうと、図４は、ダイレーション角（ダイレイタンシー角）と無次元かさ密度ρ＊との例示的な相関（一点鎖線）の図であり、これは上記の例示的な方程式に基づく仮想データを使って構成したものである。ダイレーション角は、剪断力を受けた固体が体積を増す傾向を示す尺度である。図４において、水和物飽和度は、凡例に示したデータ点に対応している。実際には、以下の手順で図４のデータを得ることができる。水和物を含んだ堆積物の試料は、対象となる地層をコアリングすることにより、または実験室で合成代用物質を製造することにより得られる。試料の特性ρｂおよびρｇは実験室で測定され、ρ＊が計算される。次に、試料を、種々の荷重体系下で三軸圧縮試験装置にかける。そして、その試料について応力−ひずみ曲線をプロットする。ダイレーション角（ダイレイタンシー角）は、これらの応力−ひずみ曲線から推定する。これによりダイレーション角とρ＊のプロットを作成し、そのプロットの点に対し曲線フィッティングを行う。この曲線を記述する数学関数が、ダイレーション角およびρ＊の相関関数である。この例ではダイレーション角およびρ＊に良好な相関があると仮定したが、実際には一般に調査を行って、最良の相関が得られる式（１６）および（１７）の変数の組み合わせを見つける。例えば、ダイレーション角は、最終的に積ρ＊γ、または３つの変数ｓｈｙｄ、ρ＊、およびφに最もよく相関する可能性がある。当業者であれば、多種の相関が存在することが理解されるであろう。図４に示した相関は、単に例示目的で提供するものであり、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程の範囲を限定するものと解釈すべきではない。
【００８４】
ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する実施形態では、多成分粒子系（例えば、クラスレート水和物を含んだ堆積物など）の機械特性を正確に予測することができる。１若しくはそれ以上の実施形態では、そのような系における機械的な作用を説明することにより、過度に保守的な機械的破壊モデルおよび弾性−脆性モデルに伴う非効率な作業手順を避けることができる。例えば、それらの実施形態を使用すると、岩石の変形または破壊を予測でき、その予測に基づいて現場作業を調整できる。より一般には、前記実施形態により現場作業を調整して、多成分粒子系の掘削に伴う問題を軽減することができる。
【００８５】
１若しくはそれ以上の実施形態では、上記の式（１）〜（１７）を参照して説明し米国特許第７，４７２，０２２号で説明された方法を使って、次式が導出される。
【数２１】
【００８６】
ここで、式（１８）および（１９）でアスタリスクが付いている項は、すべて無次元である。具体的にいうと、
【数２２】
【００８７】
は、ガス水和物飽和度がそれぞれ低い場合および高い場合の無次元の機械特性であり、
【数２３】
【００８８】
は、それぞれ無次元の封圧、疎密波（圧縮波）速度、および剪断波速度であり、ｓｈｙｄは水和物飽和度であり、φは多孔率であり、左辺の機械特性は、圧力単位の値を封圧ｐｃでスケーリングしたものである。無次元化した上記の各項は、
【数２４】
【００８９】
で定義され、ここで、ｃは圧力の次元を伴った任意のスケール変数、ＶｐおよびＶｓは、それぞれ疎密波および剪断波の速度、そしてρｇは密度である。１若しくはそれ以上の実施形態では、これら独立変数の一部が利用できない。この場合、相関は、変数の全補集合サブセットにわたり探索を行うことにより導出できる。例えば、音速ＶｐおよびＶｓがわからない場合は、式（１８）および（１９）を使って、無次元の静的排水ヤング率（Ｅ／ｐｃ）と独立変数ｓｈｙｄ、α＊、およびφとの相関を導出することができる。一実施形態では、これらの独立変数ｓｈｙｄ、α＊、およびφを、地球物理学的データおよびスケーリング定数ｐｃから測定または推定することができる。便宜上、無次元変数の上付き位置に付けるアスタリスクは、次式には含めていない。
【００９０】
本発明の１若しくはそれ以上の一実施形態において、式（１９）は次のように展開される。
【数２５】
【００９１】
この例では、音速Ｖｐ、Ｖｓが利用できない場合に式（２０）が使用される。式（１９）の展開形態を使うと、独立変数（ｓｈｙｄ、α、φ）間の２次結合が可能になる。当業者であれば、それに代わる展開が可能であり、例えばより高次の代数項、または非代数項を使用できることが理解されるであろう。さらに、展開には、音速に関する有次元の変数や項、または式（８）、（９）、および／または（１６）の右辺に見られる項の種々の組み合わせを含めることもできる。次いで、現場から得られるデータ（または試料データ）とともに数値探索スキームを使用すると、展開式から冗長項を排除し、残りの係数について最適な値を決定することができる。この場合、冗長項は、検査値に感度のない項（相関能力の改善による検査値へのフィッティングに効果を示さない項）に対応したものであってよい。１若しくはそれ以上の実施形態では、ベイズ法に基づく逆解析（インバージョン解析）が反復されて、検査値との最適マッチをもたらす係数の値が得られる。不確定性が最も大きい（データへの感度が最も小さい）係数は展開から排除され、逆解析の別の反復が行われる。この工程は、さらに排除できる項がなくなるまで反復できる。
【００９２】
１若しくはそれ以上の実施形態では、上述のベイズ法に基づいた探索スキームを使って以下の相関が決定される。
【数２６】
【００９３】
ここで、Ａ＝９０．５８、Ｂ＝７８．９０、Ｃ＝０．５８３１、Ｄ＝８００．４、Ｅ＝１．３７１、およびＦ＝１．０２２である。この例では、上記の相関が、αおよびｓｈｙｄの２次結合を説明する。一般化した展開手法を使うと、そのような結合が発見可能になる。さらに、一般化した展開を使うと、相関の所定形態を強制することに伴うバイアスなしで、相関を生成することができる。
【００９４】
前記一般化した展開手法とその結果（式（２１）など）には、水和物を含有する堆積物の機械特性に関する比較的最近の実験測定値が利用されることが、当業者であれば理解されるであろう。さらに、ベイズ法に基づく逆解析による展開式からの項の逐次排除には、著しい計算資源が必要とされる。
【００９５】
１若しくはそれ以上の実施形態では、式（２１）の相関を適用して得られた結果により、クラスレート水和物の飽和度が高い場合のヤング率は、実際の値より低く予測される可能性があると示されている。この場合、以下の条件が満たされれば上述の相関（式（２１））を使用することができる。（ｉ）クラスレート水和物の飽和度が５０％を実質的に超えず、（ｉｉ）封圧が
【数２７】
【００９６】
の範囲内である。当業者であれば、これらの条件が満たされなくとも、前記相関を使用できることが理解されるであろう。さらに当業者であれば、砂質土に存在するガス水和物および／または数値探索スキームのデータ源である砂質土若しくは環境の他の物理特性に応じて、上記の相関における定数または当該相関の形態が異なる可能性のあることが理解されるであろう。
【００９７】
実施形態は、使用するプラットフォームにかかわらず、ほぼ全タイプのコンピュータで実施できる。例えば図５に示すように、コンピュータシステム（５００）には、プロセッサ（５０２）と、それに伴うメモリ（５０４）と、記憶装置（５０６）と、今日のコンピュータに一般に備えられた他の多くの要素および機能（図示せず）とが含まれる。前記コンピュータ（５００）には、キーボード（５０８）およびマウス（５１０）などの入力手段と、モニター装置（５１２）などの出力手段とが含まれることもある。前記コンピュータシステム（５００）は、ネットワーク（５１４）（例えば、ローカルエリアネットワーク（ＬＡＮ）、インターネットなどの広域ネットワーク（ＷＡＮ）、または他の任意の類似タイプのネットワーク）に、ネットワークインターフェース接続（図示せず）を通じて接続できる。当業者であれば、これらの入出力手段が他の形態を取ることも可能なことが理解されるであろう。
【００９８】
さらに、当業者であれば、上述のコンピュータシステム（５００）の１若しくはそれ以上の要素を遠隔位置に配置することも、またネットワーク経由で他の要素に接続することも可能なことが理解されるであろう。さらに、１若しくはそれ以上の実施形態は、複数のノードを有した分散システムで実施でき、その場合、本システムの各部分（現場作業構成要素、データ解析システム、測定機構など）を前記分散システム内の異なるノードに配置することができる。１若しくはそれ以上の実施形態において、前記ノードはコンピュータシステムに対応する。あるいは、前記ノードは、物理メモリを伴ったプロセッサに対応する。あるいは、前記ノードは、共有メモリまたは共有資源（リソース）を伴ったプロセッサに対応する。さらに、前記方法の実施形態を実行するソフトウェア命令は、コンパクトディスク（ＣＤ）、ディスケット、テープ、または他の任意のコンピュータ可読記憶装置といったコンピュータ可読媒体に格納できる。
【００９９】
以上、限定された数の実施形態を参照して、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程を説明したが、当業者であれば、本明細書に開示した、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する工程の範囲を逸脱しない範囲で、本開示に基づき他の実施形態も考案できることが理解されるであろう。そのため、ガス水和物を含む堆積物のｉｎ ｓｉｔｕ機械特性を推定する当該工程の範囲は、添付の特許請求の範囲によってのみ限定されるべきである。
【特許請求の範囲】
【請求項１】
多成分粒子系の弾塑性特性相関をとる方法であって、
堆積物の水和物系の地球物理学的データから、ｓｈｙｄおよびｐｃを有するパラメータを取得する工程であって、ｓｈｙｄは水和物飽和度で、ｐｃは封圧である、前記取得する工程と、
前記パラメータおよび相関を使って前記堆積物の水和物系のヤング率（Ｅ）を決定する工程であって、前記相関は
【数２８】
であり、α＝ｐｃ／ｃで、ｃは圧力の次元を伴ったスケール変数である、前記決定する工程と、
前記堆積物の水和物系のヤング率に基づいて現場作業を調整する工程と
を有する方法。
【請求項２】
請求項１記載の方法において、前記堆積物の水和物系におけるクラスレート（包接化合物）水和物の飽和度は、５０％を実質的に超えないものである。
【請求項３】
請求項２記載の方法において、ｐｃは、
【数２９】
の範囲内である。
【請求項４】
請求項１記載の方法において、前記堆積物の水和物系は、メタン水和物を含むものである。
【請求項５】
請求項１記載の方法において、前記堆積物の水和物系は、ＴＨＦ（テトラヒドロフラン）水和物を含むものである。
【請求項６】
請求項１記載の方法において、ｃは１ＭＰａに設定されるものである。
【請求項７】
請求項１記載の方法において、前記現場作業は、坑井の完成に関連するものである。
【請求項８】
請求項１記載の方法において、前記現場作業は、坑井での生産に関連するものである。
【請求項９】
多成分粒子系における現場作業を管理する方法であって、
前記多成分粒子系に関連付けられた複数の地球物理学的パラメータおよび石油物理学的パラメータのうち少なくとも１つを取得する工程と、
前記複数の地球物理学的パラメータおよび石油物理学的パラメータに対応する少なくとも１つの特性を計算する工程と、
前記多成分粒子系の弾塑性特性の少なくとも１つの測定値を取得する工程と、
前記多成分粒子系の弾塑性特性の前記少なくとも１つの測定値と前記少なくとも１つの特性とを使って、前記多成分粒子系の少なくとも１つの弾塑性特性を予測する工程と、
前記少なくとも１つの弾塑性特性を前記少なくとも１つの特性に関係付ける一般展開式に基づいて、前記少なくとも１つの弾塑性特性について相関をとる工程と、
数値探索スキームを使って、前記一般展開式から冗長項を排除する工程と、
前記相関を使って、対象となる地層に含まれる堆積物の少なくとも１つの弾塑性特性を推定する工程と、
前記堆積物の前記少なくとも１つの弾塑性特性を使って、前記対象となる地層の機械的完全性を評価する工程と、
前記評価に基づいて前記現場作業を調整する工程と
を有する方法。
【請求項１０】
請求項９記載の方法において、前記複数の地球物理学的パラメータおよび石油物理学的パラメータは、前記多成分粒子系の種々の粒子の体積割合と、当該多成分粒子系の多孔率と、当該多成分粒子系において測定された音速と、当該多成分粒子系の有効封圧とからなる群から選択される少なくとも１つを有するものである。
【請求項１１】
請求項９記載の方法において、前記少なくとも１つの弾塑性特性はヤング率であり、前記相関は
【数３０】
であり、Ａ、Ｂ、Ｃ、Ｄ、Ｅ、およびＦは定数である。
【請求項１２】
請求項９記載の方法において、前記少なくとも１つの弾塑性特性および前記少なくとも１つの特性は、それぞれ無次元である。
【請求項１３】
請求項１１記載の方法において、Ａ＝９０．５８、Ｂ＝７８．９０、Ｃ＝０．５８３１、Ｄ＝８００．４、Ｅ＝１．３７１、およびＦ＝１．０２２である。
【請求項１４】
請求項９記載の方法において、前記多成分粒子系は、堆積物の水和物系である。
【請求項１５】
多成分粒子系の弾塑性特性相関をとるシステムであって、
堆積物の水和物系の地球物理学的データからパラメータを取得するように構成された少なくとも１つの測定機構であって、前記パラメータはｓｈｙｄおよびｐｃを有し、ｓｈｙｄは水和物飽和度で、ｐｃは封圧である、前記少なくとも１つの測定機構と、
データ解析システムであって、
ヤング率（Ｅ）と、前記地球物理学的データの少なくとも１つの特性との相関をとるように構成され、ここで、当該相関は
【数３１】
であり、α＝ｐｃ／ｃで、ｃは圧力の次元を伴ったスケール変数であり、
前記パラメータおよび前記相関を使って前記堆積物の水和物系の前記ヤング率を決定するように構成された
前記データ解析システムと、
前記堆積物の水和物系のヤング率に基づいて現場作業を調整するように構成された少なくとも１つの現場作業構成要素と
を有するシステム。
【請求項１６】
請求項１５記載のシステムにおいて、前記データ解析システムは、前記堆積物の水和物系におけるクラスレート水和物の飽和度が５０％を実質的に超えない場合、前記ヤング率を決定するようにさらに構成されており、ｐｃは、
【数３２】
の範囲内である。
【請求項１７】
請求項１５記載のシステムにおいて、前記堆積物の水和物系は、メタン水和物を含むものである。
【請求項１８】
請求項１５記載のシステムにおいて、前記堆積物の水和物系は、ＴＨＦ水和物を含むものである。
【請求項１９】
請求項１５記載のシステムにおいて、前記少なくとも１つの現場作業構成要素は、坑井の完成に関連するものである。
【請求項２０】
請求項１５記載のシステムにおいて、前記少なくとも１つの現場作業構成要素は、坑井での生産に関連するものである。
【請求項１】
多成分粒子系の弾塑性特性相関をとる方法であって、
堆積物の水和物系の地球物理学的データから、ｓｈｙｄおよびｐｃを有するパラメータを取得する工程であって、ｓｈｙｄは水和物飽和度で、ｐｃは封圧である、前記取得する工程と、
前記パラメータおよび相関を使って前記堆積物の水和物系のヤング率（Ｅ）を決定する工程であって、前記相関は
【数２８】
であり、α＝ｐｃ／ｃで、ｃは圧力の次元を伴ったスケール変数である、前記決定する工程と、
前記堆積物の水和物系のヤング率に基づいて現場作業を調整する工程と
を有する方法。
【請求項２】
請求項１記載の方法において、前記堆積物の水和物系におけるクラスレート（包接化合物）水和物の飽和度は、５０％を実質的に超えないものである。
【請求項３】
請求項２記載の方法において、ｐｃは、
【数２９】
の範囲内である。
【請求項４】
請求項１記載の方法において、前記堆積物の水和物系は、メタン水和物を含むものである。
【請求項５】
請求項１記載の方法において、前記堆積物の水和物系は、ＴＨＦ（テトラヒドロフラン）水和物を含むものである。
【請求項６】
請求項１記載の方法において、ｃは１ＭＰａに設定されるものである。
【請求項７】
請求項１記載の方法において、前記現場作業は、坑井の完成に関連するものである。
【請求項８】
請求項１記載の方法において、前記現場作業は、坑井での生産に関連するものである。
【請求項９】
多成分粒子系における現場作業を管理する方法であって、
前記多成分粒子系に関連付けられた複数の地球物理学的パラメータおよび石油物理学的パラメータのうち少なくとも１つを取得する工程と、
前記複数の地球物理学的パラメータおよび石油物理学的パラメータに対応する少なくとも１つの特性を計算する工程と、
前記多成分粒子系の弾塑性特性の少なくとも１つの測定値を取得する工程と、
前記多成分粒子系の弾塑性特性の前記少なくとも１つの測定値と前記少なくとも１つの特性とを使って、前記多成分粒子系の少なくとも１つの弾塑性特性を予測する工程と、
前記少なくとも１つの弾塑性特性を前記少なくとも１つの特性に関係付ける一般展開式に基づいて、前記少なくとも１つの弾塑性特性について相関をとる工程と、
数値探索スキームを使って、前記一般展開式から冗長項を排除する工程と、
前記相関を使って、対象となる地層に含まれる堆積物の少なくとも１つの弾塑性特性を推定する工程と、
前記堆積物の前記少なくとも１つの弾塑性特性を使って、前記対象となる地層の機械的完全性を評価する工程と、
前記評価に基づいて前記現場作業を調整する工程と
を有する方法。
【請求項１０】
請求項９記載の方法において、前記複数の地球物理学的パラメータおよび石油物理学的パラメータは、前記多成分粒子系の種々の粒子の体積割合と、当該多成分粒子系の多孔率と、当該多成分粒子系において測定された音速と、当該多成分粒子系の有効封圧とからなる群から選択される少なくとも１つを有するものである。
【請求項１１】
請求項９記載の方法において、前記少なくとも１つの弾塑性特性はヤング率であり、前記相関は
【数３０】
であり、Ａ、Ｂ、Ｃ、Ｄ、Ｅ、およびＦは定数である。
【請求項１２】
請求項９記載の方法において、前記少なくとも１つの弾塑性特性および前記少なくとも１つの特性は、それぞれ無次元である。
【請求項１３】
請求項１１記載の方法において、Ａ＝９０．５８、Ｂ＝７８．９０、Ｃ＝０．５８３１、Ｄ＝８００．４、Ｅ＝１．３７１、およびＦ＝１．０２２である。
【請求項１４】
請求項９記載の方法において、前記多成分粒子系は、堆積物の水和物系である。
【請求項１５】
多成分粒子系の弾塑性特性相関をとるシステムであって、
堆積物の水和物系の地球物理学的データからパラメータを取得するように構成された少なくとも１つの測定機構であって、前記パラメータはｓｈｙｄおよびｐｃを有し、ｓｈｙｄは水和物飽和度で、ｐｃは封圧である、前記少なくとも１つの測定機構と、
データ解析システムであって、
ヤング率（Ｅ）と、前記地球物理学的データの少なくとも１つの特性との相関をとるように構成され、ここで、当該相関は
【数３１】
であり、α＝ｐｃ／ｃで、ｃは圧力の次元を伴ったスケール変数であり、
前記パラメータおよび前記相関を使って前記堆積物の水和物系の前記ヤング率を決定するように構成された
前記データ解析システムと、
前記堆積物の水和物系のヤング率に基づいて現場作業を調整するように構成された少なくとも１つの現場作業構成要素と
を有するシステム。
【請求項１６】
請求項１５記載のシステムにおいて、前記データ解析システムは、前記堆積物の水和物系におけるクラスレート水和物の飽和度が５０％を実質的に超えない場合、前記ヤング率を決定するようにさらに構成されており、ｐｃは、
【数３２】
の範囲内である。
【請求項１７】
請求項１５記載のシステムにおいて、前記堆積物の水和物系は、メタン水和物を含むものである。
【請求項１８】
請求項１５記載のシステムにおいて、前記堆積物の水和物系は、ＴＨＦ水和物を含むものである。
【請求項１９】
請求項１５記載のシステムにおいて、前記少なくとも１つの現場作業構成要素は、坑井の完成に関連するものである。
【請求項２０】
請求項１５記載のシステムにおいて、前記少なくとも１つの現場作業構成要素は、坑井での生産に関連するものである。
【図１】
【図２】
【図３】
【図４】
【図５】
【図２】
【図３】
【図４】
【図５】
【公開番号】特開２００９−２９３３６８（Ｐ２００９−２９３３６８Ａ）
【公開日】平成２１年１２月１７日（２００９．１２．１７）
【国際特許分類】
【外国語出願】
【出願番号】特願２００９−１３２６１５（Ｐ２００９−１３２６１５）
【出願日】平成２１年６月２日（２００９．６．２）
【出願人】（５０９０５８２２１）シュルムバーガー ホールディングス リミッテッド (3)
【Ｆターム（参考）】
【公開日】平成２１年１２月１７日（２００９．１２．１７）
【国際特許分類】
【出願番号】特願２００９−１３２６１５（Ｐ２００９−１３２６１５）
【出願日】平成２１年６月２日（２００９．６．２）
【出願人】（５０９０５８２２１）シュルムバーガー ホールディングス リミッテッド (3)
【Ｆターム（参考）】
[ Back to top ]