考虑饱和-非饱和渗流的水封石油洞库围岩水位变化分析.pdfVIP

DOI编码：10.13646/j.cnki.42-1395/u.2025.04.043

考虑饱和-非饱和渗流的水封石油洞库围岩水位变化分析

蒋文健

（南昌市水利规划设计院有限公司，江西南昌330038）

摘要：考虑围岩饱和-非饱和渗流过程，采用数值分析方法研究洞库计算模型范围及复杂覆盖层渗透系数对水封石

油洞库地下水位降深、影响半径以及涌水量的影响。研究结果表明：（1）围岩渗透系数越大，洞库围岩水位降幅越

-3

大，同时降低速度越快，当围岩渗透系数小于10m/d量级时，围岩水位可形成较好的天然水封条件；（2）计算

模型范围对预测洞库上方水盖层厚度的计算影响不大，但对洞库围岩水位分布形态影响较大，并且其影响程度主要

取决于围岩渗透系数，渗透系数越小，其影响越小，反之，越大；（3）计算模型范围一定条件下，水位影响半径随

围岩渗透系数的增大而增大，只有计算模型范围小于影响半径条件下才能获得较准确地渗流场计算结果；（4）地下

洞库涌水量受洞库围岩渗透性的显著影响，洞库涌水量随时间推移而逐渐减小，且在0.5年以内涌水量降幅最大。

关键词：地下水封洞库；渗流场；渗透系数；模型范围；数值分析

中图分类号：TU443文献标识码：A文章编号：1006-7973（2025）04-0117-03

地下水封石油是利用水位以下洞室群，同时借助天然水满足达西定律。

幕或者人工水幕系统封闭油气外泄的石油储存技术[1-2]。洞库为了反映洞周围岩疏干条件下岩体的非饱和渗流过程，

储油的关键在于围岩拥有稳定且满足水封条件的地下水位，通过考虑岩体非饱和渗透系数实现。岩体渗透系数与饱和度

准确预测围岩水位变化对水封措施的设计极为重要。密切相关，饱和度越大渗透系数也越大，反之越小。同时，

国内外学者针对水封油库围岩水位变化开展了系列研饱和度与非饱和岩土体的基质吸力有关[7]，饱和度越低，基

究。Kim等[3]结合水封石油洞库实测数据，对围岩渗流场的质吸力越大，反之越小。由此，岩体饱和度与基质吸力的关

影响因素展开研究，发现石油油气压力对围岩渗流场的形成系可以通过中间变量孔隙比转化为体积含水率与基质吸力的

影响很大；李术才等[4]以我国首个大型地下水封石油为背景，关系，进而可建立非饱和渗透系数与体积含水率、基质吸力

建立洞室三维数值模型，分析了洞室开挖过程中围岩水位变的关系模型。其中较为著名的关系模型有Fredlund模型[8]、

化情况；平洋等[5]依托黄岛地下水封石油工程，建立三维数VanGenuchten（即V-G模型）[9]。

值模型对库区初始渗流场进行了反演分析；施雪松[6]等考虑二、数值研究方案

岩体渗透各项异性，建立单洞库数值模型，分析了洞室围岩为研究不同模型范围及不同渗透系数对洞库渗流场分布

渗流场及涌水量变化。的影响，建立不同模型尺寸的计算模型，考虑洞库上覆层不

文中拟考虑硐壁围岩饱和-非饱和渗流过程，采用数值分同渗透系数，研究不同模型范围及不同渗透系数对施工期洞

析方法研究硐室上覆岩层不同渗透系数及洞库模型范围对水库地下水水位、影响半径、降深以及涌水量的影响。

封石油洞库施工期地下水位、影响半径、降深以及涌水量的1．数值模型几何尺度方案

影响，以期为地下水封洞库围岩水位变化预测及水幕系统的为简化计算，计算模型仅考虑单个洞库，洞库大小宽

设计提供