CN110318742B - 基于压裂井生产数据确定裂缝闭合长度的方法和系统 - Google Patents

基于压裂井生产数据确定裂缝闭合长度的方法和系统 Download PDF

Info

Publication number
CN110318742B
CN110318742B CN201810275660.9A CN201810275660A CN110318742B CN 110318742 B CN110318742 B CN 110318742B CN 201810275660 A CN201810275660 A CN 201810275660A CN 110318742 B CN110318742 B CN 110318742B
Authority
CN
China
Prior art keywords
yield
dimensionless
data
fracture
well
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.)
Active
Application number
CN201810275660.9A
Other languages
English (en)
Other versions
CN110318742A (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.)
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
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 China Petroleum and Chemical Corp, Sinopec Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201810275660.9A priority Critical patent/CN110318742B/zh
Publication of CN110318742A publication Critical patent/CN110318742A/zh
Application granted granted Critical
Publication of CN110318742B publication Critical patent/CN110318742B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/25—Methods for stimulating production
    • E21B43/26—Methods for stimulating production by forming crevices or fractures
    • 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
    • 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

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开一种基于压裂井生产数据确定裂缝闭合长度的方法，包括以下步骤：对原始生产数据进行处理以获取产量‑裂缝闭合数据体；对产量‑裂缝闭合数据体进行处理以获取实际产量变化数据；对实际产量变化数据进行处理以获取实际产量图版数据；将实际产量图版数据和产量‑裂缝闭合图版数据在同一双对数图版中耦合处理以获取产量‑裂缝闭合特征值数据；基于实际产量图版数据、产量‑裂缝闭合图版数据与产量‑裂缝闭合特征值数据确定裂缝闭合位置，从而确定裂缝闭合长度。本发明还公开了与方法相对应的系统。上述方法和系统能够快速、准确地获取压裂井裂缝部分闭合的长度和位置，进而提高对压裂井生产动态分析的精度。

Description

