🔬 physics

Voronoi-entropy characterisation of geometric disorder in two-dimensional porous networks

本論文は、ボロノイ細胞面積分布のシャノンエントロピーが、二次元多孔質ネットワークにおける幾何学的無秩序さを定量化するための、空隙率に依存しないスカラー指標として機能することを示し、それによって、細孔空間の再編成が流体輸送にどのように影響するかを捉えることで、従来の空隙率や浸透率の尺度を補完するものである。

原著者： Vladivostok Suxo, Alexsandro Kirch, Caetano Rodrigues Miranda

公開日 2026-08-24

📖 1 分で読めます☕ さくっと読める

CC BY 4.0

原著者： Vladivostok Suxo, Alexsandro Kirch, Caetano Rodrigues Miranda

原論文は CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) でライセンスされています。 ✨ これは以下の論文のAI生成解説です。著者が執筆または承認したものではありません。技術的な正確性については原論文を参照してください。免責事項の全文を読む

地層の中を移動する流体――帯水層へと染み込む水、岩石の中を漂う石油、あるいは紙に染み込むインクなど――は、隠された、ねじれた迷路を通り抜けなければなりません。これらの液体がどれほど容易に通過できるかは、2つの要素に依存しています。それは、材料の中にどれだけの空隙が存在するか、そしてその空間がどのように配置されているかです。科学者たちは、この地下の世界を記述するために、長年2つの主要な道具を使用してきました。1つ目は「空隙率」であり、これは岩石のどれだけの割合が固体ではなく空隙であるかを示す単純な尺度です。2つ目は「浸透率（透過率）」であり、これは圧力がかかった際に流体がその空間をどれほど容易に流れることができるかを記述するものです。これら2つの数値はエンジニアや地質学者にとって不可欠ですが、物語の半分しか語っていません。2つの岩石サンプルが全く同じ量の空隙を持っていたとしても、孔（あな）の形状や接続性が異なるために、流体が通過する速度は劇的に変わることがあります。この「幾何学的無秩序」として知られる隠れた変異は、測定が困難ですが、地球の地下構造を理解する上で極めて重要です。

サンパウロ大学の研究チームは、実際の岩石の複雑さに迷い込むことなく、この無秩序さを定量化する方法を見つけ出そうと試みました。彼らは、多孔質媒体の簡略化された二次元モデルを作成しました。それは、テーブルの上に散らばったコインのように、円形の固体障害物が詰まった平らなシートを想像したものです。これらのコインの間の空間が、流体が流れる孔となります。その配置を理解するために、彼らは「ボロノイ図（ボロノイ・テセレーション）」と呼ばれる数学的手法を用いました。この手法は、空隙を個別の領域に分割するもので、各領域は最も近い固体障害物に属することになります。これらの領域の面積を測定することで、研究者たちは「シャノン・エントロピー」として知られる単一の数値を算出することができました。これは、配置がいかに無秩序であるかを示すスコアとして機能します。障害物が完璧に整列した格子状であればスコアはゼロになりますが、完全にランダムに散らばっている場合はより高いスコアになります。その後、彼らは強力なコンピュータ・シミュレーションを用いて、特定の流体、具体的には塩水溶液が、これらの異なる配置の中をどのように移動するかを観察しました。

研究者たちは、他の要因を一定に保ちながら、一度に一つの要因を変更することでモデルのテストを行いました。まず、障害物の位置が流れにどのように影響するかを調べました。彼らは整然とした格子の状態から始め、障害物の位置を徐々にバラバラにすることで、無秩序さを高めていきました。その結果、配置がより混沌とするにつれて、流体が通る経路は「屈曲度（トーチオシティ）」と呼ばれる性質を持つ、やや蛇行したものになることが分かりました。しかし、流体が通過する全体的な能力である浸透率は、ほとんど変化しませんでした。障害物がランダムに散らばっていても、流体は整然とした格子状の時とほぼ同等の容易さで通り抜ける道を見つけ出したのです。これは、無秩序さが一部の経路を狭める一方で、同時に他の経路を広げるため、流れのバランスが実質的に保たれるからです。

