主权项 |
一种基于阵列电成像测井的储层油水识别方法,其特征在于:包括如下步骤:(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<r<sub>xo</sub>过渡带:R<sub>i</sub>=ar<sup>2</sup>+br+c,r<sub>xo</sub><r<r<sub>i</sub>原状地层:R<sub>t</sub>,r>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>∂</mo><mi>f</mi></mrow><mn>1</mn></msub><msub><mrow><mo>∂</mo><mi>x</mi></mrow><mn>1</mn></msub></mfrac></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mfrac><msub><mrow><mo>∂</mo><mi>f</mi></mrow><mn>1</mn></msub><msub><mrow><mo>∂</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>∂</mo><mi>f</mi></mrow><mn>5</mn></msub><msub><mrow><mo>∂</mo><mi>x</mi></mrow><mn>1</mn></msub></mfrac></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mfrac><msub><mrow><mo>∂</mo><mi>f</mi></mrow><mn>5</mn></msub><msub><mrow><mo>∂</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>Φ</mi><mo>=</mo><munderover><mi>Σ</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),直到目标函数Φ为最小为止,这时计算的电阻率和测井的电阻率基本一致,得到真实的地层参数,从而利用真实的地层参数对三段式泥浆侵入剖面重构,进行油水识别。 |