CN112901144A - 一种特殊结构井井筒溢流动态特性分析方法 - Google Patents

一种特殊结构井井筒溢流动态特性分析方法 Download PDF

Info

Publication number
CN112901144A
CN112901144A CN202110266163.4A CN202110266163A CN112901144A CN 112901144 A CN112901144 A CN 112901144A CN 202110266163 A CN202110266163 A CN 202110266163A CN 112901144 A CN112901144 A CN 112901144A
Authority
CN
China
Prior art keywords
well
equation
target
gas
fluid flow
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.)
Pending
Application number
CN202110266163.4A
Other languages
English (en)
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 CN202110266163.4A priority Critical patent/CN112901144A/zh
Publication of CN112901144A publication Critical patent/CN112901144A/zh
Pending legal-status Critical Current

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
    • E21B47/00—Survey of boreholes or wells
    • 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/02—Determining slope or direction
    • E21B47/022—Determining slope or direction of the borehole, e.g. using geomagnetism
    • 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/06—Measuring temperature or pressure
    • 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
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Animal Husbandry (AREA)
  • Economics (AREA)
  • Agronomy & Crop Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Marine Sciences & Fisheries (AREA)
  • General Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种特殊结构井井筒溢流动态特性分析方法，涉及油气开发技术领域，包括以下步骤：S1：测井获得目标井的井深数据和井斜角数据，并根据插值计算方法获得目标井上任一点的井深和井斜角，获得连续平滑的井眼轨迹；S2：将目标井的井身分为若干节段，根据目标井内流体的质量守恒定律和动量守恒定律，建立每一节段的流体流动数值模型；S3：设定目标井的初始条件和边界条件，带入目标井的井眼轨迹，离散求解每一节段的流体流动数值模型，得到目标井井筒内的流体流动参数。本发明在特殊结构井的井筒溢流模型中引入了井斜角，使使其得到的井筒溢流模型更符合特殊结构井的实际溢流，提高井筒溢流模型对特殊结构井现场操作的指导参考价值。

Description

