随机三重孔隙热弹性中的波传播：基于Lord–Shulman理论的不确定性解析框架

《Scientific Reports》：Wave propagation in stochastic triple-porosity thermoelasticity

　　

本文提出了一个在Lord–Shulman（L–S）理论框架内、考虑随机参数不确定性的广义热弹性介质波传播综合解析框架，该介质具有三重孔隙结构。该模型在经典双重孔隙理论基础上增加了独立的第三孔隙网络，从而能够更真实地刻画裂缝性油藏、孔洞碳酸盐岩和生物组织等复杂多孔材料。研究人员考虑一个承受热-力学载荷的均质各向同性半空间，且材料属性具有随机变化。利用简正模态分析（normal mode analysis）解析求解了由位移、温度和三个孔隙压力场组成的耦合方程组。在变换域中推导出物理量的精确表达式，并通过Durbin傅里叶级数反演算法将各场量还原至物理时间域，且采用Korrektur方法和ε-算法加速，相关参数已在第4.3节中明确给出。以硅作为数值算例材料——仅用作表征良好的数值试验台——并采用层级有序的三重孔隙参数。研究中引入一种新的随机分析方法，采用蒙特卡洛模拟（Monte Carlo simulation）和摄动法量化从材料参数到场变量的不确定性传播；随机材料属性被表示为具有给定相关长度的均匀随机场，在与本文半空间问题相关的长相关极限下退化为相关随机变量。文中定量比较了确定性、一阶摄动、二阶摄动和蒙特卡洛预测结果，并作为特例恢复了正确的单孔隙和双孔隙极限。图形对比研究了第三孔隙网络对波传播和场分布的影响，在不同时刻与双孔隙及单孔隙情形进行了对照。分析表明，三重孔隙和参数不确定性会引起显著的定量差异，尤其是在应力集中、波衰减和热扩散模式方面。

一、研究背景、现存问题与研究意义

经典热弹性理论基于Fourier热传导定律，预言无限大热波速度，这一非物理结果促使Lord–Shulman（L–S）理论引入单弛豫时间，Green–Lindsay和Green–Naghdi等理论进一步建立了有限速度热传导与热弹性波传播框架。与此同时，孔隙介质理论从含孔洞线弹性、双重孔隙模型逐步发展到三重孔隙模型，用于刻画裂缝性储层、孔洞碳酸盐岩和生物组织等具有多尺度孔隙网络的材料。然而，已有三重孔隙热弹性研究几乎都是确定性的，未量化材料参数随机性带来的误差；随机孔隙力学也主要停留在双重孔隙层面，无法回答加入更慢的第三孔隙网络后不确定性如何在各传播支路间重新分配；此外，低成本摄动或代理模型在广义热弹性孔隙波问题中很少与蒙特卡洛模拟进行系统基准比较。针对上述空白，研究人员建立了同时包含三重孔隙耦合和随机参数不确定性的L–S广义热弹性半空间解析框架，并通过单孔隙→双孔隙→三孔隙的受控过渡，分离第三孔隙网络的增量贡献。该研究的意义在于为多孔介质波传播提供可复现的随机解析解，并给出低成本不确定性的可靠替代方案。

二、研究内容与主要结论

研究人员推导了均质各向同性三重孔隙热弹性半空间的控制方程，三个孔隙网络分别对应宏观孔隙、微观孔隙和纳米孔隙，通过孔隙压力、流体交换和温度场相互耦合，得到10次特征方程并给出全部系数闭式代数表达式。随机模型将弹性、热、渗透率、存储和交叉耦合参数设为截断正态或对数正态分布，采用二阶摄动解、蒙特卡洛模拟和Sobol敏感性分析量化不确定性传播。结果显示：三重孔隙显著改变波衰减、应力松弛和热扩散；随机不确定性随深度累积，孔隙压力场的变异系数高于位移场；热弛豫时间τ₀和孔弹性耦合系数b、d₁是主要不确定性来源；二阶摄动解能以约0.6%的计算成本复现蒙特卡洛统计量，且在输入变异系数约6%以内保持可靠。该工作将随机波传播理论与三重孔隙热弹性结合，为裂缝性储层和生物组织的波场解释提供了可量化不确定性的建模工具。

三、关键技术与方法

