基于投影的三维离散裂缝网络低维参数化生成对抗网络方法

《Journal of Rock Mechanics and Geotechnical Engineering》：Projection-based low-dimensional parameterization method for three-dimensional discrete fracture networks through generative adversarial network

　　

研究人员提出了一种基于深度学习辅助的投影式低维参数化方法，用于从高保真的三维离散裂缝网络（Discrete Fracture Network, DFN）生成中降低参数空间维度。由于直接生成三维DFN计算成本极高，研究将任务分解为二维DFN生成及二维与三维DFN之间的投影关系建立两部分。研究人员将二维DFN视为三维DFN在二维平面上的投影，引入裂缝倾角（dip angle）与投影距离（projection distance）两个参数描述该关系。研究采用带梯度惩罚的Wasserstein生成对抗网络（WGAN-GP）分别生成二维DFN、倾角图及投影距离图，并通过反向投影构建三维DFN。实验表明，该方法可从40维潜空间生成含50条裂缝的三维DFN（传统方法需300个参数），从72维潜空间生成含100条裂缝的三维DFN（传统方法需600个参数）。生成的DFN在裂缝连通性及流体流动特征方面均表现出与训练样本一致的高保真特性，同时保留了足够的多样性以支持后续DFN反演。

该研究发表于《Journal of Rock Mechanics and Geotechnical Engineering》，针对地下能源开发、储存及废物处置中的裂缝网络表征难题展开。当前，地下裂缝几何形态复杂且观测数据稀缺，导致反演问题病态性强，传统三维离散裂缝网络（DFN）参数化方法维度过高，难以高效求解。现有低维参数化方法多局限于二维DFN，无法捕捉三维空间的拓扑与连通性差异，因此亟需发展高效的三维DFN低维参数化框架。

研究人员提出了一种基于投影的生成式方法，将三维DFN构建转化为二维DFN生成及投影参数估计两个子任务。研究假设裂缝为圆形（Baecher盘模型），通过正交投影将三维裂缝映射为二维椭圆，并以长轴作为投影线。引入倾角图和距离图分别记录每条裂缝的倾角和投影距离，结合二维DFN实现三维重建。为降低维度，研究人员训练三个独立的WGAN-GP模型，分别从低维潜空间生成二维DFN、倾角图和距离图。在50条裂缝场景下，潜空间维度为40（32维用于二维DFN，4维用于倾角图，4维用于距离图）；100条裂缝场景下，潜空间维度为72。生成的DFN在统计分布、连通性、主干结构及流体流动特征上均与训练样本高度一致，且潜空间连续变化可使DFN几何平滑过渡，有利于反演优化。

关键技术方法

研究人员采用带梯度惩罚的Wasserstein生成对抗网络（WGAN-GP）作为核心生成模型，分别构建二维DFN、倾角图和距离图的生成器与判别器。二维DFN训练样本基于幂律分布（fracture length）和von Mises分布（fracture orientation）合成，倾角图与距离图则通过高斯滤波与直方图投影法生成，以确保空间自相关性。三维DFN重建通过反向投影实现，流体流动模拟采用嵌入式离散裂缝方法（EDFM）在GEOS数值模拟器中完成，连通性与拓扑分析使用ADFNE软件。

研究结果

4.1 案例配置

研究设置8种案例，涵盖50条与100条裂缝场景，并调整倾角均值（μ_dip）、浓度参数（κ_dip）及投影距离分布，验证方法的通用性。

4.2 三维DFN生成

随机生成的DFN在空间分布上呈现显著差异，裂缝可聚集于域中心或边界，表明生成模型具有高多样性。

4.3 潜参数对三维DFN的影响

潜参数的连续变化可使DFN几何平滑演变，这种特性源于WGAN-GP的Lipschitz连续性及倾角/距离图的空间自相关性，有利于反演中的参数微调。

4.4 统计分析

4.4.1 连通性

100条裂缝场景的连通概率（P_c）显著高于50条裂缝场景，距离分布集中于5–15 m时连通性更高，倾角均值50°时连通性降低，浓度参数影响较小。

4.4.2 拓扑性质

生成DFN的主干长度（l_bbn）与交点线总长度（l_inter）与训练样本一致，主干覆盖度（Cov_bbn）在多数案例中超过80%，表明生成模型能有效复现三维裂缝网络的拓扑特征。

4.4.3 流动特征

等效渗透率（κ_eff）与流速相关维数（D_flow）的分布与训练样本吻合，100条裂缝场景的κ_eff与D_flow显著更高，且拓扑指标与流动指标呈正相关。

讨论与结论

该方法可扩展至矩形或椭圆形裂缝，只需增加投影角图与轴比图。当前局限在于二维DFN、倾角图与距离图相互独立，未融入物探与水力测试等先验约束，未来可通过条件生成对抗网络（cGAN）引入此类信息。研究实现了三维DFN参数化维度的显著降低（约90%），生成结果在统计与物理特征上均符合预期，为地下裂缝表征与反演提供了高效工具。