发明名称 基于阵列电成像测井的储层油水识别方法及识别系统
摘要 本发明涉及一种基于阵列电成像测井的储层油水识别方法及识别系统,针对泥浆侵入地层电阻率分布特点,建立“五参数地层模型”,适应了模拟不同泥浆类型、不同侵入地层电阻率分布特点要求,在阵列电成像测井正反演模型和处理方法方面具有开创性;基于阵列电成像测井“五参数反演”方法,进行泥浆侵入地层电阻率分布剖面重构,反映了泥浆侵入油层和水层电阻率分布特征或变化规律差异的实际不同,该方法具有较强的实用性,在储层流体性质定性识别和定量计算的实际应用方面具有开创性。
申请公布号 CN102678102B 申请公布日期 2015.02.25
申请号 CN201210090672.7 申请日期 2012.03.31
申请人 中国石油大学(华东) 发明人 邓少贵;范宜仁;陈华;李虎;李智强
分类号 E21B47/002(2012.01)I;E21B47/00(2012.01)I;E21B49/00(2006.01)I 主分类号 E21B47/002(2012.01)I
代理机构 青岛联信知识产权代理事务所 37227 代理人 王月玲;王中云
主权项 一种基于阵列电成像测井的储层油水识别方法,其特征在于:包括如下步骤:(1)处理井段多地层精细划分选取,确定径向深度r;(2)对地层参数初值预设:冲洗带半径r<sub>xo</sub>、过渡带半径r<sub>i</sub>、冲洗带电阻率R<sub>xo</sub>、过渡带中部地层电阻率R<sub>i0.5</sub>以及原状地层电阻率R<sub>t</sub>;形成参数向量形式:x=(x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>,x<sub>4</sub>,x<sub>5</sub>)<sup>T</sup>   (1)(3)建立三段式地层电阻率分布模型:冲洗带:R<sub>xo</sub>,r&lt;r<sub>xo</sub>过渡带:R<sub>i</sub>=ar<sup>2</sup>+br+c,r<sub>xo</sub>&lt;r&lt;r<sub>i</sub>原状地层:R<sub>t</sub>,r&gt;r<sub>i</sub>其中:过渡带电阻率系数确定如下式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='(' close=')'><mtable><mtr><mtd><msubsup><mi>r</mi><mi>xo</mi><mn>2</mn></msubsup></mtd><mtd><msub><mi>r</mi><mi>xo</mi></msub></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><msubsup><mi>r</mi><mrow><mi>i</mi><mn>0.5</mn></mrow><mn>2</mn></msubsup></mtd><mtd><msub><mi>r</mi><mrow><mi>i</mi><mn>0.5</mn></mrow></msub></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><msup><msub><mi>r</mi><mi>i</mi></msub><mn>2</mn></msup></mtd><mtd><msub><mi>r</mi><mi>i</mi></msub></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mfenced open='(' close=')'><mtable><mtr><mtd><mi>a</mi></mtd></mtr><mtr><mtd><mi>b</mi></mtd></mtr><mtr><mtd><mi>c</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>R</mi><mi>xo</mi></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mrow><mi>i</mi><mn>0.5</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mi>t</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000633444870000011.GIF" wi="574" he="244" /></maths>(4)根据地层电阻率分布,计算阵列电阻率测井响应f<sub>i</sub>=f<sub>i</sub>(x),i=1,......,5   (2)(5)确定5×5的Frechet导数阵<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mfrac><msub><mrow><mo>&PartialD;</mo><mi>f</mi></mrow><mn>1</mn></msub><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mn>1</mn></msub></mfrac></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mfrac><msub><mrow><mo>&PartialD;</mo><mi>f</mi></mrow><mn>1</mn></msub><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mn>5</mn></msub></mfrac></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><mfrac><msub><mrow><mo>&PartialD;</mo><mi>f</mi></mrow><mn>5</mn></msub><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mn>1</mn></msub></mfrac></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mfrac><msub><mrow><mo>&PartialD;</mo><mi>f</mi></mrow><mn>5</mn></msub><msub><mrow><mo>&PartialD;</mo><mi>x</mi></mrow><mn>5</mn></msub></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000633444870000012.GIF" wi="1259" he="375" /></maths>(6)确定实际测井值y<sub>i</sub>与地层参数向量(x)计算的测井值f(x)之差构成的差值向量B=(y<sub>1</sub>‑f<sub>1</sub>,y<sub>2</sub>‑f<sub>2</sub>,y<sub>3</sub>‑f<sub>3</sub>,y<sub>4</sub>‑f<sub>4</sub>,y<sub>5</sub>‑f<sub>5</sub>)<sup>T</sup>   (3)(7)求解线性方程组AΔx=B   (3)(8)修正地层参数x=x+Δx   (4)(9)确定测井值与计算值平方差<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&Phi;</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>5</mn></munderover><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>f</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000633444870000013.GIF" wi="340" he="142" /></maths>将(4)式回代,重复步骤(3)~(9),直到目标函数Φ为最小为止,这时计算的电阻率和测井的电阻率基本一致,得到真实的地层参数,从而利用真实的地层参数对三段式泥浆侵入剖面重构,进行油水识别。
地址 266555 山东省青岛经济技术开发区长江西路66号