CN115422859A - 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 - Google Patents

一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 Download PDF

Info

Publication number
CN115422859A
CN115422859A CN202211385246.6A CN202211385246A CN115422859A CN 115422859 A CN115422859 A CN 115422859A CN 202211385246 A CN202211385246 A CN 202211385246A CN 115422859 A CN115422859 A CN 115422859A
Authority
CN
China
Prior art keywords
steam
calculating
oil
vertical
thick
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211385246.6A
Other languages
English (en)
Other versions
CN115422859B (zh
Inventor
刘�东
刘永辉
赖南君
胡廷惠
唐雷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN202211385246.6A priority Critical patent/CN115422859B/zh
Publication of CN115422859A publication Critical patent/CN115422859A/zh
Application granted granted Critical
Publication of CN115422859B publication Critical patent/CN115422859B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/24—Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection
    • E21B43/2406—Steam assisted gravity drainage [SAGD]
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，对于目标区块处，包括：S1收集地质油藏基本参数；S2计算蒸汽微元在油层垂直向的受力值；S3计算蒸汽微元在油层水平向的受力值；S4计算蒸汽微元在油藏垂直方向的合力和瞬时渗流速度；S5计算蒸汽微元在油藏水平向的合力和瞬时渗流速度；S6计算垂向分流系数，绘制垂向分流系数与井间距离的图版；S7计算注入蒸汽的超覆程度，绘制超覆程度与井间距离的图版；S8定量计算倾斜厚层稠油注蒸汽纵向波及系数。本发明达到的有益效果是：实现了在海上大井距热采且无观察井的条件下，定量预测倾斜厚层稠油注蒸汽吞吐或者蒸汽驱纵向波及系数，从而有助于稠油热采方案设计。

Description

一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法
技术领域
本发明涉及油井勘探开采技术领域，特别是一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法。
背景技术
室内研究和矿场实践表明，热力采油技术（如蒸汽吞吐、蒸汽驱、蒸汽辅助重力泄油等）是一项可以大幅度提高稠油油田采收率的有效技术，是目前陆上稠油开发的有效手段。渤海目前稠油热采处于规模化推广阶段，稠油热采油藏的粘度不断增加，油藏厚度不断扩大（地层条件原油粘度从50mPa•s到50000mPa•s，油藏单砂体厚度从6m扩展到40m），定量获取厚层稠油油藏注蒸汽吞吐或蒸汽驱纵向波及系数，有助于评价注蒸汽热采（包括蒸汽吞吐、蒸汽驱等）的有效加热范围、热采渗流模式，从而指导热采方案设计和注采参数优化。
目前关于注蒸汽的纵向波及系数，陆地油田通常在井间设置观察井，通过观察井适时监控地层中的温度，从而通过地层中的温度判断纵向上的动用情况。而海上稠油热采，由于开发成本高，通常采用大井距热采，且井间没有观察井，因此定量获取厚层稠油油藏注入蒸汽纵向波及系数成为难题。
发明内容
本发明的目的在于克服现有技术的缺点，提供一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，解决了定量获取厚层稠油油藏注入蒸汽纵向波及系数的难题。
本发明的目的通过以下技术方案来实现：一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，包括以下步骤：
S1、收集目标区块地质油藏基本参数；
S2、计算目标区块蒸汽微元在油层垂直方向上的受力值；
S3、计算目标区块蒸汽微元在油层水平方向上的受力值；
S4、计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度；
S5、计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度；
S6、计算垂向分流系数，并绘制垂向分流系数与井间距离的图版；
S7、计算注入蒸汽的超覆程度，并绘制超覆程度与井间距离的图版；蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值，用D表示；
S8、定量计算倾斜厚层稠油注蒸汽纵向波及系数。
进一步地，所述S1中的基本参数：是指收集目标区块的基础数据，包括油层厚度、地层倾角、注采井间距离、油的密度、水的密度、气体粘度、界面张力、接触角、孔隙半径、垂直方向渗透率、水平方向渗透率。为通过受力分析和渗流规律计算注入地层中蒸汽的分布做好数据准备。
进一步地，所述S2中计算目标区块蒸汽微元在油层垂直方向上的受力值，包括重力、浮力；
S21、重力的计算
根据蒸汽微元的密度
Figure 963421DEST_PATH_IMAGE001
为、蒸汽微元的体积
Figure 349403DEST_PATH_IMAGE002
，再根据重力加速度g，计算得蒸汽微元的重力F重力；
由于地层有倾角
Figure 597982DEST_PATH_IMAGE003
，以油层水平方向为X轴、油层垂直方向为Y轴，建立直角坐标系，计算蒸汽微元的重力在X方向的分量Fx重力、Y方向的分量FY重力；
S22、浮力的计算；
在注入蒸汽突破之前，蒸汽微元受到的浮力F浮力看作是恒定的，方向为竖直向上；
浮力F浮力，根据油相的密度
Figure 963497DEST_PATH_IMAGE004
、水相的密度
Figure 246711DEST_PATH_IMAGE005
、蒸汽微元的体积
Figure 538015DEST_PATH_IMAGE006
，重力加速度g，计算得到；
由于地层有倾角
Figure 70628DEST_PATH_IMAGE007
，以油层水平方向为X轴、油层垂直方向为Y轴，建立直角坐标系，计算浮力F浮力在X方向的分量FX浮力、Y方向的分量FY浮力。
在步骤S2中，通过计算的重力和浮力值是决定注入蒸汽微元在垂向上运动的动力，由于考虑了浮力和地层倾角的影响，可使描述倾斜油藏垂直油层方向的蒸汽运动更加精确。
进一步地，所述S3中计算目标区块蒸汽微元在油层水平方向上的受力值，包括毛管力和驱替压力：
S31、计算驱替压力，即计算注采井间主流线上的注采压力梯度，公式为
Figure 738369DEST_PATH_IMAGE008
其中，
Figure 876090DEST_PATH_IMAGE009
为生产井产量；
Figure 338295DEST_PATH_IMAGE010
为蒸汽注入量；
Figure 561466DEST_PATH_IMAGE011
为油层厚度；r为任意点到生产井的距离；p(r)为r处的压力；
Figure 767319DEST_PATH_IMAGE012
为r处的压力梯度；L为注采井间距离；
Figure 258081DEST_PATH_IMAGE013
为渗透率；
Figure 156767DEST_PATH_IMAGE014
为原油粘度；
S32、计算毛管力
对于蒸汽，受到毛管力为
Figure 867234DEST_PATH_IMAGE015
式中，
Figure 345620DEST_PATH_IMAGE016
为界面张力；
Figure 457933DEST_PATH_IMAGE017
为接触角；
Figure 261941DEST_PATH_IMAGE018
为孔隙半径；
在水湿情况下，毛细管力是注入水运移的动力，取正值；在油湿情况下，毛管力是注入水运移的阻力，取负值；
由于地层有倾角
Figure 459704DEST_PATH_IMAGE019
，根据该倾角，计算平行油藏方向的水平毛管力
Figure 804097DEST_PATH_IMAGE020
和垂直油藏方向的垂直毛管力
Figure 770916DEST_PATH_IMAGE021
。
在步骤S3中，通过注入蒸汽在水平方向上受力主要包括驱替力和毛管力，是决定蒸汽在水平方向上运动的动力，由于考虑了毛管力和地层倾角的影响，可使描述倾斜油藏油层水平方向的蒸汽运动更加精确。
进一步地，所述步骤S4计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度；
S41、计算目标区块蒸汽微元在油藏垂直方向的合力
根据蒸汽微元在Y方向的重力分量FY重力，蒸汽微元在Y方向的重力分量FY浮力，垂直油藏方向的垂直毛细管力FY毛管力，计算得到蒸汽微元在油藏垂直方向的合力FY合力；
S42、计算目标区块蒸汽微元在油藏垂直方向的瞬时渗流速度
注入蒸汽沿油层斜上方渗流，可将速度分解为沿油层X方向的速度和Y方向的速度；
蒸汽在重力、浮力及毛管力的共同作用下，根据达西定律，计算得到Y方向瞬时渗流速度，即在油藏垂直方向的瞬时渗流速度VX。
在步骤S4中，由于是倾斜油藏，并且S2、S3计算的各种力在油藏垂直方向均有分力，且作用方向不一致，因此需要求取S2和S3在垂直方向的合力，并根据合力和渗流理论，计算得到油藏垂直方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响，让垂直反向的瞬时渗流速度计算更加精确。
进一步地，所述S5计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度；
S51、计算目标区块蒸汽微元在油藏水平方向的合力
根据蒸汽微元在X方向的重力分量FX重力，蒸汽微元在X方向的重力分量FX浮力，驱替压力，垂直油藏方向的垂直毛细管力FX毛管力，计算得到蒸汽微元在油藏垂直方向的合力FX合力；
S52、计算目标区块蒸汽微元在油藏水平方向的瞬时渗流速度
根据达西定律，油滴在浮力及毛管力的共同作用下，计算得到X方向瞬时渗流速度为VY。
在步骤S5中，由于是倾斜油藏，第S2和S3计算的各种力在油藏水平方向均有分力，且作用方向不一致，因此需要求取S2和S3在水平方向的合力，并根据合力和渗流理论，计算得到油藏水平方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响，让水平方向瞬时渗流速度计算更加精确。
进一步地，所述S6计算垂向分流系数，并绘制垂向分流系数与井间距离的图版；
S61、计算垂向分流系数
垂向分流系数为垂向瞬时渗流速度与总瞬时渗流速度的比值，计算方程为
Figure 512870DEST_PATH_IMAGE022
；
S62、
Figure 932350DEST_PATH_IMAGE023
为纵坐标，井间距离L为横坐标，绘制垂向分流系数与井间距离L的图版。
在步骤S6中，根据S4、S5计算得到的垂直方向和水平方向瞬时渗流速度的计算结果，可以通过计算垂向分流系数，定量计算出注入蒸汽在垂向上的突进情况，解决了海上油田由于没有观察井，难以掌握注入蒸汽油层内分布情况的难题。
进一步地，所述S7计算注入蒸汽的超覆程度，并绘制超覆程度与井间距离的图版；
S71、计算注入蒸汽的超覆程度
蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值，用D表示；
将注采井间距离划分为N块，当蒸汽前缘运移到注采井间 i 块( 0 ＜ i ＜ N)时，利用此处的径向蒸汽运移量与垂向流量比的乘积和第i-1块的超覆程度迭代，计算出此时的蒸汽超覆程度，计算公式为
Figure 283697DEST_PATH_IMAGE024
；
S72、绘制超覆程度与井间距离的图版
以
Figure 105022DEST_PATH_IMAGE025
为纵坐标，井间距离为横坐标，绘制超覆程度
Figure 250833DEST_PATH_IMAGE026
与注采井间距离L的图版。
由于S6计算的垂向分流量值是一个变化的量，在注采井间的不同位置具有不同的值，通过S7，可以描绘出注采井间不同位置处的蒸汽分布情况
进一步地，所述S8定量计算倾斜厚层稠油注蒸汽纵向波及系数；
S81、根据S7的图版，计算得到厚层稠油蒸汽吞吐初期不同分块的蒸汽超覆程度，计算公式为
Figure 157609DEST_PATH_IMAGE024
，
其中，
Figure 47067DEST_PATH_IMAGE027
表示第
Figure 988479DEST_PATH_IMAGE028
块的垂向分流量；
S82、通过
Figure 69305DEST_PATH_IMAGE029
，计算每个区块纵向波及无因次厚度
Figure 260115DEST_PATH_IMAGE030
；
S83、用每个区块的长度乘以无因次厚度
Figure 953264DEST_PATH_IMAGE031
得到每个区块的面积
Figure 483603DEST_PATH_IMAGE032
；
S84、用
Figure 971216DEST_PATH_IMAGE033
计算得到每个区块的每个区块面积厚度乘积；
S85、计算得到面积加权后的纵向波及厚度
Figure 852584DEST_PATH_IMAGE034
，计算公式为
Figure 349425DEST_PATH_IMAGE035
；
计算得到面积加权后的纵向波及系数V k ，计算公式为
Figure 796586DEST_PATH_IMAGE036
。
S7计算的蒸汽超覆程度Di，在注采井间的不同位置具有不同的值，通过S8，可以描绘出注采井间任何剖面的蒸汽分布面积，根据面积加权，可以定量计算出注采井间的蒸汽波及厚度。这样得到的波及厚度和波及系数值，考虑了注采井间不同位置的超覆程度差异，使得计算结果更加贴近矿场实际情况。
本发明具有以下优点：
（1）基于注入蒸汽在地层中的受力分析，提出了通过计算注采井间距离处的纵向分流系数和不同位置的蒸汽超覆程度来预测倾斜厚层稠油油藏的纵向波及系数，能有效解决海上大井距热采且无观察井的条件下，定量预测倾斜厚层稠油注蒸汽吞吐或者蒸汽驱纵向波及系数的难题，为海上稠油油田的热采方案设计提供了支撑；
（2）传统的注蒸汽的纵向波及系数，陆地油田通常在井间设置观察井，通过观察井适时监控地层中的温度，从而通过地层中的温度判断纵向上的动用情况；
通俗地说，传统的注蒸汽的纵向波及系数，是通过经验情况进行实验的，最终获得相应系数；根据经验的方式，使得系数的结果非常不准确；若地质情况复杂，以经验的方式测得的数据不可靠；
本发明提供的测量方式，是在对注入蒸汽在地层中受力分析的基础上，定量表征出倾斜厚层稠油油藏注蒸汽开发不同方向上的瞬时渗流速度，根据距离注汽井点的不同距离计算出垂向分流系数，并计算得到注采井间不同位置的蒸汽超覆程度，从而定量计算出倾斜厚层稠油注入蒸汽的纵向波及系数；
本发明提供的测量方式，通过计算的方式获得系数，不仅中间数据计算非常精确，而且最后经计算获得的结果也非常精确；为确定海上厚层稠油油藏注蒸汽吞吐和蒸汽驱方案设计、注采参数优化、水平井纵向位置优化等提供了非常可靠的依据；
（3）本方案的测量方式，不凭借经验，而是实际计算出来的，操作性强，应用广泛。
附图说明
图1为定量预测倾斜厚层稠油注蒸汽纵向波及系数的主要步骤；
图2为实施例1倾斜厚层稠油油藏（润湿性：油湿）地层中蒸汽微元受力示意图；
图3为实施例1注采井间不同距离处的垂向分流量图版；
图4为实施例1注采井间不同距离处的注入蒸汽超覆程度图版；
图5为实施例2中典型油田注采井间不同距离处的垂向分流量图版；
图6为实施例2中典型油田注采井间不同距离处的注入蒸汽超覆程度图版；
表1 为实施例1中A目标区块的基础数据；
表2 为实施例1中A目标区注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表；
表3 为实施例2中B目标区块的基础数据；
表4 为实施2中B目标区块注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表。
具体实施方式
下面结合附图对本发明做进一步的描述，但本发明的保护范围不局限于以下所述。
(实施例1)
一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，主要包括以下步骤（图1）：
S1，收集目标区块地质油藏基本参数等静态资料。收集目标区块的基础数据，包括油层厚度、地层原油粘度、渗透率等参数（表1），为通过受力分析和渗流规律计算注入地层中蒸汽的分布做好数据准备。表1中收集的为A目标区的基础数据。
表1 A目标区块的基础数据
Figure 455101DEST_PATH_IMAGE037
S2，计算目标区块蒸汽微元在油层垂直方向上的重力和浮力值。计算的重力和浮力值是决定注入蒸汽微元在垂向上运动的动力，由于考虑了浮力和地层倾角的影响，可使描述倾斜油藏垂直油层方向的蒸汽运动更加精确。如图2所示，厚层稠油油藏地层中某蒸汽微元主要受力有浮力、毛管力、重力和驱替压力。
对于蒸汽微元，受到的重力为：
Figure 823765DEST_PATH_IMAGE038
（A-1）
式中：
Figure 348464DEST_PATH_IMAGE039
为蒸汽的密度，g /cm3；
Figure 118973DEST_PATH_IMAGE040
为蒸汽微元的体积，cm3；g 为重力加速度，m/s2。
由于地层具有一定倾角
Figure 682810DEST_PATH_IMAGE041
，以油层水平方向为X轴，油层垂直方向为Y轴，建立直角坐标系统，计算重力在X方向和Y方向的分量：
Figure 538770DEST_PATH_IMAGE042
（A-2）
Figure 377414DEST_PATH_IMAGE043
（A-3）
在注入蒸汽突破以前，蒸汽微元受到的浮力可以近似看作是恒定的，方向为垂直向上，表达式为：
Figure 736851DEST_PATH_IMAGE044
（A-4）
式中：
Figure 2747DEST_PATH_IMAGE045
为油相的密度，g /cm3；
Figure 80424DEST_PATH_IMAGE046
为水相的密度，g /cm3；
Figure 519496DEST_PATH_IMAGE047
为蒸汽微元的体积，cm3；g 为重力加速度，m/s2。
由于地层具有一定倾角
Figure 497554DEST_PATH_IMAGE048
，以油层水平方向为X轴，油层垂直方向为Y轴，建立直角坐标系统，计算浮力在X方向和Y方向的分量：
Figure 934352DEST_PATH_IMAGE049
（A-5）
Figure 499325DEST_PATH_IMAGE050
（A-6）
根据表1的数据，油的密度取值0.95g /cm3，水的密度取值1.0g /cm3，重力加速度取值9.8 m/s2，地层倾角取值30°，蒸汽微元的体积取值1.0cm3。根据式子（A-2）、（A-3）、（A-4）、（A-5）计算得到单位体积（1.0cm3，）蒸汽所受到的X方向重力和浮力分别为：-4.655mN和-0.245mN；Y方向的重力和浮力分别为8.063mN和-0.424mN。
S3，计算目标区块蒸汽微元在油层水平方向上的毛管力和驱替压力梯度值。注入蒸汽在水平方向上受力主要包括驱替力和毛管力，是决定蒸汽在水平方向上运动的动力，由于考虑了毛管力和地层倾角的影响，可使描述倾斜油藏油层水平方向的蒸汽运动更加精确。
用（A-7）计算注采井间主流线上的注采压力梯度。
Figure 679771DEST_PATH_IMAGE051
（A-7）
式中：
Figure 13800DEST_PATH_IMAGE052
为生产井产量，cm3 /s；
Figure 621499DEST_PATH_IMAGE053
为蒸汽注入量，cm3 /s；
Figure 673769DEST_PATH_IMAGE054
为油层厚度，m；r 为任意点到生产井的距离，cm；
Figure 126747DEST_PATH_IMAGE055
为 r 处的压力梯度，10-1MPa/cm； L 为注采井间距离，cm；K 为渗透率，
Figure 82327DEST_PATH_IMAGE056
；
Figure 860927DEST_PATH_IMAGE057
为原油粘度，mPa·s。
对于蒸汽，受到的毛管力为：
Figure 400492DEST_PATH_IMAGE058
（A-8）
式中：
Figure 922741DEST_PATH_IMAGE059
为界面张力，mN/m；
Figure 28100DEST_PATH_IMAGE060
为接触角，°；
Figure 977601DEST_PATH_IMAGE061
为孔隙半径，mm。
在水湿情况下，毛细管力是注入水运移的动力，取正值；在油湿情况下，毛管力是注入水运移的阻力，取负值。毛管力的方向与注入水运移方向同向或反向，为便于计算与研究，分解为平行油藏方向的水平毛管力
Figure 738884DEST_PATH_IMAGE062
和垂直油藏方向的垂直毛管力
Figure 64823DEST_PATH_IMAGE063
。
Figure 227951DEST_PATH_IMAGE064
（A-9）
Figure 846889DEST_PATH_IMAGE065
（A-10）
根据表1的数据，界面张力取值0.015 N/m，接触角取值150°，孔隙半径取值0.00005m。由于接触角较大，可以判定为油湿，毛管力为注入水运移的阻力，取负值。根据式子（A-9）、（A-10）计算得到蒸汽所受到的X和Y方向的毛管力分别为0.0015Pa和-0.003Pa。
S4，计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度。由于是倾斜油藏，S2和S3计算的各种力在油藏垂直方向均有分力，且作用方向不一致，因此需要求取S2和S3在垂直方向的合力，并根据合力和渗流理论，计算得到油藏垂直方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响，让垂直反向的瞬时渗流速度计算更加精确。
联立式子（A-3）、（A-6）、（A-10），得到Y方向合力计算方程：
Figure 361047DEST_PATH_IMAGE066
（A-11）
根据表1的数据，联立（A-3）、（A-6）、（A-10），计算得到
Figure 959519DEST_PATH_IMAGE067
为7.635Pa。
注入蒸汽沿油层斜上方渗流，可将速度分解为沿油层X方向的速度和Y方向的速度。根据达西定律，蒸汽在重力、浮力及毛管力的共同作用下，Y方向瞬时渗流速度为：
Figure 977153DEST_PATH_IMAGE068
（A-12）
根据表1的数据，垂直方向渗透率取值500ⅹ10-3µm2，气体粘度取值5 mPa·s，根据式子（A-11）、（A-12）计算得到蒸汽Y方向的瞬时渗流速度为763.534 µm2/s。
S5，计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度。由于是倾斜油藏，S2和S3计算的各种力在油藏水平方向均有分力，且作用方向不一致，因此需要求取S2和S3在水平方向的合力，并根据合力和渗流理论，计算得到油藏水平方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响，让水平方向瞬时渗流速度计算更加精确。
联立式子（A-2）、（A-5）、（A-7）、（A-9），得到X方向合力计算方程：
Figure 534037DEST_PATH_IMAGE069
（A-13）
将（A-2）、（A-5）、（A-5）计算结果代入（A-13），得到：
Figure 269911DEST_PATH_IMAGE070
（A-14）
根据达西定律，油滴在浮力及毛管力的共同作用下，X方向瞬时渗流速度为：
Figure 734391DEST_PATH_IMAGE071
（A-15）
根据表1的数据，水平方向渗透率取值2000ⅹ10-3µm2，气体粘度取值5 mPa·s，根据式子（A-13）、（A-14）计算得到蒸汽X方向的瞬时渗流速度计算方程为：
Figure 606532DEST_PATH_IMAGE072
（A-16）
S6，计算垂向分流系数，并绘制垂向分流系数与井间距离的图版。其特点在于，根据S4和S5计算得到的垂直方向和水平方向瞬时渗流速度的计算结果，可以通过计算垂向分流系数，定量计算出注入蒸汽在垂向上的突进情况，解决了海上油田由于没有观察井，难以掌握注入蒸汽油层内分布情况的难题。
垂向分流系数为垂向瞬时渗流速度与总瞬时渗流速度的比值，计算方程为（A-16）；以
Figure 334316DEST_PATH_IMAGE073
为纵坐标，井间距离L为横坐标，绘制垂向分流系数与井间距离L的图版（图3）。
Figure 58952DEST_PATH_IMAGE022
（A-17）
根据表1的数据和A区块的实际情况，
Figure 264806DEST_PATH_IMAGE074
为生产井产量，取值34722.22cm3 /s（300m3 /d）；
Figure 991453DEST_PATH_IMAGE075
为蒸汽注入量，取值34722.22cm3 /s（300m3 /d）；
Figure 890139DEST_PATH_IMAGE076
为油层厚度，取值40m；L 为注采井间距离，取值15000cm；K 为渗透率，取值
Figure 335027DEST_PATH_IMAGE077
；
Figure 344571DEST_PATH_IMAGE078
为蒸汽粘度，取值5mPa·s。将上述取值代入（A-17），计算方程变为：
Figure 456884DEST_PATH_IMAGE079
（A-18）
S7，计算注入蒸汽的超覆程度，并绘制超覆程度与井间距离的图版。S6计算的垂向分流量值是一个变化的量，在注采井间的不同位置具有不同的值，通过S7，可以描绘出注采井间不同位置处的蒸汽分布情况。其特点在于，蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值，用 D表示（表2）。具体计算方法为：将注采井间的距离划分为N块，当蒸汽前缘运移到注采井间第 i 块( 0 ＜ i ＜ N) 时，利用此处的径向蒸汽运移量与垂向流量比的乘积和第i-1块的超覆程度Di-1迭代，并用（A-19）计算出此时的蒸汽超覆程度Di。并以
Figure 260892DEST_PATH_IMAGE080
为纵坐标，井间距离为横坐标，绘制超覆程度
Figure 957190DEST_PATH_IMAGE081
与注采井间距离L的图版（图4）。其中，
Figure 301584DEST_PATH_IMAGE027
表示第
Figure 268403DEST_PATH_IMAGE028
块的垂向分流量。
Figure 977733DEST_PATH_IMAGE082
（A-19）
表2 A目标区注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表
Figure 928371DEST_PATH_IMAGE083
S8，定量计算倾斜厚层稠油注蒸汽纵向波及系数。S7计算的蒸汽超覆程度Di，在注采井间的不同位置具有不同的值，通过S8，可以描绘出注采井间任何剖面的蒸汽分布面积，根据面积加权，可以定量计算出注采井间的蒸汽波及厚度。这样得到的波及厚度和波及系数值，考虑了注采井间不同位置的超覆程度差异，使得计算结果更加贴近矿场实际情况。
根据S7的图版，用（A-19）计算厚层稠油蒸汽吞吐初期不同分块的蒸汽超覆程度（表2第5列）；通过（1-Di）计算每个区块纵向波及无因次厚度hi（表2第6列）；用每个区块的长度（10m）乘以无因次厚度1得到每个区块的面积Ai（表2第7列）；用Ai hi计算得到每个区块的每个区块面积厚度乘积（表2第8列）。
用（A-20）、（A-21），分别计算得到面积加权后的纵向波及厚度
Figure 14139DEST_PATH_IMAGE084
和纵向波及系数Vk。
Figure 835464DEST_PATH_IMAGE085
（A-20）
Figure 43592DEST_PATH_IMAGE086
（A-21）
根据表2的数据，根据（A-20）计算得到面积加权后的无因次厚度为65.3（表2第8列求和）/150（表2第7列求和）=0.435；无因次厚度h取值为1，根据（A-20）计算得到纵向厚度波及系数为43.5%。
（实施例2）
一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，对某厚层倾斜稠油油藏B进行纵向波及系数预测，过程如下：
S1，收集目标区块地质油藏基本参数等静态资料。收集目标区块的基础数据，包括油层厚度、地层原油粘度、渗透率等参数（表3）。表3中收集的为B目标区的基础数据。
表3 B目标区的基础数据
Figure 950368DEST_PATH_IMAGE037
S2，计算目标区块蒸汽微元在油层垂直方向上的受力值。对于蒸汽微元，受到的重力为：
Figure 839826DEST_PATH_IMAGE087
（A-1）
式中：
Figure 282702DEST_PATH_IMAGE088
为蒸汽的密度，g /cm3；
Figure 599414DEST_PATH_IMAGE089
为蒸汽微元的体积，cm3；g 为重力加速度，m/s2。
由于地层具有一定倾角
Figure 993486DEST_PATH_IMAGE090
，以油层水平方向为X轴，油层垂直方向为Y轴，建立直角坐标系统，计算重力在X方向和Y方向的分量：
Figure 421057DEST_PATH_IMAGE091
（A-2）
Figure 216974DEST_PATH_IMAGE092
（A-3）
在注入蒸汽突破以前，蒸汽微元受到的浮力可以近似看作是恒定的，方向为垂直向上，表达式为：
Figure 970167DEST_PATH_IMAGE093
（A-4）
式中：
Figure 851535DEST_PATH_IMAGE094
为油相的密度，g /cm3；
Figure 82796DEST_PATH_IMAGE095
为水相的密度，g /cm3；
Figure 231756DEST_PATH_IMAGE096
为蒸汽微元的体积，cm3；g 为重力加速度，m/s2。
由于地层具有一定倾角
Figure 952587DEST_PATH_IMAGE097
，以油层水平方向为X轴，油层垂直方向为Y轴，建立直角坐标系统，计算浮力在X方向和Y方向的分量：
Figure 321251DEST_PATH_IMAGE098
（A-5）
Figure 356204DEST_PATH_IMAGE099
（A-6）
根据表3的数据，油的密度取值0.95g /cm3，水的密度取值1.0g /cm3，重力加速度取值9.8 m/s2，地层倾角取值20°，蒸汽微元的体积取值1.0cm3。根据式子（A-2）、（A-3）、（A-4）、（A-5）计算得到单位体积（1.0cm3，）蒸汽所受到的X方向重力和浮力分别为：-3.184mN和-0.168mN；Y方向的重力和浮力分别为8.749mN和-0.460mN。
S3，计算目标区块蒸汽微元在油层水平方向上的受力值。主要包括驱替力和毛管力，用（A-7）计算注采井间主流线上的注采压力梯度。
Figure 861134DEST_PATH_IMAGE100
（A-7）
式中：
Figure 956129DEST_PATH_IMAGE101
为生产井产量，cm3 /s；
Figure 812090DEST_PATH_IMAGE102
为蒸汽注入量，cm3 /s；
Figure 119574DEST_PATH_IMAGE103
为油层厚度，m；r 为任意点到生产井的距离，cm；
Figure 744591DEST_PATH_IMAGE104
为 r 处的压力梯度，10-1MPa/cm； L 为注采井间距离，cm；K 为渗透率，
Figure 523671DEST_PATH_IMAGE105
；
Figure 601348DEST_PATH_IMAGE106
为原油粘度，mPa·s。
对于蒸汽，受到的毛管力为：
Figure 243682DEST_PATH_IMAGE107
（A-8）
式中：
Figure 519942DEST_PATH_IMAGE108
为界面张力，mN/m；
Figure 956740DEST_PATH_IMAGE109
为接触角，°；
Figure 521714DEST_PATH_IMAGE110
为孔隙半径，mm。
在水湿情况下，毛细管力是注入水运移的动力，取正值；在油湿情况下，毛管力是注入水运移的阻力，取负值。毛管力的方向与注入水运移方向同向或反向，为便于计算与研究，分解为平行油藏方向的水平毛管力
Figure 436580DEST_PATH_IMAGE111
和垂直油藏方向的垂直毛管力
Figure 36189DEST_PATH_IMAGE112
。
Figure 643887DEST_PATH_IMAGE113
（A-9）
Figure 194692DEST_PATH_IMAGE114
（A-10）
根据表3的数据，界面张力取值0.02 N/m，接触角取值100°，孔隙半径取值0.00005m。由于接触角较大，可以判定为油湿，毛管力为注入水运移的阻力，取负值。根据式子（A-9）、（A-10）计算得到蒸汽所受到的X和Y方向的毛管力分别为0.001Pa和-0.001Pa。
S4，计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度。联立式子（A-3）、（A-6）、（A-10），得到Y方向合力计算方程：
Figure 913249DEST_PATH_IMAGE115
（A-11）
根据表3的数据，联立（A-3）、（A-6）、（A-10），计算得到
Figure 367365DEST_PATH_IMAGE116
为8.287Pa。
注入蒸汽沿油层斜上方渗流，可将速度分解为沿油层X方向的速度和Y方向的速度。根据达西定律，蒸汽在重力、浮力及毛管力的共同作用下，Y方向瞬时渗流速度为：
Figure 145965DEST_PATH_IMAGE117
（A-12）
根据表3的数据，垂直方向渗透率取值600ⅹ10-3µm2，气体粘度取值5 mPa·s，根据式子（A-11）、（A-12）计算得到蒸汽Y方向的瞬时渗流速度为994.40 µm2/s。
S5，计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度。联立式子（A-2）、（A-5）、（A-7）、（A-9），得到X方向合力计算方程：
Figure 685530DEST_PATH_IMAGE118
（A-13）
将（A-2）、（A-5）、（A-5）计算结果代入（A-13），得到：
Figure 942199DEST_PATH_IMAGE119
（A-14）
根据达西定律，油滴在浮力及毛管力的共同作用下，X方向瞬时渗流速度为：
Figure 250821DEST_PATH_IMAGE120
（A-15）
根据表3的数据，水平方向渗透率取值3000ⅹ10-3µm2，气体粘度取值5 mPa·s，根据式子（A-13）、（A-14）计算得到蒸汽X方向的瞬时渗流速度计算方程为：
Figure 200323DEST_PATH_IMAGE121
（A-16）
S6，计算垂向分流系数，并绘制垂向分流系数与井间距离的图版。垂向分流系数为垂向瞬时渗流速度与总瞬时渗流速度的比值，计算方程为（A-16）；以
Figure 463070DEST_PATH_IMAGE122
为纵坐标，井间距离L为横坐标，绘制垂向分流系数与井间距离L的图版（图5）。
Figure 789009DEST_PATH_IMAGE022
（A-17）
根据表3的数据和A区块的实际情况，
Figure 748875DEST_PATH_IMAGE123
为生产井产量，取值34722.22cm3 /s（300m3 /d）；
Figure 869278DEST_PATH_IMAGE124
为蒸汽注入量，取值34722.22cm3 /s（300m3 /d）；
Figure 117856DEST_PATH_IMAGE125
为油层厚度，取值50m；L为注采井间距离，取值20000cm；K 为渗透率，取值
Figure 981907DEST_PATH_IMAGE126
；
Figure 999542DEST_PATH_IMAGE127
为蒸汽粘度，取值5mPa·s。将上述取值代入（A-17），计算方程变为：
Figure 556425DEST_PATH_IMAGE128
（A-18）
S7，计算注入蒸汽的超覆程度，并绘制超覆程度与井间距离的图版。蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值，用 D表示（表4）。将注采井间距离划分为N块，当蒸汽前缘运移到注采井间 i 块( 0 ＜ i ＜ N) 时，利用此处的径向蒸汽运移量与垂向流量比的乘积和第i-1块的超覆程度迭代，用（A-19）计算出此时的蒸汽超覆程度。并以
Figure 292300DEST_PATH_IMAGE129
为纵坐标，井间距离为横坐标，绘制超覆程度
Figure 491200DEST_PATH_IMAGE129
与注采井间距离L的图版（图6）。其中，
Figure 127455DEST_PATH_IMAGE027
表示第
Figure 855240DEST_PATH_IMAGE028
块的垂向分流量；。
Figure 547252DEST_PATH_IMAGE130
（A-19）
表4 B目标区块注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表
Figure 18685DEST_PATH_IMAGE131
S8，定量计算倾斜厚层稠油注蒸汽纵向波及系数。根据S7的图版，用（A-19）计算厚层稠油蒸汽吞吐初期不同分块的蒸汽超覆程度（表4第5列）；通过
Figure 10912DEST_PATH_IMAGE132
计算每个区块纵向波及无因次厚度
Figure 644019DEST_PATH_IMAGE133
（表4第6列）；用每个区块的长度（10m）乘以无因次厚度1得到每个区块的面积
Figure 354486DEST_PATH_IMAGE134
（表4第7列）；用
Figure 160768DEST_PATH_IMAGE135
计算得到每个区块的每个区块面积厚度乘积（表4第8列）。
用（A-20）、（A-21），分别计算得到面积加权后的纵向波及厚度
Figure 7501DEST_PATH_IMAGE136
和纵向波及系数
Figure 578553DEST_PATH_IMAGE137
。
Figure 510737DEST_PATH_IMAGE138
（A-20）
Figure 323972DEST_PATH_IMAGE139
（A-21）
根据表4的数据，根据（A-20）计算得到面积加权后的无因次厚度为36.5（表4第8列求和）/200（表4第7列求和）=0.1825；无因次厚度h取值为1，根据（A-20）计算得到纵向厚度波及系数为18.25%。
上述实施例仅表达了较为优选的实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。

Claims (9)

1.一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，其特征在于：包括以下步骤：
S1、收集目标区块地质油藏基本参数；
S2、计算目标区块蒸汽微元在油层垂直方向上的受力值；
S3、计算目标区块蒸汽微元在油层水平方向上的受力值；
S4、计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度；
S5、计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度；
S6、计算垂向分流系数，并绘制垂向分流系数与井间距离的图版；
S7、计算注入蒸汽的超覆程度，并绘制超覆程度与井间距离的图版；蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值，用D表示；
S8、定量计算倾斜厚层稠油注蒸汽纵向波及系数。
2.根据权利要求1所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，其特征在于：所述S1中的基本参数：
是指收集目标区块的基础数据，包括油层厚度、地层倾角、注采井间距离、油的密度、水的密度、气体粘度、界面张力、接触角、孔隙半径、垂直方向渗透率、水平方向渗透率。
3.根据权利要求2所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，其特征在于：所述S2中计算目标区块蒸汽微元在油层垂直方向上的受力值，包括重力、浮力；
S21、重力的计算
根据蒸汽微元的密度为
Figure 797767DEST_PATH_IMAGE001
、蒸汽微元的体积
Figure 978737DEST_PATH_IMAGE002
，再根据重力加速度g，计算得蒸汽微元的重力F重力；
由于地层有倾角
Figure 943150DEST_PATH_IMAGE003
，以油层水平方向为X轴、油层垂直方向为Y轴，建立直角坐标系，计算蒸汽微元的重力在X方向的分量Fx重力、Y方向的分量FY重力；
S22、浮力的计算
在注入蒸汽突破之前，蒸汽微元受到的浮力F浮力看作是恒定的，方向为竖直向上；
浮力F浮力，根据油相的密度
Figure 486127DEST_PATH_IMAGE004
、水相的密度
Figure 168781DEST_PATH_IMAGE005
、蒸汽微元的体积
Figure 822617DEST_PATH_IMAGE006
，重力加速度g，计算得到；
由于地层有倾角
Figure 638607DEST_PATH_IMAGE007
，以油层水平方向为X轴、油层垂直方向为Y轴，建立直角坐标系，计算浮力F浮力在X方向的分量FX浮力、Y方向的分量FY浮力。
4.根据权利要求3所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，其特征在于：所述S3中计算目标区块蒸汽微元在油层水平方向上的受力值，包括毛管力和驱替压力：
S31、计算驱替压力，即计算注采井间主流线上的注采压力梯度，公式为
Figure 618065DEST_PATH_IMAGE008
其中，
Figure 460119DEST_PATH_IMAGE009
为生产井产量；
Figure 979962DEST_PATH_IMAGE010
为蒸汽注入量；
Figure 856651DEST_PATH_IMAGE011
为油层厚度；r为任意点到生产井的距离；
Figure 741430DEST_PATH_IMAGE012
为r处的压力梯度；L为注采井间距离；
Figure 136027DEST_PATH_IMAGE013
为渗透率；
Figure 397244DEST_PATH_IMAGE014
为原油粘度；
S32、计算毛管力
对于蒸汽，受到毛管力为
Figure 190756DEST_PATH_IMAGE015
式中，
Figure 246437DEST_PATH_IMAGE016
为界面张力；
Figure 125400DEST_PATH_IMAGE017
为接触角；
Figure 924729DEST_PATH_IMAGE018
为孔隙半径；
在水湿情况下，毛细管力是注入水运移的动力，取正值；在油湿情况下，毛管力是注入水运移的阻力，取负值；
由于地层有倾角
Figure 507502DEST_PATH_IMAGE019
，根据该倾角，计算平行油藏方向的水平毛管力
Figure 999663DEST_PATH_IMAGE020
和垂直油藏方向的垂直毛管力
Figure 365922DEST_PATH_IMAGE021
。
5.根据权利要求4所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，其特征在于：所述步骤S4计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度；
S41、计算目标区块蒸汽微元在油藏垂直方向的合力
根据蒸汽微元在Y方向的重力分量FY重力，蒸汽微元在Y方向的重力分量FY浮力，垂直油藏方向的垂直毛细管力FY毛管力，计算得到蒸汽微元在油藏垂直方向的合力FY合力；
S42、计算目标区块蒸汽微元在油藏垂直方向的瞬时渗流速度
注入蒸汽沿油层斜上方渗流，将速度分解为沿油层X方向的速度和Y方向的速度；
蒸汽在重力、浮力及毛管力的共同作用下，根据达西定律，计算得到Y方向瞬时渗流速度，即在油藏垂直方向的瞬时渗流速度VX。
6.根据权利要求5所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，其特征在于：所述S5计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度；
S51、计算目标区块蒸汽微元在油藏水平方向的合力
根据蒸汽微元在X方向的重力分量FX重力，蒸汽微元在X方向的重力分量FX浮力，驱替压力，垂直油藏方向的垂直毛细管力FX毛管力，计算得到蒸汽微元在油藏垂直方向的合力FX合力；
S52、计算目标区块蒸汽微元在油藏水平方向的瞬时渗流速度
根据达西定律，油滴在浮力及毛管力的共同作用下，计算得到X方向瞬时渗流速度为VY。
7.根据权利要求6所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，其特征在于：所述S6计算垂向分流系数，并绘制垂向分流系数与井间距离的图版；
S61、计算垂向分流系数
垂向分流系数为垂向瞬时渗流速度与总瞬时渗流速度的比值，计算方程为
Figure 703362DEST_PATH_IMAGE022
S62、
Figure 409150DEST_PATH_IMAGE023
为纵坐标，井间距离L为横坐标，绘制垂向分流系数与井间距离L的图版。
8.根据权利要求7所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，其特征在于：所述S7计算注入蒸汽的超覆程度，并绘制超覆程度与井间距离的图版；
S71、计算注入蒸汽的超覆程度
蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值，用D表示；
将注采井间距离划分为N块，当蒸汽前缘运移到注采井间 i 块( 0 ＜ i ＜ N) 时,利用此处的径向蒸汽运移量与垂向流量比的乘积和第i-1块的超覆程度迭代，计算出此时的蒸汽超覆程度，计算公式为
Figure 868950DEST_PATH_IMAGE024
，
其中，
Figure 663119DEST_PATH_IMAGE025
表示第
Figure 804250DEST_PATH_IMAGE026
块的垂向分流量；
S72、绘制超覆程度与井间距离的图版
以
Figure 364544DEST_PATH_IMAGE027
为纵坐标，井间距离为横坐标，绘制超覆程度
Figure 198508DEST_PATH_IMAGE028
与注采井间距离L的图版。
9.根据权利要求8所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，其特征在于：所述S8定量计算倾斜厚层稠油注蒸汽纵向波及系数；
S81、根据S7的图版，计算得到厚层稠油蒸汽吞吐初期不同分块的蒸汽超覆程度，计算公式为
Figure 211464DEST_PATH_IMAGE024
；
S82、通过
Figure 156286DEST_PATH_IMAGE029
，计算每个区块纵向波及无因次厚度
Figure 571087DEST_PATH_IMAGE030
；
S83、用每个区块的长度乘以无因次厚度
Figure 575952DEST_PATH_IMAGE031
得到每个区块的面积
Figure 73273DEST_PATH_IMAGE032
；
S84、用
Figure 556207DEST_PATH_IMAGE033
计算得到每个区块的每个区块面积厚度乘积；
S85、计算得到面积加权后的纵向波及厚度
Figure 91094DEST_PATH_IMAGE034
，计算公式为
Figure 266860DEST_PATH_IMAGE035
；
计算得到面积加权后的纵向波及系数V k ，计算公式为
Figure 254408DEST_PATH_IMAGE036
。
CN202211385246.6A 2022-11-07 2022-11-07 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 Active CN115422859B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211385246.6A CN115422859B (zh) 2022-11-07 2022-11-07 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211385246.6A CN115422859B (zh) 2022-11-07 2022-11-07 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法

Publications (2)

Publication Number Publication Date
CN115422859A true CN115422859A (zh) 2022-12-02
CN115422859B CN115422859B (zh) 2023-01-24

Family

ID=84207638

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211385246.6A Active CN115422859B (zh) 2022-11-07 2022-11-07 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法

Country Status (1)

Country Link
CN (1) CN115422859B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115859530A (zh) * 2023-01-18 2023-03-28 西南石油大学 一种基于无因次泵效的抽油泵健康状态评价方法
CN116341406A (zh) * 2023-03-03 2023-06-27 常州大学 一种确定倾斜油藏注采干扰时机的计算方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4637466A (en) * 1986-04-03 1987-01-20 Texaco Inc. Method of improving conformance in steam floods with carboxylate steam foaming agents
US5214384A (en) * 1991-07-24 1993-05-25 Mobil Oil Corporation Method including electrical self potential measurements for detecting multiphase flow in a cased hole
CN104318052A (zh) * 2014-09-16 2015-01-28 中国石油大学(北京) 一种低渗透油藏水驱波及系数的评价方法
US20150198022A1 (en) * 2014-01-13 2015-07-16 Conocophillips Company Oil recovery with fishbone wells and steam
CN106150466A (zh) * 2016-08-23 2016-11-23 中国石油化工股份有限公司 凝胶泡沫抑制底水锥进的稠油热采方法
CN107435535A (zh) * 2016-05-26 2017-12-05 中国石油大学(北京) 一种采用平面重力驱开采高倾角稠油油藏的方法
CN112036033A (zh) * 2020-08-31 2020-12-04 中国海洋石油集团有限公司 一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法
US20210131244A1 (en) * 2018-05-14 2021-05-06 China University Of Petroleum Super heavy oil development method for strengthening sagd steam chamber to break through low physical property reservoir

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4637466A (en) * 1986-04-03 1987-01-20 Texaco Inc. Method of improving conformance in steam floods with carboxylate steam foaming agents
US5214384A (en) * 1991-07-24 1993-05-25 Mobil Oil Corporation Method including electrical self potential measurements for detecting multiphase flow in a cased hole
US20150198022A1 (en) * 2014-01-13 2015-07-16 Conocophillips Company Oil recovery with fishbone wells and steam
CN104318052A (zh) * 2014-09-16 2015-01-28 中国石油大学(北京) 一种低渗透油藏水驱波及系数的评价方法
CN107435535A (zh) * 2016-05-26 2017-12-05 中国石油大学(北京) 一种采用平面重力驱开采高倾角稠油油藏的方法
CN106150466A (zh) * 2016-08-23 2016-11-23 中国石油化工股份有限公司 凝胶泡沫抑制底水锥进的稠油热采方法
US20210131244A1 (en) * 2018-05-14 2021-05-06 China University Of Petroleum Super heavy oil development method for strengthening sagd steam chamber to break through low physical property reservoir
CN112036033A (zh) * 2020-08-31 2020-12-04 中国海洋石油集团有限公司 一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DONG LIU等: ""A New Method for Calculating the Inflection Point Temperature of Heavy-Oil Rheology Transforming From Non-Newton Fluid into the Newton Fluid"", 《JOURNAL OF ENERGY RESOURCES TECHNOLOGY》 *
W.SHUHONG等: ""Sequential Multiwell Steam Huff n Puff in Heavy-Oil Development"", 《2005,SPE/PS-CIM/CHOA INTERNATIONAL THERMAL OPERATIONS AND HEAVY OIL SYMPOSIUM》 *
周阴国: "蒸汽-泡沫复合驱提高稠油采收率机理研究", 《中国石油和化工标准与质量》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115859530A (zh) * 2023-01-18 2023-03-28 西南石油大学 一种基于无因次泵效的抽油泵健康状态评价方法
CN115859530B (zh) * 2023-01-18 2023-08-11 西南石油大学 一种基于无因次泵效的抽油泵健康状态评价方法
CN116341406A (zh) * 2023-03-03 2023-06-27 常州大学 一种确定倾斜油藏注采干扰时机的计算方法
CN116341406B (zh) * 2023-03-03 2024-04-23 常州大学 一种确定倾斜油藏注采干扰时间的计算方法

Also Published As

Publication number Publication date
CN115422859B (zh) 2023-01-24

Similar Documents

Publication Publication Date Title
CN115422859B (zh) 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法
CN105822302A (zh) 一种基于井地电位法的油水分布识别方法
Su et al. A fractal-monte carlo approach to model oil and water two-phase seepage in low-permeability reservoirs with rough surfaces
CN109884269A (zh) 岩心自发渗吸采收率预测方法及系统
CN104989341A (zh) 一种确定低渗透油藏有效驱替注采井距的方法
Silin et al. A model of buoyancy-driven two-phase countercurrent fluid flow
CN104975827A (zh) 预测二氧化碳驱油藏指标的物质平衡方法
Huang et al. Effects of water invasion law on gas wells in high temperature and high pressure gas reservoir with a large accumulation of water-soluble gas
CN115587674A (zh) 油藏改建储气库扩容达产过程气井动态产能预测方法
Zhao et al. A new dynamic capillary-number-recovery evaluation method for tight gas reservoirs
CN105298478A (zh) 一种用于断裂构造的地层孔隙压力确定方法
Dimri et al. Reservoir geophysics: Some basic concepts
Masalmeh et al. The importance of special core analysis in modelling remaining oil saturation in carbonate fields
Wood et al. Determining effective residual oil saturation for mixed wettability reservoirs: Endicott Field, Alaska
Qi et al. Mathematical model for prediction of dynamic reserves loss due to water invasion in water-drive gas reservoir
CN107766639A (zh) 基于压力降低系数的天然气侧向运移最大距离的计算方法
CN114676652A (zh) 一种欠饱和储层煤层气井生产边界定量预测的计算方法
Fu et al. Predicted model of relative permeability considering water distribution characteristics in tight sandstone gas reservoirs
Xu et al. Experimental research of remaining oil distribution for inverted nine-sport pattern and infill well pattern
Ran et al. Experiments and effective permeability model for multiphase flow in rock fractures with variable apertures
Gambolati Numerical models in land subsidence control
Awag et al. Comparison between downdip and updip groundwater flow on early CO2 migration in dipping storage aquifers
Yoneda et al. Deformation of seabed due to exploitation of methane hydrate reservoir
CN114136838B (zh) 确定海上稠油不同含水阶段注热水驱原油粘度界限的方法
Azin et al. Basics of Oil and Gas Flow in Reservoirs

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant