含随机系数的三维准静态电-孔弹性方程的系综蒙特卡洛有限元方法

《Applied Numerical Mathematics》：Ensemble-based Monte Carlo finite element methods for three-dimensional quasi-static electroporoelasticity equations with random coefficients

　　

电?孔弹性方程(electroporoelasticity equations)描述流体饱和多孔介质中的震电耦合(seismoelectric coupling)效应。通常该模型会受到测量误差及外部环境不确定性等随机因素的影响。本文研究具有不确定输入的电?孔弹性方程的系综数值逼近。研究人员首先建立了适定性(well?posedness)，随后应用系综算法(ensemble algorithm)、半隐式Euler方法及有限元法(finite element method, FEM)构造全离散数值格式，并证明其稳定性。研究人员在均方(mean?square)范数下推导了误差估计，即精确解与系综数值解之差的估计，获得时间与空间方向的一阶收敛性。研究人员还研究了精确解期望与系综数值解样本均值之差的误差估计。最后，研究人员通过数值实验验证了理论分析。

论文解读：含随机系数的三维准静态电?孔弹性方程的系综蒙特卡洛有限元方法（Ensemble?based Monte Carlo Finite Element Methods for Three?Dimensional Quasi?Static Electroporoelasticity Equations with Random Coefficients）——发表于《Applied Numerical Mathematics》

---

一、研究背景与意义

电?孔弹性方程（electroporoelasticity equations，亦称震电耦合方程）由Maxwell方程组与Biot孔弹性方程经电动耦合系数耦合而成，用于刻画流体饱和多孔介质中地震波与电磁场的相互作用（seismoelectric coupling phenomenon），在油气勘探和地震预警等领域有重要应用。然而实际多孔介质的渗透率、水力传导率(hydraulic conductivity)等参数因孔隙结构复杂、测量误差及温压等外部因素随机波动而具有不确定性，导致需考虑含随机系数(random coefficients)的随机偏微分方程(random PDE, RPDE)。传统求解随机PDE的方法如多项式混沌展开(polynomial chaos)、随机Galerkin(stochastic Galerkin)法等在处理多样本(multiple realizations)时计算成本高昂。系综算法(ensemble algorithm, 最初由Jiang和Layton提出用于Navier?Stokes流系综，后由Luo和Wang发展用于含随机扩散系数的抛物方程)通过将变化系数用系综平均近似，使多组右端向量共用同一系数矩阵，可大幅降低存储与计算开销。本文首次将系综?Monte Carlo?有限元方法(ensemble?based Monte Carlo finite element method)推广至三维拟静态(quasi?static)电?孔弹性方程，并给出严格的适定性与误差分析。

---

二、主要关键技术方法

研究人员考虑三维有界凸多面体域D?R³上含随机空间变化水力传导率κ(ω,x)的拟静态电?孔弹性初边值问题。关键技术包括：(1)假设随机系数满足有界性、可测性等条件(H1–H4)，在概率空间建立强解(mean?square)的适定性——通过变分形式与能量估计证明几乎必然(a.s.)唯一解及均方有界性；(2)时间方向采用半隐式向后Euler离散，空间方向电场E用Nédélec棱元(curl?conforming edge elements)、磁场H用L²元、位移u和压力p用Lagrange元(P₁)，随机系数κ(ω)在各样本中以系综平均E[κ]代入刚度矩阵，使J个样本共享同一系数矩阵仅需一次LU分解、多次回代不同右端项——即系综Monte Carlo FEM格式；(3)在均方范数下推导全离散格式的存在唯一性、稳定性及误差估计：时间步长Δt和单元尺寸h下获O(Δt+h)一阶收敛，并结合Monte Carlo采样误差O(N^?1/2)给出期望值与样本均值误差界；(4)数值验证取单位立方体D=[0,1]³，κ=1+(1+ξ(ω))sin(xyz)其中ξ~U[2.7,3.3]，用FEniCS平台在不同网格与样本数下检验收敛阶。

---

三、研究结果

Preliminaries（预备知识）

引入随机拟静态电?孔弹性变分形式（含电场E(ω,t,x)、磁场H(ω,t,x)、固相位移u(ω,t,x)、孔隙压力p(ω,t,x)及电动耦合项），给出随机系数κ(ω,x)的四点基本假设(H1–H4)：几乎必然正且有界、可测、均方有界及其逆的均方有界。明确解函数在Sobolev空间H₀(curl)×L²(D)×H¹₀×H¹₀中取值。

Regularity estimates and well?posedness（正则性估计与适定性）

研究人员通过引理3.1–3.4推导解的能量不等式，证明对任意固定ω∈Ω几乎必然存在唯一解且连续依赖于初值与参数；进一步取期望并利用It?等距思想得均方(mean?square)先验估计‖E(ω)‖_{L²ω(H(curl))}及‖u(ω)‖_L²ω(H¹)等有界，从而确立随机电?孔弹性方程在均方意义下的适定性(well?posedness)。

Ensemble?based Monte Carlo finite element approximation（系综蒙特卡洛有限元逼近）

将时间区间[0,T]作均匀剖分，用半隐式Euler离散时间导数，空间用Nédélec元离散E、Lagrange P₁元离散u,p，磁场H显式更新。系综格式以E[κ_h]（随机水力传导率的有限元插值之期望）替代各样本变化的κ_h^(j)构造共同刚度矩阵A_h=A(E[κ_h])，对每个样本j仅右端向量F^(j)_n不同。研究人员证明该全离散系综格式存在唯一解且满足离散稳定性（能量不增长）。在均方范数下得误差估计‖(E,u,p)?(E_h^(j),u_h^(j),p_h^(j))‖_L²ω(·)≤C(Δt+h)；结合Monte Carlo样本均值1/J Σ_j=1^J(E_h^(j),u_h^(j),p_h^(j))与精确解期望E[(E,u,p)]之差的误差界为O(Δt+h+J^?1/2)。

Numerical experiments（数值实验）

在单位立方体域D=[0,1]³取已知精确解构造强制项，随机系数κ(ω,x)=1+(1+ξ(ω))sin(xyz)，ξ服从[2.7,3.3]上均匀分布。分别取网格尺寸h=1/4,1/8,1/12,1/16及不同Monte Carlo样本数J，计算L²误差。数值结果表明电场E、位移u及压力p的时间与空间收敛阶均接近1，与理论预测吻合；样本均值随J增大趋近期望值，验证系综Monte Carlo FEM的有效性及误差估计理论。

---

四、讨论与结论（翻译结论部分）

本研究首先建立了含随机系数的三维拟静态电?孔弹性方程的适定性理论。随后研究人员推导了系综有限元逼近格式，证明了数值格式的适定性与收敛阶。研究人员还研究了精确解期望与系综数值解样本均值之间的差异，获得了误差估计。数值实验结果支持了理论分析。该系综?Monte Carlo?FEM框架为具随机参数的多物理场耦合问题提供了一种高效、理论保证的数值求解途径，可推广至其他含随机系数的耦合PDE系统。

（注：CRediT声明—Xuan Liu: 初稿撰写、方法论、调研、形式分析；Yongkui Zou, Shimin Chai, Amnon J. Meir: 审阅修改、方法论、调研。作者声明无利益冲突。）