次に、チームは障害物のサイズや全空隙量を変更したときに何が起こるかを検証しました。無秩序のレベルを同じに保ったまま障害物を大きくした場合、浸透率は劇的に増加し、7倍近くに達しました。これは、障害物が大きくなると、空間を埋めるために必要な数が少なくなり、その結果、流体が駆け抜けるための隙間がより広くなったためです。同様に、障害物のサイズを一定に保ったまま空隙量を増やした場合、浸透率は70倍以上に跳ね上がりました。これらの結果は、隙間の大きさ、すなわち「スロート（喉部）」のサイズこそが、流体が移動する速さを決定づける主要な要因であることを裏付けました。無秩序の度合いは、迷路の形を変えることはあっても、流れの速度を制御するような方法ではなかったのです。

本研究は、新しいエントロピー・スコアは価値のあるツールであるが、それは流速を予測するためのものではないと結論付けています。むしろ、それは幾何学そのものの精密な記述として機能します。それは、孔の空間における局所的な再編成を捉え、幅広な通路と狭い通路の複雑な混在を単一の変数へと要約します。これにより、科学者は、空隙率と浸透率において見た目が同一であっても、内部構造が全く異なる2つの岩石を区別できるようになります。本研究で使用されたモデルは、簡略化された二次元表現であり、三次元の岩石の完全な複雑さや複数の流体の相互作用をまだ考慮に入れてはいませんが、無秩序の役割を分離することには成功しています。この研究は、地球における流体輸送を完全に理解するためには、単にどれだけの空間があるかだけでなく、その空間がどのように配置されているかという特定の統計的特徴を測定する必要があることを示唆しています。

技術要約：二次元多孔質ネットワークにおける幾何学的無秩序性のボロノイ・エントロピーによる特性評価

問題提起
多孔質媒体における流体輸送の定量化は、依然として課題である。なぜなら、標準的な記述子である空隙率（空隙率）や浸透率（流れやすさ）だけでは不十分だからである。同一の空隙率を持つサンプルであっても、孔隙空間の幾何学的な配置の違いによって、浸透率は劇的に変化し得る。空隙率は「どれだけの空間が存在するか」を測定するが、「その空間がどのように配置されているか」を定量化するものではない。既存の記述子は、空隙率や孔径から幾何学的な無秩序性を独立して分離できていないことが多い。本研究は、情報理論から導出されたスカラー量を用いて、この幾何学的無秩序性を特徴付ける手法を検討し、ボロノイ細胞面積分布のシャノン・エントロピーが、空隙率に依存しない無秩序性の尺度として機能するかどうか、およびそれが単相流特性とどのように相関するかを検証することで、この空白を埋めるものである。

手法
著者らは、多孔質媒体の固相を表す円形障害物からなる、理想化された二次元モデルを用いた。本研究では、分子動力学から得られた流体特性と細孔スケールのシミュレーションを組み合わせた階層的アプローチを採用した。

ネットワーク生成： 二次元ドメイン（横方向の長さ $5.68 \, \mu\text{m}$ ）内に、円形障害物を配置した。本研究では、位置の無秩序性、障害物のサイズ、空隙率、およびそれに伴うボロノイ・エントロピーの4つの要因を系統的に変化させた。変数を分離するため、各構成に対して1,000個の独立した実現系を生成した。これらは、規則的な正方格子から完全にランダムな配置まで多岐にわたるが、連結性を確保するために最小限の表面間距離を設けている。
ボロノイ解析： 幾何学的無秩序性は、ボロノイ細胞面積の分布を用いたシャノン・エントロピー ( $H$ ) を用いて定量化した。ボロノイ・テセレーションは障害物の中心に基づいて構築された。エントロピーは、スケーリングされた細胞面積 ( $A/\langle A \rangle$ ) の分布をビン（bin）に分割することで計算された。この指標は、細胞の辺の数に基づくトポロジー的ボロノイ・エントロピーとは異なり、細胞サイズの広がりに対して敏感である。
流体シミュレーション： 単相流のシミュレーションには、Taxilaパッケージを用いた格子ボルツマン法（LBM）を用いた。流体の特性（密度、動粘度）および流体・固体相互作用パラメータ（接触角）は、モンモリロナイト・マトリックス中のAPI塩水に関する先行の分子動力学研究からインポートされた。シミュレーションは、線形応答を保証するため、クリーピング流領域（レイノルズ数 $< 1$ ）で行われた。
輸送特性： 定常状態の速度場から、水力学的屈曲度（tortuosity）およびダルシー浸透率を抽出した。

主な貢献