研究人员采用简正模态分析与Helmholtz分解将耦合偏微分方程组化为常微分方程组，推导10次特征方程并利用边界条件确定待定常数；通过Durbin傅里叶级数反演算法配合Korrektur方法和ε-算法将变换域解还原至时间域。随机分析方面，采用截断正态分布和对数正态分布描述随机参数，以Karhunen–Loève展开联系随机场与随机变量，并实施一阶、二阶摄动展开和N=10,000独立样本的蒙特卡洛模拟，最后用Sobol指数进行全局敏感性分析。研究未涉及实验样本队列，硅仅作为数值试验台材料。

四、研究结果

验证单、双孔隙极限。研究人员通过将第三网络解耦系数置零，使特征方程分别退化为既有双孔隙和单孔隙理论，数值结果与独立实现的理论对达到机器精度，确认新增孔隙网络的引入在极限情形下完全一致。

时间演化。与单、双孔隙相比，三重孔隙引入更高频振荡和额外模式耦合；垂直位移峰值比单孔隙约高40%，表明活化额外孔隙网络降低等效刚度。三种孔隙压力按宏观、微观、纳米顺序依次衰减，1/e衰减时间之比为1:2.4:6.1，与存储系数和渗透率之比的理论预期一致，验证了参数标定的内部一致性。

热弛豫时间效应。增大τ₀使时间曲线更光滑、峰值更低、响应持续时间更长；温度峰值约下降18%，脉冲半宽约增大32%，而热能基本守恒，体现了L–S模型有限热波速度对热-力耦合的调控作用。

三维空间分布。三维云图显示位移随深度快速衰减并形成干涉分层，三重孔隙增强热扩散并使温度场更缓慢衰减；不同尺度孔隙压力在空间上呈现分离模式，纳米孔隙压力向更深部扩展，产生驱动跨孔隙流体交换的三维压力梯度。

随机分析。样本路径和方差分析表明，场变量变异系数随深度增加，孔隙压力场比位移场具有更高不确定性，变异系数约15%–25%；位移与温度样本呈负相关，反映热弹性耦合项的作用；蒙特卡洛估计的均值相对标准误差低于0.3%，标准差相对误差约0.7%。Sobol指数识别出τ₀、b、d₁为主导参数，τ₀约占温度方差35%；最大应力95%分位数比确定性预测约高15%。

摄动解与蒙特卡洛模拟对比。在输入变异系数5%时，一阶摄动低估标准差7.5%–17.4%，二阶摄动误差降至1.2%–4.9%；二阶摄动约以0.6%计算成本获得蒙特卡洛统计量，误差随支路扩散性增强而增大，因而在ε约6%以内可采用摄动方案。

敏感性排序稳健性。改变分布族、输入变异系数和参数间相关性后，排序τ₀>b>d₁>κ>d₂保持不变，τ₀的一阶Sobol指数为0.35±0.03，表明主导参数来自三重孔隙物理本身而非统计假设。

第三孔隙网络的增量贡献。通过连续变化第三网络的耦合系数d₂、渗透率γ₂和存储系数a₃，发现d₂控制力学软化，γ₂控制波衰减，a₃控制长时孔隙压力记忆；三种效应作用于不同可观测物理量，因此无法通过重调双孔隙模型同时复现。

五、讨论与结论

讨论部分明确指出，该研究使用硅仅作为数值试验台，未声称实验验证，均匀半空间假设限制在相关长度远大于波长的情形，理论为线性小变形且采用弱耦合流体交换，二阶摄动仅在输入变异系数约6%以内定量可靠。结论部分可概括为：研究建立了L–S理论下三重孔隙热弹性波传播的随机解析框架；10,000样本蒙特卡洛模拟统计收敛；τ₀、b、d₁是主导不确定来源；孔隙压力场不确定性高于位移场；三重孔隙产生增强波衰减、降低相速度、延长孔隙压力扩散和复杂应力分布等多尺度效应；概率超越界限为设计提供约15%的应力裕度提示；二阶摄动解可低成本复现蒙特卡洛结果；敏感性排序对分布族、变异系数和相关符号保持稳健；公式能精确退化到单、双孔隙极限；第三孔隙网络的耦合、渗透率和存储效应可分离且不可被双孔隙模型替代；实验鉴定三重孔隙交叉耦合系数仍是开放问题。