一种特殊结构井井筒溢流动态特性分析方法
技术领域
本发明涉及石油和天然气开发技术领域，特别涉及一种特殊结构井井筒溢流动态特性分析方法。
背景技术
特殊结构井包括定向井，大位移井，大斜度井，水平井等井身结构不是直井的井。这些特殊结构的井，无论从设计钻井工艺、完井措施等方面，要求远大于普通井的钻井。
井筒溢流是指在钻井过程中，地层气体侵入井筒后与钻井液形成气液两相流，并从井口溢出。特殊结构井发生溢流时，由于其井身结构特点，井筒流动特征与常规直井井筒流动有很大不同，研究特殊结构井井筒溢流特征，必须考虑其井身结构的影响，然而现有井筒溢流模型并未考虑特殊结构井的实际井眼轨迹，来建立特殊结构井环空流动模型。对特殊结构井发生溢流，现有研究对指导现场操作存在局限性。
发明内容
本申请的目的在于克服现有技术中特殊结构井在构建井筒溢流模型时未考虑实际井眼轨迹的问题，提供一种特殊结构井井筒溢流动态特性分析方法，在特殊结构井井筒的流体流动数值模型中引入了井斜角，使其得到的井筒溢流模型更符合特殊结构井的实际溢流，提高井筒流体流动数值模型对特殊结构井现场操作的指导参考价值。
为了实现上述发明目的，本申请提供了以下技术方案：一种特殊结构井井筒溢流动态特性分析方法，包括以下步骤：
S1：测井获得若干目标井的井深数据和井斜角数据，并根据插值计算方法获得所述目标井上任一点的井斜角，获得所述目标井的井眼轨迹；
S2：将所述目标井的井身分为若干节段，根据所述目标井内流体的质量守恒定律和动量守恒定律，建立每一节段的流体流动的数值模型；
S3：设定所述目标井的初始条件和边界条件，带入所述目标井的井眼轨迹，离散求解每一所述节段的流体的流动数值模型，得到所述目标井井筒内的流体流动参数。
进一步地，所述S1具体包括以下步骤：
S11：沿所述目标井选取若干个测点，测井获得若干个所述目标井的井深数据、井斜角数据；
S12：以井深为自变量建立井斜角的计算式，根据插值计算方法构建所述目标井两个测点之间的井斜角样条函数α(s)的表达式；
S13：求解井斜角样条函数α(s)的表达式，得到所述目标井任一井深处的井斜角，获得所述目标井的井眼轨迹。
进一步地，所述井斜角样条函数α(s)的表达式为：
Figure BDA0002972361620000021
其中：
Figure BDA0002972361620000022
Mk＝αn(sk)，Mk-1＝αn(sk-1)；
k——测点序号；
Lk——节段长度，Lk＝sk-sk-1，m；
s——测井深，m；
sk-1——k-1点处的井深，m；
sk——k点处的井深，m；
N——节段个数。
进一步地，所述井斜角样条函数α(s)中Mk通过以下方法获得：
S121：设所述目标井的井口和井底的井斜角的二阶倒数均为常数，建立含有N-1个Mk的线性方程组；
S122：采用追赶法求解线性方程组，得到Mk的值。
进一步地，流体流动数值模型为每一节段的气相连续性方程、液相连续性方程和、总动量方程联立的方程组，其方程组为：
Figure BDA0002972361620000023
其中，ρg和ρl分别为气体和液体密度，kg/m3；vg和vl分别为气体和液体速度，m/s；Eg为含气率；El为持液率；A为截面积，m2；α为井斜角，°；Qg为侵入气体流量，m3/min；
Figure BDA0002972361620000031
为摩阻压降，MPa/m。
进一步地，流体流动数值模型通过以下步骤建立：
S21：将所述目标井的井身分为若干节段，根据质量守恒定律，获得一所述节段内流体的气相连续性方程和液相连续性方程；其气相连续性方程为：
Figure BDA0002972361620000032
其液相连续性方程为：
Figure BDA0002972361620000033
S22：根据动量守恒定律，获得该节段的气相动量方程和液相动量方程；
其中该节段的气相动量方程为：
Figure BDA0002972361620000034
其中，该节段的液相动量方程为：
Figure BDA0002972361620000035
S23：代入含气率和持液率的关系式、摩阻压降的计算式，建立该节段的总动量方程，所述总动量方程为：
Figure BDA0002972361620000036
其中，所述含气率和持液率的关系式为：
Eg+El＝1 (14)；
其中，所述摩阻压降的计算式为：
Figure BDA0002972361620000037
其中，ρg和ρl分别为气体和液体密度，kg/m3；vg和vl分别为气体和液体速度，m/s；Eg为含气率；El为持液率；A为节段的横截面积，m2；α为井斜角，°；Qg为侵入气体流量，m3/min；
Figure BDA0002972361620000038
为摩阻压降，MPa/m。
S24：联立式(8)、式(9)、式(13)，获得流体流动数值模型。
进一步地，所述目标井的初始条件为所述目标井在气体未侵入井底时得到的所述目标井环空内的压力和流体速度；所述目标井的初始条件为：
Figure BDA0002972361620000041
其中，h为井深，m；Ql为泥浆排量，L/s；Pb为井底压力，MPa。
进一步地，所述目标井的边界条件为所述目标井内气体达到井口时的压力和井筒内总流量；所述目标井的边界条件为：
Figure BDA0002972361620000042
其中，P0是井口大气压，MPa。
进一步地，步骤S3中，流体流动数值模型的离散求解包括以下步骤：
S31：选取所述目标井的一节段L，将节段L进行时间离散和空间离散，建立流体流动数值模型的通用偏微分方程；
S32：将流体流动数值模型的通用偏微分方程在节段L的区域进行积分，并根据格林定理得到节段L的边界曲线积分方程并简化，得到简化的边界积分方程；
S33：将流体流动数值模型的每一方程分别带入所述简化的边界积分方程，得到流体流动数值模型中每一方程的差分方程；
S34：将所述目标井的初始条件和边界条件带入流体流动数值模型中每一方程的差分方程中，得到所述目标井的井筒溢流动态模型。
进一步地，在所述节段L进行时间离散和空间离散中，进行时间离散的时间步长与进行空间离散的步长满足以下关系：
Figure BDA0002972361620000043
其中，Δt为时间步长；Δz为空间步长；vg为气体速度。
与现有技术相比，本发明的具有以下有益效果：本发明公开的特殊结构井井筒溢流动态特性分析方法，在建立特殊结构井井筒溢流动态模型时，考虑特殊结构井的井眼轨迹对溢流的影响，在井筒溢流动态模型的建立过程中，引入了井斜角这一参数，充分考虑了特殊结构井的井眼轨迹，得到井筒内考虑了井眼轨迹的流体流动参数，为特殊结构井的溢流现场操作提供更适用的理论指导。
附图说明
图1本发明公开的特殊结构井井筒溢流动态特性分析方法的流程示意图；
图2本发明公开的特殊结构井井筒溢流动态特性分析方法中一节段L的网格划分和离散单元区域示意图；
图3本发明公开的实施例中得到的直井、大斜度井、水平井的井眼轨迹；
图4本发明公开的实施例中得到直井、大斜度井、水平井的环空压力变化曲线；
图5本发明公开的实施例中得到直井、大斜度井、水平井的井底压力变化曲线；
图6本发明公开的实施例中得到直井、大斜度井、水平井的泥浆增量变化曲线；
图7本发明公开的一实施例采用特殊结构井井筒溢流动态特性分析方法与现有分析方法得到的泥浆增量的变化曲线对比图；
图8本发明公开的一实施例采用特殊结构井井筒溢流动态特性分析方法与现有分析方法得到的井底压力的变化曲线对比图；
具体实施方式
下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例，凡基于本发明内容所实现的技术均属于本发明的范围。
特殊结构井发生溢流时，由于其井身结构特点，井筒流动特征与常规直井井筒流动有很大不同，研究特殊结构井井筒溢流特征，必须考虑其井身结构的影响，然而现有井筒溢流模型并未考虑特殊结构井的实际井眼轨迹，来建立特殊结构井环空流动模型。对特殊结构井发生溢流，现有研究对指导现场操作存在局限性。
基于上述技术问题的考虑，参阅图1，本申请公开了一种特殊结构井井筒溢流动态特性分析方法，包括以下步骤：
S1：测井获得若干目标井的井深数据和井斜角数据，并根据插值计算方法获得所述目标井上任一点的井斜角，获得所述目标井的井眼轨迹；
S2：将所述目标井的井身分为若干节段，根据所述目标井内流体的质量守恒定律和动量守恒定律，建立每一节段的流体流动数值模型；
S3：设定所述目标井的初始条件和边界条件，带入所述目标井的井眼轨迹，离散求解每一所述节段的流体流动数值模型，得到所述目标井井筒内的流体流动参数。
需要说明的是，本申请文件建立的所述目标井井筒的流体流动数值模型，基于以下假设：
(1)井筒内的气体和钻井液流动是一维的；
(2)钻井液不可压缩；
(3)每一节段内气、液相连续；
(4)忽略环空偏心的影响；
(5)不考虑气体在钻井液中的溶解。
需要说明的是，所述S1具体包括以下步骤：
S11：沿所述目标井选取N+1个测点，测井获得N+1个所述目标井的井深数据、井斜角数据；
S12：以井深为自变量建立井斜角的计算式，根据插值计算方法构建所述目标井两个测点之间的井斜角样条函数α(s)的表达式；
所述井斜角样条函数α(s)中Mk通过以下方法获得：
S121：设所述目标井的井口和井底的井斜角的二阶倒数均为常数，建立含有N-1个Mk的线性方程组；
S122：采用追赶法求解线性方程组，得到Mk的值。
其中，线性方程组的建立步骤如下：
S1211：得到所述目标井N+1个井深数据和对应的N+1个井斜角数据，并根据测点位置(即井深)将井深数据依次表示为S0、S1……SN；对应的井斜角数据依次表示为：α0、α1……αN；得到所述目标井的井眼轨迹数据：
Figure BDA0002972361620000061
S1212：任意取两个测点之间的区间，[sk-1,sk](k＝1,2,3,......,N)，根据三次样条函数，以井深为自变量，以井斜角表达为井深的函数，建立在[sk-1,sk](k＝1,2,3,......,N)区间内井斜角样条函数α(s)的表达式；
所述井斜角样条函数α(s)的表达式为：
Figure BDA0002972361620000062
其中：
Figure BDA0002972361620000063
Mk＝αn(sk)，Mk-1＝αn(sk-1)；
k——测点序号；
Lk——节段长度，Lk＝sk-sk-1，m；
s——测井深，m；
sk-1——k-1点处的井深，m；
sk——k点处的井深，m；
N——节段个数。
需要说明的是，此处的井深s是指所述目标井的井深。
S1213：由于Mk既与井斜角的测量值有关，又与所述目标井的井口和井底的边界条件有关；因此设井口和井底的井斜角的二阶倒数为常数，则有：
M0＝MN＝0 (3)；
因此，对于所述目标井全井的N+1个测点，得到含有N-1个未知数Mk的线性方程组：
Figure BDA0002972361620000071
其中，
Figure BDA0002972361620000072
λ0＝1，
Figure BDA0002972361620000073
由于，式(4)为对角线方程组，即可采用追赶法求解，得到Mk的值。
S13：求解井斜角样条函数α(s)的表达式，即可得到所述目标井两个测点之间任一位置的井斜角。
然后将Mk的值带入α(s)的表达式中，即可求出在[sk-1,sk](k＝1,2,3,......,N)区间上任意井深处的井斜角，进而得到所述目标井的井眼轨迹。
需要说明的是，所述井眼轨迹就是一系列的井深和该井深数据点对应的井斜角数据，其井深数据和井斜角数据构成了井眼轨迹的坐标数据。但测井获得的井眼轨迹的坐标数据是不连续的，难以进行分析计算，因此为了得到所述目标井的井眼轨迹，在本申请中采用插值计算方法对井深和井斜角进行曲线拟合。
在步骤S2中，流体流动数值模型通过以下步骤建立：
S21：将所述目标井的井身分为若干节段，然后通过所述目标井的井眼轨迹计算出对应节段的垂深h。
所述目标井任一节段的垂深通过以下方法获得：
在所述目标井的垂直段，井斜角α为0，每一节段的垂深等于该节段的长度，因此：
Figure BDA0002972361620000081
在所述目标井的造斜段，每一节段的垂深与该节段的长度为关系为：
Figure BDA0002972361620000082
需要说明的是，此处的每一节段的垂深dz是指该节段在垂直方向上的长度。
根据质量守恒定律，获得一所述节段内流体的气相连续性方程和液相连续性方程；其具体包括：
根据质量守恒定律，每一节段的质量变化量为该节段内进入的流体质量减去流出的流体质量。本申请文件中所述结构井内的流体视为气液两相组成的多相流流体，因此，以每一节段内的气相为例，进入该节段内的气体质量为：
Figure BDA0002972361620000083
流出该节段的气体质量为：
Figure BDA0002972361620000084
当地层气体侵入井筒，节段内气体发生膨胀，则该节段出现含气率变化，因此该节段流动的气体质量发生变化，因此该节段内流动的气体质量变化为：
Figure BDA0002972361620000085
因此，通过式(5)和式(7)即可获得该节段内流体的气相连续性方程，其气相连续性方程具体为：
Figure BDA0002972361620000086
同理的，可以得到该节段内流体的液相连续性方程，其液相连续性方程具体为：
Figure BDA0002972361620000087
S22：根据动量守恒定律，物体的动量相对于时间的变化率等于施加在节段或流体上的外力之和。物体的动量相对于时间的变化率有两部分:局部导数和对流导数，因此可以得到动量方程:
Figure BDA0002972361620000091
因此，可以获得该节段的气相动量方程和液相动量方程；
其中该节段的气相动量方程为：
Figure BDA0002972361620000092
其中，该节段的液相动量方程为：
Figure BDA0002972361620000093
S23：代入含气率和持液率的关系式、摩阻压降的计算式，建立该节段的总动量方程，所述总动量方程为：
Figure BDA0002972361620000094
其中，所述含气率和持液率的关系式为：
Eg+El＝1 (14)；
其中，所述摩阻压降的计算式为：
Figure BDA0002972361620000095
其中，ρg和ρl分别为气体和液体密度，kg/m3；vg和vl分别为气体和液体速度，m/s；Eg为含气率；El为持液率；A为节段的横截面积，m2；α为井斜角，°；Qg为侵入气体流量，m3/min；
Figure BDA0002972361620000096
为摩阻压降，MPa/m；i为l或者g。
S24：联立式(8)、式(9)、式(12)，获得流体流动数值模型。
流体流动数值模型为：为每一节段的气相连续性方程、液相连续性方程和、总动量方程联立的方程组，其方程组为：
Figure BDA0002972361620000101
在所述步骤S3中，所述数值方程的约束条件包括边界条件和初始条件。当所述结构井井底未溢流时，即地层气体还没有侵入所述目标井的井筒，气体还没有进入环空，钻井液填充整个环空，可以得到的所述目标井环空内的压力和流体速度；因此，所述目标井的初始条件为：
Figure BDA0002972361620000102
其中，h为井深，m；Ql为泥浆排量，L/s；Pb为井底压力，MPa。
其中，所述目标井的边界条件为所述目标井内气体达到井口时的压力和井筒内总流量；此时所述目标井的井口压力为大气压力，井筒内总流量等于钻井液流量和侵入的地层气体流量之和，因此，所述目标井的边界条件为：
Figure BDA0002972361620000103
其中，P0是井口大气压，MPa。
流体流动数值模型的离散求解包括以下步骤：
S31：选取所述目标井的一节段L，并对节段L进行网格划分，将节段L进行时间离散和空间离散，建立流体流动数值模型的通用偏微分方程；
在所述节段L进行时间离散和空间离散中，进行时间离散的时间步长与进行空间离散的步长满足以下关系：
Figure BDA0002972361620000104
其中，Δt为时间步长；Δz为空间步长；vg为气体速度。
所述节段L的网格划分和进行离散的单元网格区域示意图如图2所示。
流体流动数值模型的通用偏微分方程可以写成：
Figure BDA0002972361620000105
S32：将流体流动数值模型的通用偏微分方程在节段L的区域进行积分，并根据格林定理得到节段L的边界曲线积分方程，其边界曲线积分方程为：
Figure BDA0002972361620000111
再将式(20)简化，得到简化的边界积分方程：
Figure BDA0002972361620000112
S33：将流体流动数值模型的每一方程分别带入所述简化的边界积分方程，得到流体流动数值模型中每一方程的差分方程；
在气相连续性方程中，令：
Figure BDA0002972361620000113
因此，根据式(8)和式(9)、式(21)可得气相差分方程：
Figure BDA0002972361620000114
同理，在液相连续性方程中，令：
Figure BDA0002972361620000115
得到液相差分方程：
Figure BDA0002972361620000116
在总动量方程中，令
Figure BDA0002972361620000117
根据式(21)即可得到混合动量差分方程：
Figure BDA0002972361620000118
S34：将所述目标井的初始条件和边界条件带入流体流动数值模型中每一方程的差分方程中，即可获得所述目标井内的流体流动参数。
其中，节段L的初始条件为：
Figure BDA0002972361620000121
边界条件为：
Figure BDA0002972361620000122
具体实施例：
以下以某油藏的三口试验井为例，其中试验井A(以下简称为A井)为直井，其总井深为4923m；其中试验井B(以下简称为B井)为大斜度井，其垂深为4923m，井深为5325m，，其中，直井段垂深3923m，造斜段垂深1000m，造斜段长1246m；其中试验井C(以下简称为C井)为水平井，其总垂深为4923m，其中，直井段垂深3923m，造斜段垂深1000m，造斜段长1246m，水平段365m。
分别在三个试验井均匀每30m选取一个个测点，进行测井，获得井深数据和井斜角数据，采用三次样条函数进行插值计算，对井深数据和井斜角数据进行拟合，得到三个试验井的井眼轨迹，其得到井眼轨迹如图3所示。
三个试验井的主要参数如表1所示，分别将三个试验井的主要参数带入数值模型中，得到三个试验井的溢流参数；其中每一试验井的(具体说环空压力，井底压力变化，泥浆增量等参数)溢流参数如图4～6所示。
从图4～6可以看出，井筒发生溢流后，大斜度井和水平井的环空压力更低，泥浆池泥浆增加更多，在溢流初期，直井井底压力就开始不断降低，泥浆池泥浆开始不断增加，而大斜度井和水平井井底压力变化较小，泥浆池泥浆变化不明显；溢流后期，直井、大斜度井和水平井井底压力开始快速降低，泥浆池泥浆快速增加，其中大斜度井和水平井井底压力降低更快，泥浆池泥浆增加更多。这是因为在溢流初期，由于大斜度井和水平井的井身结构特点，进入井筒的气体不会像直井那样在环空中迅速向上运动，在水平段或造斜段，环空压力变化较小，气体与钻井液之间没有滑脱或滑脱较小，气体移动缓慢，此时进入的气体对井底压力影响较小；当气体运移到直井段，气体向上运动速度加快，造成井底压力大幅下降，泥浆池泥浆快速增加。因此，相对于直井，大斜度井和水平井发生溢流，井底压力变化滞后，溢流后期，环空气体膨胀剧烈，井底压力快速降低，大大增加了大斜度井和水平井溢流风险。
为验证本文模型的有效性，选取酒东一口发生溢流定向井的基础数据进行动态溢流模拟。在22:14时用密度1.79g/cm3钻井液钻进至井深3725m时发生了溢流，此时井底压力为65.41Mpa，泥浆池泥浆增加0.45m3，22:36时发生了剧烈的井喷事故，现场每3分钟记录一次泥浆池泥浆增量，并在22:37溢流达到5.9m3时关井。以该井的基础参数用本文建立的流体流动数值模型进行模拟，将本文的计算结果与现场数据和未考虑实际井眼轨迹的模型(土林，《特殊结构井井控技术研究》，西南石油大学，2009)计算结果进行比较，如图7、图8所示。
从图7、图8可知，采用本发明公开的分析方法建立的流体流动数值模型与现场数据温和较好，具有较高的准确度。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种特殊结构井井筒溢流动态特性分析方法，其特征在于，包括以下步骤：
S1：测井获得若干目标井的井深数据和井斜角数据，并根据插值计算方法获得所述目标井上任一点的井深和井斜角，获得所述目标井连续平滑的井眼轨迹；
S2：将所述目标井的井身分为若干节段，根据所述目标井内流体的质量守恒定律和动量守恒定律，建立每一节段的流体流动数值模型；
S3：设定所述目标井的初始条件和边界条件，带入所述目标井的井眼轨迹，离散求解每一所述节段的流体流动数值模型，得到所述目标井井筒内的流体流动参数。
2.根据权利要求1所述的特殊结构井井筒溢流动态特性分析方法，其特征在于，所述S1具体包括以下步骤：
S11：沿所述目标井选取若干个测点，测井获得若干个所述目标井的井深数据、井斜角数据；
S12：以井深为自变量建立井斜角的计算式，根据插值计算方法构建所述目标井两个测点之间的井斜角样条函数α(s)的表达式；
S13：求解井斜角样条函数α(s)的表达式，得到所述目标井任一井深处的井斜角，获得所述目标井的井眼轨迹。
3.根据权利要求2所述的特殊结构井井筒溢流动态特性分析方法，其特征在于，所述井斜角样条函数α(s)的表达式为：
Figure FDA0002972361610000011
其中：
Figure FDA0002972361610000012
Mk＝αn(sk)，Mk-1＝αn(sk-1)；
k——测点序号；
Lk——节段长度，Lk＝sk-sk-1，m；
s——测井深，m；
sk-1——k-1点处的井深，m；
sk——k点处的井深，m；
N——节段个数。
4.根据权利要求3所述的特殊结构井井筒溢流动态特性分析方法，其特征在于，所述井斜角样条函数α(s)中Mk通过以下方法获得：
S121：设所述目标井的井口和井底的井斜角的二阶倒数均为常数，建立含有N-1个Mk的线性方程组；
S122：采用追赶法求解线性方程组，得到Mk的值。
5.根据权利要求1所述的特殊结构井井筒溢流动态特性分析方法，其特征在于，所述流体流动数值模型为每一节段的气相连续性方程、液相连续性方程、总动量方程联立的方程组，其方程组为：
Figure FDA0002972361610000021
其中，ρg和ρl分别为气体和液体密度，kg/m3；vg和vl分别为气体和液体速度，m/s；Eg为含气率；El为持液率；A为截面积，m2；α为井斜角，°；Qg为侵入气体流量，m3/min；
Figure FDA0002972361610000022
为摩阻压降，MPa/m。
6.根据权利要求5所述的特殊结构井井筒溢流动态特性分析方法，其特征在于，所述流体流动数值模型通过以下步骤建立：
S21：将所述目标井的井身分为若干节段，根据质量守恒定律，获得一所述节段内流体的气相连续性方程和液相连续性方程；其气相连续性方程为：
Figure FDA0002972361610000023
其液相连续性方程为：
Figure FDA0002972361610000031
S22：根据动量守恒定律，获得该节段的气相动量方程和液相动量方程；
其中该节段的气相动量方程为：
Figure FDA0002972361610000032
其中，该节段的液相动量方程为：
Figure FDA0002972361610000033
S23：代入含气率和持液率的关系式、摩阻压降的计算式，建立该节段的总动量方程，所述总动量方程为：
Figure FDA0002972361610000034
其中，所述含气率和持液率的关系式为：
Eg+El＝1 (14)；
其中，所述摩阻压降的计算式为：
Figure FDA0002972361610000035
其中，ρg和ρl分别为气体和液体密度，kg/m3；vg和vl分别为气体和液体速度，m/s；Eg为含气率；El为持液率；A为节段的横截面积，m2；α为井斜角，°；Qg为侵入气体流量，m3/min；
Figure FDA0002972361610000036
为摩阻压降，MPa/m。
S24：联立气相连续性方程、液相连续性方程、总动量方程，获得所述流体流动数值模型。
7.根据权利要求1所述的特殊结构井井筒溢流动态特性分析方法，其特征在于，所述目标井的初始条件为所述目标井在气体未侵入时得到的所述目标井环空内的压力和流体速度；所述目标井的初始条件为：
Figure FDA0002972361610000041
其中，Eg为含气率；El为持液率；vl为液体速度，m/s；A为节段的横截面积，m2；h为井深，m；Ql为泥浆排量，L/s；Pb为井底压力，MPa。
8.根据权利要求1所述的特殊结构井井筒溢流动态特性分析方法，其特征在于，所述目标井的边界条件为所述目标井内气体达到井口时的压力和井筒内总流量；所述目标井的边界条件为：
Figure FDA0002972361610000042
其中，P0是井口大气压，MPa。
9.根据权利要求1所述的特殊结构井井筒溢流动态特性分析方法，其特征在于，步骤S3中，所述流体流动数值模型的离散求解包括以下步骤：
S31：选取所述目标井的一节段L，将节段L进行时间离散和空间离散，建立所述流体流动数值模型的通用偏微分方程；
S32：将所述流体流动数值模型的通用偏微分方程在节段L的区域进行积分，并根据格林定理得到节段L的边界曲线积分方程并简化，得到简化的边界积分方程；
S33：将所述流体流动数值模型的每一方程分别带入所述简化的边界积分方程，得到所述流体流动数值模型中每一方程的差分方程；
S34：将所述目标井的初始条件和边界条件带入流体流动数值模型中每一方程的差分方程中，得到所述目标井的井筒溢流动态模型。
10.根据权利要求9所述的特殊结构井井筒溢流动态特性分析方法，其特征在于，在所述节段L进行时间离散和空间离散中，进行时间离散的时间步长与进行空间离散的步长满足以下关系：
Figure FDA0002972361610000043
其中，Δt为时间步长；Δz为空间步长；vg为气体速度。
CN202110266163.4A 2021-03-11 2021-03-11 一种特殊结构井井筒溢流动态特性分析方法 Pending CN112901144A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110266163.4A CN112901144A (zh) 2021-03-11 2021-03-11 一种特殊结构井井筒溢流动态特性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110266163.4A CN112901144A (zh) 2021-03-11 2021-03-11 一种特殊结构井井筒溢流动态特性分析方法