基于压裂井生产数据确定裂缝闭合长度的方法和系统
技术领域
本发明涉及油气勘探开发领域，具体涉及一种基于压裂井生产数据确定裂缝闭合长度方法和系统。
背景技术
压裂井在增加可采储量、提高单井产量、减少开采成本以及复杂油气藏的开发等方面具有重要的战略意义和经济效益，因此对压裂井的产量监测、效果评价以及增产改造等已成为当今油气藏开采领域的热点问题。监测油井的产量递减过程是现代油气藏开发过程中了解地层特征的一种技术。
在实际生产中，压裂井沿着裂缝的产油量或产气量是不均匀的，或者仅有部分裂缝段产油或产气，而其他裂缝段的产油量或产气量几乎为零。导致上述情况发生的原因如下：储集层的非均质性；地层压力降低，裂缝重新闭合；压裂液反排不充分，使裂缝堵塞；其他污染因素等。
目前，压裂井的产量递减分析方法是以整个压裂裂缝段全部投入生产为依据，若采用传统的产量递减分析方法解释部分裂缝段不产油或不产气的压裂井将导致解释参数(地层渗透率K、表皮系数S等)失真。传统的产量递减分析方法假设各段裂缝的流量相等且不为零，该假设与实际生产中沿裂缝段的流量不均匀 (例如：只有部分裂缝段产油或产气；或全部裂缝段均产油或产气，但是各段裂缝的产油量或产气量不相等情况)的压裂井在流动形态上有很大差异，不能准确描述各段裂缝的产油量或产气量不均匀的压裂井的流体流动状态，因此有必要研制一种确定压裂井压裂裂缝闭合长度的方法，以适应各裂缝段产油量或产气量不均匀现象的监测和控制。
授权公告号为CN101560879B的发明专利公开了一种用于低渗透气藏的试井分析方法，包括如下步骤：测量气井在开井一段时间内的井底压力值，得到压力差与时间的关系值；判断关系值与预设关系值是否相等，若相等，则进行下一步骤，否则修正预设关系值对应的试井解释参数；将预设关系值对应的试井解释参数作为气井的试井解释参数；根据气井的试井解释参数分析该气井的气体渗流情况，以控制气井开采。该发明专利还公开了一种试井分析控制系统，包括气井压力测量单元、判断单元、试井解释参数处理单元和试井分析控制单元。虽然该发明专利的试井分析方法及系统能够描述气藏的气体渗流状态，但是无法满足该产量递减分析与气藏实际渗流情况相符合的要求。
申请公布号为CN 105298479A的发明专利申请公开了一种压裂井产油(气) 位置的诊断方法，该方法通过压力与时间的变化来寻找裂缝闭合的位置，在缺少压力数据或者压力数据不准确时无法进行。
申请公布号为CN104612658A的发明专利申请公开了一种用于水平井分段产液识别的试井分析控制方法，该方法的步骤与申请公布号为CN 105298479A的发明专利申请中的步骤相同，只有预设关系值与对应的试井解释参数的关系式不同。同时该发明专利申请还公开了一种试井分析控制系统，与申请公布号为CN 105298479A的发明专利申请中所公开的系统相同。因此，该发明专利同样无法满足产量递减分析与油藏实际渗流情况相符合的要求。
在《低渗透储层产量递减模型的渗流机理及应用》一文中公布了一种使用 Arps广义产量递减模型来解释油藏实际渗流情况的方法，该文刊登于《石油学报》 2011年01期。该方法步骤简单，但对于较复杂的地层情况，比如压裂井不均匀产气等情况，该方法的准确性较低，需要改进。
在《压裂气井产能预测新方法研究》一文中公布了压裂气井产能预测的方法，利用Duhamel褶积通过数值反演得到了Fetkovich/Arps无限导流垂直裂缝井联合产量递减图版，通过图版拟合方法得到产能预测参数。该文刊登在《石油天然气学报》2012年02期。该方法的缺点是只能解决无限导流裂缝的相关问题，与现场实际参数相比，误差较大。
在《不对称垂直裂缝井产量递减规律》一文中公布了一种描述垂直裂缝的不对称性和裂缝导流能力的时变性的方法，该文刊登在《西安石油大学学报》2014 年02期。该方法可以相对准确的描述裂缝的不对称性，但针对不对称性的原因没有进行深入的探究，对于裂缝部分闭合现象导致的不对称性也没有相应的说明，改进空间较大。
发明内容
为了解决现有技术中存在的技术问题，本发明提供一种基于压裂井生产数据确定裂缝闭合长度的方法，包括以下步骤：
步骤一：对原始生产数据进行处理以获取产量-裂缝闭合数据体；
步骤二：对产量-裂缝闭合数据体进行处理以获取实际产量变化数据；
步骤三：对实际产量变化数据进行处理以获取实际产量图版数据；
步骤四：将实际产量图版数据和产量-裂缝闭合图版数据在同一双对数图版中耦合处理以获取产量-裂缝闭合特征值数据；
步骤五：基于实际产量图版数据、产量-裂缝闭合图版数据与产量-裂缝闭合特征值数据确定裂缝闭合位置，从而确定裂缝闭合长度。
在一个优选实施例中，在步骤一中，所述的产量-裂缝闭合数据体包括生产时间、井底流压和日产量数据。
在一个优选实施例中，在步骤二中，所述实际产量变化数据包括：不同物质平衡时间下的规整化产量，规整化产量积分和规整化产量积分导数数据。
在一个优选实施例中，在步骤三中，所述实际产量图版数据包括：物质平衡时间-规整化产量、物质平衡时间-规整化产量积分和物质平衡时间-规整化产量积分导数三条曲线数据。
在一个优选实施例中，在步骤四中，产量-裂缝闭合图版数据是以考虑裂缝流量不均的压裂井模型为基础，调整产量-裂缝闭合图版数据的特征值，获得的不同特征值下的曲线图版数据。
在一个优选实施例中，所述产量-裂缝闭合图版数据的特征值包括无量纲边界尺寸、无量纲近井段裂缝长度和无量纲近井段裂缝流量，分别有下式获得：
reD＝re/xf
xf1D＝xf1/xf
q1D＝q1/q
式中：reD为无量纲边界尺寸，无量纲；xf1D为无量纲近井段裂缝长度，无量纲；q1D为无量纲近井段裂缝流量，无量纲；re为圆形边界尺寸，单位为m； xf为裂缝总长度，单位为m；xf1为近井段裂缝长度，单位为m；q1为近井段裂缝流量，单位为m3；q为裂缝总流量，单位为m3。
在一个优选实施例中，在步骤五中，根据无量纲边界尺寸计算产量-裂缝闭合系数，
Figure BDA0001613582950000031
其中，reD为无量纲边界尺寸，无量纲；bDpss为产量-裂缝闭合系数，无量纲。
在一个优选实施例中，在步骤五中，在实际产量图版数据中选择一个特征数据点(tc，q/dp)M，以及选择对应的产量-裂缝闭合图版数据中的点(tcDd，qDd)M，从而获得参数：储层渗透率、表皮系数、井控面积和储量，进而获取裂缝闭合长度，
其中实际产量图版数据中的特征数据点中：tc为物质平衡时间，单位为D； q/dp为规整化产量，单位为m3/D/MPa；产量-裂缝闭合图版数据中的点中：tcDd为无量纲物质平衡时间，无量纲；qDd为无量纲规整化产量，无量纲；M为下标，表示这两个点是对应的。
在一个优选实施例中，其中储层渗透率由下式获得：
Figure BDA0001613582950000041
(油)
Figure BDA0001613582950000042
(气)
式中：K为油井/气井储层渗透率，单位为mD；q/△p为油井规整化产量，单位为m3/D/MPa；q/△pp为气井规整化产量，单位为m3/D/MPa；qDd为无量纲规整化产量，无量纲；μ为粘度，单位为mPa.s；B为体积系数，单位为m3/m3； h为地层厚度，单位为m；
其中储量由下式获得：
Figure BDA0001613582950000043
(油)
Figure BDA0001613582950000044
(气)
式中：N为油藏储量，单位为m3；G为气藏储量，单位为m3；tc为物质平衡时间，单位为D；tcDd为无量纲物质平衡时间，无量纲；tca为物质平衡拟时间，单位为D；tcaDd为无量纲物质平衡拟时间，无量纲；q/△p为油井规整化产量，单位为m3/D/MPa；q/△pp为气井规整化产量，单位为m3/D/MPa；qDd为无量纲规整化产量，无量纲；Ct为综合压缩系数，单位为MPa-1；Sw为水相饱和度，无量纲；
其中井控面积由下式获得：
Figure BDA0001613582950000051
(油)
Figure BDA0001613582950000052
(气)
式中：A为油/气井空面积，单位为m2；N为油藏储量，单位为m3；G为气藏储量，单位为m3；Bo为油相体积系数，单位为m3/m3；Bgi为油相体积系数，单位为m3/m3；φ为渗透率，无量纲；h为地层厚度，单位为m；Swi为水相饱和度，无量纲；
其中表皮系数由下式获得：
Figure BDA0001613582950000053
Figure BDA0001613582950000054
Figure BDA0001613582950000055
式中：re为圆形边界尺寸，单位为m；A为油/气井空面积，单位为m2；xf为裂缝总长度，单位为m；reD为无量纲边界尺寸，无量纲；s为表皮系数，无量纲； rw为井筒半径，单位为m；
其中裂缝闭合长度由下式获得：
xf1＝xf1Dxf
xf2＝xf-xf1
式中：xf1为近井段裂缝长度，单位为m；xf1D无量纲近井段裂缝长度，无量纲；xf为裂缝总长度，单位为m；xf2为裂缝闭合长度，单位为m。
本发明的另一方面还公开一种基于压裂井生产数据确定裂缝闭合长度的系统，包括以下单元：
原始数据处理单元，所述数据处理单元对原始生产数据进行处理以获取产量 -裂缝闭合数据体；
第一中间处理单元，所述第一中间处理单元对所述原始数据处理单元得到的产量-裂缝闭合数据体进行处理以获取实际产量变化数据；
第二中间处理单元，所述第二中间处理单元对所述第一中间处理单元得到的实际产量变化数据进行处理以获取实际产量图版数据；
第三中间处理单元，所述第三中间处理单元将所述第二中间处理单元得到的实际产量图版数据和产量-裂缝闭合图版数据在同一双对数图版中耦合处理以获取产量-裂缝闭合特征值数据；以及
裂缝闭合长度确定单元：所述裂缝闭合长度确定单元基于所述第二中间处理单元得到的实际产量图版数据、所述产量-裂缝闭合图版数据与所述第三中间处理单元得到的产量-裂缝闭合特征值数据确定裂缝闭合位置，从而确定裂缝闭合长度。
本发明公开的基于压裂井生产数据确定裂缝闭合长度的方法和系统，针对压裂井裂缝部分闭合导致的产气不均匀的现象，依据生产数据能够快速、准确地获取压裂井裂缝部分闭合的长度和位置，以适应各裂缝段产油量或产气量不均匀现象的监测和控制，进而提高对压裂井生产动态分析的精度。总之，本方法和系统可以确定压裂缝部分闭合后每段裂缝长度，可以准确预测产量变化趋势。
附图说明
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中：
图1是本发明实施例一中的基于压裂井生产数据确定裂缝闭合长度的方法的流程图；
图2是本发明实施例一中的步骤四中的特定特征值下的产量-裂缝闭合图版数据；
图3是本发明实施例一中的步骤四中的产量-裂缝闭合特征数据曲线；
图4是本发明实施例二中的基于压裂井生产数据确定裂缝闭合长度的系统的结构示意图。
在附图中，相同的部件使用相同的附图标记。附图并未按照实际的比例。
具体实施方式
在本发明中涉及很多单位，下面做具体的说明，其中m表示单位米，m2表示单位平方米，m3表示单位立方米，D表示时间单位天，m3/D/Mpa表示单位立方米每天每兆帕，mPa.s表示单位毫帕秒，m3/m3表示立方米每立方米，MPa-1表示单位每兆帕，m3/D表示单位立方米每天，mD表示的单位毫达西。
下面将结合附图对本发明作进一步说明。
实施例一
如图1所示，为本实施例中的基于压裂井生产数据确定裂缝闭合长度的方法的流程图；该方法包括一下步骤：
步骤S1：对原始生产数据进行处理以获取产量-裂缝闭合数据体。
步骤S2：对产量-裂缝闭合数据体进行处理以获取实际产量变化数据；
步骤S3：对实际产量变化数据进行处理以获取实际产量图版数据；
步骤S4：将实际产量图版数据和产量-裂缝闭合图版数据在同一双对数图版中耦合处理以获取产量-裂缝闭合特征值数据；
步骤S5：基于实际产量图版数据、产量-裂缝闭合图版数据与产量-裂缝闭合特征值数据确定裂缝闭合位置，从而确定裂缝闭合长度。
下面在图1的基础上，结合生产中的实例对本发明进行说明。
目前，生产动态分析需要获取基础参数。这些基础参数包括油气藏类型、井型、有效厚度、孔隙度、地层温度、原始地层压力和气体比重等。
本实施例以大牛地气藏中的一口压裂井作为一个实施例，井号为#D1-1-24，基础参数为：井筒半径0.108m，产层厚度8.6m，孔隙度0.11，平均残余水饱和度23％，天然气相对密度为0.7，原始条件下的体积系数是0.0056m3/m3，原始条件下气体粘度为0.0231mPa.s，原始条件下的天然气压缩系数为0.0035MPa-1，原始地层压力为23.55MPa。
在步骤S1中，对原始生产数据进行处理以获取产量-裂缝闭合数据体。
原始生产数据包括每次对油气井进行测试的时间、油气井的日产量和井底流压或井口压力。将原始生产数据按时间顺序插入来获取产量-裂缝闭合数据体。优选的，产量-裂缝闭合数据体包括生产时间、井底压力和日产量等数据。在本实施例中，#D1-1-24的产量-裂缝闭合数据体就是#D1-1-24的测试时间、#D1-1-24 的日产量和#D1-1-24的井底流压或井口压力数据。
在步骤S2中，对产量-裂缝闭合数据体进行处理以获取实际产量变化数据。优选的，所述实际产量变化数据包括：不同物质平衡时间下的规整化产量，规整化产量积分和规整化产量积分导数数据。
在该步骤中，首先计算物质平衡时间：
对于油井来说，物质平衡时间tc：
Figure BDA0001613582950000071
式中，Np：可采储量，108m3；
q：日产量，m3/D。
对于气井来说，物质平衡时间为tca：
Figure BDA0001613582950000081
式中，μ：气体粘度，mPa.s；
Cg：气体压缩系数，MPa-1。
其次，计算规整化产量、规整化产量积分和规整化产量积分导数。油井和气井的规整化产量分别
Figure BDA0001613582950000082
和
Figure BDA0001613582950000083
为：
Figure BDA0001613582950000084
(油)
Figure BDA0001613582950000085
(气)
式中，q：日产量，m3/D；
pi：原始地层压力，M Pa；
pwf：井底流压，M Pa；
ppi：气体拟原始地层压力，M Pa；
ppwf：气体拟井底流压，M Pa；
油井和气井的规整化产量积分分别为：
Figure BDA0001613582950000086
(油)
Figure BDA0001613582950000087
(气)
油井和气井的规整化产量积分导数分别为：
Figure BDA0001613582950000088
(油)
Figure BDA0001613582950000089
(气)
在本实施例中，由于测试井为气井，#D1-1-24的实际产量变化数据就是 #D1-1-24的规整化产量
Figure BDA00016135829500000810
规整化产量积分
Figure BDA00016135829500000811
和规整化产量积分导数
Figure BDA00016135829500000812
在步骤S3中，对实际产量变化数据进行处理以获取实际产量图版数据。优选的，其中实际产量图版数据包括物质平衡时间-规整化产量、物质平衡时间-规整化产量积分和物质平衡时间-规整化产量积分导数三条曲线数据。在本实施例中，#D1-1-24的实际产量图版数据就是#D1-1-24的物质平衡时间-规整化产量、物质平衡时间-规整化产量积分和物质平衡时间-规整化产量积分导数三条曲线数据。
在步骤S4中，将实际产量图版数据和产量-裂缝闭合图版数据在同一双对数图版中耦合处理以获取产量-裂缝闭合特征值数据。优选的，产量-裂缝闭合图版数据是以考虑裂缝流量不均的压裂井模型为基础，调整产量-裂缝闭合图版数据的特征值，获得的不同特征值下的曲线图版数据。
在该步骤中，图2即是特征值分别为xf1D＝0.4、q1D＝0.6和reD＝10时的产量- 裂缝闭合图版。本实施例中，#D1-1-24的实际产量图版数据与特征值为xf1D＝0.4、 q1D＝0.6和reD＝5的产量-裂缝闭合数据的耦合数据曲线拟合效果较好，如图3所示。
在步骤S5中，基于实际产量图版数据、产量-裂缝闭合图版数据与产量-裂缝闭合特征值数据确定裂缝闭合位置，从而确定裂缝闭合长度。
优选的，在该步骤中：
首先根据无量纲边界尺寸reD计算产量-裂缝闭合系数bDpss：
Figure BDA0001613582950000091
其中，reD为无量纲边界尺寸，无量纲；bDpss为产量-裂缝闭合系数，无量纲。
其次，在实际产量图版数据中选择一个特征数据点(tc,q/dp)M以及相应产量- 裂缝闭合图版数据中的点(tcDd,qDd)M，其中实际产量图版数据中的特征数据点中： tc为物质平衡时间，单位为D；q/dp为规整化产量，单位为m3/D/MPa；产量-裂缝闭合图版数据中的点中：tcDd为无量纲物质平衡时间，无量纲；qDd为无量纲规整化产量，无量纲；M为下标，表示这两个点是对应的。
优选的，储层渗透率K可以由下式计算获得：
Figure BDA0001613582950000092
(油)
Figure BDA0001613582950000093
(气)
式中，K为油井/气井储层渗透率，单位为mD；q/△p为油井规整化产量，单位为m3/D/MPa；q/△pp为气井规整化产量，单位为m3/D/MPa；qDd为无量纲规整化产量，无量纲；μ为粘度，单位为mPa.s；B为体积系数，单位为m3/m3； h为地层厚度，单位为m。
优选的，油藏储量N、气藏储量G：
Figure BDA0001613582950000101
(油)
Figure BDA0001613582950000102
(气)
式中，式中：N为油藏储量，单位为m3；G为气藏储量，单位为m3；tc为物质平衡时间，单位为D；tcDd为无量纲物质平衡时间，无量纲；tca为物质平衡拟时间，单位为D；tcaDd为无量纲物质平衡拟时间，无量纲；q/△p为油井规整化产量，单位为m3/D/MPa；q/△pp为气井规整化产量，单位为m3/D/MPa；qDd为无量纲规整化产量，无量纲；Ct为综合压缩系数，单位为MPa-1；Sw为水相饱和度，无量纲。
优选的，井控面积A由下式可以得到：
Figure BDA0001613582950000103
(油)
Figure BDA0001613582950000104
(气)
式中：A为油/气井空面积，单位为m2；N为油藏储量，单位为m3；G为气藏储量，单位为m3；Bo为油相体积系数，单位为m3/m3；Bgi为油相体积系数，单位为m3/m3；φ为渗透率，无量纲；h为地层厚度，单位为m；Swi为水相饱和度，无量纲。
优选的，表皮系数s由下式可以得到：
Figure BDA0001613582950000105
Figure BDA0001613582950000106
Figure BDA0001613582950000107
式中，式中：re为圆形边界尺寸，单位为m；A为油/气井空面积，单位为 m2；xf为裂缝总长度，单位为m；reD为无量纲边界尺寸，无量纲；s为表皮系数，无量纲；rw为井筒半径，单位为m。
进而计算裂缝剩余长度xf1与裂缝闭合长度xf2：
xf1＝xf1Dxf
xf2＝xf-xf1
式中：xf1为近井段裂缝长度，单位为m；xf1D无量纲近井段裂缝长度，无量纲；xf为裂缝总长度，单位为m；xf2为裂缝闭合长度，单位为m。
本实施例最终结果参数如下：K＝0.34mD，G＝0.58×108m3，A＝0.83×106m2， re＝515m，s＝-6.64，xf1＝41.2m，xf2＝61.8m。总裂缝半长达到103m，裂缝闭合长度为61.8m。有以上计算结果可知，本发明可以准确地确定压裂缝部分闭合长度，为准确预测产量变化趋势奠定基础。
实施例二
如图4所示，为本实施例中的基于压裂井生产数据确定裂缝闭合长度的系统的结构示意图。该系统，包括：
原始数据处理单元1，所述数据处理单元1对原始生产数据进行处理以获取产量-裂缝闭合数据体；
第一中间处理单元2，所述第一中间处理单元2对所述原始数据处理单元1 得到的产量-裂缝闭合数据体进行处理以获取实际产量变化数据；
第二中间处理单元3，所述第二中间处理单元3对所述第一中间处理单元2 得到的实际产量变化数据进行处理以获取实际产量图版数据；
第三中间处理单元4，所述第三中间处理单元4将所述第二中间处理单元3 得到的实际产量图版数据和产量-裂缝闭合图版数据在同一双对数图版中耦合处理以获取产量-裂缝闭合特征值数据；以及
裂缝闭合长度确定单元5：所述裂缝闭合长度确定单元5基于所述第二中间处理单元3得到的实际产量图版数据、所述产量-裂缝闭合图版数据与所述第三中间处理单元4得到的产量-裂缝闭合特征值数据确定裂缝闭合位置，从而确定裂缝闭合长度。
本发明中的系统与本发明中的方法为相互对应的关系，本领域技术人员跟具前面所述的方法可以得到相应的系统具体的结构和操作过程，在此不再赘述。
虽然已经参考优选实施例对本发明进行了描述，但在不脱离本发明的范围的情况下，可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是，只要不存在结构冲突，各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文中公开的特定实施例，而是包括落入权利要求的范围内的所有技术方案。

