吉林大学学报(地球科学版)

• 地质工程与环境工程 • 上一篇    下一篇

深部咸水层CO2地质封存数值模拟参数的全局敏感性分析:以苏北盆地盐城组为例

郑菲,施小清,吴吉春,赵良,陈旸   

  1. 南京大学地球科学与工程学院,南京210093
  • 收稿日期:2013-05-04 出版日期:2014-01-26 发布日期:2014-01-26
  • 通讯作者: 施小清(1979-),男,副教授,博士,主要从事反应溶质运移模拟研究 E-mail:shixq@nju.edu.cn
  • 作者简介:郑菲(1987-),女,博士研究生,主要从事多相流数值模拟方面的研究,E-mail:xiaomi2008happy@163.com
  • 基金资助:

    国家自然科学基金项目(41172206);江苏省自然科学基金项目(BK2012313);江苏省产学研联合创新资金计划项目(BY2010136)

Global Parametric Sensitivity Analysis of Numerical Simulation for CO2 Geological Sequestration in Saline Aquifers: A Case Study of Yancheng Formation in Subei Basin

Zheng Fei, Shi Xiaoqing, Wu Jichun, Zhao Liang, Chen Yang   

  1. School of Earth Sciences and Engineering, Nanjing University, Nanjing210093, China
  • Received:2013-05-04 Online:2014-01-26 Published:2014-01-26

摘要:

基于对苏北盆地盐城组下段砂岩储层概化建立二维径向剖面模型,运用TOUGH2/ECO2N程序模拟深部咸水层中CO2迁移分布过程。采用定性(Morris法)和定量(Sobol和EFAST法)全局敏感性分析方法,以储层中注入井处的压力、气相CO2总量及CO2气相羽扩散距离作为模型响应变量,对储层的水平渗透率、孔隙度、残余液体饱和度、孔隙分布指数、压缩系数、进气压力的倒数和盐度7个参数进行全局敏感性分析,讨论了咸水层的水文地质参数、盐度及其他模型参数对CO2封存运移泄漏过程的影响。Morris和Sobol法的敏感性分析结果都表明,对于不同的响应变量,参数的敏感性排序不同:以注入井处的压力为响应变量时,孔隙度的敏感性最高;以气相CO2总量和CO2羽扩散距离为响应变量时,水平渗透率的敏感性最高。EFAST 1阶及总敏感度分析结果与Sobol分析结果基本一致,但EFAST法相对Sobol法计算更高效稳健,需要的样本数较少。

关键词: CO2地质封存, 全局敏感性分析, Morris法, Sobol法, EFAST法, 苏北盆地

Abstract:

TOUGH2/ECO2N code was used to simulate the complex coupled processes taking place during and after CO2 injection in the deep saline sandstone aquifer of Lower Yancheng Formation in Subei basin. Morris, Sobol and EFAST methods were used to analyze the global sensitivity of parameters, i.e., kx, n, Slr,λ, C, 1/p0 and S, on three response variables: 1)the pressure located in the injection well; 2)the amount of gas-phase CO2 and 3)the spread distance of CO2 gas plume. Results from the Morris and Sobol methods show that n is the most sensitive parameter for the injection well pressure, kx is the most sensitive parameter both for the amount of gas-phase CO2 and the spread distance of CO2 gas plume. In summary, the sensitivity orders of these parameters are totally different for different corresponding variables. Comparative results from EFAST and Sobol show that the main and total effect obtained from the two quantitative sensitivity analysis methods are basically identical, however, EFAST is computationally more efficient than Sobol in numerical experiment.

Key words: CO2 geological sequestration, global sensitivity analysis, Morris, Sobol, EFAST, Subei basin

中图分类号: 

  • P641
[1] 蔡潇, 夏威, 马晓东, 臧素华, 丁安徐, 刘玉霞, 花彩霞, 朱一川, 李辉. 苏北盆地溱潼凹陷古近系阜宁组页岩油储集空间特征及其影响因素[J]. 吉林大学学报(地球科学版), 2024, 54(5): 1482-1493.
[2] 莫绍星, 龙星皎, 李瀛, 郑菲, 施小清, 张可霓, 赵良. 基于TOUGHREACT-MP的苏北盆地盐城组咸水层CO2矿物封存数值模拟[J]. 吉林大学学报(地球科学版), 2014, 44(5): 1647-1658.
[3] 罗建男,卢文喜,康柱,宋晶,张一博. 表面活性剂强化含水层修复过程的全局灵敏度[J]. 吉林大学学报(地球科学版), 2013, 43(2): 558-0563.
[4] 杨峰田, 庞忠和, 王彩会, 段忠丰, 罗璐, 李义曼. 苏北盆地老子山地热田成因模式[J]. J4, 2012, 42(2): 468-475.
[5] 刘长利, 朱筱敏, 王兴明, 任红民, 李树静. 苏北盆地海安凹陷泰州组层序地层及沉积特征[J]. J4, 2010, 40(3): 519-526.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 尤敏鑫,刘建民. 同位素地球化学在峨眉山大火成岩省研究中的应用现状与进展[J]. 吉林大学学报(地球科学版), 2014, 44(4): 1231 -1243 .
[2] 杨春梅, 李洪奇,陆大卫,张方礼,高 原,邵英超. 不同驱替方式下岩石电阻率与饱和度的关系[J]. J4, 2005, 35(05): 667 -671 .
[3] 祝洪臣,张炯飞,权 恒. 大兴安岭中生代两期成岩成矿作用的元素、同位素特征及其形成环境[J]. J4, 2005, 35(04): 436 -0442 .
[4] 朱建伟, 赵刚, 刘博, 郭巍, 成俊. 油页岩测井识别技术及应用[J]. J4, 2012, 42(2): 289 -295 .
[5] 陈力,梁海安,张文娟,荣帆. 模糊数学方法在城市工程地质环境区划中的应用--以抚顺市城区为例[J]. J4, 2008, 38(5): 837 -0840 .
[6] 高桂梅,苏 克,王文颖,甘树才,刘招君. 吉林省桦甸油页岩中稀土元素和微量元素的研究[J]. J4, 2006, 36(6): 974 -0979 .
[7] 吴孔运,蒋忠诚,叶 晔. 不同植物群落对灰岩试块溶蚀速率的影响[J]. J4, 2007, 37(5): 967 -0971 .
[8] 周彦章,迟宝明,刘中培. 山东夏甸金矿床充水机理构造控制模式[J]. J4, 2008, 38(2): 255 -0260 .
[9] 张渊,刘连登,孙景贵,陈国华,张洪喜,闫复传,杨开春. 胶东西北部黄埠岭金矿床两期次叠加成矿[J]. J4, 2008, 38(1): 21 -0026 .
[10] 鲁程鹏, 束龙仓, 苑利波, 张蓉蓉, 黄币娟, 王彬彬. 基于示踪试验求解岩溶含水层水文地质参数[J]. J4, 2009, 39(4): 717 -721 .