ダルシーの法則は、多孔質媒体およびヘレ・ショーセルを通る流体の流れを記述する方程式です。この法則は、ヘンリー・ダルシーが砂層を通る水の流れに関する実験結果[ 1 ]に基づいて定式化され、地球科学の一分野である水文地質学の基礎を形成しました。これは、静電気学におけるオームの法則に類似しており、流体の体積流量を水頭差(多くの場合、圧力差に比例する)と水理伝導率を介して線形的に関連付けます。実際、ダルシーの法則は、運動量フラックスに関するストークス方程式の特殊な場合であり、運動量ナビエ・ストークス方程式から導出されます。
ダルシーの法則は、熱伝導の分野におけるフーリエの法則、電気回路の分野におけるオームの法則、拡散理論におけるフィックの法則に類似している。

1830年代以降、フランスの水理技師アンリ・ダルシーは、フランスのディジョン市の水道システムを改善するために、砂ろ過器内の水の流量を研究した。彼は経験的に法則を確立し、1856年に建設報告書で発表した。[ 2 ]

積分形式では、重力が存在せず、均質な透過性媒体においては、ダルシーの法則は体積流量と体積流量の間の単純な比例関係によって与えられる。圧力低下多孔質媒体 を通して。比例定数は透過率と関連している。媒体の流体の動粘度与えられた距離圧力降下が計算される領域と断面積形式は以下のとおりです。
比率:
ダルシーの法則による水力抵抗とも定義できます。
ダルシーの法則は、局所的な形式に一般化することができる。
どこ水理勾配とは体積流量であり、ここでは表面速度とも呼ばれます。比率に注意してください。
これは、ダルシーの法則の水理伝導率と考えることができます。
(より一般的な形式ではない)積分形式では、体積流量と圧力勾配は次の比率に対応する。
。
異方性多孔質媒体の場合、透過率は2階テンソルであり、テンソル表記ではより一般的な法則を記述することができる。
量に注目してください(しばしばダルシー流束またはダルシー速度と呼ばれる)は、流体が細孔を通過する速度ではありません。これは体積流束、つまり単位面積あたりの流量です(「Flux Flummoxed」、 Ground Water、44(2):125–128、2006を参照)。流速(u )は、多孔度(φ)によって、次の式で流束(q)と関連付けられます。
単相(流体)流に対するダルシーの構成方程式は、絶対透過率(単相透過率)を定義する方程式である。
右図を参照すると、流速はSI単位で表されている。多孔度φは無次元数であるため、ダルシー流束は、または単位面積当たりの排出量も、単位で定義されます。透過性単位動粘度単位水理勾配の単位は。
積分形式では、全圧力降下単位は、 そしてはサンプルの長さ(単位)ですダルシーの体積流量、または排出量は、単位でも定義されます。断面積単位これらのパラメータのいくつかは、以下の代替定義で使用されます。 流体は高圧領域から低圧領域へ流れるという標準的な物理学の慣例に従い、フラックスの定義では負の符号が使用されます。入口と出口が異なる高さにある場合は、標高ヘッドを考慮する必要があることに注意してください。 圧力の変化が負の場合、流れは正のx方向になります。 絶対透過率の構成方程式についてはいくつかの提案があり、最も有名なのはおそらくKozeny 方程式( Kozeny–Carman 方程式とも呼ばれる) です。
静水圧の関係(ステビンの法則)を考慮すると、次のようになります。
積分形式を次の式に展開することもできます。 ここでνは動粘度である。したがって、対応する透水係数は次のようになる。
ダルシーの法則は、帯水層を流れる地下水が示すいくつかのよく知られた性質を簡潔にまとめた単純な数学的表現であり、以下のような性質が含まれます。
定常状態の地下水流方程式(ダルシーの法則と質量保存則に基づく)の使用例として、ダムの下を流れる地下水の量を定量化するための流線網の構築が挙げられる。
ダルシーの法則は、低速で粘性のある流れにのみ有効ですが、ほとんどの地下水流はこのカテゴリーに該当します。通常、レイノルズ数が1未満の流れは明らかに層流であり、ダルシーの法則を適用しても問題ありません。実験では、地下水流の場合のように、レイノルズ数が10までの流れでもダルシーの法則が成り立つ場合があることが示されています。多孔質媒体の流れのレイノルズ数(無次元パラメータ)は、通常次のように表されます。
ここで、νは水の動粘度、qは比流量(細孔速度ではなく、単位は長さ/時間)、 dは多孔質媒体の代表的な粒径(標準的な選択はmath| d30で、これはふるいを使用した粒度分析からの30 %通過サイズであり、単位は長さ)です。
定常、クリープ、非圧縮性の流れ、すなわち D ( ρu i ) / Dt ≈ 0の場合、ナビエ・ストークス方程式はストークス方程式に簡略化され、バルク項を無視すると次のようになります。
ここで、μは粘度、u iはi方向の速度、pは圧力である。粘性抵抗力が速度に比例すると仮定すると、次のように書ける。
ここで、φは多孔度、k ijは2次透過率テンソルである。これにより、 n方向の速度が得られる。
これはn方向の体積流量密度に関するダルシーの法則を与える。
等方性多孔質媒体では、透過率テンソルの非対角要素はゼロ、k ij = 0 ( i ≠ j)であり、対角要素は同一、k ii = kとなり、以下の共通形式が得られます。これにより、与えられた領域で一連の方程式を解くことにより、液体の流れ速度を決定できます。[ 3 ]
上記の式は、多孔質媒体における単相流体の流れを支配する方程式である。
ダルシーの法則の別の導出は、石油工学において透過性媒体を通る流れを決定するために広く用いられており、その中で最も単純なものは、単一の流体相と一定の流体粘度を持つ一次元均質岩石層の場合である。
ほぼすべての油田には油層の下に水層があり、油層の上にガス層がある油田もあります。石油生産によって油田の圧力が低下すると、下から水が油層に流れ込み、上からガスが油層に流れ込みます(ガス層が存在する場合)。その結果、油層内のすべての流体相が同時に流れ、混ざり合わない混合状態になります。油田の操業者は、石油生産量を増やすために水(またはガス)を注入することもあります。そのため、石油業界では、Morris Muskat らが開発した多相流の一般化された Darcy 方程式を使用しています。 Darcy の名前は広く知られており、多孔質媒体の流れと強く結びついているため、多相流の方程式は、多相流の Darcy の法則、一般化された Darcy 方程式(または法則)、あるいは単に Darcy 方程式(または法則)、または文脈上、本文がMuskatの多相流方程式について論じているとわかる場合は流れ方程式と呼ばれます。石油・ガス貯留層における多相流は包括的なテーマであり、このテーマに関する多くの記事の1つに、多相流に関するダルシーの法則があります。
ダルシーの法則は、モカポットでのコーヒー抽出の物理現象、特に加圧下で熱湯がコーヒー粉を通過する様子をモデル化するために、数多くの論文で利用されてきました。2001年のA. VarlamovとG. Balestrinoの論文[ 4 ]を皮切りに、2007年のConcetto Gianinoの論文[ 5 ] 、 2008年のL. Navariniらの論文[ 6 ] 、そして2008年のW. Kingの論文[ 7 ]へと続いています。これらの論文では、簡略化のためにコーヒーの透過率を一定とするか、抽出過程を通して変化を測定するかのいずれかになります。
ダルシーの法則は、非常に一般的には次のように表すことができます。
ここで、 qはその媒体内の特定点における流体の体積流量ベクトル、hはその点における全水頭、Kはその点における水理伝導率テンソルである。水理伝導率はしばしばスカラーで近似できる。(静電気学におけるオームの法則との類似性に注意。流量ベクトルは電流密度に、水頭は電圧に、水理伝導率は電気伝導率にそれぞれ対応する。)
レイノルズ数が約 1 ~ 10 より大きい多孔質媒体内の流れでは、慣性効果も重要になることがあります。通常、ダルシーの式に慣性項が追加され、これはフォルヒハイマー項として知られています。この項は、圧力差と流量データの非線形挙動を説明することができます。 [ 8 ] [ 9 ]
ここで、追加項k 1は慣性透磁率として知られており、長さの単位で表されます。。
レイノルズ数が1 より大きい場合、ダルシーの法則から遷移した後、慣性流領域の開始が観察されます。ダルシーの法則はもはや有効ではなく、非線形慣性補正が加えられます。この挙動の起源は、細孔スケールでの流れの流線の形状の変化の結果です。この変化は、流れの構造の変化に重要な役割を果たす再循環領域によって誘発されます。[ 9 ]
砂岩貯留層の中央部における流れは非常に遅いため、通常はフォルヒハイマーの式は不要ですが、ガス生産井へのガス流入量は、フォルヒハイマーの式を使用するのに十分なほど高い場合があります。この場合、3Dモデルのグリッドセルではなく、井戸全体の流入性能計算はフォルヒハイマーの式に基づいて行われます。この結果、流入性能の式に、流量に依存するスキン効果が追加されます。
炭酸塩岩貯留層には多くの亀裂が存在するため、多相流に関するダルシーの式は、亀裂内の流れとマトリックス(すなわち、従来の多孔質岩)内の流れの両方を記述できるように一般化される。亀裂壁の不規則な表面と亀裂内の高い流量は、フォルヒハイマーの式の使用を正当化する可能性がある。
特性寸法が小さいガス流(例えば、非常に細かい砂、ナノ多孔質構造など)では、粒子と壁の相互作用がより頻繁になり、追加の壁摩擦(クヌーセン摩擦)が発生します。粘性摩擦とクヌーセン摩擦の両方が存在するこの領域の流れには、新しい定式化を使用する必要があります。Martin Knudsen は、小さな毛細管での実験に基づいて、遷移領域の流れの半経験的モデルを発表しました。[ 10 ] [ 11 ]多孔質媒体の場合、クヌーセン方程式は次のように与えられます。 [ 11 ]
ここで、Nはモル流束、Rは気体定数、Tは温度、D eff Kは多孔質媒体の有効クヌッセン拡散率です。このモデルは、第一原理に基づく二成分摩擦モデル (BFM) からも導出できます。[ 12 ] [ 13 ] BFM に基づく多孔質媒体の遷移流の微分方程式は次のように与えられます。[ 12 ]
この式は毛細管と多孔質媒体の両方に有効です。クヌーセン効果とクヌーセン拡散率という用語は、機械工学と化学工学でより一般的です。地質工学と石油化学工学では、この効果はクリンケンベルグ効果として知られています。モル流束の定義を使用すると、上記の式は次のように書き換えることができます。
この方程式は次のように変形できます。
この方程式を従来のダルシーの法則と比較すると、新しい定式化は次のようになる。
どこ
これは、クリンケンベルクによって提案された有効透過率の定式化と同等である。[ 14 ]
ここで、bはクリンケンベルグパラメータとして知られており、気体と多孔質媒体の構造に依存します。これは、上記の式を比較すれば明らかです。クリンケンベルグパラメータbは、透過率、クヌッセン拡散率、粘度(つまり、気体と多孔質媒体の両方の特性)に依存します。
非常に短い時間スケールでは、ダルシーの法則にフラックスの時間微分を加えることができ、その結果、非常に短い時間で有効な解が得られます(熱伝達では、これはフーリエの法則の修正形と呼ばれます)。
ここでτは非常に小さな時定数であり、この時定数によってこの方程式は「通常の」時間(>ナノ秒)でダルシーの法則の通常の形式に簡略化されます。これを行う主な理由は、通常の地下水流方程式(拡散方程式)が非常に短い時間で定水頭境界で特異点を引き起こすためです。この形式は数学的に厳密ですが、双曲線型の地下水流方程式につながり、解くのがより難しく、非常に短い時間でのみ有用であり、通常は実用的な範囲外です。
ダルシーの法則の伝統的な形式に対するもう1つの拡張は、境界間の遷移流を説明するために使用されるブリンクマン項です(1949年にアンリ・ブリンクマンによって導入されました[ 15 ])。
ここでβは有効粘度項である。この補正項は、媒体の粒子自体が多孔質である媒体を通過する流れを考慮するものであるが、使用が難しく、通常は無視される。
1936年、同僚のラルフ・デューイ・ワイコフとホルブルック・ゴーラム・ボセットの実験結果に基づき、モリス・マスカットとミラン・W・メレスは、石油貯留層の多孔質媒体における水、油、ガスの多相流をカバーするようにダルシーの法則を一般化した。[ 16 ]マスカットらによる多相流のダルシーの法則は、今日まで続く貯留層工学の分析的基礎を提供している。
ダルシーの法則は、堆積物中の層流に対して有効です。細粒堆積物では、間隙の寸法が小さいため、流れは層流になります。粗粒堆積物も同様の挙動を示しますが、非常に粗粒の堆積物では、流れが乱流になる場合があります。[ 17 ]したがって、ダルシーの法則は、このような堆積物では常に有効とは限りません。市販の円形パイプ内の流れでは、レイノルズ数が2000未満の場合は流れが層流になり、4000を超える場合は乱流になりますが、一部の堆積物では、レイノルズ数が1未満の場合に流れが層流になることがわかっています。[ 18 ]