发明名称 光透法定量孔隙介质内两相流中流体饱和度的方法
摘要 本发明公开了一种光透法定量孔隙介质内两相流中流体饱和度的方法,即光强‑饱和度(LIS)模型。该方法首先利用光透系统中CCD监测确定孔隙介质内饱水条件时的光强值;然后利用相同的方法实测目标时刻孔隙介质内两相流时的光强值;再利用LIS模型定量求解孔隙介质内各相流体的饱和度;最后对各像素位置单一相流体体积进行累加求取单一相流体的体积总量。通过两组密封砂箱室内试验对LIS模型进行了验证,发现2个NAPL/水系统和2个水/气系统的LIS模型都得到与实测资料相吻合的结果。该方法最大的优点是利用光透系统获得的实验数据获得不可直接测量的饱和度,为定量研究孔隙介质内流体运移规律奠定了基础。
申请公布号 CN103913429B 申请公布日期 2016.11.23
申请号 CN201410108700.2 申请日期 2014.03.21
申请人 南京大学 发明人 叶淑君;吴吉春;章艳红
分类号 G01N21/31(2006.01)I 主分类号 G01N21/31(2006.01)I
代理机构 江苏圣典律师事务所 32237 代理人 贺翔
主权项 一种光透法定量孔隙介质内两相流中流体饱和度的方法,其特征在于,包括以下步骤:(1)首先,利用光透系统中CCD监测确定孔隙介质内饱水条件时的光强值I<sub>w</sub>;(2)再利用相同的方法实测目标时刻孔隙介质内两相流的光强值I;(3)然后再利用四个LIS模型,即光强‑饱和度,定量求解孔隙介质内流体的饱和度S<sub>g</sub>(S<sub>o</sub>);水/气两相系统中的模型WG‑A和模型WG‑B,具体表达式如下:S<sub>g</sub>=ln(I/I<sub>w</sub>)/ln(C<sub>1</sub>)S<sub>g</sub>=(1‑I/I<sub>w</sub>)/(1‑C<sub>1</sub>)C<sub>1</sub>=I<sub>res</sub>/I<sub>w</sub>其中,S<sub>g</sub>为气相饱和度,I<sub>res</sub>是孔隙介质内只有残余水时的光强值,I<sub>w</sub>是孔隙介质内饱水条件时的光强值;NAPL/水两相系统中的模型NW‑A和模型NW‑B以C<sub>2</sub>的形式分别可以表达为:S<sub>o</sub>=ln(I/I<sub>w</sub>)/ln(C<sub>2</sub>)S<sub>o</sub>=(1‑I/I<sub>w</sub>)/(1‑C<sub>2</sub>)C<sub>2</sub>=I<sub>o</sub>/I<sub>w</sub>其中,S<sub>o</sub>为NAPL相饱和度,I<sub>o</sub>是孔隙介质内NAPL完全饱和时的光强值;以上饱和度均是指有效饱和度,其中水相的绝对饱和度S<sub>w‑abs</sub>计算如下所示:S<sub>w‑abs</sub>=S<sub>w</sub>(1‑S<sub>w‑res</sub>)+S<sub>w‑res</sub>其中S<sub>w</sub>为水相饱和度,S<sub>w‑res</sub>是孔隙介质的残余水饱和度;其它相X的饱和度,则表达为:S<sub>X‑abs</sub>=S<sub>X</sub>(1‑S<sub>w‑res</sub>)其中,S<sub>X</sub>是NAPL相或气相的饱和度,S<sub>X‑abs</sub>是NAPL相或气相的绝对饱和度;(4)最后对各像素位置单一相流体体积进行累加求取单一相流体的体积总量<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>V</mi><mi>X</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mn>2</mn></msub></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mn>1</mn></msub></munderover><msub><mi>S</mi><mrow><mi>X</mi><mo>-</mo><mi>a</mi><mi>b</mi><mi>s</mi></mrow></msub><mo>&times;</mo><mi>A</mi><mo>&times;</mo><mi>T</mi><mo>&times;</mo><mi>&theta;</mi></mrow>]]></math><img file="FDA0001014645660000011.GIF" wi="556" he="131" /></maths>其中,i,j是像素单元横向,纵向坐标,n<sub>1</sub>,n<sub>2</sub>分别为横向,纵向像素单元的数量,A是像素单元的面积,T是砂箱的厚度,θ是孔隙度。
地址 210000 江苏省南京市栖霞区仙林大道163号南大环境学院