储层流体高压物性参数计算方法(专题).ppt
储层流体高压物性参数计算方法（专题）,什么是储层流体？,储层是指具有孔隙性和渗透性、油气能在其中流动的岩层叫储集层，简称储层。,储层流体是指油气藏中存在的主要流体，它们就是通常所指的油、气、水。,储藏有石油的储集层叫储油层，简称油层；储藏有天然气的储集层叫储气层，简称气层；同时储藏有石油和天然气的储集层叫储油气层，简称油气层；同时储藏有石油、天然气和水的储集层叫储油气水层，简称油气水层；等等。,什么是储层流体高压物性？,储层流体物性是指储层内流体的物理化学性质及其在地层条件下的相态和体积特征。,储层流体高压物性是指储层内流体在地层条件下（高温、高压条件下）的物理化学性质。由于原油、天然气以及地层水都不是单一物质，而是混合物。因此，它们都不可以采用固定的模式去评价。所以，只有“具体问题具体解决”。,哪些是储层流体高压物性参数？,它们主要是指流体的粘度、相对密度、体积系数、压缩系数、分子量、天然气的偏差因子、原油的溶解油气比和两相体积系数等等。,为什么要进行储层流体高压物性参数计算？,对储层流体物性的评价是油气藏工程研究中的首要环节，也是最重要的环节。由于储层流体物性参数是油气藏的重要参数，因此，在可能的情况下，应当在实验室中进行测定。,然而，在实际油田开发和生产过程中不易获得更多的实际测定数值，尤其是新近开发的油气藏，因此采用以“最少的、最容易收集的资料来较为准确地估算储层流体物性参数”就显得十分必要了。,依据储层流体物性的参数是压力、温度、油气相对密度以及其组成组分的函数，在对比分析研究的基础上，从国内外的许多相关经验公式中，筛选出了一套最佳的经验公式，用来计算储层流体的高压物性参数。,1.经验关系式法,2.非线性神经网络法,主要方法:,第一节 地层天然气的物性参数计算,地层天然气主要是指干气气藏气体、凝析气藏气体和煤层气气体，其高压物性参数包括天然气的偏差因子、压缩系数、体积系数和粘度。,1、拟临界压力ppc和拟临界温度 Tpc的计算,计算方法一：组分分析方法,表3 天然气中常见组分主要物理化学性质,计算方法二：相关经验公式方法,干气:,凝析气:,注意：上式是对于纯天然气适用，而对于含非烃CO2、H2S 等可以用Wichert和Aziz修正。,Standing(1941)提供的相关经验公式,注意：上式是对于纯天然气适用，而对于含非烃CO2、H2S 等可以用Wichert和Aziz修正。修正常数的计算公式为:,修正后的拟临界压力和温度公式为:,若再考虑N2的影响，根据Smith等人的研究成果，将上面两式改为：,2、拟对比压力PPr和拟对比温度TPr的计算,对比参数就是指某一参数与其应对应的临界参数之比：即,3、天然气偏差因子Z的计算,定义：天然气偏差因子Z的计算是指在某一压力和温度条件下，同一质量气体的真实体积与理想体积之比值。,1974年,Dranchuk和Purvis等人在拟合Standing-Katz图版的基础上,提出了计算偏差因子Z的如下牛顿迭代公式。,在已知ppr和Tpr的情况下，由(54)式求解Z时，采用迭代法。即首先给定Z的一个初值Zo(例如Zo=1.0)，由(54)式求出ro，作为(51)式迭代的初值。比较ro与用(51)式计算所得的r之值，如两者相差较小(一般小于0.00001)，则可认为求得了r值，再由(54)式求得Z值。否则，用新的r代替原ro进行重复计算，直到满足精度要求为止。,(51),(54),4、天然气的压缩系数,天然气的压缩系数就是指在恒温条件下，随压力变化的单位体积变化量，即,经过一系列的推导，可以获得如下的表达式：,由真实气体的状态方程，得下式：,式中：,5、天然气的体积系数,天然气的体积系数就是指在地层条件下，某一摩尔气体占有的实际体积，除以在地面标准条件下同样摩尔量气体占有的体积，由下式表示：,在实际计算时，通常取Zsc=1.0，而当Psc=0.101MPa，Tsc=293K时，由上式得：,6、天然气的粘度,流体的粘度被定义为流体中任一点上单位面积的剪应力与速度梯度的比值，它是流体(气体或液体)内摩擦而引起的阻力。因此，天然气的粘度即为天然气内部某一部分相对于另一部分流动时摩擦阻力的量度。作为一种混合气体的天然气，其粘度与压力、温度和组成有关。天然气粘度的大小对于它在地层中或管路中的流动计算是很重要和必备的参数。粘度的高低表明流体流动的难易，粘度愈大，流动阻力愈大，流动越难。,Lee和Gonzalez等人根据四个石油公司提供的8个天然气样品，在温度为37.8171.2和压力为0.10155.16MPa的条件下，进行粘度和密度的实验测定，利用测定的结果得到了如下的相关经验公式：,其中：,7.凝析气藏露点压力的计算,在原始条件下，确定凝析气藏露点压力的大小，是一项重要的工作内容。它可以通过井下高压物性取样和PVT分析确定，也可以利用Nemeth和Kennedy提供的如下相关经验公式确定：,式中：,该方法计算结果的平均偏差为7.4%。,第二节 地层原油的物性参数计算,地层原油的高压物性参数包括原油饱和压力、溶解气油比、压缩系数、体积系数和粘度。,1、原油饱和压力,原油饱和压力系指在地层条件下，原油中的溶解气开始分离出来时的压力。,饱和压力Pb（泡点压力）是原油中第一个气泡出现时的最高压力，该项重要参数可通过等组成膨胀实验得到。,1947年，Standing根据加利福尼亚油田的22种烃类系统的105个实验测量的泡点压力，确定了图解相关关系，1981年用下述数学表达式表述了图解相关关系：,其中：,1)Standing法,其平均相对误差为4.8%。,Glaso于1980年根据北海6个油藏的26个和其他的19个流体的PVT分析样品，按照Standing的研究方法，获得：,其中：,2)Glaso方法,其在1.03448.263MPa压力范围内的标准差为6.89%，而在13.78948.263MPa压力范围内的标准差为3.84%。对挥发油，Glaso建议把方程中的温度指数修改为0.130。,1988年，Marhoun对来自中东的69种烃类混合物进行PVT分析后，得到160个实验确定的泡点压力，并得出了如下估算Pb的经验公式：,3)Marhoun方法,当计算值与实测值比较时，其平均相对误差为3.66%。,1993年，Petroksy和Farshad依据墨西哥海湾原油系统进行了81次实验分析，并用非线性多元回归软件得到如下估算Pb的经验公式：,4)Petroksy-Farshad方法,当计算值与实测值比较时，其平均相对误差为3.66%。,其中：,2.溶解气油比,在地层条件下的原油溶解有天然气，单位体积原油中天然气溶解量称为天然气溶解度，也称为溶解气油比。,1947年Standing采用加利福尼亚的22种原油和天然气烃类体系的105个实验数据，以压力、气体相对密度、原油相对密度和体系温度为关联参数，提出一种估算Rs的图版，1981年Standing利用该图版，得到一种数学计算方法。,1)Standing方法,(1)当ppb时,(2)当ppb时,Beggs给出了一种根据油藏压力、油藏温度、脱气原油相对密度和采出气的相对密度来计算溶解气油比的相关经验公式。,2)Beggs方法,3.原油压缩系数,原油压缩系数是油藏弹性能量的一个量度。它定义为在地层条件下每变化1MPa压力，单位体积原油的体积变化率。,Vaquez和Beggs提出了一种用于估算泡点压力以上压力条件下的原油压缩系数的相关经验公式：,Villena-Lanzi给出了一种估算重质油压缩系数的方法，这一关系式在压力低于泡点压力时精度较高。,4.原油地层体积系数,原油地层体积系数为采出地面条件下1m3 的脱气原油体积所占有的地层原油体积量，即：,1）Beggs-Standing的计算公式(1)Standing给出了一种根据溶解汽油比、溶解气的相对密度、脱气原油的相对密度以及油藏温度等估算单相原油地层体积系数的关系式。Beggs以方程的形式给出了Standing的估算原油地层体积系数的计算公式。,当ppb时,当p=pb时,当ppb时,(2)原油的两相体积系数（总体积系数）,原油的总体积系数是指当油层压力低于饱和压力时，地层中原油和析出气体的总体积与它在标准状态下的体积之比。,经过一系列推导得原油的两相体积系数为：,Glaso利用北海油田和其他地区的pVT分析资料，由回归分析法获得的相关经验公式：,2）Glaso的经验公式,5.原油粘度的计算,(1)Egbogah给出了计算压力小于或等于饱和压力的脱气原油的粘度：,(2)Beggs和Robinson给出了含溶解气的原油与脱气原油粘度之间的关系：,根据对2073个油样分析的结果表明：计算的平均绝对误差为1.83。,(3)当压力高于饱和压力时，Vazques和Beggs给出了下式：,根据对3143 个油样的分析表明，上式的平均绝对误差为7.54。,第三节 地层水的物性参数计算,地层水的高压物性参数包括溶解气水比、压缩系数、体积系数和粘度。,1.溶解气水比,溶解气水比系指标准状况下1m3的水在气藏温度、压力条件下所溶解的天然气在标准状态下的体积数。McCain提出了一个用来估算溶解气水比的关系式：,关于系数A、B、C值：,2.地层水的等温压缩系数,地层水的等温压缩系数取决于压力、温度、溶解气水比以及地层水的矿化度。,地层水的压缩系数是指：单位体积地层水在单位压力改变时的体积变化值，即：,当已知气藏的地层压力、地层温度和天然气在地层水中的溶解度之后，可由如下的相关经验公式计算地层水的压缩系数：,其中：,3.地层水的体积系数,McCain提出了一个用来估算地层水的体积系数的关系式：,地层水的体积系数定义为气层条件下水的体积与标准状况下纯水(即脱气水)的体积之比，即,其中：,4.地层水的粘度,(1)McCain提出了计算在大气压和油（气）藏温度下水的粘度的关系式：,上式的精度度为95%。,(2)在油（气）藏条件下的地层水的粘度为：,当地层压力小于68.95MPa时，其精度为96%,5.地层水的密度,(1)在地层条件下，确定纯水密度的相关经验公式为：,(2)若考虑矿化度的地层水的密度：,6.气井的水气比,气井凝析水的产量可由如下的经验公式计算：,第四节 BP神经网络法,利用BP神经网络模型建立油藏流体高压物性参数的预测模型，来对原油高压物性参数进行预测，是油藏高压物性参数研究的一种较为新的方法，其优势在于，不管参数之间的关系多么复杂，它都能通过自学习和自适应过程，把这种复杂的关系隐含到神经元之间的连接权重，因此它对复杂的高压物性参数问题，特别是非线性问题有较强的解决能力。虽然这种方法获得的解释模型是一个二维矩阵，不能直观地反映参数之间的相互关系，但使用起来却极为方便。,图 BP网络的结构,1.BP网络的结构,BP网络是应用最广泛的一类网络，它由输入层、隐含层、输出层构成，输入、输出层的节点数目可由输入、输出参数的数目来决定；隐含层可以是多层，一般而言，隐含层较多，则较少次数的权数的权值调整网络就能学到样本知识，并以权值分布的形式贮存起来：不过隐含层多，需要调整的权值个数就大幅度增加，网络的学习就会占用更多的时间。,2.激活函数的选择,构成BP网络的神经元所用的激活函数必须是处处可导的。一般地，多数设计者都使用s型函数。对一个神经元来说，取它的网络输入,该神经元的输出为,3.BP算法,BP算法主要包含4步，这4步被分为两个阶段：,A 向前传播阶段（1）从样本集中取一个样本（Xp,Yp），将Xp输入网络；（2）计算相应的实际输出Op。在此阶段，信息从输入层经过逐级的变换，传送到输出层。这个过程也是网络在完成训练后正常运行时执行的过程。在此过程中，网络执行的是下列运算：,B 向后传播阶段（1）计算实际输出Op与相应的理想输出Yp的差；（2）按极小化误差的方式调整权矩阵。,这两个阶段的工作一般应受到精度要求的控制，取,作为网络关于第p个样本的误差测度。将网络关于整个样本集的误差测度定义为,4.高压物性参数预测网络模型的建立过程,(1)网络模型的建立 神经网络方法用于油藏流体物性参数预测的基本思路是，首先应用已知油藏区块的取样数据(样本)对网络进行训练，然后再将未知井的取样数据输入训练后的网络，最后预测出该区块未知井段的流体性质。,图PVT BP网络预测模型拓扑结构,(2)数据的归一化预处理,由于各种取样数据的量纲不同，其数值相差很大。就像泡点压力和原油粘度这种情况，但是如果输出数据的变化范围很大的话，神经网络是不稳定的。因此，必须在ANN进行学习训练阶段前，消除这些困难。,在运用神经网络预测油藏流体高压物性参数之前，无论是学习样本还是预测数据，都需要进行归一化处理，将它们统一刻度在01之间。在规一化过程中，各取样数据归一化极值的选取至关重要。在实际工作中，为了尽可能减小人为性和经验性选取归一化极值带来的误差，最好在一个地区对每种取样数据均采用统一的归一化标准，在统一标准的基础上，也可对特殊地层段归一化极值进行微调。,(2)数据的归一化预处理,对于具有近似线性特征的输入信息，采用线性归一化公式：,而对于地层原油粘度等具有非线性特征的曲线，则采用对数归一化公式：,(3)学习样本的选择,经过大量学者的研究和实践，发现在样本进入学习之前，需按以下原则做适当的取舍:(1)剔除个别明显异常的奇异点；(2)适当减少相同特征点数；(3)适当补充特征明显的典型样本；(4)适当控制训练样本的总量；,(4)BP神经网络学习和训练程序框图,5实际预测及效果分析,(1)饱和压力预测分析,图 某地区饱和压力预测值与实测值(ANN模型)比较,图 某地区饱和压力预测值与实测值(Standig模型)比较,(2)体积系数预测分析,图1某地区Bo预测值与实测值(ANN模型)比较,图 某地区Bo预测值与实测值(Standing模型)比较,Good luck every day!,Thank You!,