Claims (3)

1.一种基于压裂井生产数据确定裂缝闭合长度的方法，其特征在于，包括以下步骤：
步骤一：对原始生产数据进行处理以获取产量-裂缝闭合数据体，所述的产量-裂缝闭合数据体包括生产时间、井底流压和日产量数据；
步骤二：对产量-裂缝闭合数据体进行处理以获取实际产量变化数据；
步骤三：对实际产量变化数据进行处理以获取实际产量图版数据，所述实际产量图版数据包括：物质平衡时间-规整化产量、物质平衡时间-规整化产量积分和物质平衡时间-规整化产量积分导数三条曲线数据；
步骤四：将实际产量图版数据和产量-裂缝闭合图版数据在同一双对数图版中耦合处理以获取产量-裂缝闭合图版数据的特征值，所述产量-裂缝闭合图版数据是以考虑裂缝流量不均的压裂井模型为基础，调整产量-裂缝闭合图版数据的特征值，获得的不同特征值下的曲线图版数据；
步骤五：基于实际产量图版数据、产量-裂缝闭合图版数据与产量-裂缝闭合图版数据的特征值确定裂缝闭合位置，从而确定裂缝闭合长度；
在步骤五中，根据无量纲边界尺寸计算产量-裂缝闭合系数，
Figure FDA0003632405580000011
其中，reD为无量纲边界尺寸，无量纲；bDpss为产量-裂缝闭合系数，无量纲；
在步骤五中，在实际产量图版数据中选择一个特征数据点(tc，q/dp)M，以及选择对应的产量-裂缝闭合图版数据中的点(tcDd，qDd)M，从而获得参数：储层渗透率、表皮系数、井控面积和储量，进而获取裂缝闭合长度，
其中实际产量图版数据中的特征数据点中：tc为物质平衡时间，单位为D；q/dp为规整化产量，单位为m3/D/MPa；产量-裂缝闭合图版数据中的点中：tcDd为无量纲物质平衡时间，无量纲；qDd为无量纲规整化产量，无量纲；M为下标，表示这两个点是对应的；
其中储层渗透率由下式获得：
Figure FDA0003632405580000012
Figure FDA0003632405580000021
式中：K为油井/气井储层渗透率，单位为mD；q/Δp为油井规整化产量，单位为m3/D/MPa；q/Δpp为气井规整化产量，单位为m3/D/MPa；qDd为无量纲规整化产量，无量纲；μ为粘度，单位为mPa.s；B为体积系数，单位为m3/m3；h为地层厚度，单位为m；
其中储量由下式获得：
Figure FDA0003632405580000022
Figure FDA0003632405580000023
式中：N为油藏储量，单位为m3；G为气藏储量，单位为m3；tc为物质平衡时间，单位为D；tcDd为无量纲物质平衡时间，无量纲；tca为物质平衡拟时间，单位为D；tcaDd为无量纲物质平衡拟时间，无量纲；q/Δp为油井规整化产量，单位为m3/D/MPa；q/Δpp为气井规整化产量，单位为m3/D/MPa；qDd为无量纲规整化产量，无量纲；Ct为综合压缩系数，单位为MPa-1；Sw为水相饱和度，无量纲；
其中井空面积由下式获得：
Figure FDA0003632405580000024
Figure FDA0003632405580000025
式中：A为油/气井空面积，单位为m2；N为油藏储量，单位为m3；G为气藏储量，单位为m3；Bo为油相体积系数，单位为m3/m3；Bgi为气相体积系数，单位为m3/m3；φ为渗透率，无量纲；h为地层厚度，单位为m；Swi为水相饱和度，无量纲；
其中表皮系数由下式获得：
Figure FDA0003632405580000026
Figure FDA0003632405580000027
Figure FDA0003632405580000028
式中：re为圆形边界尺寸，单位为m；A为油/气井空面积，单位为m2；xf为裂缝总长度，单位为m；reD为无量纲边界尺寸，无量纲；s为表皮系数，无量纲；rw为井筒半径，单位为m；
其中裂缝闭合长度由下式获得：
xf1＝xf1Dxf
xf2＝xf-xf1
式中：xf1为近井段裂缝长度，单位为m；xf1D无量纲近井段裂缝长度，无量纲；xf为裂缝总长度，单位为m；xf2为裂缝闭合长度，单位为m；
所述产量-裂缝闭合图版数据的特征值包括无量纲边界尺寸、无量纲近井段裂缝长度和无量纲近井段裂缝流量，分别有下式获得：
reD＝re/xf
xf1D＝xf1/xf
q1D＝q1/q
式中：reD为无量纲边界尺寸，无量纲；xf1D为无量纲近井段裂缝长度，无量纲；q1D为无量纲近井段裂缝流量，无量纲；re为圆形边界尺寸，单位为m；xf为裂缝总长度，单位为m；xf1为近井段裂缝长度，单位为m；q1为近井段裂缝流量，单位为m3；q为裂缝总流量，单位为m3。
2.根据权利要求1所述的基于压裂井生产数据确定裂缝闭合长度的方法，其特征在于，在步骤二中，所述实际产量变化数据包括：不同物质平衡时间下的规整化产量，规整化产量积分和规整化产量积分导数数据。
3.基于压裂井生产数据确定裂缝闭合长度的系统，其特征在于，包括：
原始数据处理单元，所述数据处理单元对原始生产数据进行处理以获取产量-裂缝闭合数据体，所述的产量-裂缝闭合数据体包括生产时间、井底流压和日产量数据；
第一中间处理单元，所述第一中间处理单元对所述原始数据处理单元得到的产量-裂缝闭合数据体进行处理以获取实际产量变化数据；
第二中间处理单元，所述第二中间处理单元对所述第一中间处理单元得到的实际产量变化数据进行处理以获取实际产量图版数据，所述实际产量图版数据包括：物质平衡时间-规整化产量、物质平衡时间-规整化产量积分和物质平衡时间-规整化产量积分导数三条曲线数据；
第三中间处理单元，所述第三中间处理单元将所述第二中间处理单元得到的实际产量图版数据和产量-裂缝闭合图版数据在同一双对数图版中耦合处理以获取产量-裂缝闭合图版数据的特征值，所述产量-裂缝闭合图版数据是以考虑裂缝流量不均的压裂井模型为基础，调整产量-裂缝闭合图版数据的特征值，获得的不同特征值下的曲线图版数据；以及
裂缝闭合长度确定单元：所述裂缝闭合长度确定单元基于所述第二中间处理单元得到的实际产量图版数据、所述产量-裂缝闭合图版数据与所述第三中间处理单元得到的产量-裂缝闭合图版数据的特征值确定裂缝闭合位置，从而确定裂缝闭合长度；
根据无量纲边界尺寸计算产量-裂缝闭合系数，
Figure FDA0003632405580000041
其中，reD为无量纲边界尺寸，无量纲；bDpss为产量-裂缝闭合系数，无量纲；
在实际产量图版数据中选择一个特征数据点(tc，q/dp)M，以及选择对应的产量-裂缝闭合图版数据中的点(tcDd，qDd)M，从而获得参数：储层渗透率、表皮系数、井控面积和储量，进而获取裂缝闭合长度，
其中实际产量图版数据中的特征数据点中：tc为物质平衡时间，单位为D；q/dp为规整化产量，单位为m3/D/MPa；产量-裂缝闭合图版数据中的点中：tcDd为无量纲物质平衡时间，无量纲；qDd为无量纲规整化产量，无量纲；M为下标，表示这两个点是对应的；
其中储层渗透率由下式获得：
Figure FDA0003632405580000042
Figure FDA0003632405580000043
式中：K为油井/气井储层渗透率，单位为mD；q/Δp为油井规整化产量，单位为m3/D/MPa；q/Δpp为气井规整化产量，单位为m3/D/MPa；qDd为无量纲规整化产量，无量纲；μ为粘度，单位为mPa.s；B为体积系数，单位为m3/m3；h为地层厚度，单位为m；
其中储量由下式获得：
Figure FDA0003632405580000051
Figure FDA0003632405580000052
式中：N为油藏储量，单位为m3；G为气藏储量，单位为m3；tc为物质平衡时间，单位为D；tcDd为无量纲物质平衡时间，无量纲；tca为物质平衡拟时间，单位为D；tcaDd为无量纲物质平衡拟时间，无量纲；q/Δp为油井规整化产量，单位为m3/D/MPa；q/Δpp为气井规整化产量，单位为m3/D/MPa；qDd为无量纲规整化产量，无量纲；Ct为综合压缩系数，单位为MPa-1；Sw为水相饱和度，无量纲；
其中井空面积由下式获得：
Figure FDA0003632405580000053
Figure FDA0003632405580000054
式中：A为油/气井空面积，单位为m2；N为油藏储量，单位为m3；G为气藏储量，单位为m3；Bo为油相体积系数，单位为m3/m3；Bgi为气相体积系数，单位为m3/m3；φ为渗透率，无量纲；h为地层厚度，单位为m；Swi为水相饱和度，无量纲；
其中表皮系数由下式获得：
Figure FDA0003632405580000055
Figure FDA0003632405580000056
Figure FDA0003632405580000057
式中：re为圆形边界尺寸，单位为m；A为油/气井空面积，单位为m2；xf为裂缝总长度，单位为m；reD为无量纲边界尺寸，无量纲；s为表皮系数，无量纲；rw为井筒半径，单位为m；
其中裂缝闭合长度由下式获得：
xf1＝xf1Dxf
xf2＝xf-xf1
式中：xf1为近井段裂缝长度，单位为m；xf1D无量纲近井段裂缝长度，无量纲；xf为裂缝总长度，单位为m；xf2为裂缝闭合长度，单位为m；
所述产量-裂缝闭合图版数据的特征值包括无量纲边界尺寸、无量纲近井段裂缝长度和无量纲近井段裂缝流量，分别有下式获得：
reD＝re/xf
xf1D＝xf1/xf
q1D＝q1/q
式中：reD为无量纲边界尺寸，无量纲；xf1D为无量纲近井段裂缝长度，无量纲；q1D为无量纲近井段裂缝流量，无量纲；re为圆形边界尺寸，单位为m；xf为裂缝总长度，单位为m；xf1为近井段裂缝长度，单位为m；q1为近井段裂缝流量，单位为m3；q为裂缝总流量，单位为m3。
CN201810275660.9A 2018-03-30 2018-03-30 基于压裂井生产数据确定裂缝闭合长度的方法和系统 Active CN110318742B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810275660.9A CN110318742B (zh) 2018-03-30 2018-03-30 基于压裂井生产数据确定裂缝闭合长度的方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810275660.9A CN110318742B (zh) 2018-03-30 2018-03-30 基于压裂井生产数据确定裂缝闭合长度的方法和系统

