CN113803054A - 油水界面深度确定方法及防止生产井水窜的预警方法 - Google Patents
油水界面深度确定方法及防止生产井水窜的预警方法 Download PDFInfo
- Publication number
- CN113803054A CN113803054A CN202010531436.9A CN202010531436A CN113803054A CN 113803054 A CN113803054 A CN 113803054A CN 202010531436 A CN202010531436 A CN 202010531436A CN 113803054 A CN113803054 A CN 113803054A
- Authority
- CN
- China
- Prior art keywords
- water
- production well
- reservoir
- oil
- water cone
- 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
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 358
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 226
- 238000000034 method Methods 0.000 title claims abstract description 69
- 230000005465 channeling Effects 0.000 title claims abstract description 11
- 239000003921 oil Substances 0.000 claims description 49
- 239000010779 crude oil Substances 0.000 claims description 28
- 230000005484 gravity Effects 0.000 claims description 13
- 230000001133 acceleration Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 4
- 238000005553 drilling Methods 0.000 description 10
- 239000012530 fluid Substances 0.000 description 9
- 238000011161 development Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000009472 formulation Methods 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000012407 engineering method Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000009916 joint effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 210000005239 tubule Anatomy 0.000 description 1
- 235000020681 well water Nutrition 0.000 description 1
- 239000002349 well water Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/04—Measuring depth or liquid level
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/10—Locating fluid leaks, intrusions or movements
-
- 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
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Theoretical Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Pipeline Systems (AREA)
Abstract
本发明提供一种油水界面深度确定方法及防止生产井水窜的预警方法。油水界面深度确定方法包括：基于水锥界面的压力平衡、裂缝型生产井的产油量方程、水锥高度和储集体厚度以及含水率三者之间的关系、裂缝宽度与储集体厚度之间的关系，构建裂缝型储集体顶部到裂缝型生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型；检测裂缝型生产井的含水率和储集体顶部距离地表的深度；根据测得的含水率，利用所述关系模型，确定储集体顶部到该生产井下方的水锥顶部之间的距离；根据储集体顶部距离地表的深度和储集体顶部到该生产井下方的水锥顶部之间的距离，确定该生产井的油水界面深度。该方法能够准确确定油水界面深度。
Description
技术领域
本发明属于缝洞型油藏开发技术领域，具体涉及一种钻遇裂缝型储集体生产井动态油水界面的计算方法。
背景技术
缝洞型油藏大部分生产井钻遇奥陶系顶部完钻，缺少储层完整的测井数据及产液剖面测试数据，难以准确确定动态油水界面深度，而动态油水界面深度的确定对于生产井防止水窜措施的制定和开发调整措施的制定具有重要意义。
目前试图确定动态油水界面的方法主要包括测井法和经验公式法两大类。
测井法：基于储层完整的饱和度测井数据，利用含水饱和度来判别油水界面的方法。如：郭振彬，曹文丽，吴剑峰等利用测井资料和观察井油水界面资料综合确定了裂缝性油藏动态油水界面(郭振彬，曹文丽，吴剑峰.一种确定裂缝性油藏动态油水界面的综合方法.测井技术,1997,21(4):269～271)。
基于目的层深度的经验公式法：根据动静态资料确定生产井的动态油水界面，然后拟合确定动态油水界面的经验公式。如：郭分乔等综合利用静态和单井动态资料综合分析缝洞型碳酸盐岩油藏油水界面，并拟合得到了油水界面深度与目的层顶面深度的经验公式(郭分乔，宋付英，李远钦等.塔河缝洞型碳酸盐岩油藏油水界面确定方法探讨.内蒙古石油化工,2008,15:74～76)。
基于压力的经验公式法：运用原始地层压力和油层中部深度压力参数，推导地层油水界面深度的计算公式。如闫晓芳，邹伟宏，陈戈等通过对塔里木油田轮古13区块单井资料的分析，利用各井动静态资料确定了油水界面位置与油层压力的经验关系式(闫晓芳，邹伟宏，陈戈等.碳酸盐岩缝洞型油藏油水界面计算方法—以塔里木油田轮古15区块为例.石油地质与工程,2012,26(5):67～69)。
现有技术存在的问题主要是：利用测井法计算动态油水界面深度，需要完整的饱和度测井数据，而实际生产过程中很少进行测井，现场产液剖面数据也比较稀缺，因此通过饱和度测井和产液剖面测试难以对所有钻遇裂缝型生产井的动态油水界面深度进行确定。而经验公式法只适用于某些特定区块，对未知区块则难以准确确定动态油水界面深度，而且经验公式中的某些特定的参数难以准确确定，从而影响了经验公式法的推广应用。
因此，现有的测井法和经验公式法，很难准确确定缝洞型油藏钻遇裂缝型储集体生产井动态油水界面的深度。有必要探索一种针对缝洞型油藏裂缝型生产井动态油水界面深度的确定方法。
发明内容
本发明提供一种油水界面深度确定方法及防止生产井水窜的预警方法，以准确确定裂缝型生产井的油水界面深度。
第一方面，本申请的实施方式提供一种油水界面深度确定方法，应用于裂缝型储集体生产井，包括以下步骤：基于水锥界面的压力平衡、裂缝型生产井的产油量方程、水锥高度和储集体厚度以及含水率三者之间的关系、裂缝宽度与储集体厚度之间的关系，构建裂缝型储集体顶部到裂缝型生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型；检测裂缝型生产井的含水率和储集体顶部距离地表的深度；根据测得的含水率，利用所述关系模型，确定储集体顶部到该生产井下方的水锥顶部之间的距离；根据储集体顶部距离地表的深度和储集体顶部到该生产井下方的水锥顶部之间的距离，确定该生产井的油水界面深度。
在一个实施例中，基于水锥界面的压力平衡、裂缝型生产井的产油量方程、水锥高度和储集体厚度以及含水率三者之间的关系、裂缝宽度与储集体厚度之间的关系，构建裂缝型储集体顶部到裂缝型生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型，包括：基于水锥界面的压力平衡确定裂缝型生产井下方的水锥界面的压力平衡方程，确定裂缝型生产井的产油量方程，通过联立裂缝型生产井下方的水锥界面的压力平衡方程和该生产井的产油量方程建立该生产井下方与储集体厚度相关的水锥界面到裂缝边缘的垂向水锥高度计算式；确定裂缝型生产井下方的水锥高度、储集体厚度以及该生产井的含水率三者之间的第一函数关系；确定裂缝宽度与储集体厚度之间的第二函数关系；根据所述水锥界面到裂缝边缘的垂向水锥高度计算式、所述第一函数关系和所述第二函数关系确定该生产井下方的水锥高度与所述含水率之间的第三函数关系；确定该生产井下方的水锥高度、储集体厚度以及储集体顶部到该生产井下方的水锥顶部之间的距离的第四函数关系；根据所述第一函数关系和所述第四函数关系，确定所述水锥高度、所述含水率以及储集体顶部到该生产井下方的水锥顶部之间的距离的第五函数关系；根据所述第三函数关系和所述第五函数关系构建储集体顶部到该生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型。
在一个实施例中，裂缝型生产井的产油量方程为：
其中，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，ω为裂缝宽度，f(x)为水锥界面到裂缝边缘的垂向水锥高度，p为距离井筒中心x处的裂缝内水锥界面的压力，L为裂缝长度方向的距离。
在一个实施例中，该生产井下方与储集体厚度相关的水锥界面到裂缝边缘的垂向水锥高度计算式为：
其中，f(x)为水锥界面到裂缝边缘的垂向水锥高度，ω为裂缝宽度，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，Δρow表示油水密度差，g为重力加速度，单位为m/s2，θ为裂缝方向与竖直方向的夹角，Lf为裂缝长度。
在一个实施例中，裂缝型生产井下方的水锥高度、储集体厚度以及该生产井的含水率三者之间的第一函数关系为：
其中，fw为裂缝型生产井的含水率，h表示裂缝型生产井下方的水锥高度，H表示储集体厚度。
在一个实施例中，裂缝宽度与储集体厚度之间的第二函数关系为：
ω＝H·cosθ
其中，ω为裂缝宽度，H表示储集体厚度，θ为裂缝方向与竖直方向的夹角。
在一个实施例中，该生产井下方的水锥高度与所述含水率之间的第三函数关系为：
其中，h表示裂缝型生产井下方的水锥高度，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，Δρow表示油水密度差，g为重力加速度，单位为m/s2，Lf为裂缝长度，θ为裂缝方向与竖直方向的夹角，rw为井筒半径，fw为裂缝型生产井的含水率。
在一个实施例中，该生产井下方的水锥高度、储集体厚度以及储集体顶部到该生产井下方的水锥顶部之间的距离的第四函数关系为：
how＝H-h
其中，how表示所述储集体顶部到该生产井下方的水锥顶部之间的距离，h表示裂缝型生产井下方的水锥高度，H表示储集体厚度。
在一个实施例中，所述水锥高度、所述含水率以及储集体顶部到该生产井下方的水锥顶部之间的距离的第五函数关系为：
其中，how表示所述储集体顶部到该生产井下方的水锥顶部之间的距离，h表示裂缝型生产井下方的水锥高度，fw为裂缝型生产井的含水率。
在一个实施例中，储集体顶部到该生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型为：
其中，how表示储集体顶部到该生产井下方的水锥顶部之间的距离，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，Δρow表示油水密度差，g为重力加速度，单位为m/s2，Lf为裂缝长度，θ为裂缝方向与竖直方向的夹角，rw为井筒半径，fw为裂缝型生产井的含水率。
第二方面，本申请的实施方式提供一种防止生产井水窜的预警方法，包括以下步骤：检测生产井钻遇的储集体的类型；当钻遇的储集体为裂缝型储集体时，采用如上文所述的油水界面深度确定方法来确定生产井的动态油水界面深度；根据生产井的动态油水界面深度确定生产井下的水锥位置；根据生产井下的水锥位置判断生产井发生水淹事故的风险，并由此判断是否进行水淹预警。
第三方面，本申请的实施方式提供一种存储介质，存储有计算机程序，所述计算机程序被处理器执行时，实现如上文所述的油水界面深度确定方法的步骤或防止生产井水窜的预警方法的步骤。
第四方面，本申请的实施方式提供一种电子设备，包括处理器和存储有程序代码的存储介质，所述程序代码被所述处理器执行时，实现如上文所述的油水界面深度确定方法的步骤或防止生产井水窜的预警方法的步骤。
本发明针对现有油水界面确定方法难以准确确定动态油水界面深度的问题，提出一种针对缝洞型油藏钻遇裂缝类储集体生产井的动态油水界面深度的确定方法。该方法基于塔河缝洞型油藏生产井含水率上升的主要原因为井控范围内出水储集体的不断增加的认识，将生产井的含水率表示为水锥高度与储集体总厚度的比值，考虑裂缝中的流体流动特征及重力的影响，基于长直细管流动模型理论，建立了钻遇裂缝型储集体生产井的动态油水界面的确定方法。本发明利用生产井日常的产油产水数据、流体参数及裂缝物性参数，基于油藏工程方法确定的水锥高度和含水率的表达式，求解得到了裂缝型生产井动态油水界面深度的表达式，能够方便快捷的确定裂缝型生产井的动态油水界面深度。
本发明的方法能够克服常规方法无法准确确定动态油水界面深度的缺陷，快速有效地确定缝洞型油藏裂缝型生产井目前的动态油水界面深度，为缝洞型油藏钻遇裂缝型储集体生产井的水窜预警和开发调整方案的制定提供重要的参考依据。
附图说明
构成本申请的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定，在附图中：
图1为根据本申请一实施方式的油水界面深度确定方法的流程图；
图2为根据本申请一实施方式的构建裂缝型储集体顶部到裂缝型生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型的流程图；
图3为根据本申请一实施方式的裂缝型生产井下方水锥形态的示意图。
具体实施方式
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
实施例一
图1为根据本申请一实施方式的油水界面深度确定方法的流程图。如图1所示，本申请的实施方式提供一种油水界面深度确定方法，应用于裂缝型储集体生产井，包括以下步骤：
S100：基于水锥界面的压力平衡、裂缝型生产井的产油量方程、水锥高度和储集体厚度以及含水率三者之间的关系、裂缝宽度与储集体厚度之间的关系，构建裂缝型储集体顶部到裂缝型生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型。
图2为根据本申请一实施方式的构建裂缝型储集体顶部到裂缝型生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型的流程图。如图2所示，S100可以包括以下步骤：
S110：基于水锥界面的压力平衡确定裂缝型生产井下方的水锥界面的压力平衡方程，确定裂缝型生产井的产油量方程，通过联立裂缝型生产井下方的水锥界面的压力平衡方程和该生产井的产油量方程建立该生产井下方与储集体厚度相关的水锥界面到裂缝边缘的垂向水锥高度计算式。
基于水锥界面的压力平衡确定裂缝型生产井下方的水锥界面的压力平衡方程可以按照如下方式实施：
油井射孔后沟通高角度裂缝，流体经裂缝流入井底，流体流动为稳态N-S流，油、水流入裂缝后开始形成油-水两相界面。裂缝远井端压力恒定为pL，生产井井底压力为pwf，裂缝长度为Lf，裂缝宽度为ω，Lf>>ω，流体在裂缝中流动可以按无限长直细管流动进行简化。
裂缝型生产井下方的水锥形态如图3所示，其中，H为储集体厚度，h为裂缝型生产井下方的水锥高度，how为T7 4界面到水锥最高点的距离，Qo为产油量，裂缝中油的流速为v，裂缝与竖直方向的夹角为θ。在水锥界面上任取A、B两点：A点处，水锥半径(即距离井筒中心的水平距离)为x1，水锥界面到裂缝边缘的垂向水锥高度(即水锥界面距离裂缝边缘的垂向距离)为y1，原始油水界面(OWC，the Oil-Water Contact)上对应压力为PA，水锥界面上对应压力为P1；B点处，水锥半径为x2，水锥界面到裂缝边缘的垂向水锥高度为y2，原始油水界面上对应压力为PB，水锥界面上等高处对应压力为P2。
在裂缝内忽略毛管力的差异，根据静力学原理，有pA＝pB，从而得到裂缝型生产井下方的水锥界面的压力平衡方程如式(1)所示：
其中，ρo表示油的密度，ρw表示水的密度，g表示重力加速度。
将式(1)变形后得到：
其中，Δρow＝ρo-ρw，Δρow表示油水密度差。
将式(1)的两边同除以(x1-x2)，得到：
当(x1-x2)→0时，水锥界面方程f(x)＝y，则式(3)经微分变化得到：
其中，p为距离井筒中心的水平距离为x处的裂缝内油水界面的压力。
裂缝边界无滑移，裂缝内沿流体流动方向的压力梯度为：
将裂缝内沿流体流动方向的压力梯度分解后得到水平方向的压力梯度为：
依据哈根-泊肃叶定律，水锥形成时，裂缝型生产井的产油量方程可以表示为：
其中，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，ω为裂缝宽度，f(x)为水锥界面到裂缝边缘的垂向水锥高度，p为距离井筒中心x处的裂缝内水锥界面的压力，L为裂缝长度方向的距离。
将式(7)变形后得到：
联立式(4)、式(6)和式(8)得到：
在式(9)中代入边界条件：x＝Lf·sinθ，f(x)＝0，并对式9进行积分得到该生产井下方与储集体厚度相关的水锥界面到裂缝边缘的垂向水锥高度计算式为：
其中，f(x)为水锥界面到裂缝边缘的垂向水锥高度，ω为裂缝宽度，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，Δρow表示油水密度差，g为重力加速度，单位为m/s2，θ为裂缝方向与竖直方向的夹角，Lf为裂缝长度。
S120：确定裂缝型生产井下方的水锥高度、储集体厚度以及该生产井的含水率三者之间的第一函数关系。
将式(10)中的x代入井筒半径rw，得到裂缝型生产井下方与储集体厚度相关的水锥界面到裂缝边缘的垂向水锥高度表达式为：
其中，f(x)为水锥界面到裂缝边缘的垂向水锥高度，ω为裂缝宽度，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，Δρow表示油水密度差，g为重力加速度，单位为m/s2，θ为裂缝方向与竖直方向的夹角，Lf为裂缝长度，rw为井筒半径。
由于塔河缝洞型油藏高角度缝发育，不同部位储集体出水时间不同，生产井含水率的上升主要由井控范围内出水储集体的不断增加所导致的，当生产井钻遇裂缝型储集体时，含水率fw可以近似为井筒附近水锥高度h与储集体总厚度H的比值，即裂缝型生产井下方的水锥高度、储集体厚度以及该生产井的含水率三者之间的第一函数关系为：
其中，fw为裂缝型生产井的含水率，h表示裂缝型生产井下方的水锥高度，H表示储集体厚度。
变换式(12)，得到储集体厚度H的表达式为：
S130：确定裂缝宽度与储集体厚度之间的第二函数关系。
对于沟通井底的裂缝，裂缝宽度ω与储集体厚度H存在如下第二函数关系：
ω＝H·cosθ (14)
将式(13)和(14)代入式(11)，得到：
将式(15)依次变形，得到：
S140：根据所述水锥界面到裂缝边缘的垂向水锥高度计算式、所述第一函数关系和所述第二函数关系确定该生产井下方的水锥高度与所述含水率之间的第三函数关系。
将式(18)变形后，得到裂缝型生产井下方的水锥高度与所述含水率之间的第三函数关系为：
其中，h表示裂缝型生产井下方的水锥高度，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，Δρow表示油水密度差，g为重力加速度，单位为m/s2，Lf为裂缝长度，θ为裂缝方向与竖直方向的夹角，rw为井筒半径，fw为裂缝型生产井的含水率。
在实际钻井过程中，钻遇储集体顶部完钻，则无法了解到储集体厚度到底是多少。在本申请中，将储集体厚度与含水率联系起来，则可以通过生产井的含水率来避免直接求储集体厚度的问题，从而确定水锥高度。
S150：确定该生产井下方的水锥高度、储集体厚度以及储集体顶部到该生产井下方的水锥顶部之间的距离的第四函数关系。
根据裂缝型生产井下方的水锥高度、储集体厚度以及储集体顶部与该生产井下方的水锥顶部之间的位置关系，确定所述第四函数关系为：
how＝H-h (20)
其中，how表示所述储集体顶部到该生产井下方的水锥顶部之间的距离，h表示裂缝型生产井下方的水锥高度，H表示储集体厚度。
需要说明的是，在实际的裂缝型储集体中，存在大量的裂缝，水锥顶部指的是在所有裂缝的共同作用下所形成的水锥的顶部。而图3中仅示意性的画出了诸多裂缝中的一条裂缝，不能将图3中所画出的一条裂缝所形成的水锥的顶部作为在所有裂缝的共同作用下所形成的水锥的顶部。
S160：根据所述第一函数关系和所述第四函数关系，确定所述水锥高度、所述含水率以及储集体顶部到该生产井下方的水锥顶部之间的距离的第五函数关系。
将式(13)代入式(20)，得到所述水锥高度、所述含水率以及储集体顶部到该生产井下方的水锥顶部之间的距离的第五函数关系为：
其中，how表示所述储集体顶部到该生产井下方的水锥顶部之间的距离，h表示裂缝型生产井下方的水锥高度，fw为裂缝型生产井的含水率。
S170：根据所述第三函数关系和所述第五函数关系构建储集体顶部到该生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型。
将水锥高度表达式(19)代入式(21)，得到储集体顶部到该生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型为：
其中，how表示储集体顶部到该生产井下方的水锥顶部之间的距离，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，Δρow表示油水密度差，g为重力加速度，单位为m/s2，Lf为裂缝长度，θ为裂缝方向与竖直方向的夹角，rw为井筒半径，fw为裂缝型生产井的含水率。
S200：检测裂缝型生产井的含水率和储集体顶部距离地表的深度。
裂缝型生产井的含水率可以实时测得，可以将生产井在某一时刻的产水量与产液量的比值作为该时刻该生产井的含水率。
储集体顶部距离地表的深度即T7 4界面距离地表的深度，可以为钻遇储集体顶部时的钻入深度。
S300：根据测得的含水率，利用所述关系模型，确定储集体顶部到该生产井下方的水锥顶部之间的距离。
将测得的含水率代入所述关系模型(即式(22))，得到储集体顶部到该生产井下方的水锥顶部之间的距离。
S400：根据储集体顶部距离地表的深度和储集体顶部到该生产井下方的水锥顶部之间的距离，确定该生产井的油水界面深度。
该生产井的油水界面深度可以如式(23)所示：
其中，hd-ow表示裂缝型生产井的油水界面深度，
表示储集体顶部距离地表的深度。
本发明针对现有油水界面确定方法难以准确确定动态油水界面深度的问题，提出一种针对缝洞型油藏钻遇裂缝类储集体生产井的动态油水界面深度的确定方法。该方法基于塔河缝洞型油藏生产井含水率上升的主要原因为井控范围内出水储集体的不断增加的认识，将生产井的含水率表示为水锥高度与储集体总厚度的比值，考虑裂缝中的流体流动特征及重力的影响，基于长直细管流动模型理论，建立了钻遇裂缝型储集体生产井的动态油水界面的确定方法。本发明利用生产井日常的产油产水数据、流体参数及裂缝物性参数，基于油藏工程方法确定的水锥高度和含水率的表达式，求解得到了裂缝型生产井动态油水界面深度的表达式，能够方便快捷的确定裂缝型生产井的动态油水界面深度。
本发明的方法能够克服常规方法无法准确确定动态油水界面深度的缺陷，快速有效地确定缝洞型油藏裂缝型生产井目前的动态油水界面深度，为缝洞型油藏钻遇裂缝型储集体生产井的水窜预警和开发调整方案的制定提供重要的参考依据。
实施例二
W-1井是某油田的一口生产井，为典型的缝洞型油藏钻遇裂缝型储集体的生产井。投产初期产量达到42t/d，不含水，无水采油期长，目前日产油量为13t/d，含水率为69％。该井的T7 4界面的深度为5630m。W-1井原油性质参数和计算所用到的油藏参数见表1。
表1W-1井原油物性参数及计算所用油藏参数表
以上参数代入公式(23)，得到W-1井目前的动态油水界面深度为5651.5m，由此可以判断，W-1井动态油水界面距离井底还有21.5m。
实施例三
本申请的实施方式提供一种防止生产井水窜的预警方法，包括以下步骤：
检测生产井钻遇的储集体的类型；
当钻遇的储集体为裂缝型储集体时，采用如上文所述的油水界面深度确定方法来确定生产井的动态油水界面深度；
根据生产井的动态油水界面深度确定生产井下的水锥位置；
根据生产井下的水锥位置判断生产井发生水淹事故的风险，并由此判断是否进行水淹预警。
本发明的防止生产井水窜的预警方法能够根据钻遇储集体的类型准确确定生产井的动态油水界面深度，准确的动态油水界面深度能够有效帮助施工人员正确认识目前的水锥位置，为缝洞型油藏裂缝型生产井水淹预警及开发调整措施的制定提供重要的参考依据。
实施例四
本申请的实施方式提供一种存储介质，存储有计算机程序，所述计算机程序被处理器执行时，实现如上文所述的油水界面深度确定方法的步骤或防止生产井水窜的预警方法的步骤。
存储介质包括但不限于RAM、ROM、EEPROM、闪存或其他存储器技术、CD-ROM、数字多功能盘(DVD)或其他光盘存储、磁盒、磁带、磁盘存储或其他磁存储装置、或者可以用于存储期望的信息并且可以被计算机访问的任何其他的介质。
实施例五
本申请的实施方式提供一种电子设备，包括处理器和存储有程序代码的存储介质，所述程序代码被所述处理器执行时，实现如上文所述的油水界面深度确定方法的步骤或防止生产井水窜的预警方法的步骤。
需要注意的是，这里所使用的的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
应当理解的是，本说明书中的示例性实施方式可以由多种不同的形式来实施，并且不应当被解释为只限于这里所阐述的实施方式。提供这些实施方式是为了使得本申请的公开彻底且完整，并且将这些示例性实施方式的构思充分传达给本领域普通技术人员，而不应当理解为对本发明的限制。
Claims (13)
1.一种油水界面深度确定方法，应用于裂缝型储集体生产井，其特征在于，包括以下步骤：
基于水锥界面的压力平衡、裂缝型生产井的产油量方程、水锥高度和储集体厚度以及含水率三者之间的关系、裂缝宽度与储集体厚度之间的关系，构建裂缝型储集体顶部到裂缝型生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型；
检测裂缝型生产井的含水率和储集体顶部距离地表的深度；
根据测得的含水率，利用所述关系模型，确定储集体顶部到该生产井下方的水锥顶部之间的距离；
根据储集体顶部距离地表的深度和储集体顶部到该生产井下方的水锥顶部之间的距离，确定该生产井的油水界面深度。
2.根据权利要求1所述的油水界面深度确定方法，其特征在于，基于水锥界面的压力平衡、裂缝型生产井的产油量方程、水锥高度和储集体厚度以及含水率三者之间的关系、裂缝宽度与储集体厚度之间的关系，构建裂缝型储集体顶部到裂缝型生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型，包括：
基于水锥界面的压力平衡确定裂缝型生产井下方的水锥界面的压力平衡方程，确定裂缝型生产井的产油量方程，通过联立裂缝型生产井下方的水锥界面的压力平衡方程和该生产井的产油量方程建立该生产井下方与储集体厚度相关的水锥界面到裂缝边缘的垂向水锥高度计算式；
确定裂缝型生产井下方的水锥高度、储集体厚度以及该生产井的含水率三者之间的第一函数关系；
确定裂缝宽度与储集体厚度之间的第二函数关系；
根据所述水锥界面到裂缝边缘的垂向水锥高度计算式、所述第一函数关系和所述第二函数关系确定该生产井下方的水锥高度与所述含水率之间的第三函数关系；
确定该生产井下方的水锥高度、储集体厚度以及储集体顶部到该生产井下方的水锥顶部之间的距离的第四函数关系；
根据所述第一函数关系和所述第四函数关系，确定所述水锥高度、所述含水率以及储集体顶部到该生产井下方的水锥顶部之间的距离的第五函数关系；
根据所述第三函数关系和所述第五函数关系构建储集体顶部到该生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型。
3.根据权利要求2所述的油水界面深度确定方法，其特征在于，裂缝型生产井的产油量方程为：
其中，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，ω为裂缝宽度，f(x)为水锥界面到裂缝边缘的垂向水锥高度，p为距离井筒中心x处的裂缝内水锥界面的压力，L为裂缝长度方向的距离。
4.根据权利要求2所述的油水界面深度确定方法，其特征在于，该生产井下方与储集体厚度相关的水锥界面到裂缝边缘的垂向水锥高度计算式为：
其中，f(x)为水锥界面到裂缝边缘的垂向水锥高度，ω为裂缝宽度，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，Δρow表示油水密度差，g为重力加速度，单位为m/s2，θ为裂缝方向与竖直方向的夹角，Lf为裂缝长度。
5.根据权利要求2所述的油水界面深度确定方法，其特征在于，裂缝型生产井下方的水锥高度、储集体厚度以及该生产井的含水率三者之间的第一函数关系为：
其中，fw为裂缝型生产井的含水率，h表示裂缝型生产井下方的水锥高度，H表示储集体厚度。
6.根据权利要求2所述的油水界面深度确定方法，其特征在于，裂缝宽度与储集体厚度之间的第二函数关系为：
ω＝H·cosθ
其中，ω为裂缝宽度，H表示储集体厚度，θ为裂缝方向与竖直方向的夹角。
7.根据权利要求2所述的油水界面深度确定方法，其特征在于，该生产井下方的水锥高度与所述含水率之间的第三函数关系为：
其中，h表示裂缝型生产井下方的水锥高度，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，Δρow表示油水密度差，g为重力加速度，单位为m/s2，Lf为裂缝长度，θ为裂缝方向与竖直方向的夹角，rw为井筒半径，fw为裂缝型生产井的含水率。
8.根据权利要求2所述的油水界面深度确定方法，其特征在于，该生产井下方的水锥高度、储集体厚度以及储集体顶部到该生产井下方的水锥顶部之间的距离的第四函数关系为：
how＝H-h
其中，how表示所述储集体顶部到该生产井下方的水锥顶部之间的距离，h表示裂缝型生产井下方的水锥高度，H表示储集体厚度。
9.根据权利要求2所述的油水界面深度确定方法，其特征在于，所述水锥高度、所述含水率以及储集体顶部到该生产井下方的水锥顶部之间的距离的第五函数关系为：
其中，how表示所述储集体顶部到该生产井下方的水锥顶部之间的距离，h表示裂缝型生产井下方的水锥高度，fw为裂缝型生产井的含水率。
10.根据权利要求2所述的油水界面深度确定方法，其特征在于，储集体顶部到该生产井下方的水锥顶部之间的距离与该生产井的含水率之间的关系模型为：
其中，how表示储集体顶部到该生产井下方的水锥顶部之间的距离，Qo为产油量，单位为m3，Bo为原油体积系数，μo为原油粘度，单位为mpa·s，Δρow表示油水密度差，g为重力加速度，单位为m/s2，Lf为裂缝长度，θ为裂缝方向与竖直方向的夹角，rw为井筒半径，fw为裂缝型生产井的含水率。
11.一种防止生产井水窜的预警方法，其特征在于，包括以下步骤：
检测生产井钻遇的储集体的类型；
当钻遇的储集体为裂缝型储集体时，采用如权利要求1至10中任一项所述的油水界面深度确定方法来确定生产井的动态油水界面深度；
根据生产井的动态油水界面深度确定生产井下的水锥位置；
根据生产井下的水锥位置判断生产井发生水淹事故的风险，并由此判断是否进行水淹预警。
12.一种存储介质，存储有计算机程序，其特征在于，所述计算机程序被处理器执行时，实现如权利要求1-10中任一项所述的油水界面深度确定方法的步骤或如权利要求11所述的防止生产井水窜的预警方法的步骤。
13.一种电子设备，包括处理器和存储有程序代码的存储介质，所述程序代码被所述处理器执行时，实现如权利要求1-10中任一项所述的油水界面深度确定方法的步骤或如权利要求11所述的防止生产井水窜的预警方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010531436.9A CN113803054B (zh) | 2020-06-11 | 2020-06-11 | 油水界面深度确定方法及防止生产井水窜的预警方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010531436.9A CN113803054B (zh) | 2020-06-11 | 2020-06-11 | 油水界面深度确定方法及防止生产井水窜的预警方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113803054A true CN113803054A (zh) | 2021-12-17 |
| CN113803054B CN113803054B (zh) | 2024-04-09 |
Family
ID=78943760
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010531436.9A Active CN113803054B (zh) | 2020-06-11 | 2020-06-11 | 油水界面深度确定方法及防止生产井水窜的预警方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113803054B (zh) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2213853C2 (ru) * | 2001-08-22 | 2003-10-10 | Общество с ограниченной ответственностью "ЛУКОЙЛ-ВолгоградНИПИморнефть" | Способ разработки массивной нефтяной залежи |
| CN104463345A (zh) * | 2014-10-31 | 2015-03-25 | 中国石油天然气股份有限公司 | 一种预测地层油水界面位置的方法及装置 |
| CN105550773A (zh) * | 2015-12-10 | 2016-05-04 | 中国石油天然气股份有限公司 | 一种预测油水界面位置的方法及装置 |
| CN107246900A (zh) * | 2017-05-27 | 2017-10-13 | 中国石油天然气股份有限公司 | 一种确定油水界面位置的方法及装置 |
| CN107288618A (zh) * | 2016-03-31 | 2017-10-24 | 中国石油化工股份有限公司 | 一种缝洞型油藏油水界面确定方法及系统 |
| CN107975358A (zh) * | 2017-11-07 | 2018-05-01 | 中国石油天然气股份有限公司 | 一种油藏开采方法及装置 |
| CN109558609A (zh) * | 2017-09-26 | 2019-04-02 | 中国石油天然气股份有限公司 | 裂缝性地层油基钻井液的性能参数的确定方法和装置 |
| CN110056346A (zh) * | 2019-04-17 | 2019-07-26 | 中国石油大港油田勘探开发研究院 | 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 |
| US20190292884A1 (en) * | 2018-03-21 | 2019-09-26 | ResFrac Corporation | Systems and methods for hydraulic fracture and reservoir simulation |
-
2020
- 2020-06-11 CN CN202010531436.9A patent/CN113803054B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2213853C2 (ru) * | 2001-08-22 | 2003-10-10 | Общество с ограниченной ответственностью "ЛУКОЙЛ-ВолгоградНИПИморнефть" | Способ разработки массивной нефтяной залежи |
| CN104463345A (zh) * | 2014-10-31 | 2015-03-25 | 中国石油天然气股份有限公司 | 一种预测地层油水界面位置的方法及装置 |
| CN105550773A (zh) * | 2015-12-10 | 2016-05-04 | 中国石油天然气股份有限公司 | 一种预测油水界面位置的方法及装置 |
| CN107288618A (zh) * | 2016-03-31 | 2017-10-24 | 中国石油化工股份有限公司 | 一种缝洞型油藏油水界面确定方法及系统 |
| CN107246900A (zh) * | 2017-05-27 | 2017-10-13 | 中国石油天然气股份有限公司 | 一种确定油水界面位置的方法及装置 |
| CN109558609A (zh) * | 2017-09-26 | 2019-04-02 | 中国石油天然气股份有限公司 | 裂缝性地层油基钻井液的性能参数的确定方法和装置 |
| CN107975358A (zh) * | 2017-11-07 | 2018-05-01 | 中国石油天然气股份有限公司 | 一种油藏开采方法及装置 |
| US20190292884A1 (en) * | 2018-03-21 | 2019-09-26 | ResFrac Corporation | Systems and methods for hydraulic fracture and reservoir simulation |
| CN110056346A (zh) * | 2019-04-17 | 2019-07-26 | 中国石油大港油田勘探开发研究院 | 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113803054B (zh) | 2024-04-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108301825B (zh) | 一种高温高压储层孔隙压力分析方法 | |
| US7653488B2 (en) | Determination of point of sand production initiation in wellbores using residual deformation characteristics and real time monitoring of sand production | |
| EP2904530B1 (en) | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties | |
| MX2008015192A (es) | Metodo para determinar dimensiones de una fractura hidraulica de formacion. | |
| CN104863577A (zh) | 利用地震纵波传播时间预测地层孔隙压力的方法 | |
| CN101634620A (zh) | 稠油松散岩心饱和度参数校正方法 | |
| CN115422859B (zh) | 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 | |
| CN109798108B (zh) | 深层低渗-致密储层复杂流体性质的识别方法和装置 | |
| CN112576240B (zh) | 基于封闭井筒压力波动进行水力压裂裂缝监测的方法 | |
| US9556724B2 (en) | Method for determining parameters of a bottomhole and a near-bottomhole zone of a wellbore | |
| CN110007344B (zh) | 一种断溶体储层连通性的地震识别方法及装置 | |
| EP3486427B1 (en) | Method of determining maximum stress in a well drilled in a reservoir | |
| US10753203B2 (en) | Systems and methods to identify and inhibit spider web borehole failure in hydrocarbon wells | |
| CN112560246A (zh) | 一种目标井散点地层压力系数的预测方法 | |
| CN113803054B (zh) | 油水界面深度确定方法及防止生产井水窜的预警方法 | |
| Akram et al. | A model to predict wireline formation tester sample contamination | |
| US20080230221A1 (en) | Methods and systems for monitoring near-wellbore and far-field reservoir properties using formation-embedded pressure sensors | |
| CN113803055B (zh) | 油水界面深度确定方法及防止生产井水窜的预警方法 | |
| CN113738344B (zh) | 油水界面深度确定方法及防止生产井水窜的预警方法 | |
| CN114427444B (zh) | 一种自喷期采油井的井底压力预测方法 | |
| CN114638179A (zh) | 一种微量气侵条件下钻井安全关井周期的确定方法 | |
| CN110764161B (zh) | 一种通过钻探资料判识碳酸盐岩断裂破碎带的综合方法 | |
| CN108930534B (zh) | 随钻核磁共振地层温度影响的校正方法与装置 | |
| Tjomsland et al. | Application of Multirate Well Tests to Scale Management: Part 2—Interpretation of MRTs With Known Produced-Water Origin | |
| RU2768341C1 (ru) | Способ прогнозирования дебита скважин с учетом анизотропии проницаемости карбонатных горных пород |
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 |