高压物性参数预测方法研究 - 油气田地面工程.doc

Prediction Method of High Pressure Physical property 高压物性参数预测方法 李瑞琪 （大庆油田勘探开发研究院） 摘要 地层原油高压物性参数是确定油藏类型、制定开发方案和进行油藏工程计算不可缺少的重要参数，是研究油田驱动类型、确定油田开采方式、计算油田储量、选择油井工作制度的基础。 关键词 溶解气油比 饱和压力 压缩系数 体积系数 粘度 地层原油高压物性参数是确定油藏类型、制定开发方案和进行油藏工程计算不可缺少的重要参数，是研究油田驱动类型、确定油田开采方式、计算油田储量、选择油井工作制度的基础。目前取得高压物性参数的方法有3种： （1）对于饱和压力低于原始地层压力的未饱和油藏，试油、试采过程中，在井底流压高于饱和压力的条件下，通过井底取样器，取得有代表性的地层油样，然后在实验室内测定； （2）对于饱和压力等于或接近于原始地层压力的饱和油藏，难以通过井底取样获得有代表性的油样，一般采用地面取样进行重组复配的方法。在地面同时取得油样和气样，然后模拟地层条件的PVT组合配样，确定地层原油高压物性； （3）当不具备取样和PVT分析条件，或因地层含水较高而取不到有代表性的样品时，可以借用一些国内外已发表的经验公式预测高压物性。 经验公式的优点是实用方便。一般只要将所得的经验公式用在与原始实验数据相类似的情况下，都可以得到相当满意的结果。但如果不加选择地将经验关系式用于其它条件下，将产生不可靠的结果。例如当我们应用下式预测B10井溶解气油比时，预测值为79.1m3/t，实测值为45.4m3/t，误差非常大，必须对经验公式进行改进。 （1） 改进方法主要有两种： （1）根据理论、专业经验或实验数据的散点图确定经验公式的形式，然后根据实测数据确定公式中的系数； （2）仍采用目前已有公式的形式，再根据实测数据确定公式中的系数。 第一种方法实现起来比较困难，这里采用第二种方法，即公式形式不变，例如对上面的溶解气油比，其形式为： （2） 然后根据实测数据确定公式中的4个系数。确定公式系数时使用了MATLAB数学软件中的非线性回归函数nlinfit（）。 1、非线性回归函数nlinfit（） 在MATLAB数学软件中，非线性回归可由nlinfit（）来实现，调用格式为 [betar，j] = nlinfit（x，y，model’，beta0 其中x为自变量数据矩阵，维数为n×m；y为因变量数据向量n ；beta0为模型系数初值，model是非线性函数beta是估计出的回归系数，r是残差，j是矩阵，它们是估计预测误差需要的数据。 中的4个系数赋初值，最后调用nlinfit（）函数。 x = [1.0763 0.8369 72.8 18.56; …………………………………… 0.6809 0.8401 52.5 14.14]; y = [45.4 ……………………………………13.0]; beta0 = [1.1661 1.187 6.3967 0.036595]; [beta,r,j] = nlinfit(x, y, Rs, beta0) （3）执行Rs_HLE.m文件，返回能够满足精度要求的回归系数。 beta = 0.0025 1.8920 14.4156 0.0002 2、公式系数的确定 这里主要研究溶解气油比、饱和压力、压缩系数、体积系数和粘度5个高压物性参数经验公式系数的确定。这5个高压物性参数经验公式来源于国内外文献，形式见表1。 表1 高压物性参数经验公式形式 参数 公式 溶解气油比 饱和压力 压缩系数 体积系数 粘度 利用实测数据采用非线性回归 （3） 其中 （2）平均相对误差 （4） 其中 从表4可以看出对于原系数，5个高压物性参数的平均误差和平均相对误差均较大，平均相对误差为40.8%，必须对公式中的系数进行修正，从而显示出本文工作的意义。与原系数计算结果相比，新系数计算结果的误差则要小得多，平均相对误差为11.7%，应用新系数预测高压物性参数，能够得到比较可靠的结果。 表2 高压物性参数预测公式系数 表3 15口井高压物性参数计算结果对比 表3 15口井高压物性参数计算结果对比（续） 表4 15口井高压物性参数平均误差及平均相对误差 符号说明 Bo — 体积系数，m3/ m3； Ct — 综合压缩系数，1/MPa； Cf — 岩石孔隙体积压缩系数，1 /MPa； Co — 地层油压缩系数，1 /MPa； Cw — 地层水压缩系数，1 /MPa； h — 有效厚度，m； J — 采油指数，m3/(d·MPa)； K — 渗