Publications (2)

Publication Number Publication Date
CN110318742A CN110318742A (zh) 2019-10-11
CN110318742B true CN110318742B (zh) 2022-07-15

Family

ID=68111354

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810275660.9A Active CN110318742B (zh) 2018-03-30 2018-03-30 基于压裂井生产数据确定裂缝闭合长度的方法和系统

Country Status (1)

Country Link
CN (1) CN110318742B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112835124B (zh) * 2021-03-10 2024-01-26 长江大学 基于成像测井和阵列声波测井资料的裂缝有效性评价方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0456424A2 (en) * 1990-05-07 1991-11-13 Halliburton Company Method of determining fracture characteristics of subsurface formations
EP0589591A1 (en) * 1992-09-11 1994-03-30 Halliburton Company Downhole fracture test and analysis
US6101447A (en) * 1998-02-12 2000-08-08 Schlumberger Technology Corporation Oil and gas reservoir production analysis apparatus and method
WO2003102371A1 (en) * 2002-05-31 2003-12-11 Schlumberger Canada Limited Method and apparatus for effective well and reservoir evaluation without the need for well pressure history
CN101560879A (zh) * 2008-04-15 2009-10-21 中国石油大学(北京) 用于低渗透气藏的试井分析控制系统及方法
CN104594872A (zh) * 2015-01-04 2015-05-06 西南石油大学 一种优化致密气藏压裂水平井裂缝导流能力的方法
CN104612658A (zh) * 2014-11-25 2015-05-13 中国石油大学(北京) 用于水平井分段产液识别的试井分析控制系统及方法
CN105298479A (zh) * 2015-09-09 2016-02-03 中国石油大学(北京) 压裂直井产油(气)位置的诊断方法及其系统
CN106353809A (zh) * 2016-08-11 2017-01-25 中国石油大学(华东) 一种压裂裂缝网络的反演表征方法
WO2017041074A1 (en) * 2015-09-03 2017-03-09 Schlumberger Technology Corporation Method of integrating fracture, production, and reservoir operations into geomechanical operations of a wellsite
CN107480383A (zh) * 2017-08-21 2017-12-15 中国石油大学(北京) 一种通过测压数据监测注水动态裂缝的方法
CN107578471A (zh) * 2017-09-18 2018-01-12 西南石油大学 一种自支撑裂缝初始形态构建方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7788037B2 (en) * 2005-01-08 2010-08-31 Halliburton Energy Services, Inc. Method and system for determining formation properties based on fracture treatment
RU2324810C2 (ru) * 2006-05-31 2008-05-20 Шлюмберже Текнолоджи Б.В. Способ определения размеров трещины гидроразрыва пласта
CA2808858C (en) * 2012-03-16 2016-01-26 Weatherford/Lamb, Inc. Wellbore real-time monitoring and analysis of fracture contribution
CN105350960B (zh) * 2015-12-07 2018-11-23 西南石油大学 确定低渗透非均质气藏压裂水平井裂缝参数的方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0456424A2 (en) * 1990-05-07 1991-11-13 Halliburton Company Method of determining fracture characteristics of subsurface formations
EP0589591A1 (en) * 1992-09-11 1994-03-30 Halliburton Company Downhole fracture test and analysis
US6101447A (en) * 1998-02-12 2000-08-08 Schlumberger Technology Corporation Oil and gas reservoir production analysis apparatus and method
WO2003102371A1 (en) * 2002-05-31 2003-12-11 Schlumberger Canada Limited Method and apparatus for effective well and reservoir evaluation without the need for well pressure history
CN101560879A (zh) * 2008-04-15 2009-10-21 中国石油大学(北京) 用于低渗透气藏的试井分析控制系统及方法
CN104612658A (zh) * 2014-11-25 2015-05-13 中国石油大学(北京) 用于水平井分段产液识别的试井分析控制系统及方法
CN104594872A (zh) * 2015-01-04 2015-05-06 西南石油大学 一种优化致密气藏压裂水平井裂缝导流能力的方法
WO2017041074A1 (en) * 2015-09-03 2017-03-09 Schlumberger Technology Corporation Method of integrating fracture, production, and reservoir operations into geomechanical operations of a wellsite
CN105298479A (zh) * 2015-09-09 2016-02-03 中国石油大学(北京) 压裂直井产油(气)位置的诊断方法及其系统
CN106353809A (zh) * 2016-08-11 2017-01-25 中国石油大学(华东) 一种压裂裂缝网络的反演表征方法
CN107480383A (zh) * 2017-08-21 2017-12-15 中国石油大学(北京) 一种通过测压数据监测注水动态裂缝的方法
CN107578471A (zh) * 2017-09-18 2018-01-12 西南石油大学 一种自支撑裂缝初始形态构建方法

