数值模拟技术简介 数值模拟主要过程和步骤
\u6570\u503c\u6a21\u62df\u7684\u6570\u503c\u6a21\u62df\u7684\u53d1\u5c55\u53f2\u6570\u503c\u6a21\u62df\u6280\u672f\u8bde\u751f\u4e8e1953\u5e74Bruce G.H\u548cPeacemanD.W\u6a21\u62df\u4e86\u4e00\u7ef4\u6c14\u76f8\u4e0d\u7a33\u5b9a\u5f84\u5411\u548c\u7ebf\u5f62\u6d41\u3002\u53d7\u5f53\u65f6\u8ba1\u7b97\u673a\u80fd\u529b\u53ca\u89e3\u6cd5\u9650\u5236\uff0c\u6570\u503c\u6a21\u62df\u6280\u672f\u53ea\u662f\u521d\u6b65\u5e94\u7528\u4e8e\u89e3\u4e00\u7ef4\u4e00\u76f8\u95ee\u9898\u3002\u4e24\u76f8\u6d41\u52a8\u6a21\u62df\u8bde\u751f\u4e8e1954\u5e74\uff0cWest W J\u548cGarvin W.W\u6a21\u62df\u4e86\u6cb9\u85cf\u4e0d\u7a33\u5b9a\u4e24\u76f8\u6d41\u30021955\u5e74Peaceman\u4e0eRachford\u7814\u53d1\u7684\u4ea4\u66ff\u9690\u5f0f\u89e3\u6cd5(ADI)\u662f\u6570\u503c\u6a21\u62df\u6280\u672f\u7684\u91cd\u5927\u7a81\u7834\u3002\u8be5\u89e3\u6cd5\u975e\u5e38\u7a33\u5b9a,\u800c\u4e14\u901f\u5ea6\u5feb\uff0c\u6240\u4ee5\u8fc5\u901f\u5728\u5305\u62ec\u77f3\u6cb9,\u6838\u7269\u7406\uff0c\u70ed\u4f20\u5bfc\u7b49\u9886\u57df\u5f97\u5230\u5e7f\u6cdb\u5e94\u7528\u30021958\u5e74Douglas,Jim\u548cBlair,P.M\u7b2c\u4e00\u6b21\u8fdb\u884c\u4e86\u8003\u8651\u6bdb\u7ba1\u538b\u529b\u6548\u679c\u7684\u6c34\u9a71\u6a21\u62df\u30021959\u5e74\uff0cDouglas Jim\u548cPeaceman D.W\u7b2c\u4e00\u6b21\u8fdb\u884c\u4e86\u4e24\u7ef4\u4e24\u76f8\u6a21\u62df\uff0c\u8fd9\u6807\u5fd7\u7740\u73b0\u4ee3\u6570\u503c\u6a21\u62df\u6280\u672f\u7684\u5f00\u59cb\u3002\u5728\u4ed6\u4eec\u7684\u6a21\u62df\u5668\u91cc\u5168\u9762\u8003\u8651\u4e86\u76f8\u5bf9\u6e17\u900f\u7387\uff0c\u7c98\u5ea6\uff0c\u5bc6\u5ea6\uff0c\u91cd\u529b\u53ca\u6bdb\u7ba1\u538b\u529b\u7684\u5f71\u54cd\u300260\u5e74\u4ee3\u6570\u503c\u6a21\u62df\u6280\u672f\u7684\u53d1\u5c55\u4e3b\u8981\u5728\u6570\u503c\u89e3\u6cd5\uff0c\u7b2c\u4e00\u4e2a\u6709\u6548\u7684\u6570\u503c\u6a21\u62df\u89e3\u6cd5\u5668\u662f1968\u5e74Stone\u63a8\u51fa\u7684SIP(Strong Implicit Procedure)\u3002\u8be5\u89e3\u6cd5\u53ef\u4ee5\u5f88\u597d\u5730\u7528\u6765\u6a21\u62df\u975e\u5747\u8d28\u6cb9\u85cf\u548c\u5f62\u72b6\u4e0d\u89c4\u5219\u6cb9\u85cf\u3002\u53e6\u4e00\u4e2a\u7a81\u7834\u662f\u65f6\u95f4\u9690\u5f0f\u6cd5\uff0c\u8be5\u65b9\u6cd5\u53ef\u4ee5\u7528\u6765\u6709\u6548\u7684\u89e3\u9ad8\u6d41\u901f\u95ee\u9898\uff0c\u6bd4\u5982\u9525\u8fdb\u95ee\u9898\u300260\u5e74\u4ee3\u5176\u4ed6\u65b9\u9762\u7684\u53d1\u5c55\u8fd8\u67091967\u5e74Coats K.H\u548cNielsen R.L\u9996\u6b21\u8fdb\u884c\u4e86\u4e09\u7ef4\u4e24\u76f8\u6a21\u62df\uff0c\u800c\u4e14\u63d0\u51fa\u4e86\u5782\u76f4\u5e73\u8861\u548c\u62df\u76f8\u5bf9\u6e17\u900f\u7387\u53ca\u6bdb\u7ba1\u538b\u529b\u65b9\u6cd5\u30021968\u5e74Breitenbach E.A\u53d1\u8868\u4e86\u4e09\u7ef4\u4e09\u76f8\u6a21\u62df\u89e3\u6cd5\u3002Stone\u572870\u5e74\u4ee3\u53d1\u8868\u4e86\u4e09\u76f8\u76f8\u5bf9\u6e17\u900f\u7387\u6a21\u578b\uff0c\u7531\u6cb9\u6c34\u548c\u6cb9\u6c14\u4e24\u76f8\u76f8\u5bf9\u6e17\u900f\u7387\u8ba1\u7b97\u6cb9\u3001\u6c14\u3001\u6c34\u4e09\u76f8\u6d41\u52a8\u65f6\u7684\u76f8\u5bf9\u6e17\u900f\u7387\uff0c\u8be5\u6280\u672f\u73b0\u5728\u8fd8\u5e7f\u4e3a\u5e94\u7528\u300270\u5e74\u4ee3\u53e6\u4e00\u9879\u4e3b\u8981\u6210\u5c31\u662fPeaceman\u63d0\u51fa\u7684\u4ece\u7f51\u683c\u538b\u529b\u6765\u786e\u5b9a\u4e95\u5e95\u6d41\u538b\u7684\u6821\u6b63\u65b9\u6cd5\uff0c\u53ca\u73b0\u5728\u901a\u7528\u7684Peaceman\u65b9\u7a0b\u3002\u5728\u89e3\u6cd5\u65b9\u9762\u7684\u53d1\u5c55\u662f\u91c7\u7528\u4e86\u6b63\u4ea4\u52a0\u901f\u7684\u8fd1\u4f3c\u5206\u89e3\u6cd5\u300270\u5e74\u4ee3\u5728\u7ec4\u5206\u548c\u70ed\u91c7\u6a21\u62df\u65b9\u9762\u4e5f\u53d6\u5f97\u5f88\u5927\u8fdb\u5c55\uff0c1973\u5e74Nolen J.S\u63cf\u8ff0\u4e86\u8003\u8651\u6cb9\u6c14\u4e2d\u95f4\u7ec4\u5206\u5206\u5e03\u7684\u7ec4\u5206\u6a21\u62df\uff0cCook\u63d0\u51fa\u53d8\u9ed1\u6cb9\u6a21\u62df\u6765\u8fdb\u884c\u7ec4\u5206\u6a21\u62df\u3002Shutler\u57281970\u5e74\u53d1\u8868\u4e86\u5bf9\u4e24\u7ef4\u4e09\u76f8\u6a21\u578b\u7684\u84b8\u6c14\u6ce8\u5165\u6a21\u62df\u300270\u5e74\u4ee3\u5728EOR\u65b9\u9762\u4e5f\u53d6\u5f97\u4e86\u6781\u5927\u8fdb\u5c55\u300280\u5e74\u4ee3\u6700\u5927\u7684\u6210\u5c31\u662fAppleyyard J R\u548cCheshire I.M\u53d1\u8868\u4e86\u5d4c\u5957\u56e0\u5f0f\u5206\u89e3\u6cd5\uff0c\u8be5\u89e3\u6cd5\u975e\u5e38\u7a33\u5b9a\u800c\u4e14\u901f\u5ea6\u5feb\uff0c\u662f\u76ee\u524d\u6700\u4e3a\u5e7f\u6cdb\u5e94\u7528\u7684\u89e3\u6cd5\u3002\u6b63\u662f\u57fa\u4e8e\u8be5\u89e3\u6cd5\uff0cCheshire I.M\u4e8e1981\u5e74\u540cJohn Appleyard\u548cJon Holmes\u6210\u7acbECL\u516c\u53f8\uff0c\u5f00\u59cb\u7814\u53d1\u540e\u6765\u4e3b\u5bfc\u6570\u503c\u6a21\u62df\u8f6f\u4ef6\u5e02\u573a\u7684ECLIPSE\u8f6f\u4ef6\u300280\u5e74\u4ee3\u89c1\u8bc1\u7684\u53e6\u4e00\u4e2a\u4e3b\u8981\u53d1\u5c55\u662f\u7ec4\u5206\u6a21\u578b\uff0c\u867d\u7136\u7ec4\u5206\u6a21\u578b\u572860\u5e74\u4ee3\u5c31\u5df2\u7ecf\u63a8\u51fa\uff0c\u4f46\u5f88\u4e0d\u7a33\u5b9a\u300280\u5e74\u4ee3\u63d0\u51fa\u7684\u4f53\u79ef\u5e73\u8861\u548cYong-Stephenson\u65b9\u7a0b\u89e3\u51b3\u4e86\u7ec4\u5206\u6a21\u578b\u7a33\u5b9a\u95ee\u9898\uff0c\u4f7f\u7ec4\u5206\u6a21\u578b\u53ef\u4ee5\u5e7f\u4e3a\u5e94\u7528\u3002Ponting D.K\u63d0\u51fa\u4e86\u89d2\u70b9\u7f51\u683c\u6765\u6a21\u62df\u6a21\u578b\uff0c\u8fd9\u6837\u53ef\u4ee5\u771f\u5b9e\u5730\u63cf\u8ff0\u6cb9\u85cf\u300290\u5e74\u4ee3\u6570\u503c\u6a21\u62df\u7684\u8fdb\u5c55\u4e3b\u8981\u5728\u7c97\u5316\u6280\u672f\uff0c\u5e76\u884c\u8ba1\u7b97\uff0cPEBI\u7f51\u683c\u7b49\u65b9\u9762\u3002Zoltan E.Heinemann\u63d0\u51fa\u4e86PEBI\u7f51\u683c\uff0cPEBI\u7f51\u683c\u7ed3\u5408\u4e86\u6b63\u4ea4\u7f51\u683c\u548c\u89d2\u70b9\u7f51\u683c\u7684\u4f18\u70b9\uff0c\u73b0\u5728\u6b63\u9010\u6e10\u6210\u4e3a\u4e3b\u6d41\u6570\u503c\u6a21\u62df\u7f51\u683c\u4f53\u7cfb\u3002VIP\u4e8e1994\u5e74\u63a8\u51fa\u5e76\u884c\u7b97\u6cd5\uff0cECLIPSE\u4e8e1996\u5e74\u63a8\u51fa\u5e76\u884c\u7b97\u6cd5\u3002CMG\u4e8e2001\u5e74\u63a8\u51fa\u5e76\u884c\u7b97\u6cd5\u3002\u7c97\u5316\u6280\u672f\u7684\u96be\u70b9\u5728\u4e8e\u6e17\u900f\u7387\u7684\u7c97\u5316\uff0c\u57fa\u4e8e\u6d41\u52a8\u8ba1\u7b97\u8fdb\u884c\u7684\u6e17\u900f\u7387\u7c97\u5316\u53ef\u4ee5\u8f83\u771f\u5b9e\u7684\u7b26\u5408\u5730\u8d28\u6a21\u578b\uff0c\u73b0\u5728\u65b0\u7684\u7c97\u5316\u6280\u672f\u8fd8\u5728\u53d1\u5c55\u3002\u3000 21\u4e16\u7eaa\u6570\u503c\u6a21\u62df\u6280\u672f\u53d1\u5c55\u4f53\u73b0\u5728\u4e24\u65b9\u9762\uff0c\u4e00\u65b9\u9762\u662f\u4e00\u4f53\u5316\u6a21\u62df\u6280\u672f\uff0c\u6570\u503c\u6a21\u62df\u5c06\u4e0d\u53ea\u662f\u5bf9\u6cb9\u85cf\u7684\u6a21\u62df\uff0c\u6570\u503c\u6a21\u62df\u5c06\u5bf9\u6cb9\u85cf\uff0c\u4e95\u7b52\uff0c\u5730\u9762\u8bbe\u5907\uff0c\u7ba1\u7f51\u4ee5\u53ca\u6cb9\u6c14\u5904\u7406\u5382\u8fdb\u884c\u4e00\u4f53\u5316\u6a21\u62df\uff0c\u4ece\u800c\u6700\u4f18\u5316\u7ba1\u7406\u6cb9\u7530\u3002\u53e6\u4e00\u65b9\u9762\u662f\u5b9a\u91cf\u8fdb\u884c\u5c5e\u6027\u4e0d\u786e\u5b9a\u6027\u5206\u6790\uff0c\u5b9a\u91cf\u5206\u6790\u5c5e\u6027\u4e0d\u786e\u5b9a\u6027\u5bf9\u8ba1\u7b97\u7ed3\u679c\u7684\u5f71\u54cd\u3002
1\u3001\u9996\u5148\u8981\u5efa\u7acb\u53cd\u6620\u95ee\u9898\uff08\u5de5\u7a0b\u95ee\u9898\u3001\u7269\u7406\u95ee\u9898\u7b49\uff09\u672c\u8d28\u7684\u6570\u5b66\u6a21\u578b\u3002
\u5177\u4f53\u8bf4\u5c31\u662f\u8981\u5efa\u7acb\u53cd\u6620\u95ee\u9898\u5404\u91cf\u4e4b\u95f4\u7684\u5fae\u5206\u65b9\u7a0b\u53ca\u76f8\u5e94\u7684\u5b9a\u89e3\u6761\u4ef6\u3002\u8fd9\u662f\u6570\u503c\u6a21\u62df\u7684\u51fa\u53d1\u70b9\u3002\u6ca1\u6709\u6b63\u786e\u5b8c\u5584\u7684\u6570\u5b66\u6a21\u578b\uff0c\u6570\u503c\u6a21\u62df\u5c31\u65e0\u4ece\u8c08\u8d77\u3002\u725b\u987f\u578b\u6d41\u4f53\u6d41\u52a8\u7684\u6570\u5b66\u6a21\u578b\u5c31\u662f\u8457\u540d\u7684\u7eb3\u7ef4\u2014\u65af\u6258\u514b\u65af\u65b9\u7a0b\uff08\u7b80\u79f0\u65b9\u7a0b\uff09\u53ca\u5176\u76f8\u5e94\u7684\u5b9a\u89e3\u6761\u4ef6\u3002
2\u3001\u5bfb\u6c42\u9ad8\u6548\u7387\u3001\u9ad8\u51c6\u786e\u5ea6\u7684\u8ba1\u7b97\u65b9\u6cd5
\u7531\u4e8e\u4eba\u4eec\u7684\u52aa\u529b\uff0c\u76ee\u524d\u5df2\u53d1\u5c55\u4e86\u8bb8\u591a\u6570\u503c\u8ba1\u7b97\u65b9\u6cd5\u3002\u8ba1\u7b97\u65b9\u6cd5\u4e0d\u4ec5\u5305\u62ec\u5fae\u5206\u65b9\u7a0b\u7684\u79bb\u6563\u5316\u65b9\u6cd5\u53ca\u6c42\u89e3\u65b9\u6cd5\uff0c\u8fd8\u5305\u62ec\u8d34\u4f53\u5750\u6807\u7684\u5efa\u7acb\uff0c\u8fb9\u754c\u6761\u4ef6\u7684\u5904\u7406\u7b49\u3002\u8fd9\u4e9b\u8fc7\u53bb\u88ab\u4eba\u4eec\u5ffd\u7565\u6216\u56de\u907f\u7684\u95ee\u9898\uff0c\u73b0\u5728\u53d7\u5230\u8d8a\u6765\u8d8a\u591a\u7684\u91cd\u89c6\u548c\u7814\u7a76\u3002
3\u3001\u5f00\u59cb\u7f16\u5236\u7a0b\u5e8f\u548c\u8fdb\u884c\u8ba1\u7b97
\u5b9e\u8df5\u8868\u660e\u8fd9\u4e00\u90e8\u5206\u5de5\u4f5c\u662f\u6574\u4e2a\u5de5\u4f5c\u7684\u4e3b\u4f53\uff0c\u5360\u7edd\u5927\u90e8\u5206\u65f6\u95f4\u3002\u7531\u4e8e\u6c42\u89e3\u7684\u95ee\u9898\u6bd4\u8f83\u590d\u6742\uff0c\u6bd4\u5982\u65b9\u7a0b\u5c31\u662f\u4e00\u4e2a\u975e\u7ebf\u6027\u7684\u5341\u5206\u590d\u6742\u7684\u65b9\u7a0b\uff0c\u5b83\u7684\u6570\u503c\u6c42\u89e3\u65b9\u6cd5\u5728\u7406\u8bba\u4e0a\u4e0d\u591f\u5b8c\u5584\uff0c\u6240\u4ee5\u9700\u8981\u901a\u8fc7\u5b9e\u9a8c\u6765\u52a0\u4ee5\u9a8c\u8bc1\u3002\u6b63\u662f\u5728\u8fd9\u4e2a\u610f\u4e49\u4e0a\u8bb2\uff0c\u6570\u503c\u6a21\u62df\u53c8\u53eb\u6570\u503c\u8bd5\u9a8c\u3002\u5e94\u8be5\u6307\u51fa\u8fd9\u90e8\u5206\u5de5\u4f5c\u51b3\u4e0d\u662f\u8f7b\u800c\u6613\u4e3e\u7684\u3002
\u6269\u5c55\u8d44\u6599\uff1a
\u6570\u503c\u6a21\u62df\u7684\u53d1\u5c55\u53f2\uff1a
1955\u5e74Peaceman\u4e0eRachford\u7814\u53d1\u7684\u4ea4\u66ff\u9690\u5f0f\u89e3\u6cd5(ADI)\u662f\u6570\u503c\u6a21\u62df\u6280\u672f\u7684\u91cd\u5927\u7a81\u7834\u3002\u8be5\u89e3\u6cd5\u975e\u5e38\u7a33\u5b9a,\u800c\u4e14\u901f\u5ea6\u5feb\uff0c\u6240\u4ee5\u8fc5\u901f\u5728\u5305\u62ec\u77f3\u6cb9,\u6838\u7269\u7406\uff0c\u70ed\u4f20\u5bfc\u7b49\u9886\u57df\u5f97\u5230\u5e7f\u6cdb\u5e94\u7528\u30021958\u5e74Douglas,Jim\u548cBlair,P.M\u7b2c\u4e00\u6b21\u8fdb\u884c\u4e86\u8003\u8651\u6bdb\u7ba1\u538b\u529b\u6548\u679c\u7684\u6c34\u9a71\u6a21\u62df\u3002
60\u5e74\u4ee3\u6570\u503c\u6a21\u62df\u6280\u672f\u7684\u53d1\u5c55\u4e3b\u8981\u5728\u6570\u503c\u89e3\u6cd5\uff0c\u7b2c\u4e00\u4e2a\u6709\u6548\u7684\u6570\u503c\u6a21\u62df\u89e3\u6cd5\u5668\u662f1968\u5e74Stone\u63a8\u51fa\u7684SIP(Strong Implicit Procedure)\u3002\u8be5\u89e3\u6cd5\u53ef\u4ee5\u5f88\u597d\u5730\u7528\u6765\u6a21\u62df\u975e\u5747\u8d28\u6cb9\u85cf\u548c\u5f62\u72b6\u4e0d\u89c4\u5219\u6cb9\u85cf\u3002
Stone\u572870\u5e74\u4ee3\u53d1\u8868\u4e86\u4e09\u76f8\u76f8\u5bf9\u6e17\u900f\u7387\u6a21\u578b\uff0c\u7531\u6cb9\u6c34\u548c\u6cb9\u6c14\u4e24\u76f8\u76f8\u5bf9\u6e17\u900f\u7387\u8ba1\u7b97\u6cb9\u3001\u6c14\u3001\u6c34\u4e09\u76f8\u6d41\u52a8\u65f6\u7684\u76f8\u5bf9\u6e17\u900f\u7387\uff0c\u8be5\u6280\u672f\u73b0\u5728\u8fd8\u5e7f\u4e3a\u5e94\u7528\u300270\u5e74\u4ee3\u53e6\u4e00\u9879\u4e3b\u8981\u6210\u5c31\u662fPeaceman\u63d0\u51fa\u7684\u4ece\u7f51\u683c\u538b\u529b\u6765\u786e\u5b9a\u4e95\u5e95\u6d41\u538b\u7684\u6821\u6b63\u65b9\u6cd5\u3002
\u53c2\u8003\u8d44\u6599\u6765\u6e90\uff1a\u767e\u5ea6\u767e\u79d1\u2014\u6570\u503c\u6a21\u62df
（一）研究现状
地下多相、多组分流体运移数值模拟是在质量和能量守恒的基础上，建立的多相流体运动以及反映地球化学运移扩散的数学模型，通过离散建立大量的线形或非线形方程组，然后利用计算机计算求解，再通过图像显示模拟结果，达到对工程问题和物理问题乃至相关其他问题研究的目的。CO2地质封存数值模拟就是利用计算机模拟的方法，来解决CO2进入地质封存系统后运移、转化、水-岩-气之间的相互反应、CO2泄漏对浅部含水层影响及诱发的储盖层物性变化等一系列问题，从而指导CO2地质封存工程的实施。
目前，国内外已开展的关于CO2地质封存数值模拟的研究工作包括以下几个方面：
1.超临界CO2-水多相流体运动模拟
Pruess等（2003）模拟了均质各向同性咸水含水层中以恒定流量灌注CO2条件下，灌注井井周非等温径向流情况。当忽略重力和惯性力效应时，模拟结果中存在相似变量ζ＝R2/t（其中，R为径向流动距离，t为时间），CO2饱和度、溶解CO2质量分数、沉淀盐的体积分数及流体压力都是相似变量的函数。这与O' Sullivan（1981）及Doughty和Pruess（1992）的结果一致。两相流的模拟考虑了CO2和水的相对渗透率及毛细管力作用问题（Van Genuchten，1980），考虑了流体密度、黏度和CO2溶解性随压力、温度和含盐度的变化，以及盐的沉淀导致含水层渗透率的减小等因素。
Doughty和Pruess（2004）利用Fro咸水含水层封存CO2监测资料，反推了CO2灌注后发生的物理和化学过程。他们采用TOUGH2数值模拟软件对两相（液、气）三组分（CO2、水和溶解NaCl）系统进行模拟。考虑15MPa和65℃条件下，超临界CO2在咸水中为非混溶流体，并能部分溶解于咸水的情况，分析了多相流系统边界设定的影响及相对渗透率取值问题，即模拟中对侧向边界的设置为均开（或均闭），结果导致压力的模拟结果与实际相比过低（或过高）。研究表明，由于上覆断层对CO2的封堵作用，侧向边界对CO2弥散羽的影响不大。模拟结果还显示，相对渗透率函数对CO2弥散羽的演化有很强的影响。如何确定一个合适的相对渗透率以表征CO2灌注咸水含水层的变化，仍是亟待解决的问题。Doughty和Pruess模拟了两种CO2残余饱和度条件下的CO2羽扩展随时间的变化，发现存在较大差异。残余饱和度较大的情况下，CO2羽表现出紧缩状，在浮力作用下运移较慢；相反，在残余饱和度较小的情况下，CO2羽流弥散很快，溶解性显著提高。
2.多组分反应地球化学运移模拟
水-砂岩-CO2相互作用往往形成一系列次生矿物，或次生矿物组合。Worden et al.（2006）通过岩石学和CO2灌注长石砂岩的地球化学模拟表明，北海Magnus油田上侏罗统浊积亚长石砂岩中的铁白云石、高岭石和石英可能具有成因联系。其中，铁白云石中的碳来自有机成因的CO2。Watson et al.（2004）通过CO2气与CH4气储集砂岩的比较岩石学研究，证实澳大利亚Otway盆地Ladbroke Grove CO2气储集砂岩中与CO2气灌注有关的次生矿物组合为铁白云石-高岭石-次生石英。
Xu et al.（2005）利用一维砂岩-页岩系统模型模拟了储层中灌注的CO2与矿物发生的化学反应过程，以及对储层环境的影响。模拟显示，在砂岩环境下，CO2主要被方解石所固定，而方解石的沉淀导致孔隙度减小，进而导致渗透率相应减小。10万年间，砂岩封存能力达到90kg/m3的封存能力，这些被矿物固定下来的CO2可以永久封存。Zwingmann等运用地球化学模拟软件EQ3/6进行的水-矿物-CO2相互作用模拟也表明，若将CO2灌注到日本本州岛中北部新潟盆地更新世灰爪组砂岩，CO2以溶于水和形成碳酸盐矿物两种形式封存，其中后者封存量最大为21.3mol/kgH2O，可达总封存量的90％，形成的碳酸盐矿物中也出现了片钠铝石。
3.耦合岩石力学模拟
从目前发表的论文及各国研究计划的综合报告上看，在CO2咸水含水层封存研究方面，对于CO2运移机制的分析和模拟很少考虑应力场的耦合作用。事实上，CO2灌注压力和超临界CO2的浮力作用将改变地层应力状态，即CO2在上浮运移和侧向扩散过程中，孔隙压力可能会对原始裂隙和断裂产生影响；CO2在咸水含水层中的长时期（千年级尺度以上）的封存，将改变含水层的地球化学状态，CO2-咸水-含水层矿物的化学作用将可能导致岩体力学和水力学性质发生变化。
日本因位于4大板块交界处与环太平洋构造带中，活断层密集发育，地震频繁，地应力分布复杂，在CO2地质封存评价方面，非常重视CO2地质封存的力学稳定性研究（李琦等，2002；李小春等，2003）。李琦等（2002；2004；2006）提出了一个考虑初始地应力场、灌注压力、CO2浮力及含热传导作用的热-水-力（THM）耦合模拟框架，考虑盖层底部附近存在不同倾角断层的二维平面应变地质封存问题。采用有限元算法，对灌注CO2流体对断层稳定性的影响进行模拟分析。计算结果表明，为了避免断层位移需要特别注意对灌注压力的控制，因为CO2灌注压力对断层滑动的影响远大于CO2羽流浮力带来的影响。停止灌注CO2后，CO2羽流的上升则成为应力场扰动的主要因素。
（二）主要软件介绍
近年来，计算机模拟技术在许多研究领域得到了广泛的应用，开发出了许多优秀的模拟软件和程序。同样，可用于研究CO2地质封存的数值模拟软件也很多，主要有PHREEQC、GEM、ECLIPSE、TOUGHREACT、PetroMod、MUFTE-UG和NUFT等，它们都有各自的特点和适用性。在进行数值模拟之前，需对这些数值模拟软件进行评价分析，选择适用于所研究问题的模拟软件。现对国际上常用的几款软件简介如下。
1.PHREEQC
PHREEQC是一款用于计算多种低温水文地球化学反应的计算机软件。以离子缔合水模型为基础，PHREEQC可完成以下任务：(1)计算物质形成种类与矿物的溶解饱和指数；(2)模拟地球化学反演过程；(3)计算批反应与一维运移反应。另外，与多组分溶质-运移模型耦合的PHREEQC可生成PHAST，一个用于模拟地下水流系统的三维反应-运移模拟器。但由于PHREEQC是在单相水流的基础上建立的模型，因而不能模拟超临界CO2-水的两相流运动。
PHREEQC最简单的应用就是计算溶液中各种化学物质的分布，以及溶液中矿物与气体的饱和状态。反演模拟功能可推导和量化在流动过程中，能够反应化学物质变化的化学反应方程。PHREEQC可处理的反应方程包括建立矿物、表面配合物、阳离子交换剂、土壤溶液、气体组分单位分压、给定压力或给定体积气相间平衡的物质运移反应。在模拟这些均衡反应的同时，PHREEQC还可以模拟动力化学与生物反应，以及模拟从简单的线性衰变（代谢物降解或放射性衰变）到复杂的依赖于溶液化学组成和微生物数量确定的反应速度。这些反应处理功能可在批反应模拟或一维对流、弥散、反应型运移模拟中使用。
2.GEM
GEM v.2009.13（Nghiem et al.，2004）是一款用来模拟利用CO2和酸性气体提高石油采收率的模拟器，该模拟器完全耦合了地球化学组成状态方程。GEM采用一步求解法进行状态方程的求解。GEM可以用来模拟：对流和弥散流体、油（或超临界CO2）、气和咸水间的平衡、水相物种间的化学平衡，以及矿物的动态溶解和沉淀。该模拟器采用自适应的隐式离散技术利用一维、二维或者三维模型来模拟孔隙介质中溶质的运移。油相和气相用一个状态方程来模拟，气体在水相的溶解度采用亨利定律模型来计算。水通过蒸发进入到气相、盖层的穿透、热效应和裂隙的封闭作用也可以利用GEM来模拟。
3.ECLIPSE
ECLIPSE是一个并行化的可以模拟黑油、组分、热采等问题的成熟软件。1994年，胜利石油管理局引进了ECLIPSE油藏数值模拟串行软件，广泛开展了从油藏到气藏，从常用油田到特殊油气田、从常规模拟研究到特殊模拟研究等多方面的应用。主要模块有主模型、黑油、组分、热采、流线法、运行平台和ECLIPSE Office等。
ECLIPSE是一个商业软件，在使用中其内核部分是封闭的，使用者只能将其作为一个“黑箱”来操作。其不足之处有：不可能免费的获得和随意地使用和修改；无法耦合最前沿的地质流体热力学模型；无法加入更多影响因素来研究具体问题。因此，ECLIPSE不适宜用于前沿科学研究。
4.TOUGH2/TOUGHREACT
TOUGH2是Transport of Unsaturated Groundwater and Heat（非饱和地下水流及热流传输）的英文缩写，是一个模拟一维、二维和三维孔隙或裂隙介质中，多相流、多组分及非等温的水流及热量运移的数值模拟程序。TOUGH2使用积分有限差分（Integral Finite Differences，IFD）（图3-8）的方法来解决多相流动和多组分化学运移模拟中的空间离散化问题（Pruess et al.，1999s；Xu et al.，2004）。为了满足大规模计算需要，Zhang et al.（2008）开发了TOUGH2的平行计算版本，即TOUGH2-MP。
该方法在对地质介质的离散化上是比较灵活的，允许使用不规则的网格，十分适合对多区域非均质系统和裂隙岩石系统中流体流动、运移和水岩相互作用的模拟。而对于规则的网格剖分，积分有限差分方法相当于传统的有限差分法。其中，对于任意区域Vn，它的质量（对于水、气体和其他化学组分）和能量（对于热）守恒方程可以用积分的方式（式3-5）表达：
图3-8 积分有限差分法中的空间离散化和几何参数数据构成图
中国二氧化碳地质封存选址指南研究
式中：下角标n为表示一个单元格；下角标m为表示和单元格n互相连接的网格m；Δt为时间步长；Mn为单元格n的平均质量或能量密度；Anm为单元网格n和m交界的面段；Fnm为通过面积为Anm的质量或能量通量；qn为单元格n内单位体积的平均源汇率。
许天福等（1998）在TOUGH2的框架基础之上，增加了多组分溶质运移和地球化学反应的模拟功能，形成了一套较为完善的可变饱和地质介质中非等温多相流体反应地球化学运移模拟软件——TOUGHREACT。该软件不仅包括了TOUGH2的全部功能，而且适用于不同温度、压力、水饱和度、离子强度、pH值和氧化还原电位（Eh）等水文地质和地球化学条件下的热-物理-化学过程。还可以应用于一维、二维或三维非均质（物理和化学的）孔隙或裂隙介质中的相关数值模拟研究。在理论上可以容纳任意数量的以固相、液相或气相存在的化学组分（但是在实际模拟中会受到计算能力和计算时间等硬件条件的限制），并且考虑了一系列化学平衡反应，如溶液中的配合反应、气体的溶解或脱溶、离子吸附作用、阳离子交换及受平衡控制或反应动力学控制的矿物溶解或沉淀反应等。可以说TOUGHREACT、是TOUGH2的升级版，近年来在世界范围内CO2地质封存研究和工程实践中得到了广泛的应用。
除包含TOUGH2所有的功能外，TOUGHREACT还可以应用于一系列的反应性流体和地球化学迁移问题。比如：(1)伴随Kd线性吸附和放射性衰变的污染物迁移问题；(2)在周围环境条件下，自然界中地下水的化学演变；(3)核废料处置地点评估；(4)深部岩层的沉积成岩作用；(5)CO2地质处置。多相流体运动，多组分反应地球化学，各种封存形式封存量以及随时间空间变化；(6)矿物沉积（如表生铜矿富集）；(7)自然和补给环境下热水系统中的矿物变化。
通过最近几年相关研究者的不懈努力，TOUGHREACT在实际应用中得到了进一步完善和提高，增加了部分新功能，如水相内部反应动力学和生物降解作用，改进了矿物-水反应表面积计算方法，以及气-水反应中气的活度系数的修正等。
5.PetroMod
由德国IES（Integrated Exploration System）公司研究开发的PetroMod多组分、多相态的多维含油气系统模拟软件综合平台已被世界石油业所公认。该软件融入了断层活动性、盐丘上涌和刺穿、火山岩的侵入、气体扩散效应、油气水三相运移和油气吸附模型等相关技术。
该模拟软件平台推出和采用的油气运移组合模拟算法（Hybird）是当今最先进的油气运移模拟算法，既可以保证模拟的精度，又可以极大地提高模拟的运算速度。其中的PetroFlow3D用于油气运移、聚集、圈闭和散失等情况的模拟，同时PetroCharge Express为我们提供了基于图件的油气运移和圈闭模拟的快速分析工具。
6.MUFTE-UG
MUFTE-UG是MUFTE和UG.MUFTE的结合。MUFTE即多相流（Muliphase Flow）、运移（Transport）和能量（Energy）模型。该软件包主要包括物理模型概念和孔隙裂隙介质中等温和非等温多相多组分流动和运移过程的离散方法（Helmig，1997；Helmig et al.，1998）。它能对裂隙孔隙介质进行离散性描述（Dietrich et al.，2005）。UG是非结构性网格（Unstructured Grid）的缩写，它提供的数据结构能快速解算以平行、自适应多网格法为基础的离散型偏微分方程。具有模块化结构的MUFTE-UG很容易解决各种有特殊要求的问题。
模块化结构的MUFTE-UG具有许多不同的环境与技术应用。例如，在环境应用领域，MUFTE-UG能够模拟如下两个问题。
（1）NAPL（非液相流体）向饱和与非饱和土壤的渗流。优化改进的修复技术在MUFTE中具有广泛的研究和发展空间。
（2）地下CO2的消散。CO2以高温高压灌注地表以下几百米的地层中，MUFTE-UG可用于非均质含水层（对流和弥散运移）中羽状体演化评价，伴随温度效应（由于膨胀和压缩）和组分间相互溶解（卤水和CO2）。
7.NUFT
NUFT（Nonisothermal Unsaturated-Saturated Flowand Transport model）是一套用来解决在多孔介质中多相、多组分非等温流动和溶质运移过程中地下污染物运移的数值解法器。此软件利用简单的代码来利用通用的实用程序和输入文件的格式。最近，此代码在Unix和DOS系统下运行成功。
该程序利用一套完整的有限差分空间离散法求解平衡方程组。每一个时间步长内利用Newton-Raphson方法求解非线性方程组，而在每一步迭代过程中利用直接解法和预共轭梯度法求解线性方程组。该模型可以解决一、二和三维水流及溶质运移问题。将来该模型会耦合进毛细滞后、非正交网格离散、有限单元剖分和固体非线性等温吸附等功能。
（三）研究方法
通常情况下，CO2地质封存数值模拟包括以下主要过程。
（1）建立概念模型：根据各种方法获取的实际资料来概化和建立CO2地质封存概念模型，包括边界范围、地层或储盖层高程、储盖层确定、参数及分区、源汇项、主要物理化学过程以及模型维度（一维、二维和三维）。
（2）建立数学模型：建立一套描述深部咸水层中多相流动和多组分反应性溶质运移的偏微分方程组，包括初始条件和边界条件问题。
（3）模型离散化：把概念模型中的各种信息通过网格剖分进行离散，形成大量的网格单元，然后通过有限差分、有限单元和积分有限差分等方法转化成单元的质量和能量守恒方程组，再用多种方法将非线性方程组线性化，形成线性代数方程组，然后求解方程组。
（4）模型识别和校正：根据模型计算结果和实际监测数据进行对比拟合，适度合理调整参数，使模型能够综合反映实际情况。在历史拟合过程中出现较大误差，应重新检查概念模型，修正概念模型。对所建模型进行参数敏感性分析，对于较敏感的参数应该慎重选取，甚至需要做大量的试验来确定。
（5）模型预测：建立了可靠的模型后，便可以进行模拟预测。
数值模拟的关键是地质模型概化、计算精度和计算速度。由于计算的精度取决于离散的程度，而离散的程度又决定了计算的速度，这是一对矛盾，要根据解决问题的需要来选择离散化的程度和计算速度。
CO2在储层中的运移、溶解以及与围岩的化学反应形成了一个多相、多组分的反应体系，涉及的主要数学方程有超临界CO2-水的两相流体运动控制方程、溶质运移控制方程和化学反应方程等。建立数值模型时，通常采用有限差分法、有限元法和积分有限差分法等。
由于实际应用时多采用已有的数值模拟软件对CO2地质封存的某一过程进行模拟，不涉及软件的开发及程序代码的编写，只需根据研究的需要选择合适的软件进行模拟预测，而软件一旦选定，数学模型和数值模型基本上已经确定。
绛旓細閲戝睘浣撶Н鎴愬舰宸ヨ壓涓鏁板兼ā鎷熸妧鏈姒傝 1. 宸ヨ壓涓庣壒鐐癸細閲戝睘浣撶Н鎴愬舰锛屽鑷敱閿汇佹ā閿诲拰鍥炶浆鎴愬舰绛夛紝浠ュ叾鐙壒鐨勫鎬у彉褰㈣繃绋嬪閫犻噾灞為浂浠躲傚伐鑹哄洜绱犲娓╁害銆佸帇鍔涘拰妯″叿璁捐瀵归浂浠惰川閲忔湁閲嶅ぇ褰卞搷锛屽父瑙佺殑缂洪櫡闇瑕佸瘑鍒囧叧娉ㄣ2. 鐞嗚鍩虹锛氭暟鍊兼ā鎷熶腑锛屽垰/绮樺鎬ф湁闄愬厓娉曟槸鏍稿績锛屽熀浜庡鎬у姏瀛﹀師鐞嗭紝鑰冭檻鏉愭枡鐨勫彉鍒...
绛旓細鎺ョ潃锛屾枃绔犳繁鍏ュ墫鏋愪簡涓ょ鍏抽敭鐨鏁板兼ā鎷熸妧鏈锛氬熀浜庣撼缁-鏂墭鍏嬫柉鏂圭▼鐨勬ā鎷燂紝浠ュ強鍩轰簬姘村姩鍔涙柟绋嬬殑妯℃嫙銆傚湪杩欎簺鏂规硶涓紝浣滆呴噸鐐规彮绀轰簡鍏朵腑鐨勫叧閿棶棰樺拰鎸戞垬锛屽绋冲畾鎬с佺簿搴﹀拰璁＄畻鏁堢巼鐨勪紭鍖栥備负浜嗙‘淇濇ā鎷熺粨鏋滅殑鍙俊搴︼紝鏂囩珷杩浠嬬粛浜嗘暟鍊兼ā鎷熺殑楠岃瘉绛栫暐鍜屾敼杩涙帾鏂斤紝杩欏浜庢彁楂樻ā鎷熺殑鍑嗙‘鎬ц嚦鍏抽噸瑕併傞氳繃杩欎簺...
绛旓細鍦ㄨ绠楁満涓婂疄鐜颁竴涓壒瀹氱殑璁＄畻锛岄潪甯哥被浼间簬灞ヨ涓涓墿鐞嗗疄楠屻傝繖鏃跺垎鏋愪汉鍛樺凡璺冲嚭浜嗘暟瀛︽柟绋 鐨勫湀瀛愭潵瀵瑰緟鐗╃悊鐜拌薄鐨勫彂鐢燂紝灏卞儚鍋氫竴娆＄墿鐞嗗疄楠屻鏁板兼ā鎷瀹為檯涓婂簲璇ョ悊瑙ｄ负鐢ㄨ绠楁満鏉ュ仛瀹為獙銆傛瘮濡傛煇涓鐗瑰畾鏈虹考鐨勭粫娴侊紝閫氳繃璁＄畻骞跺皢鍏惰绠楃粨鏋滃湪鑽у厜灞忎笂鏄剧ず锛屽彲浠ョ湅鍒版祦鍦虹殑鍚勭缁嗚妭锛氬婵娉㈡槸鍚﹀瓨鍦紝瀹冪殑浣嶇疆...
绛旓細姘斿鏁板兼ā鎷熸妧鏈鏄竴绉嶅己澶х殑宸ュ叿锛屽畠閫氳繃鏋勫缓绮剧粏鐨勬皵鍊欐ā鍨嬫潵鐞嗚В鍜岄娴嬫湭鏉ョ殑姘斿欑姸鍐点傝繖涓ā鍨嬪熀浜庡褰撳墠鍜屽巻鍙叉皵鍊欐暟鎹殑娣卞叆鍒嗘瀽锛岀壒鍒叧娉ㄦ俯瀹ゆ皵浣撴祿搴﹀鍔犲彲鑳藉甫鏉ョ殑姘斿欏彉鍖栥傚湪杩囧幓鐨勫崄骞撮噷锛屾皵鍊欐暟鍊兼ā鎷熸妧鏈殑绮剧‘鎬у拰鍙潬鎬у緱鍒颁簡鏄捐憲鎻愬崌锛岀壒鍒槸鍦ㄨ繃鍘荤殑浜斿勾涓紝闅忕潃璁＄畻鏈烘妧鏈殑椋為熷彂灞曪紝娴锋皵...
绛旓細鎺ョ潃锛岀浜岀珷娣卞叆鎺㈣闀挎睙娌冲彛鍦板尯鐨勫湡澹ゆ按鐩愯繍绉绘暟鍊兼ā鎷熷強鍏跺姩鎬侀娴嬨備綔鑰呯壒鍒叧娉ㄤ笁宄″伐绋嬪涓嬫父娌冲彛鍦熷￥姘寸洂鍔ㄦ佺殑褰卞搷锛岄氳繃瀹炰緥鍒嗘瀽锛岄槓杩颁簡濡備綍杩愮敤鏁板兼ā鎷熸妧鏈鏉ョ悊瑙ｅ拰棰勬祴杩欎竴鐜拌薄銆傝繖閮ㄥ垎鍐呭瀵圭幆澧冨鍜屾按鏂囧涓撲笟浜哄＋鍏锋湁閲嶈鍙傝冧环鍊硷紝鏈夊姪浜庝粬浠悊瑙ｅ鏉傜殑鐜鍙樺寲鍔ㄦ併
绛旓細1銆 澶嶆潅绯荤粺鐨勮璁★紝渚嬪椋炴満鍜岀墿娴佺郴缁熴2銆 鐢靛姏绯荤粺浠跨湡銆3銆 宸ヤ笟杩囩▼鐨勮璁★紝渚嬪鍖栧鍔犲伐鍘傘4銆 杩囩▼宸ョ▼浠跨湡宸ュ叿銆5銆 鐢ㄤ簬璁捐鏈哄櫒浜哄拰鏈哄櫒浜烘帶鍒剁畻娉曠殑鏈哄櫒浜妯℃嫙鍣ㄣ6銆 妯℃嫙姹借溅纰版挒浠ユ祴璇曟柊杞﹀瀷鐨勫畨鍏ㄦ満鍒躲備汉浠璁＄畻鏈轰豢鐪熺殑鍙潬鎬у拰淇′换鍙栧喅浜庝豢鐪熸ā鍨嬬殑鏈夋晥鎬э紝鍥犳楠岃瘉鍜岀‘璁ゅ浜庤绠楁満...
绛旓細闅忕潃娌圭敯鎶曚骇鍚庨潤鎬佸強鍔ㄦ佽祫鏂欏鍔,杩橀渶瑕佷慨鏀瑰師鏈夌殑鍦拌川妯″瀷,閫氳繃鍏ㄤ綋娌硅棌妯″瀷鏁板兼ā鎷熺爺绌跺姞娣卞鍦拌川妯″瀷鐨勬柊璁よ瘑,骞跺湪娌圭敯鐢熶骇鍘嗗彶鎷熷悎鍩虹涓婅繘琛岀敓浜ч娴嬨 鍥犳,娌硅棌鏁板兼ā鎷熸妧鏈鏄补钘忓伐绋嬬爺绌躲佹补鐢板姩鎬佸垎鏋愪腑鐨勪竴椤瑰崄鍒嗛噸瑕佺殑鎵嬫銆 涓浗娴锋补鐨勬补钘忔暟鍊兼ā鎷熺爺绌惰捣姝ヤ簬20涓栫邯80骞翠唬鍒濄備负浜嗗敖蹇缉鐭繖椤规妧鏈笌鍥介檯...
绛旓細鍦ㄥ昂瀵告柟闈紝涔︾睄鐨勫昂瀵告槸25鍘樼背瀹姐17.5鍘樼背楂橈紝鍘氬害浠呬负2.5鍘樼背锛屾暣浣撻噸閲忕害涓898鍏嬶紝渚夸簬鎼哄甫鍜屽瓨鏀俱傛澶栵紝瀹冭繕鏈変竴涓狝SIN缂栧彿锛屽嵆B0011CLQVS锛岃繖鏄簹椹婂湪绾垮钩鍙扮殑鍞竴鏍囪瘑绗︺傝繖鏈功娑电洊浜嗗喍閲戣繃绋嬩腑鐨勫叧閿鏁板兼ā鎷熸妧鏈锛屼粠鐞嗚鍩虹鍒板疄闄呭簲鐢紝娣卞叆娴呭嚭鍦版帰璁ㄤ簡鍏跺湪宸ヤ笟鐢熶骇涓殑浣滅敤銆傛棤璁烘槸...
绛旓細SWANTEC 鍏徃绠浠SWANTEC, 绠绉"绉戝鐒婃帴鍜鏁板兼ā鎷熸妧鏈", 璇炵敓浜1999骞6鏈1鏃ワ紝婧愪簬涓归害鎶鏈ぇ瀛(DTU)鍦ㄧ數闃荤剨鎺ュ拰閲戝睘鐮旂┒棰嗗煙鐨勬繁鍘氱Н绱傝繖瀹剁钀ユ湁闄愯矗浠诲叕鍙告簮浜庡ぇ瀛︾殑绉戠爺鎴愭灉锛屼笓娉ㄤ簬鐢甸樆鐒婃帴鍙婄浉鍏冲伐鑹虹殑绉戞妧鍒涙柊銆傚叕鍙哥殑鏍稿績涓氬姟鍦ㄤ簬寮鍙戝苟鍟嗕笟鍖栨ā鎷熶紭鍖栬蒋浠讹紝涓撲负鐢甸樆鐒婃帴杩囩▼鎻愪緵绮惧噯鐨勬ā鎷熷拰宸ヨ壓...
绛旓細瀵规湁鏉′欢杩涜瀹為獙鐨勬潗鏂欙紝灏介噺閲囩敤瀹為獙鏂规硶锛岃緟浠鏁板兼ā鎷妫楠屻傝屽湪宸 绋嬪簲鐢ㄤ腑锛屽緢澶氭儏鍐典笅鏃犳硶杩涜瀹為獙锛屽閲囩熆闂绛夛紝鏁板兼ā鎷熷唴閮ㄧ▼搴忔湁鐩稿簲鐨勮绠楁柟娉曪紝鑳芥ā鎷熻緝澶嶆潅杩囩▼銆傜洿瑙傛т笌姹傝В閫熷害锛氬疄楠岀洿瑙傛у己锛屾暟鍊兼ā鎷熺洿瑙傛т笉濡傚疄楠屾柟娉曞ソ锛岃緝鎶借薄锛屼絾鍙互 蹇熷緱鍒扮粨鏋溿傚疄楠屾搷浣滃鏉傘傛垚鏈細瀹為獙鎴愭湰...
