CN113553746A - 快速诊断缝洞油藏储层参数的方法和处理器 - Google Patents
快速诊断缝洞油藏储层参数的方法和处理器 Download PDFInfo
- Publication number
- CN113553746A CN113553746A CN202110950790.XA CN202110950790A CN113553746A CN 113553746 A CN113553746 A CN 113553746A CN 202110950790 A CN202110950790 A CN 202110950790A CN 113553746 A CN113553746 A CN 113553746A
- Authority
- CN
- China
- Prior art keywords
- fracture
- reservoir
- curve
- coincidence rate
- parameters
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 206010017076 Fracture Diseases 0.000 claims description 173
- 208000010392 Bone Fractures Diseases 0.000 claims description 166
- 238000012937 correction Methods 0.000 claims description 14
- 238000011156 evaluation Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 14
- 230000035699 permeability Effects 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000003745 diagnosis Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 11
- 238000004458 analytical method Methods 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 62
- 239000003921 oil Substances 0.000 description 47
- 238000010586 diagram Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 7
- 210000003491 skin Anatomy 0.000 description 7
- 238000004519 manufacturing process Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000005465 channeling Effects 0.000 description 5
- 239000012530 fluid Substances 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 239000000126 substance Substances 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000009530 blood pressure measurement Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 239000010779 crude oil Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 239000003129 oil well Substances 0.000 description 2
- 230000003287 optical effect Effects 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
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 210000002615 epidermis Anatomy 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009420 retrofitting Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- 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)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Agronomy & Crop Science (AREA)
- Algebra (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Fluid Mechanics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computing Systems (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
本发明公开了一种快速诊断缝洞油藏储层参数的方法、处理器和机器可读存储介质，该方法包括：获取缝洞油藏无限导流裂缝模型的第一早期近似解和第一解析解；根据第一早期近似解构建第一导流能力影响函数；根据第一解析解、第一导流能力影响函数构建有限导流裂缝解析解模型，以获取缝洞油藏压裂井井底的理论压力数据；获取缝洞油藏压裂井井底的实际压力数据；根据实际压力数据和理论压力数据确定缝洞油藏储层参数的求取值，简明易懂，操作简单，大幅度减少了计算时间，提高了缝洞油藏储层参数的解析效率。
Description
技术领域
本发明涉及油藏工程领域，具体地涉及一种快速诊断缝洞油藏储层参数的方法、处理器和机器可读存储介质。
背景技术
缝洞油藏的探明储量占全球探明储量的60％以上，水力压裂已经成为油藏开发中增加单井产量、提高开发效果的主流技术，每年油田现场有大量的井需要压裂后投产以获得最大的油气生产效益。因此，对缝洞油藏中压裂直井的压力监测、评价诊断是当前油气田开发的一个热点问题。通过对压裂井压力监测得到的数据进行分析，达到获取地层信息，诊断地层的目的，为开发方案设计，开发方案调整和增产措施的施行提供技术支撑。
试井是目前最常用于分析测压数据的技术，利用渗流力学理论建立符合油藏地质情况的模型，然后将该模型与实际压力数据进行拟合，最终达到反演储层参数的目的。然而，现有压裂井试井模型主要是半解析和数值模型，需要对裂缝和地层进行离散才能获得模型的解，这样造成模型的计算速度慢、得到的计算结果精度低。另外，在计算过程中常常出现结果不收敛的情况，非常不便于油田工程师使用，所以有必要提出一种快速、简单的通过压裂井测压数据诊断缝洞油藏储层参数的方法。
发明内容
本发明的目的是为了提供一种快速诊断缝洞油藏储层参数的方法、处理器和机器可读存储介质，该快速诊断缝洞油藏储层参数的方法、处理器和机器可读存储介质具有可操作性强，能快速、简单地诊断出缝洞油藏储层参数。
为了实现上述目的，本发明第一方面提供一种快速诊断缝洞油藏储层参数的方法，该方法包括：
获取缝洞油藏无限导流裂缝模型的第一早期近似解和第一解析解；
根据第一早期近似解构建第一导流能力影响函数；
根据第一解析解、第一导流能力影响函数构建有限导流裂缝解析解模型，以获取缝洞油藏压裂井井底的理论压力数据；
获取缝洞油藏压裂井井底的实际压力数据；
根据实际压力数据和理论压力数据确定缝洞油藏储层参数的求取值。
在本发明的实施例中，缝洞油藏储层参数至少包括井筒储集系数、表皮系数、地层渗透率、裂缝半长、裂缝储容比以及基质窜流系数。
在本发明的实施例中，根据第一早期近似解构建第一导流能力影响函数包括：
获取缝洞油藏有限导流裂缝近似解模型以及有限导流裂缝近似解模型的第二早期近似解；
基于第二早期近似解和第一早期近似解的差值构建第一导流能力影响函数。
在本发明的实施例中，根据第一解析解、第一导流能力影响函数构建有限导流裂缝解析解模型，以获取缝洞油藏压裂井井底的理论压力数据包括：
根据第一导流能力影响函数构建第二导流能力影响函数；
对第一解析解和第二导流能力影响函数进行叠加操作，以构建有限导流裂缝解析解模型；
对缝洞油藏储层参数赋予初始值，将初始值输入有限导流裂缝解析解模型中，以获取缝洞油藏压裂井井底的理论压力数据。
在本发明的实施例中，根据第一导流能力影响函数构建第二导流能力影响函数包括：
获取修正项函数，其中，修正项函数中包括拉普拉斯变量、窜流函数和无量纲裂缝导流系数；
对第一导流能力影响函数和修正项函数进行叠加操作，以构建第二导流能力影响函数。
在本发明的实施例中，根据实际压力数据和理论压力数据确定缝洞油藏储层参数的求取值包括：
根据实际压力数据绘制实际压力双对数曲线，其中，实际压力双对数曲线包括实际压力值曲线和实际压力值导数曲线；
根据理论压力数据绘制理论压力双对数曲线，其中，理论压力双对数曲线包括理论压力值曲线和理论压力值导数曲线；
根据实际压力双对数曲线和理论压力双对数曲线确定缝洞油藏储层参数的求取值。
在本发明的实施例中，根据实际压力双对数曲线和理论压力双对数曲线确定缝洞油藏储层参数的求取值包括：
计算实际压力值曲线、理论压力值曲线之间的第一重合率和实际压力值导数曲线、理论压力值导数曲线之间的第二重合率，根据第一重合率和第二重合率确定缝洞油藏储层参数的求取值。
在本发明的实施例中，计算实际压力值曲线、理论压力值曲线之间的第一重合率和实际压力值导数曲线、理论压力值导数曲线之间的第二重合率，根据第一重合率和第二重合率确定缝洞油藏储层参数的求取值包括：
将第一重合率和第一预设重合率范围进行对比，以判断第一重合率是否在第一预设重合率范围内；
将第二重合率和第二预设重合率范围进行对比，以判断第二重合率是否在第二预设重合率范围内；
在确定第一重合率在第一预设重合率范围内，且第二重合率在第二预设重合率范围内的情况下，初始值为缝洞油藏储层参数的求取值。
本发明的第二方面提供一种处理器，该处理器被配置成执行上述的快速诊断缝洞油藏储层参数的方法。
本发明的第三方面提供一种机器可读存储介质，该机器可读存储介质上存储有指令，指令用于使得机器执行上述的快速诊断缝洞油藏储层参数的方法。
通过上述技术方案，获取缝洞油藏压裂井井底的实际压力数据，通过无限导流裂缝模型的第一早期近似解获取到第一导流能力影响函数，再由第一解析解、第一导流能力影响函数构建有限导流裂缝解析解模型，该模型公式更加简洁易懂，并计算出缝洞油藏压裂井井底的理论压力数据，再根据实际压力数据和理论压力数据确定缝洞油藏储层参数，该方法简明易懂，操作简单，大幅度减少了计算时间，提高了储层参数的解析效率，解决了现有技术中诊断缝洞油藏地层参数计算速度慢、效率低、计算精度差以及不容易收敛的问题，为之后缝洞油藏压裂井测压数据解释、开发方案设计和增产改造措施的施行提供有利技术支撑。
附图说明
附图是用来提供对本发明实施例的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明实施例，但并不构成对本发明实施例的限制。在附图中：
图1是本发明实施例中快速诊断缝洞油藏储层参数的方法的流程示意图；
图2是本发明实施例中缝洞油藏压裂井物理模型；
图3是本发明实施例中流动关系示意图；
图4是本发明实施例中理论压力值与理论压力值导数双对数曲线；
图5是本发明实施例中理论双对数曲线图与实际双对数曲线拟合结果示意图。
附图标记说明
1 理论压力值曲线 2 理论压力值导数曲线
3 实际压力值曲线 4 实际压力值导数曲线
具体实施方式
以下结合附图对本发明实施例的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明实施例，并不用于限制本发明实施例。
需要说明，若本申请实施方式中有涉及方向性指示(诸如上、下、左、右、前、后……)，则该方向性指示仅用于解释在某一特定姿态下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
另外，若本申请实施方式中有涉及“第一”、“第二”等的描述，则该“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外，各个实施方式之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本申请要求的保护范围之内。
本发明的一个实施例中提供一种新型的快速诊断缝洞油藏储层参数的方法，适用于快速诊断缝洞油藏储层参数的控制装置，具体地，如图1所示，该方法包括如下步骤：
步骤S101：构建缝洞油藏无限导流裂缝模型并获取无限导流裂缝模型的第一早期近似解和第一解析解。
可以理解，本实施例中使用经典的Warren-Root模型来表示缝洞油藏，即如图2所示，本实施例中基于经典的Warren-Root模型提出缝洞油藏压裂直井的物理模型，为之后构建缝洞油藏的有限导流裂缝解析解模型奠定基础，本实施例中用于快速诊断缝洞油藏储层参数的控制装置包括处理器和存储器，存储器和处理器电连接，上述缝洞油藏压裂直井的物理模型预存在存储器中，在需要时可通过处理器调取出来。缝洞油藏压裂直井的物理模型中涉及的流动关系阐述如下：缝洞油藏基质中的流体主要流向天然裂缝，天然裂缝中的流体主要流向人工裂缝，但是基质和天然裂缝不直接与井筒相连，储存在基质和天然裂缝中的流体必须通过人工裂缝流入井筒而被开采出来，流体流动关系示意图如图3所示。
在本发明的实施例中，步骤S101，获取缝洞油藏无限导流裂缝模型的第一早期近似解和第一解析解包括如下步骤：
步骤S201：根据以下公式获取无限导流裂缝模型的第一早期近似解；
其中，
为拉普拉斯空间下缝洞油藏无限导流裂缝模型的第一早期近似解；s为拉普拉斯变量，s＝ln(2)/tD，tD为无量纲时间；上述公式(1)预存在存储器中，在需要时可通过处理器调取出来。
步骤S202：获取积分变量、窜流函数和贝塞尔函数；
步骤S203：根据以下公式获取无限导流裂缝模型的第一解析解：
其中，
为拉普拉斯空间下缝洞油藏无限导流裂缝模型的第一解析解；K0为贝塞尔函数；f(s)为窜流函数；β为积分变量；积分变量、窜流函数、贝塞尔函数和上述公式(2)均预存在存储器中，在需要时可通过处理器调取出来。
本实施例中，获取窜流函数包括如下步骤：
获取基质窜流系数、裂缝储容比；
根据以下公式获得窜流函数：
其中，f(s)为窜流函数；ωf为裂缝储容比，单位为m；λm为基质窜流系数，上述公式(3)预存在存储器中，在需要时可通过处理器调取出来。
步骤S102：根据第一早期近似解构建第一导流能力影响函数。
可以理解，在本发明的实施例中，步骤S102，根据第一早期近似解构建第一导流能力影响函数包括如下步骤：
步骤S301：获取缝洞油藏有限导流裂缝近似解模型以及有限导流裂缝近似解模型的第二早期近似解。
可以理解，在本发明的实施例中，步骤S301，获取缝洞油藏有限导流裂缝近似解模型以及有限导流裂缝近似解模型的第二早期近似解包括如下步骤：
步骤S401：获取无量纲裂缝导流系数。
本实施例中，获取无量纲裂缝导流系数包括：
获取裂缝渗透率、裂缝宽度、地层渗透率和裂缝半长；
根据以下公式计算获得无量纲裂缝导流系数：
其中，FcD为无量纲裂缝导流系数；kf为裂缝渗透率，单位为mD；wf为裂缝宽度，单位为m；kfb为地层渗透率，单位为mD；xf为裂缝半长，单位为m。进一步地，本实施例中裂缝渗透率、裂缝宽度、地层渗透率和裂缝半长均预先存储在存储器中。
步骤S402：根据以下公式获取有限导流裂缝近似解模型的第二早期近似解：
其中，
为有限导流裂缝近似解模型的第二早期近似解，上述公式(5)预存在存储器中，在需要时可通过处理器调取出来。
步骤S302：基于第二早期近似解和第一早期近似解的差值构建第一导流能力影响函数。
可以理解，在本发明的实施例中，利用缝洞油藏有限导流裂缝模型的第二早期近似解减去无限导流裂缝模型的第一早期近似解即可得到第一导流能力影响函数，第一导流能力影响函数如以下公式所示：
其中，fFD1为第一导流能力影响函数。
步骤S103：根据第一解析解、第一导流能力影响函数构建有限导流裂缝解析解模型，以获取缝洞油藏压裂井井底的理论压力数据。
可以理解，在本发明的一个实施例中，步骤S103：根据第一解析解、第一导流能力影响函数构建有限导流裂缝解析解模型，以获取缝洞油藏压裂井井底的理论压力数据包括步骤S501-S503，步骤S501-S503如下所示：
步骤S501：根据第一导流能力影响函数构建第二导流能力影响函数。
步骤S501中，根据第一导流能力影响函数构建第二导流能力影响函数包括步骤S601-S602，S601-S602如下所示：
步骤S601：获取修正项函数，其中，修正项函数中包括拉普拉斯变量、窜流函数和无量纲裂缝导流系数；
步骤S602：对第一导流能力影响函数和修正项函数进行叠加操作，以构建第二导流能力影响函数。
具体地，处理器根据以下公式获取到修正项函数：
其中，Y为修正项函数，上述公式(7)预存在存储器中，在需要时可通过处理器调取出来。
处理器在获得修正项函数后，将第一导流能力影响函数和修正项函数进行叠加操作，即可获得第二导流能力影响函数，第二导流能力影响函数如下式所示：
其中，fFD2为第二导流能力影响函数。
步骤S502：对第一解析解和第二导流能力影响函数进行叠加操作，以构建有限导流裂缝解析解模型。
处理器获取到第二导流能力影响函数后，将第一解析解和第二导流能力影响函数进行叠加操作，即可完成有限导流裂缝解析解模型的构建，有限导流裂缝解析解模型如以下公式所示：
其中，
为有限导流裂缝解析解模型。
步骤S503：对缝洞油藏储层参数赋予初始值，将初始值输入有限导流裂缝解析解模型中，以获取缝洞油藏压裂井井底的理论压力数据。
可以理解，缝洞油藏储层参数至少包括井筒储集系数、表皮系数、地层渗透率、裂缝半长、裂缝储容比以及基质窜流系数。在通过公式(9)获取获取缝洞油藏压裂直井井底的理论压力数据前，需要先确定待求的缝洞油藏储层参数，本实施例中确定为井筒储集系数、表皮系数、地层渗透率、裂缝半长、裂缝储容比以及基质窜流系数7个参数，在确定待求的缝洞油藏储层参数后，再根据测井、完井、地震和岩芯等地质资料为上述待求的各个缝洞油藏储层参数赋予初始值，以使得公式(9)的计算结果更符合地层的实际情况，并减少后续拟合过程的修正次数，提高了解析的效率。例如，本实施例中预设井筒储集系数被赋予的初始值为第一初始值，表皮系数被赋予的初始值为第二初始值，地层渗透率被赋予的初始值为第三初始值，裂缝长度被赋予的初始值为第四初始值，无量纲裂缝导流系数被赋予的初始值为第五初始值，基质窜流系数被赋予的初始值为第六初始值，裂缝储容比被赋予的初始值为第七初始值，上述各个初始值均被预先存储在存储器中，且第一初始值为0.01m3/MPa，第二初始值为0.03，第三初始值为1mD，第四初始值为45m，第五初始值为20，第六初始值为0.03，第七初始值为0.3，之后将上述各个初始值代入公式(9)中，以计算出有限导流裂缝模型的解析解(即缝洞油藏压裂井井底的理论压力数据)。
本实施例中的理论压力数据包括理论压力值和理论压力值导数，又由于在实际生产过程中，井筒存在储集效应和井底污染，因此处理器在基于公式(9)得到的有限导流裂缝模型的解析解后，可再利用杜哈美原理获得考虑井筒储集效应和表皮系数时缝洞油藏压裂直井井底的理论压力值。具体地，理论压力值可根据如下公式获得：
其中，
为拉普拉斯空间下考虑井储表皮的压力降，单位为MPa；
为拉普拉斯空间下地层原始压力，单位为MPa；
为拉普拉斯空间下考虑井储表皮的井底压力(即本实施例中的理论压力值)，进一步的，该
为拉普拉斯空间下考虑井储表皮的井底压力随时间变化的数据，单位为MPa；q为油井产量，m3/d；μ为原油粘度，单位为mPa.s；B为原油体积系数；h为储层厚度，单位为m；C为井筒储集系数，单位为m3/MPa；S为表皮系数；φf为裂缝介质孔隙度；ctf为裂缝介质综合压缩系数，单位为MPa-1；φm为基质孔隙度；ctm为基质综合压缩系数，单位为MPa-1。进一步地，本实施例中裂缝介质孔隙度、裂缝介质综合压缩系数、基质孔隙度和基质综合压缩系数以及公式(10)均预先存储在存储器中，且ctf＝ctm＝4×10- 4MPa-1，φf+φm＝0.1。
步骤S104：获取缝洞油藏压裂井井底的实际压力数据。
可以理解，本实施例中用于快速诊断缝洞油藏储层参数的控制装置包括压力计和处理器，其中，压力计和处理器电连接。在对缝洞油藏储层参数进行诊断之前，操作人员需要先选取缝洞油藏中一口用于测试的压裂直井，关井后将压力计下入油井产层中部，并用封隔器封隔压力计与上部流体，以便实时检测关井后压裂直井井底的实际压力值，压力计在检测完成后将实际压力值传送给处理器，以便处理器基于该实际压力值获取到缝洞油藏压裂井井底的实际压力数据。本实施例中的实际压力数据包括实际压力值和实际压力值导数，处理器对实际压力值进行计算处理即可得到实际压力值导数。
步骤S105：根据实际压力数据和理论压力数据确定缝洞油藏储层参数的求取值。
可以理解，步骤S105中，根据实际压力数据和理论压力数据确定缝洞油藏储层参数的求取值包括如下步骤：
步骤S601：根据实际压力数据绘制实际压力双对数曲线，其中，实际压力双对数曲线包括实际压力值曲线3和实际压力值导数曲线4。
可以理解，处理器在获得实际压力值后，计算实际压力值对对数时间的导数，以获得实际压力值导数，之后处理器再将实际压力值和实际压力值导数同时绘制在双对数图版上，以获得缝洞油藏压裂直井试井双对数实际曲线图版，该缝洞油藏压裂井双对数实际曲线图版如图4所述。
步骤S602：根据理论压力数据绘制理论压力双对数曲线，其中，理论压力双对数曲线包括理论压力值曲线1和理论压力值导数曲线2。
可以理解，处理器在获得理论压力值后，对该理论压力值进行Stehfest数值反演，以进一步获得实空间下井底压力随时间变化的数据，再求取上述数据对对数时间的导数即可获得理论压力值导数，之后处理器再将理论压力值和理论压力值导数同时绘制在双对数图版上，以获得缝洞油藏压裂直井试井双对数理论曲线图版，该缝洞油藏压裂井双对数理论曲线图版如图4所述。
步骤S603：根据实际压力双对数曲线和理论压力双对数曲线确定缝洞油藏储层参数的求取值。
可以理解，处理器在获取到实际压力双对数曲线和理论压力双对数曲线后，对上述的两对曲线进行拟合，如图5所示，即实际压力值曲线3和理论压力值曲线1进行拟合，实际压力值导数曲线4和理论压力值导数曲线2进行拟合，以便根据拟合的结果来确定缝洞油藏储层参数的求取值，其中，拟合算法包括最小二乘法或遗传算法中的至少一种。
在本发明的一个实施例中，根据实际压力双对数曲线和理论压力双对数曲线确定缝洞油藏储层参数的求取值包括：
计算实际压力值曲线、理论压力值曲线1之间的第一重合率和实际压力值导数曲线、理论压力值导数曲线2之间的第二重合率，根据第一重合率和第二重合率确定缝洞油藏储层参数的求取值。
可以理解，处理器在对实际压力双对数曲线和理论压力双对数曲线进行拟合后，要计算实际压力双对数曲线和理论压力双对数曲线之间的重合率，即处理器要分别计算出实际压力值曲线3和理论压力值曲线1拟合后的第一重合率，以及实际压力值导数曲线4和理论压力值导数曲线2拟合后的第二重合率，以便根据第一重合率和第二重合率确定缝洞油藏储层参数的求取值。
在本发明的一个实施例中，计算实际压力值曲线3、理论压力值曲线1之间的第一重合率和实际压力值导数曲线4、理论压力值导数曲线2之间的第二重合率，根据第一重合率和第二重合率确定缝洞油藏储层参数的求取值包括：
将第一重合率和第一预设重合率范围进行对比，以判断第一重合率是否在第一预设重合率范围内；
将第二重合率和第二预设重合率范围进行对比，以判断第二重合率是否在第二预设重合率范围内；
在确定第一重合率在第一预设重合率范围内，且第二重合率在第二预设重合率范围内的情况下，初始值为缝洞油藏储层参数的求取值。
可以理解，处理器在获得第一重合率和第二重合率后，需要判断第一重合率是否在第一预设重合率范围内，以及第二重合率是否在第二预设重合率范围内，其中，第一预设重合率范围为0.7-1；第二预设重合率范围为0.7-1。若第一重合率在第一预设重合率范围内，且第二重合率在第二预设重合率范围内，则说明实际压力双对数曲线和理论压力双对数曲线达到了拟合成功的标准，进而确定缝洞油藏储层参数的各个初始值即为其对应的各个求取值；若第一重合率不在第一预设重合率范围内，和/或第二重合率不在第二预设重合率范围内，则说明实际压力双对数曲线和理论压力双对数曲线没有达到拟合成功的标准，进而确定缝洞油藏储层参数的各个初始值不是其对应的各个求取值，此时处理器需要对缝洞油藏储层参数的各个初始值重新进行赋值，并在赋值后继续执行后续步骤，直到实际压力双对数曲线和理论压力双对数曲线达到拟合成功的标准(即第一重合率在第一预设重合率范围内，且第二重合率在第二预设重合率范围内)，进而确定缝洞油藏储层参数被重新赋值的各个初始值是其对应的各个求取值。本发明的实施例中，在第一初始值重新赋值为0.03m3/MPa，第二初始值重新赋值为0.05m3/MPa，第三初始值重新赋值为2mD，第四初始值重新赋值为50m，第五初始值重新赋值为20，第六初始值赋值为0.028，第七初始值重新赋值为0.2时，实际压力双对数曲线和理论压力双对数曲线才达到了拟合成功的标准。
本发明实施例中提供的快速诊断缝洞油藏储层参数的方法，获取缝洞油藏压裂井井底的实际压力数据，通过无限导流裂缝模型的第一早期近似解获取到第一导流能力影响函数，再由第一解析解、第一导流能力影响函数构建有限导流裂缝解析解模型，该模型公式更加简洁易懂，并计算出缝洞油藏压裂井井底的理论压力数据，再根据实际压力数据和理论压力数据确定缝洞油藏储层参数，该方法简明易懂，操作简单，大幅度减少了计算时间，提高了储层参数的解析效率，解决了现有技术中诊断缝洞油藏地层参数计算速度慢、效率低、计算精度差以及不容易收敛的问题，为之后缝洞油藏压裂井测压数据解释、开发方案设计和增产改造措施的施行提供有利技术支撑。
本发明的另一个实施例中提供一种处理器，该处理器被配置成执行上述的快速诊断缝洞油藏储层参数的方法。
本发明的另一个实施例中提供一种机器可读存储介质，该机器可读存储介质上存储有指令，指令用于使得机器执行上述的快速诊断缝洞油藏储层参数的方法。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (10)
1.一种快速诊断缝洞油藏储层参数的方法，其特征在于，所述方法包括：
获取缝洞油藏无限导流裂缝模型的第一早期近似解和第一解析解；
根据所述第一早期近似解构建第一导流能力影响函数；
根据所述第一解析解、所述第一导流能力影响函数构建有限导流裂缝解析解模型，以获取缝洞油藏压裂井井底的理论压力数据；
获取所述缝洞油藏压裂井井底的实际压力数据；
根据所述实际压力数据和所述理论压力数据确定所述缝洞油藏储层参数的求取值。
2.根据权利要求1所述的快速诊断缝洞油藏储层参数的方法，其特征在于，所述缝洞油藏储层参数至少包括井筒储集系数、表皮系数、地层渗透率、裂缝半长、裂缝储容比以及基质窜流系数。
3.根据权利要求1所述的快速诊断缝洞油藏储层参数的方法，其特征在于，所述根据所述第一早期近似解构建第一导流能力影响函数包括：
获取缝洞油藏有限导流裂缝近似解模型以及所述有限导流裂缝近似解模型的第二早期近似解；
基于所述第二早期近似解和所述第一早期近似解的差值构建所述第一导流能力影响函数。
4.根据权利要求2所述的快速诊断缝洞油藏储层参数的方法，其特征在于，所述根据所述第一解析解、所述第一导流能力影响函数构建有限导流裂缝解析解模型，以获取所述缝洞油藏压裂井井底的理论压力数据包括：
根据所述第一导流能力影响函数构建第二导流能力影响函数；
对所述第一解析解和所述第二导流能力影响函数进行叠加操作，以构建所述有限导流裂缝解析解模型；
对所述缝洞油藏储层参数赋予初始值，将所述初始值输入所述有限导流裂缝解析解模型中，以获取所述缝洞油藏压裂井井底的理论压力数据。
5.根据权利要求4所述的快速诊断缝洞油藏储层参数的方法，其特征在于，所述根据所述第一导流能力影响函数构建第二导流能力影响函数包括：
获取修正项函数，其中，所述修正项函数中包括拉普拉斯变量、窜流函数和无量纲裂缝导流系数；
对所述第一导流能力影响函数和所述修正项函数进行叠加操作，以构建所述第二导流能力影响函数。
6.根据权利要求1所述的快速诊断缝洞油藏储层参数的方法，其特征在于，所述根据所述实际压力数据和所述理论压力数据确定所述缝洞油藏储层参数的求取值包括：
根据所述实际压力数据绘制实际压力双对数曲线，其中，所述实际压力双对数曲线包括实际压力值曲线和实际压力值导数曲线；
根据所述理论压力数据绘制理论压力双对数曲线，其中，所述理论压力双对数曲线包括理论压力值曲线和理论压力值导数曲线；
根据所述实际压力双对数曲线和所述理论压力双对数曲线确定所述缝洞油藏储层参数的求取值。
7.根据权利要求6所述的快速诊断缝洞油藏储层参数的方法，其特征在于，所述根据所述实际压力双对数曲线和所述理论压力双对数曲线确定所述缝洞油藏储层参数的求取值包括：
计算所述实际压力值曲线、所述理论压力值曲线之间的第一重合率和所述实际压力值导数曲线、所述理论压力值导数曲线之间的第二重合率，根据所述第一重合率和所述第二重合率确定所述缝洞油藏储层参数的求取值。
8.根据权利要求7所述的快速诊断缝洞油藏储层参数的方法，其特征在于，所述计算所述实际压力值曲线、所述理论压力值曲线之间的第一重合率和所述实际压力值导数曲线、所述理论压力值导数曲线之间的第二重合率，根据所述第一重合率和所述第二重合率确定所述缝洞油藏储层参数的求取值包括：
将所述第一重合率和所述第一预设重合率范围进行对比，以判断所述第一重合率是否在所述第一预设重合率范围内；
将所述第二重合率和所述第二预设重合率范围进行对比，以判断所述第二重合率是否在所述第二预设重合率范围内；
在确定所述第一重合率在所述第一预设重合率范围内，且所述第二重合率在所述第二预设重合率范围内的情况下，所述初始值为所述缝洞油藏储层参数的求取值。
9.一种处理器，其特征在于，所述处理器被配置成执行根据权利要求1至8中任意一项所述的快速诊断缝洞油藏储层参数的方法。
10.一种机器可读存储介质，所述机器可读存储介质上存储有指令，其特征在于，所述指令用于使得机器执行根据权利要求1至8任一项所述的快速诊断缝洞油藏储层参数的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110950790.XA CN113553746B (zh) | 2021-08-18 | 2021-08-18 | 快速诊断缝洞油藏储层参数的方法和处理器 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110950790.XA CN113553746B (zh) | 2021-08-18 | 2021-08-18 | 快速诊断缝洞油藏储层参数的方法和处理器 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113553746A true CN113553746A (zh) | 2021-10-26 |
| CN113553746B CN113553746B (zh) | 2024-03-29 |
Family
ID=78105726
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110950790.XA Active CN113553746B (zh) | 2021-08-18 | 2021-08-18 | 快速诊断缝洞油藏储层参数的方法和处理器 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113553746B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117076956A (zh) * | 2023-10-16 | 2023-11-17 | 西安石油大学 | 缝洞油藏物理模型相似准数优化方法及装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050222852A1 (en) * | 2004-03-30 | 2005-10-06 | Craig David P | Method and an apparatus for detecting fracture with significant residual width from previous treatments |
| CN103590824A (zh) * | 2013-10-21 | 2014-02-19 | 中国石油天然气股份有限公司 | 经过多段压裂改造后的致密气藏水平井的产能计算方法 |
| CN107725034A (zh) * | 2017-08-21 | 2018-02-23 | 中国石油大学(北京) | 一种用于多级压裂水平井判别来水方向的压力监测方法 |
| CN111581786A (zh) * | 2020-04-19 | 2020-08-25 | 东北石油大学 | 用于分析缝洞串联模式双孔复合储层的试井解释模型 |
| CN112377183A (zh) * | 2020-11-30 | 2021-02-19 | 中国石油天然气股份有限公司 | 多洞多缝型碳酸盐岩储层产量分析模型构建方法及装置 |
| CN113111582A (zh) * | 2021-04-12 | 2021-07-13 | 中国石油大学(北京) | 基于深度学习的页岩储层试井智能解释分析方法及装置 |
-
2021
- 2021-08-18 CN CN202110950790.XA patent/CN113553746B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050222852A1 (en) * | 2004-03-30 | 2005-10-06 | Craig David P | Method and an apparatus for detecting fracture with significant residual width from previous treatments |
| CN103590824A (zh) * | 2013-10-21 | 2014-02-19 | 中国石油天然气股份有限公司 | 经过多段压裂改造后的致密气藏水平井的产能计算方法 |
| CN107725034A (zh) * | 2017-08-21 | 2018-02-23 | 中国石油大学(北京) | 一种用于多级压裂水平井判别来水方向的压力监测方法 |
| CN111581786A (zh) * | 2020-04-19 | 2020-08-25 | 东北石油大学 | 用于分析缝洞串联模式双孔复合储层的试井解释模型 |
| CN112377183A (zh) * | 2020-11-30 | 2021-02-19 | 中国石油天然气股份有限公司 | 多洞多缝型碳酸盐岩储层产量分析模型构建方法及装置 |
| CN113111582A (zh) * | 2021-04-12 | 2021-07-13 | 中国石油大学(北京) | 基于深度学习的页岩储层试井智能解释分析方法及装置 |
Non-Patent Citations (7)
| Title |
|---|
| CAO WEI等: "A hybrid analytic solution for a well with a finite-conductivity vertical fracture", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》, vol. 188, pages 1 - 14 * |
| CAO WEI等: "A quick analysis approach for early time well test data of vertically fractured wells:Applications in Changqing oilfield,china", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》, vol. 201, pages 1 - 13 * |
| MINGXIAN WANG等: "Rate Decline Analysis for Modeling Volume Fractured WEll Production in Naturally Fractured Reservoirs", 《ENERGIES》, pages 1 - 21 * |
| WANJING LUO等: "A simple and accurate calculation method for finite-conductivity fracture", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》, vol. 161, pages 590 - 598, XP085336631, DOI: 10.1016/j.petrol.2017.11.053 * |
| 曾杨等: "存在有限导流断层的油藏压力动态分析", 《天然气与石油》, vol. 36, no. 2, pages 83 - 88 * |
| 王西强等: "导流能力不均匀的部分压开垂直裂缝试井分析方法", 《科学技术与工程》, vol. 20, no. 19, pages 7644 - 7649 * |
| 王链: "致密气藏压裂井试井分析模型研究", 《中国优秀硕士学位论文全文数据库(工程科技Ⅰ辑)》, no. 8, pages 019 - 59 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117076956A (zh) * | 2023-10-16 | 2023-11-17 | 西安石油大学 | 缝洞油藏物理模型相似准数优化方法及装置 |
| CN117076956B (zh) * | 2023-10-16 | 2024-01-26 | 西安石油大学 | 缝洞油藏物理模型相似准数优化方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113553746B (zh) | 2024-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109441422B (zh) | 一种页岩气井间距优化开采方法 | |
| US11163923B2 (en) | Automated upscaling of relative permeability and capillary pressure in multi-porosity systems | |
| Schembre et al. | A technique for measuring two-phase relative permeability in porous media via X-ray CT measurements | |
| EP3108099B1 (en) | Measuring behind casing hydraulic conductivity between reservoir layers | |
| CN103982179A (zh) | 一种油藏储层的古压力定量反演探测方法 | |
| EP3245384A1 (en) | Measuring inter-reservoir cross flow rate between adjacent reservoir layers form transient pressure tests | |
| CN104535472A (zh) | 一种煤岩动态渗透率检测方法及装置 | |
| US20150062300A1 (en) | Wormhole Structure Digital Characterization and Stimulation | |
| EP3574184B1 (en) | Evaluation of pressure-transient behavior of wells | |
| CN110735636B (zh) | 多层合采井储层等效边界距离的测算方法、系统 | |
| Zhang et al. | High-order streamline simulation and macro-scale visualization experimental studies on waterflooding under given pressure boundaries | |
| CN113553746A (zh) | 快速诊断缝洞油藏储层参数的方法和处理器 | |
| CN113255092A (zh) | 体积压裂水平井特征点拟合反演方法、电子设备及介质 | |
| CN111577264A (zh) | 裂缝孔隙型油藏水平井产能预测方法及装置 | |
| CN105844011A (zh) | 一种基于毛管模型的渗透率计算方法 | |
| CN111582532B (zh) | 应力敏感性油藏水平井产液能力预测方法及装置 | |
| Jia et al. | A pore-network-modeling approach to predict petrophysical properties of diatomaceous reservoir rock | |
| CA2649483A1 (en) | Refined analytical model for formation parameter calculation | |
| CN114153007B (zh) | 水锁效应损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统 | |
| Liu et al. | Transient transfer shape factor for fractured tight reservoirs: Effect of the dynamic threshold pressure gradient in unsteady flow | |
| Schembre et al. | Direct measurement of dynamic relative permeability from CT monitored spontaneous imbibition experiments | |
| CN113673101B (zh) | 串珠状缝洞体溶洞体积的计算方法和处理器 | |
| CN111651848B (zh) | 应力敏感性油藏直井产液能力预测方法及装置 | |
| CN113705122B (zh) | 润湿反转损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统 | |
| Kanj et al. | Taming complexities of coupled geomechanics in rock testing: from assessing reservoir compaction to analyzing stability of expandable sand screens and solid tubulars |
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 |