Non-Patent Citations (14)

* Cited by examiner, † Cited by third party
Title
Numerical simulation of fracture path and nonlinear closure for simultaneous and sequential fracturing in a horizontal well;WanCheng et al.;《Computers and Geotechnics》;20170831;第88卷;第242-255页 *
Patrick W et al..Reservoir Simulation Model for Fractured and Partially Fractured Reservoirs based on PEBI Grids .《Society of Petroleum Engineers》.2001,全文. *
The use of source and green’s functions in solving unsteady-flow problems in reservoirs;Gringarten A C et al.;《Society of Petroleum Engineers》;19731231;第1-12页 *
youwei he et al..A Semianalytical Methodology To Diagnose the Locations of Underperforming Hydraulic Fractures Through Pressure-Transient Analysis in Tight Gas Reservoir .《Society of Petroleum Engineers》.2016,第22卷(第3期),第924–939页. *
不对称垂直裂缝井产量递减规律;熊健等;《西安石油大学学报(自然科学版)》;20140331;第29卷(第2期);第74-77页 *
低渗透储层产量递减模型的渗流机理及应用;陈新彬等;《石油学报》;20110131;第32卷(第1期);第113-116页 *
低渗透储层生产数据分析技术及其应用研究;王宁博等;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20140815(第8期);第26-35页 *
压裂气井产能预测新方法研究;唐帅等;《石油天然气学报》;20120228;第34卷(第2期);第148-151页 *
基于数值模型的气井现代产量递减分析及动态预测;孙贺东等;《石油学报》;20171030;第38卷(第10期);第1194-1199页 *
多段压裂水平井不均匀产油试井模型;何佑伟等;《中国石油大学学报(自然科学版)》;20170830;第41卷(第4期);第116-123页 *
王宁博等.低渗透储层生产数据分析技术及其应用研究.《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》.2014,(第8期),第B019-132页. *
致密气藏动态分析方法及软件研制;赵学峰等;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20111115(第11期);第B019-80页 *
致密气藏水平井压裂缝不均匀产气试井分析;严谨等;《天然气地球科学》;20170630;第28卷(第6期);第839-845页 *
过量顶替液作业下压裂水平气井的产能模拟;严向阳等;《岩性油气藏》;20070228;第29卷(第1期);第140-146页 *