幾何学的無秩序性の分離： 本研究は、位置の無秩序性、障害物のサイズ、および空隙率が輸送特性に与える影響を、自然界の岩石サンプルでは困難な分離手法によって、正常に分離することに成功した。
ボロノイ細胞面積エントロピーという記述子： 本論文は、流体シミュレーションを行わずに障害物の位置のみから計算可能な、幾何学的無秩序性のコンパクトなスカラー尺度として、ボロノイ細胞面積のシャノン・エントロピーを導入し、その妥当性を検証した。
浸透率の駆動因子の明確化： 本研究は、空隙率と浸透率が標準的な指標である一方で、それらが孔隙空間の具体的な幾何学的再編成を捉えきれていないことを示している。エントロピー指標は、一定の空隙率において、ある場所ではスロート（孔隙の首）を広げ、別の場所では狭めるという、このような「局所的な再編成」を捉えるものである。

結果

浸透率への依存性： バルク浸透率は、主に孔隙スロートのスケールによって支配されていることが判明した。単一の冪乗則 ( $k \propto w_{med}^{2.9}$ ) により、空隙率や障害物のサイズの変化にわたるデータが収束した。
位置の無秩序性の影響： 空隙率と障害物のサイズを固定した場合、位置の無秩序性が増加するにつれ（規則的な格子からランダムへ）、浸透率はわずかに増加し（約 $4.75 $から$ 5.01 $mD）、屈曲度は増加した（$ 1.04 $から$ 1.16$）。浸透率の増加は、スロート幅の分布が広がり、優先的な広いチャネルが出現することで、他のチャネルの狭まりを相殺したことに起因する。
エントロピーの影響： 完全ランダムなネットワーク内において、エントロピーのみを変化させた場合（空隙率と障害物サイズを固定）、浸透率や屈曲率に系統的な傾向は見られず、変動は構成ごとの統計的なばらつきの範囲内に留まった。これは、結晶的な秩序が失われると、さらなるエントロピーの上昇が輸送特性を大きく変えることはないことを示している。
障害物のサイズと空隙率： （空隙率を固定した状態で）障害物のサイズを大きくすると、浸透率は大幅に増加し（約7倍）、屈曲度も増加した。空隙率を増加させると（障害物サイズを固定）、浸透率は急上昇し（70倍以上）、屈曲度は減少した。
エントロピーと屈曲度： エントロピーは、規則的な格子から無秩序な状態への遷移過程において、屈曲度と相関を示した。しかし、この相関は、規則的な格子が流れの方向に整列しており、経路長を最小化しているという側面があるため、部分的にアーティファクト（人工的な現象）である。完全に無秩序な状態では、この相関は弱まる。

意義と主張
本論文は、ボロノイ細胞面積のシャノン・エントロピーが、直接的な予測因子ではなく、補完的な記述子として機能することを主張している。

補完的な役割： エントロピーは、空隙率や浸透率では測定できない幾何学的な無秩序を定量化する。それは、孔隙空間の再編成の複雑さを単一の変数へと凝縮するものである。
予測の限界： 本研究は、エントロピーが孔隙スロートのスケールよりも浸透率を良く予測することはないと明示している。浸透率のエントロピーへの依存性の弱さは、臨界パス輸送理論と一致しており、これは流れがグローバルな配置統計ではなく、ボトルネックのスケール（最も広い連結スロート）によって制御されることを示唆している。
冗長性： 著者らは、観察されたガンマ分布に近い分布において、エントロピーは細胞面積の変動係数（CV）と数学的に冗長であることを指摘している。したがって、それは細胞面積の広がりの便利な単一数値の要約として機能する。
適用範囲： 結論は、円形障害物を通過する二次元、単相、クリーピング流という理想化された条件に厳密に限定されている。著者らは、これらの結果が追加の検証なしに、複雑な三次元、多相の貯留岩に直接適用できるとは主張していない。ただし、デジタルロックのワークフローにおいて、完全な孔隙ネットワークの抽出なしに、粒子の重心からエントロピーを計算することは可能である。

要約すると、本研究は、幾何学的無秩序（ボロノイ・エントロピーで測定）が孔隙空間の局所的なトポロジーを変化させる一方で、この領域におけるマクロな浸透率は孔隙スロートのスケールによって支配されていることを確立した。エントロピー指標は、媒体の「配置」を特徴付ける上では価値があるが、主要な輸送記述子として空隙率や浸透率に取って代わるものではない。

関連論文