Publications (1)

Publication Number Publication Date
CN112901144A true CN112901144A (zh) 2021-06-04

Family

ID=76105671

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110266163.4A Pending CN112901144A (zh) 2021-03-11 2021-03-11 一种特殊结构井井筒溢流动态特性分析方法

Country Status (1)

Country Link
CN (1) CN112901144A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114818516A (zh) * 2022-06-27 2022-07-29 中国石油大学(华东) 一种井筒腐蚀形态剖面智能预测方法及程序产品

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130062077A1 (en) * 2011-09-13 2013-03-14 Jim Basuki Surjaatmadja Methods and equipment to improve reliability of pinpoint stimulation operations
CN105239944A (zh) * 2015-10-14 2016-01-13 中国石油天然气股份有限公司 一种复杂井眼轨迹扶正器加装方法
CN109403907A (zh) * 2018-10-18 2019-03-01 西南石油大学 一种深水钻完井上下一体井控安全控制新方法
CN110424954A (zh) * 2019-08-09 2019-11-08 西南石油大学 基于井筒环空瞬态多相流流动特征的环空瞬态水击模型
US20200232316A1 (en) * 2019-01-23 2020-07-23 Saudi Arabian Oil Company Well kick detection
CN112459747A (zh) * 2020-12-08 2021-03-09 中国石油天然气集团有限公司 基于三次样条曲线的控压钻井限压调阀的控制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130062077A1 (en) * 2011-09-13 2013-03-14 Jim Basuki Surjaatmadja Methods and equipment to improve reliability of pinpoint stimulation operations
CN105239944A (zh) * 2015-10-14 2016-01-13 中国石油天然气股份有限公司 一种复杂井眼轨迹扶正器加装方法
CN109403907A (zh) * 2018-10-18 2019-03-01 西南石油大学 一种深水钻完井上下一体井控安全控制新方法
US20200232316A1 (en) * 2019-01-23 2020-07-23 Saudi Arabian Oil Company Well kick detection
CN110424954A (zh) * 2019-08-09 2019-11-08 西南石油大学 基于井筒环空瞬态多相流流动特征的环空瞬态水击模型
CN112459747A (zh) * 2020-12-08 2021-03-09 中国石油天然气集团有限公司 基于三次样条曲线的控压钻井限压调阀的控制方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
刘修善: "如何用样条函数模拟实际井眼轨迹", 《大庆石油学院学报》 *
孙锦山等: "《理论爆轰物理》", 31 December 1995, 国防工业出版社 *
李相方等: "欠平衡钻井期间地层流体流入规律研究", 《石油学报》 *
毛良杰等: "含H2S天然气井溢流动态模拟", 《科学技术与工程》 *
白方方等: "水平井井控井筒压力分布规律研究", 《重庆科技学院学报(自然科学版)》 *
董世明等: "《水平井有杆抽油系统设计》", 31 December 1996, 石油工业出版社 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114818516A (zh) * 2022-06-27 2022-07-29 中国石油大学(华东) 一种井筒腐蚀形态剖面智能预测方法及程序产品