Also Published As

Publication number Publication date
CN110318742A (zh) 2019-10-11

Similar Documents

Publication Publication Date Title
CN109594968B (zh) 一种页岩气多段压裂水平井压后裂缝参数评价方法及系统
CN107301306B (zh) 用于致密砂岩气藏压裂水平井的动态无阻流量预测方法
US9568642B1 (en) Permeability and inflow performance determination for horizontal wells
CN107462936B (zh) 利用压力监测资料反演低渗透储层非达西渗流规律的方法
US10018029B2 (en) Method and device using productivity index in drill guidance for drilling slanted water injection wells
US9348058B2 (en) Method for determining the profile of an inflow and the parameters of a well-surrounding area in a multipay well
RU2455482C2 (ru) Способ определения профиля притока флюидов и параметров околоскважинного пространства
CN107506948A (zh) 一种基于动态泄流体积的页岩油气综合产量分析方法
CN106194154B (zh) 一种非常规油气藏中长期产能预测方法
CN104101564A (zh) 一种非稳态高温高压测试低渗透岩心启动压力梯度的装置及方法
US10208548B2 (en) Method for detecting gain or loss of drilling fluid in a drilling installation associated calculation system and associated drilling installation
CN104975852A (zh) 多段压裂水平井产油(气)位置的诊断方法及其系统
CN111691869B (zh) 一种高含水致密气藏气井产量变化趋势的预测方法
CN104504604A (zh) 一种定性气井井筒积液的方法
RU2460878C2 (ru) Способ определения профиля притока флюидов и параметров околоскважинного пространства
Li et al. Effect of adsorption and permeability correction on transient pressures in organic rich gas reservoirs: Vertical and hydraulically fractured horizontal wells
CN110318742B (zh) 基于压裂井生产数据确定裂缝闭合长度的方法和系统
CN111734394B (zh) 一种确定致密油藏压裂井不定常流井底压力的方法
CN114595504A (zh) 一种页岩气藏压裂水平井产出剖面反演解释方法
CN110630243B (zh) 基于压裂井生产数据确定压裂裂缝参数的方法
CN111963149B (zh) 一种考虑滞地液量增压的压裂后地层压力求取方法
CN112257349B (zh) 一种判断致密砂岩可动水气藏气井是否具有开发价值的方法
US11193370B1 (en) Systems and methods for transient testing of hydrocarbon wells
CN112878987A (zh) 利用生产数据计算页岩气井控制储量的方法
RU2569522C1 (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