数值分析课设——测井解释POR程序油气校正迭代研究

1 问题的定义

多数情况下是得不到一般数学方法所需的函数表达式，或难以找到原函数，都可以利用数值方法来求解，在计算机中实现的数值方法也称为数值算法。牛顿迭代法是数值分析中一个重要的计算方法和思想。迭代法的主要功能：计算方程

时可以比较快速。牛顿迭代法的一个迭代关系式：

Atlas公司的POR程序虽然简介使用，但它采用的是含水泥质岩石模型，未考虑地层孔隙中油气的影响，因此在解释评价含油气地层，尤其是含气地层时有很大的局限性，可能使得解释结果与实际地层情况相差较大，因为气的密度和油水差别很大，而在未知流体性质的基础上一般参数卡中默认的流体密度通常是1，即为地层水的密度。

利用SAND2程序的油气迭代方法，编写了一个油气校正子例程，并把该油气校正子例程加到已修改的POR程序中，对POR程序中已做过泥质校正的CNL，DEN曲线再作油气影响校正，再用校正后的CNL，DEN曲线精细解释有关地质参数和油藏物性参数。

2 迭代方法研究

POR程序增加油气校正的原理

油气校正基本思路及步骤：

（1）第一次调用油气校正子例程时，要对油气密度赋初值； （2）由POR程序中的方法,计算SW；

（3）当SW＞0.7时，无需迭代（说明地层无油气或油气较少）； （4）当SW≤0.7时，作油气校正：

重新计算冲洗带残余油气饱和度(SRH)及冲洗带残余油气相对体积（GAS） SRH=SRHM*(1-SW) GAS=SRH*POR

（5）计算油和气的密度及油气密度DHY

DHO0.7(1x)SRH(1x)(10.72x)SRHDHG(10.72x)SRH(1x)(2.20.8x)SRH(maVsh(mash)b)NDxD(ma(maf)f)DHYmaxDHO，DHG

（6）计算油、气层的校正系数A、B值

NNVshNsh油气校正方法是采用CORIBAND程序中的校正公式。其校正系数采用以下式子计算： 对气层：

A1.07[(1.110.15PPMF)*DF1.24DHY]BDF(1PPMF)2.2DHYDF(1PPMF)

对油层：

A1.07[(1.110.15PPMF)*DF1.11DHY0.03]BDF(1PPMF)DHY0.3DF(1PPMF)

式中：DHY ――油气密度

PPMF――泥浆滤液含盐量，ppm*10-6 DF ――流体密度 （7）对中子、密度进行油气校正

DENDENA*GAS CNLCNLB*GAS*100*EN （8）油气校正中迭代控制原理

a、如果 SW>=0.7 则不作油气校正，DHY=1，水层处 若 SW<0.7 作油气校正 b、如果 DHY>=0.4则认为是油层

若 DHY<0.4 则认为是气层

c、若 |GAS-GASO|<0.01 则油气校正完毕，退出迭代；（终止迭代） 若 |GAS-GASO|>=0.01 则需继续迭代，重新计算ρh，重复上述各步迭代运算，进行油气校正。

式中：GAS―当前的冲洗带残余油气相对体积 GASO－上一次迭代的冲洗带残余油气相对体积 d、迭代次数ND控制 当ND>= 5则退出迭代 当ND< 5 则需继续迭代校正

注：迭代次数是主程序调子例程次数，即迭代不是子例程完成。

（9）油气校正完后，需计算以下参数

迭代次数ND=－1 （赋小于零的值）

上述校正思路采用了SAND2油气校正思路，基本公式采用了CORIBAND油气校正公式。

3 实际资料的验证

所选测井资料概况

处理的井：LAN34井

测井曲线：RT，RLLS，RXO，SP，GR，CAL，CNL，DEN

处理井段：1700～1900米，第三纪砂泥岩剖面，淡水泥浆，CNL对砂岩刻度，测量井底深度为1930米，地层温度为820C，地温梯度为3.50C/100米。

图-1 LAN34井1700~1750m校正前的解释成果图

图-2 LAN34井1700~1750m校正后的解释成果图

由图-1和图-2处理的结果对比可以看出，经过油气的影响校正，残余油的饱和度明显减少，即可动油的体积相对增大，这对于估算油气资源的地质储量意义重大。

4 结论

通过本次课程设计，了解了牛顿迭代法的实际运用效果，其算法简单，能快速求出其他商法求不出或者难以求出的解。并且与我的专业联系起来达到了较好的效果，更加深层次的意识到数值计算的重要性，对以后工作学习起到了的巨大帮助。