Similar Documents

Publication Publication Date Title
US11203922B1 (en) Method and equipment for optimizing hydraulic parameters of deepwater managed pressure drilling in real time
CN104213906B (zh) 一种钻井井筒压力校正方法
CN102943620B (zh) 基于钻井环空井筒多相流动计算的控压钻井方法
CN108509703B (zh) 一种气藏状态参数随钻数值反演分析方法
CN104453804A (zh) 一种气驱油藏开发动态监测与评价方法
CN111894561B (zh) 一种适用于欠平衡钻井的地层特性随钻解释方法
CN115587674B (zh) 油藏改建储气库扩容达产过程气井动态产能预测方法
CN112576248A (zh) 一种底水气藏早期产能评价预测方法
Liu et al. A unified formula for determination of wellhead pressure and bottom-hole pressure
Da Silva et al. A review of coupled dynamic well-reservoir simulation
CN112901144A (zh) 一种特殊结构井井筒溢流动态特性分析方法
CN103122756B (zh) 一种确定深水隔水管气举钻井注气量的方法
CN113836742A (zh) 一种基于移动大井法原理的动态计算矿井涌水量的方法
CN111626001B (zh) 一种提高采油井精细化注水的方法
US6999878B2 (en) Method and installation for locating the position of the boundary formed at the interface between two fluids contained in a reservoir
CN110242261A (zh) 垂向注气孔隙油气水微观渗流规律的预测方法及系统
CN114638179A (zh) 一种微量气侵条件下钻井安全关井周期的确定方法
CN114810023A (zh) 一种适应于密集井眼防碰的施工方法
Shi et al. Analytical model for transient pressure analysis in a horizontal well intercepting with multiple faults in karst carbonate reservoirs
Wang et al. Modeling the flow of carbon dioxide during the drilling of oil, gas, and geothermal energy
CN111119856B (zh) 基于实际钻井数据的地层自然造斜特性反演方法
CN111364955A (zh) 注采井间流场演化模拟方法
Ehirim et al. Modelling and simulation of a gas condensate reservoir
Al-Hajri et al. Estimating the Static Bottom-hole Pressure of Gas Wells by Top Node Calculation Using Apparent Molecular Weight Profiling
Perneder et al. Steady-state solutions of a propagating borehole: Helical trajectory

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210604