发明名称 一种油气藏水锁伤害半径的确定方法及系统
摘要 本申请提供一种油气藏水锁伤害半径的确定方法,该确定方法包括:根据毛管力P<sub>c</sub>以及井筒内液相压力与地层压力的压差ΔP对水相的驱替作用,确定水相运动速度;根据所述水相运动速度和水相与储层的接触截面积,确定水相进入储层的水量;根据所述水相进入储层的水量、储层的孔道截面积,孔隙度计算出水相侵入储层的半径。本申请还提供一种油气藏水锁伤害半径的确定系统。本申请实施例的确定方法不需要进行岩心实验和试井的方法,即能够准确预测气井水锁侵入的半径,而且还能计算各种作业时的动态侵入半径,该系统操作简单、计算误差小。能够准确描述周期长的作业措施的水锁伤害半径,弥补了常规实验方法不能得出水锁伤害半径的缺陷。
申请公布号 CN103266870B 申请公布日期 2016.03.02
申请号 CN201310206671.9 申请日期 2013.05.29
申请人 中国海洋石油总公司;中海油田服务股份有限公司 发明人 赵利昌;李敬松;黄子俊;刘子雄;王杏尊;吴英;曾鸣
分类号 E21B43/20(2006.01)I;E21B49/00(2006.01)I 主分类号 E21B43/20(2006.01)I
代理机构 北京安信方达知识产权代理有限公司 11262 代理人 苏蕾;武晶晶
主权项 一种油气藏水锁伤害半径的确定方法,其特征在于,该确定方法包括:根据岩心的毛管力实验曲线,确定岩心实验条件下的毛管力P<sub>c1</sub>,再根据如下表达式计算出储层条件下的毛管力P<sub>c</sub>:P<sub>c</sub>=1/15P<sub>c1</sub>,P<sub>c1</sub>为岩心实验条件下的含气或含水饱和度对应的毛管力;根据如下表达式确定压差△P:△P=P<sub>wf</sub>‑P<sub>i</sub>,P<sub>wf</sub>为井底压力,P<sub>i</sub>为地层静压力;根据如下表达式确定水相运动速度:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>v</mi><mi>w</mi></msub><mo>=</mo><mfrac><msub><mi>k</mi><mi>w</mi></msub><msub><mi>&mu;</mi><mi>w</mi></msub></mfrac><mo>&CenterDot;</mo><mfrac><mrow><mi>&Delta;</mi><mi>p</mi><mo>+</mo><msub><mi>p</mi><mi>c</mi></msub></mrow><mi>x</mi></mfrac></mrow>]]></math><img file="FDA0000862624990000011.GIF" wi="432" he="159" /></maths>其中,v<sub>w</sub>为水相运动速度;k<sub>w</sub>为水相渗透率;μ<sub>w</sub>为水相粘度;x为运动的距离;根据如下表达式确定含水饱和度S<sub>wf</sub>:当S<sub>wi</sub>≥Sor时,S<sub>wf</sub>为P<sub>c2</sub>‑△P相对应的含水饱和度;当S<sub>wi</sub>&lt;Sor时,S<sub>wf</sub>为P<sub>c2</sub>‑△P‑P<sub>c</sub>相对应的含水饱和度;P<sub>c2</sub>为与原始岩心中含水饱和度S<sub>wi</sub>相对应的毛管力,Sor为束缚水饱和度;根据如下表达式确定进入储层的水量:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>N</mi><mrow><mi>w</mi><mi>t</mi></mrow><mn>2</mn></msubsup><mo>=</mo><mi>A</mi><mo>&CenterDot;</mo><msub><mi>V</mi><mi>w</mi></msub><mo>&CenterDot;</mo><mi>t</mi><mo>&CenterDot;</mo><mi>A</mi><mo>&CenterDot;</mo><mi>x</mi><mo>&CenterDot;</mo><mi>&phi;</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>S</mi><mrow><mi>w</mi><mi>f</mi></mrow></msub><mo>-</mo><msub><mi>S</mi><mrow><mi>w</mi><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><msup><mi>A</mi><mn>2</mn></msup><mfrac><mrow><msub><mi>k</mi><mi>w</mi></msub><mi>&phi;</mi><mrow><mo>(</mo><msub><mi>S</mi><mrow><mi>w</mi><mi>f</mi></mrow></msub><mo>-</mo><msub><mi>S</mi><mrow><mi>w</mi><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>&Delta;</mi><mi>p</mi><mo>+</mo><msub><mi>p</mi><mi>c</mi></msub><mo>)</mo></mrow></mrow><msub><mi>&mu;</mi><mi>w</mi></msub></mfrac><mi>t</mi></mrow>]]></math><img file="FDA0000862624990000012.GIF" wi="1316" he="143" /></maths>其中,N<sub>wt</sub>为累积进入储层中的水量;A为孔道与水相的接触面积;k<sub>w</sub>为水相渗透率;φ为岩心孔隙度;S<sub>wf</sub>为水相进入地层后的含水饱和度;S<sub>wi</sub>为原始岩心中含水饱和度;μ<sub>w</sub>为水相粘度;x为运动的距离;t为作业过程中液相与储层的接触时间;根据如下表达式计算水相侵入储层的半径:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>L</mi><mi>c</mi></msub><mo>=</mo><mfrac><msub><mi>N</mi><mrow><mi>w</mi><mi>t</mi></mrow></msub><mrow><mi>A</mi><mi>&phi;</mi></mrow></mfrac><mo>=</mo><msqrt><mfrac><mrow><msub><mi>k</mi><mi>w</mi></msub><mi>t</mi><mrow><mo>(</mo><msub><mi>S</mi><mrow><mi>w</mi><mi>f</mi></mrow></msub><mo>-</mo><msub><mi>S</mi><mrow><mi>w</mi><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>&Delta;</mi><mi>p</mi><mo>+</mo><msub><mi>p</mi><mi>c</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>&mu;</mi><mi>w</mi></msub><mi>&phi;</mi></mrow></mfrac></msqrt></mrow>]]></math><img file="FDA0000862624990000021.GIF" wi="772" he="159" /></maths>其中,N<sub>wt</sub>为累积进入储层中的水量;A为孔道与水相的接触面积;k<sub>w</sub>为水相渗透率;φ为储层的岩心孔隙度;S<sub>wf</sub>为水相进入地层后的含水饱和度;S<sub>wi</sub>为原始岩心中含水饱和度;μ<sub>w</sub>为水相粘度;t为作业过程中液相与储层的接触时间。
地址 100010 北京市东城区朝阳门北大街25号