专利名称：用于推算在亚层土的间隔区域中运动的流体的方法
技术领域：
本发明涉及一种用于推算在亚层土的间隔区域中运动的流体体积 的方法，例如天然沉积物，也称为储层，其在工业上可适用于油田， 另外，也可以用于监测碳氢化合物的开采和储藏，并减少勘探阶段期 间和路上及海上储层的开采期间的采矿危险。
背景技术：
在亚层土资源勘探时，通常依赖于重力场的竖直分量和竖直梯度 的测量。根据这些数据的分析，可以推导出有关特定地点的亚层土密 度分布的信息。
基于该方法，也可以获得储层内部存在的碳氢化合物的质量变化。 该技术基于的思想是，由于储层内部的碳氢化合物运动与密度变 化相关联，所以可以借助于竖直重力梯度测量进行评估。
实际上，自从三十年代以来，重力场梯度的测量就已经成功地用
于亚层土资源勘探。自从1936年以来，人们就已经认识到利用竖直梯 度的重要性，但由于竖直梯度具有较好的分辨力并且对地域效应相对 不是那么敏感，所以竖直梯度常常具有根据重力场数据不能容易地获 得的特定结构。
竖直重力场梯度的测量可以借助于称为重力梯度仪的特定仪器进行。
作为选择，某个点的重力场竖直梯度也可以借助于差不多同期获 得的不同高度的两个重力测量值利用好的近似法进行测量。
在该第二种情况下，在利用地质条件解释在场地获得的数据之前， 布格(Bouguer)异常方面的减少很频繁，据此排除不希望的效应， 随后实施竖直梯度的计算和分析。
5要进行的最重要的校正如下
一仪器零点漂移
一潮沙校正
一綷度校正
一自由空气校正
一布格校正
一地形4交正
这些校正的描述如下。
仪器零点漂移重力计的数据读数由于形成仪器本身的材料的弹 性特性而随时间变化。仪器零点漂移可以通过对同一测量点在不同时 间，典型每l-2小时，进行重复测量而容易地确定。参照笛卡尔轴的 表示给出零点漂移曲线，对许多重力计来说，零点漂移曲线为线型。 根据测量时间，在随后的测量点进行的各个测量减去定值。 潮汐校正测量的零点漂移实际上还包含由于月亮-太阳吸引力引 起的海洋类型效应(潮汐)的进一步作用。根据理论基础，借助于对 该效应进行量化的公式，例如朗曼(Longman)公式，计算要进行的 校正。
绰度校正地球的旋转及其赤道的膨胀都使重力随绵度而增大， 当观测的重力数据减少时，必须对此进行考虑。
自由空气校正这是为了考虑测量点海拔而使用的校正。
布格校正该校正用来考虑由介于测量点与基准面之间的质量引 起的吸引力。在1749年，布格建议，该附加吸引力可以像由厚度等于 测量点从海平面的海拔的无限大的水平板引起的作用那样来计算。
地形校正所述板的近似在带有铰接地貌趋势的区域是不令人满 意的。
在这些条件下，适于增加一校正以便考虑上述板之上的质量以及 那些在布格校正中错误地减掉的作用。
在减去所列出的之后，按照如下所述的方式计算场地的竖直梯度。 亚层土密度随时间的变化可以根据竖直重力场梯度的测量进行测量和监测。
该方法早已用来测量和监测水层和地热场地。
本发明的目的是提供一种用于推算流体例如在地下的间隔区域中 运动、即例如开采、注入和/或存储的流体例如液体和/或气态碳氢化 合物的体积的方法。
在上述目的中，必须对每单个场地储层画出校准特性曲线。
依照本发明，这些及其他目的可以借助于如权利要求1所述的推 算在间隔区域中移动的流体体积的方法来实现。
该方法的进一步的特征是从属权利要求的主题。
依照本发明，推算在间隔区域中的流体体积的方法的特征和优点 将从下面直观且非模拟的描述中变得更加显而易见。
发明内容
依照本发明，推算亚层土中的流体体积的方法包括"构架/校准区 域"的第一梯度测量阶段，基于此，按照统计近似法确定模型，以便 量化在亚层土中开采、注入和/或存储的流体、例如碳氢化合物体积的 运动。
更具体地说，在某一时间段重复实施一系列测量活动，以便借助
于模型量化开釆、注入和/或存储的流体运动量。
然后通过统计分析获得的数据求出该模型的精确度。
测量作业活动在勘测区域中适当打乱的 一 系列测量点进行。
对每单个测量点，利用高精度重力仪和适当定位的三脚架实施测
量，以允许在离开地面的某一恒定距离进行测量。
对于如上所述的地球物理应用，建议使用精确度不低于nGal(微 伽)的重力计。
根据操作的观点，测量活动由下列阶段构成。
获得第一地面重力测量Gb。t，接着利用三脚架将重力计定位在离 开地面一段距离dh处，获得第二重力测量Gt。p。
在各单个测量点和整个勘测测量期间，与地面的距离dh优选保持恒定。事实上，申请人已经注意到，通过在各单个测量点和对整个 勘测使与地面的距离保持恒定，梯度异常的测量和确定更加精确。
基于所测量的Gb。t和Gt。p，然后确定对于前述的效应进行校正的 重力值G、。t、 G
top0
借助于下列公式，从校正的重力值获得重力场的竖直梯度值
借助于带有毫米精确度的激光测距仪，确定在高度水平dh的各 单个测量点的两个相对测量值之间的差。
然后随着时间变化在同 一 区域重复该获得和详细设计的程序，强 调是同一测量点。
然后获得时滞(time lapse)信号，即与亚层土中一定时间段的重 力变化有关的信号。
在不同的时间，根据同 一测量点上实施的两次梯度勘测之间的差， 计算时滞信号
1Xj产VGGi画VGGj( 2 )
这样，可以评价所检验的区域的哪些地带遇到了密度变化，从而 披露存在流体相对运动的地方。
这样，获得有关亚层土中流体体积变化的定性结果，该体积变化 与不同情况的开采、注入和/或存储作业有关。
申请人随后鉴定每单个存储层或开采储层的校准特性曲线，所述 校准特性曲线将某一时滞所检测的梯度值变化与同一时间段内移动或 开采、注入和/或存储的流体对应体积相关联。
申请人还开发、测试和证实了所述校准曲线的评价方法。
只要考虑质量随时间的运动的定性评价，申请人就认为它适于找 到一条校准曲线，所述校准曲线将流体体积变化与时滞数据相互关联， 而不是与单个评价的重力梯度数据相互关联。
为此，必须评定与开采、注入和/或存储作业相关的流体体积变化 的特征参数。已经证明，利用时滞值在检验区域的整个限定范围上的积分适用
于该目的
JJtt^o,:^*( 3a )
由此获得的值优选相对于第一时滞TLu的积分进行标准化，从而
将第一评价作为基准点。
由此可获得量纲参数Pu，赋值为体积的定量总体变化，并限定为
I JJ"7Z2(x，力^:办
可以根据测量活动编号i所获得的值计算参数plin
同样，在整个勘测区域上分布的n个测量点所测量的时滞值总和
也适于作为参数
Z7X,(") (3b )
而且在这种情况下，由此获得的值优选相对于第一时滞TLu的总 和进行标准化，以便获得量纲参数
2 )
凡三
Ub)
然后申请人根据至少三次获得活动，在已知运动(开采、注入和/ 或存储)的流体体积的情况下构造校准曲线。
根据与已知运动的体积AVj的时滞相关联的参数P1P按照近似法 重构参数Pu与相应的已知开采、注入和/或存储的体积厶Vj之间的关 系。
由此获得的结果的各种控制表明并证实由此确定的变化规律允 许运动(开采、注入和/或存储)的流体体积变化从"用于构架/校准 区域"的第一阶段之后的梯度测量活动开始进行推算。
以纯粹示例的方式，在此描述了用于确定将测量的梯度值与运动 流体即开采、注入和/或存储的流体、例如碳氢化合物的体积相互关联
的可能的规律的方法。
假定参数Pu与运动的碳氢化合物体积AVi的线性关系为
9P=a+bAV ( 5 )
该直线可构造成用于推算从向前第四勘测开始的AV。 与第一个三次测量活动和已知的运动的相对体积相关联的参数
Pn是用于限定参数a和b所必须的。
为了求出参数a和b，利用相对不确定度&a和8b,例如可以利
用使理想直线的理论值与所观测的数据之间的差最小化的最小卡方近
似法，即下列表达式
%2=，"，2 (6)
一旦这样确定了a和b,已知参数Pu的值，就可以求出AV的对 应值。
该方法的不确定度和灵敏度取决于参数a和b的不确定度以及有
关来自重力测量的参数P的不确定度8P。 用于推算运动的体积的等式是
△r = ^^ (7)
所以，相关的误差8AV通过下式给出
(8 ) 必k =
3AK
3尸
3a
5AF
的:
的
同样，可以采取合适程度的多项式变化规律和使用与近似法一样 .小卡方法(square minima method )或内插法。 本发明的方法目的的特征以及相应优点从上述描述中是显而易见
申请人增加了确定竖直重力梯度和对应时滞流程的算法，用于推 算一定时间段在储层中运动的流体体积。
最后，4艮明显，由此构思的方法可以进行进一步的改进和变化， 所有这些都包括在本发明的范围之内。
10
权利要求
1.一种用于推算在亚层土的间隔区域中运动的流体体积的方法，包括下列步骤a)在预先设定的时间间隔，对勘测场地实施一系列“用于构架/校准区域”的梯度测量活动，对于每次梯度测量活动，获悉容纳在亚层土的相应间隔区域中的流体体积的相对变化ΔVi；b)根据实施的测量，对于在前一个之后的每个“用于构架/校准区域”的测量活动，计算与整个勘察区域有关的参数P1i；c)根据已知的流体体积变化ΔVi和所计算的相应参数P1i，按照近似法，确定有关参数P1i与所述间隔区域中的流体体积变化的变化规律，基于在步骤c)所确定的变化规律，关联在那些在步骤a)下限定的“用于构架/校准区域”的测量活动之后的任何测量活动。
2. 如权利要求1所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，根据所述阶段a)中在预先设定的时间间隔进行的所述梯度测量，将在各个测量点在不同时间进行的两次梯度测量VGGj、 VGGj之间的差计算为时滞TLji。
3. 如权利要求2所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，根据时滞TLu值在勘测区域的整个限定范围上的积分计算所述参数P"。
4. 如权利要求3所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，所述参数Pu是时滞TLu值在勘测区域的整个限定范围上的积分，并相对于第一时滞TLu的积分进行标准化。
5. 如权利要求2所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，根据在勘测区域的整个限定范围上分布的n个测量点所检测的时滞TLu值的相加计算所述参数Pu。
6. 如权利要求5所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，根据在勘测区域的整个限定范围上分布的n个测量点所检测的时滞TLu值的相加，并相对于第一时滞TLu的总和标准化，来计算所述参数Pu。
7. 如权利要求4或6所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，所述阶段c)在于确定所述计算的Pu参数和已知的对应体积变化AVi之间的线性变化规律。
8. 如权利要求4或6所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，所述阶段c)在于确定所述计算的Pli参数和已知的对应体积变化△ Vi之间的多项式变化规律。
9. 如权利要求7或8所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，按照近似法，通过最小卡方法或内插法确定所述变化规律。
10. 如上述任一权利要求所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，所述阶段a)中的所述系列梯度测量活动至少为三次。
11. 如上述任一权利要求所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，借助于重力计和三脚架，对各个测量点实施用于梯度测量的所述阶段a),以便允许几乎同期获得有关不同海拔的两个重力测量值Gb。t和Gt。p。
12. 如权利要求11所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，从所述重力测量值Gb。t和Gt。p导出根据 一 系列影响因素校正的重力测量值Gb。t和Gt。p 。
13. 如权利要求12所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，按照以下规律，由校正的所述重力测量值G、。t和G 。p计算所述阶段a)的所述梯度测量，謡，。'—G(。p) (1)其中dh是单个测量点的两次测量之间的高度差。
14. 如权利要求13所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，借助于激光测距仪实现所述高度差dh。
15. 如权利要求13或14所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，对于整个勘测的所有测量点，所述高度差dh必须保持恒定。
16.如权利要求l-10任一所述的用于推算在亚层土的间隔区域中运动的流体体积的方法，其特征在于，借助于重力计，在各个测量点实施所述梯度测量活动。
全文摘要
本发明涉及一种用于推算在亚层土的间隔区域中运动的流体体积的方法，包括步骤在预先设定的时间间隔，对勘测场地实施一系列“用于构架/校准区域”的梯度测量活动，对于每次梯度测量活动，获悉容纳在亚层土的相应间隔区域中的流体体积的相对变化ΔV<sub>i</sub>；根据实施的测量，对于在前一个之后的每个“用于构架/校准区域”的测量活动，计算与整个勘察区域有关的参数P<sub>1i</sub>；根据已知的流体体积变化ΔV<sub>i</sub>和所计算的相应参数P<sub>1i</sub>，按照近似法确定有关参数P<sub>1i</sub>与间隔区域中的流体体积变化的变化规律，基于所确定的变化规律，关联在那些“用于构架/校准区域”的测量活动之后的任何测量活动。
文档编号G01V7/00GK101680966SQ200880018537
公开日2010年3月24日 申请日期2008年3月28日 优先权日2007年4月13日
发明者I·焦里, L·泰尔齐, L·萨维尼, M·安东内利, S·贾梅蒂 申请人:艾尼股份公司