CN108331574A - 一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法 - Google Patents
一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法 Download PDFInfo
- Publication number
- CN108331574A CN108331574A CN201810013090.6A CN201810013090A CN108331574A CN 108331574 A CN108331574 A CN 108331574A CN 201810013090 A CN201810013090 A CN 201810013090A CN 108331574 A CN108331574 A CN 108331574A
- Authority
- CN
- China
- Prior art keywords
- sand control
- well
- formula
- horizontal segment
- section
- 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
- 239000004576 sand Substances 0.000 title claims abstract description 170
- 230000011218 segmentation Effects 0.000 title claims abstract description 53
- 239000011435 rock Substances 0.000 claims abstract description 37
- 238000004519 manufacturing process Methods 0.000 claims abstract description 20
- 238000009826 distribution Methods 0.000 claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims description 34
- 230000003068 static effect Effects 0.000 claims description 26
- 238000000034 method Methods 0.000 claims description 15
- 201000004569 Blindness Diseases 0.000 abstract description 3
- 230000007774 longterm Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 14
- 239000003921 oil Substances 0.000 description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 4
- 239000003129 oil well Substances 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000006735 deficit Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 239000000295 fuel oil Substances 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000004080 punching Methods 0.000 description 2
- 239000009671 shengli Substances 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000007596 consolidation process Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，包括：S1：计算储层岩石力学参数沿水平井水平段的分布；S2：预测水平井水平段出砂临界生产压差分布；S3：计算水平段各点的实际井底流压，以及水平段各点的实际生产压差；标定井底流压即水平段根端位置的井底流压；S4：计算得到出砂净压差，进行无量纲化计算得到出砂强度指数；S5：绘制水平井相对出砂剖面图，划分出砂亏空程度等级，并实现出砂程度的分段；S6：决定分段分级防砂的分段分级方案。本发明通过相对出砂剖面预测指导分段分级防砂，解决水平井笼统防砂的盲目性强、防砂效果差的问题，降低了防砂作业的成本，保证水平井长期稳定生产，提高经济效益。
Description
技术领域
本发明涉及一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，属于油气勘探技术领域。
背景技术
疏松砂岩油气资源储量在全球及国内都占较大的比例，开发潜力巨大。水平井开采技术因其泄油面积大、采油速度快、有效减缓底水锥进，被广泛应用于油气藏开采，特别是疏松砂岩油气藏，是新油田开发、老油田挖潜的重要技术手段。
由于疏松砂岩储层岩石胶结程度差，在水平井开采疏松砂岩油气藏过程中常伴有不同程度的出砂现象，轻则造成抽油泵磨蚀、油嘴冲蚀；重则形成地层亏空、储层坍塌、套管挤毁、甚至油气井报废。此外，水平井水平段长度大，储层均质性差，各井段出砂程度各异，增加了防砂施工的复杂性。但由于目前对于水平井长水平段出砂位置和剖面难以预测，导致难以根据出砂严重层段进行针对性防砂，只能全井段笼统防砂，既增加成本，影响油气井产量，又难以取得理想的综合防砂效果，不利于提高经济效益。
发明内容
针对现有技术的疏松砂岩油气藏易出砂的问题，本发明提供了一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法；
本发明通过测井资料预测水平井段出砂临界生产压差分布，计算出砂净压差和出砂强度指数，进而绘制水平井相对出砂剖面图，划分出砂亏空程度等级；据此决定分段分级防砂的分段分级方法，优化防砂方案，提高了水平井防砂的针对性，提高了防砂效果,提高油井产量同时降低了防砂成本，提高了经济效益。
术语解释：
1、地层静压，全称为地层静止压力，也叫油层压力，是指油井在关井后，待压力恢复到稳定状态时所测得的油层中部压力，简称静压。
2、井底流压，全称井底流动压力，又简称流动压力或流压。是油、气井生产时的井底压力。它表示油、气从地层流到井底后剩余的压力，对自喷井来讲，也是油气从井底流到地面的起点压力。
3、流动压降计算模型，是指对于给定井段根据一端的压力和温度用于计算整个井段压力和温度分布的计算方法、公式等，常用的有Beggs-Brill模型、Aziz模型和Okasweki模型等。
4、井筒多相流模型，对于油井是指适用于油气水三相流动条件的流动压降计算模型，对于气井是指适用于气水两相流动条件的流动压降计算模型。
本发明的技术方案为：
一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，包括：
S1：计算储层岩石力学参数沿水平井水平段的分布，储层岩石力学参数包括动态泊松比、静态泊松比、动态弹性模量、单轴抗压强度、内聚力以及内摩擦角；
S2：根据储层岩石力学参数及地应力分布，预测水平井水平段出砂临界生产压差分布；地应力分布是指水平井周围岩石所处的垂向应力、切向应力和径向应力三个主应力；
S3：根据地层静压和标定井底流压，以及水平段流动压降计算模型，计算水平段各点的实际井底流压，以及水平段各点的实际生产压差(地层静压减去井底流压)；标定井底流压即水平段根端位置的井底流压；
S4：定义出砂净压差为实际生产压差与出砂临界生产压差的差值，计算得到出砂净压差，进行无量纲化计算得到出砂强度指数；
S5：绘制水平井相对出砂剖面图，划分出砂亏空程度等级，并实现不同等级出砂程度的分段；水平井相对出砂剖面图是指用来表达沿水平井井轴方向上出砂严重程度和亏空程度相对大小的剖面图。根据计算得到的出砂强度指数绘制，是出砂强度指数沿井深的分布曲线。
S6：根据出砂程度等级决定分段分级防砂的分段分级方案。
本发明通过相对出砂剖面预测指导分段分级防砂，解决水平井笼统防砂的盲目性强、防砂效果差的问题，降低了防砂作业的成本，保证水平井长期稳定生产，提高经济效益。
根据本发明优选的，所述步骤S1中，
动态泊松比的计算公式如式(Ⅰ)所示：
式(Ⅰ)中，μd—岩石的动态泊松比，无量纲；Δth—横波时差，μs/m；Δtv—纵波时差，μs/m；
静态泊松比的计算公式如式(Ⅱ)所示：
μ＝A·μd+B (Ⅱ)
式(Ⅱ)中，μ—岩石的静态泊松比，无量纲；A、B—经验系数，通过实测数据拟合得到，无量纲；A＝0.38，B＝0.082；
动态弹性模量的计算公式如式(Ⅲ)所示：
式(Ⅲ)中，Ed—岩石的动态弹性模量，MPa；ρr—地层岩石密度，kg/m3；
单轴抗压强度的计算公式如式(Ⅳ)所示：
σc＝[C·(1-VCL)+D·VCL]·Ed (Ⅳ)
式(Ⅳ)中，σc—岩石单轴抗压强度，MPa；C、D—为经验拟合系数，C＝0.00459，D＝0.00816；VCL—泥质含量，小数；
内聚力的计算公式如式(Ⅴ)所示：
式(Ⅴ)中，C0—内聚力，c/MPa；
内摩擦角的计算公式如式(Ⅵ)所示：
φf＝36.545-0.4952C0 (Ⅵ)
式(Ⅵ)中，φf—岩石内摩擦角，rad。
根据本发明优选的，所述步骤S2，出砂临界生产压差的计算公式如式(Ⅶ)所示：
式(Ⅶ)中，△Pc—出砂临界生产压差，MPa；C0—岩石内聚力，MPa；μ—岩石的静态泊松比，无量纲；σze—外边界处的垂向应力，MPa；β—比奥特常数，无量纲；α—失效角，rad；Pr—地层静压，MPa。
根据本发明优选的，所述步骤S3，水平段各点的实际井底流压的计算公式如式(Ⅷ)所示，水平段各点的实际生产压差的计算公式如式(Ⅸ)所示：
Pwf(i)＝Pwf0+ΔPw(i) (Ⅷ)
ΔP(i)＝Pr-Pwf(i) (Ⅸ)
式(Ⅷ)、(Ⅸ)中，Pwf0—标定井底流压，MPa；ΔPw(i)—第i点的压力与根端压力之差，MPa；i是水平井水平段按照长度计算的任一分段，1≤i≤I；分段长度为1-3m，分段总数I等于水平井水平段总长度除以分段长度；根据井筒多相流模型计算得到；Pwf(i)—第i点实际井底流压，MPa；Pr—地层静压，MPa；ΔP(i)—第i点的实际生产压差，MPa。
根据本发明优选的，所述步骤S4中，
定义出砂净压差为实际生产压差与出砂临界生产压差的差值，计算得到出砂净压差，如式(Ⅹ)所示：
Δp(i)＝ΔP(i)-ΔPc(i),ifΔp(i)≤0thenΔp(i)＝0 (Ⅹ)
式(Ⅹ)中，ΔPc(i)—第i点的出砂临界生产压差，MPa；ΔP(i)—第i点的出砂净压差，MPa；
计算得到出砂强度指数，计算公式如式(Ⅺ)所示：
式(Ⅺ)中，△pmax—全井段出砂净压差的最大值，Δpmax＝max{p(i)}；△pmin—全井段出砂净压差的最小值，Δpmin＝min{Δp(i)},ifΔp#in>0.25thenΔpmin＝0.25；Js(i)—第i点的出砂强度指数，无量纲。
根据本发明优选的，所述步骤S5，包括：
以5-15m为最小计量井段长度间隔划分水平井水平段，当求取的某水平段的出砂强度指数不大于0.1时，将该水平段划分为不出砂井段；当求取的某水平段的出砂强度指数为0.1-0.35时，将该水平段划分为轻微出砂井段；当求取的某水平段的出砂强度指数为0.35-0.70时，将该水平段划分为中等出砂井段；当求取的某水平段的出砂强度指数为0.70-1.0时，将该水平段划分为严重出砂井段。
进一步优选的，以10m为最小计量井段长度间隔划分水平井水平段。
根据本发明优选的，所述步骤S6，根据出砂程度等级决定分段分级防砂的分段分级方案，包括：
对于步骤S5中划为的不出砂井段，如果不出砂井段的长度(是指出砂等级被判别为不出砂的井段的长度)大于20m，则设计为零级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中防砂等级较低的防砂等级；
对于步骤S5中划为的轻微出砂井段，如果轻微出砂井段的长度大于20m，则设计为一级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中的防砂等级较低的防砂等级；
对于步骤S5中划为的中等出砂井段，如果中等出砂井段的长度大于20m，则设计为二级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中的防砂等级较低的防砂等级；
对于步骤S5中划为的严重出砂井段，如果严重出砂井段的长度大于20m，则设计为三级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中的防砂等级较低的防砂等级。
零级防砂井段、一级防砂井段、二级防砂井段、三级防砂井段是指根据出砂井段的出砂程度和需要防砂的强度所划分的防砂等级，等级越高，表示越需要重点防砂；等级越低，表示防砂的要求越低。其中，零级防砂是指不需要防砂，具体实施时零级防砂井段不需要管柱或仅使用打孔管(在光油管上打孔孔密50-120孔/m，孔径8-10mm)；一级防砂是指筛管挡砂介质厚度不低于5-8mm或砾石层厚度不低于20mm；二级防砂是指筛管挡砂介质厚度不低于10-15mm或砾石层厚度不低于30mm；三级防砂是指筛管挡砂介质厚度不低于20-25mm或砾石层厚度不低于50mm。
本发明的有益效果为：
1、本发明通过测井资料预测水平井段出砂临界生产压差分布，计算出砂净压差和出砂强度指数，进而绘制水平井相对出砂剖面图，划分出砂亏空程度等级；有助于准确了解水平井水平段的严重出砂位置，通过控制各分段的生产压差实现水平井水平段分段生产控制，利于获得均匀产出剖面和提高产量；并且降低了总体防砂附加渗流阻力和表皮系数。相比未知出砂剖面分布情况下的笼统生产，表皮系数降低0.2-0.5，平均产量提高15％左右。
2、本发明根据划分得到的水平井水平段出砂亏空和严重程度等级，决定分段分级防砂的分段分级方法，优化防砂方案，提高了水平井防砂的针对性和综合防砂效果。相比未知出砂剖面分布情况下的笼统防砂，防砂工艺设计和实施的目的性和针对性更强，防砂有效率提高10-15％。
3、本发明根据划分得到的水平井水平段出砂亏空和严重程度等级，决定分段分级防砂的分段分级方法，优化防砂方案，避免了不出砂和轻微出砂井段的过度防砂或不必要的防砂，降低施工难度和防砂成本。相比全井段笼统防砂，平均成本节约30％以上。
4、本发明提出的出砂相对剖面预测和分级分段防砂的分段分级方法，为提高疏松砂岩水平井综合防砂效果提供了一条有效途径。尤其与水平段控水结合起来，可以实现控水防砂一体化设计，同时实现高效防砂和控水，提高了含水疏松砂岩油气藏水平井开发效果。
附图说明
图1为实施例1水平井水平段出砂临界生产压差分布计算结果示意图；
图2为实施例1水平井水平段实际生产压差分布计算结果示意图；
图3为实施例1预测得到的水平井水平段出砂强度指数分布示意图；
图4为实施例1预测得到的水平井相对出砂剖面及出砂亏空程度等级划分结果示意图；
图5为实施例1水平井分级防砂多段多级划分结果示意图；
图6为实施例2预测得到大港油田某水平井水平段的出砂临界压差、生产压差、出砂净压差分布示意图；
图7为实施例2预测得到大港油田某水平井水平段的分级防砂多段多级划分结果示意图；
图8为实施例2预测得到胜利油田某稠油油藏水平井水平段的分级防砂多段多级划分结果示意图。
具体实施方式
下面结合说明书附图和实施例对本发明作进一步限定，但不限于此。
实施例1
一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，包括：
S1：根据水平井测井资料(密度测井、纵波测井、横波测井)计算储层岩石力学参数沿水平井水平段的分布，储层岩石力学参数包括动态泊松比、静态泊松比、动态弹性模量、单轴抗压强度、内聚力以及内摩擦角；
动态泊松比的计算公式如式(Ⅰ)所示：
式(Ⅰ)中，μd—岩石的动态泊松比，无量纲；Δth—横波时差，μs/m；Δtv—纵波时差，μs/m；
静态泊松比的计算公式如式(Ⅱ)所示：
μ＝A·μd+B (Ⅱ)
式(Ⅱ)中，μ—岩石的静态泊松比，无量纲；A、B—经验系数，通过实测数据拟合得到，无量纲；A＝0.38，B＝0.082；
动态弹性模量的计算公式如式(Ⅲ)所示：
式(Ⅲ)中，Ed—岩石的动态弹性模量，MPa；ρr—地层岩石密度，kg/m3；
单轴抗压强度的计算公式如式(Ⅳ)所示：
σc＝[C·(1-VCL)+D·VCL]·Ed (Ⅳ)
式(Ⅳ)中，σc—岩石单轴抗压强度，MPa；C、D—为经验拟合系数，C＝0.00459，D＝0.00816；VCL—泥质含量，小数；
内聚力的计算公式如式(Ⅴ)所示：
式(Ⅴ)中，C0—内聚力，c/MPa；
内摩擦角的计算公式如式(Ⅵ)所示：
φf＝36.545-0.4952C0 (Ⅵ)
式(Ⅵ)中，φf—岩石内摩擦角，rad。
S2：根据储层岩石力学参数及地应力分布，预测水平井水平段出砂临界生产压差分布；地应力分布是指水平井周围岩石所处的垂向应力、切向应力和径向应力三个主应力；出砂临界生产压差的计算公式如式(Ⅶ)所示：
式(Ⅶ)中，△Pc—出砂临界生产压差，MPa；C0—岩石内聚力，MPa；μ—岩石的静态泊松比，无量纲；σze—外边界处的垂向应力，MPa；β—比奥特常数，无量纲；α—失效角，rad；Pr—地层静压，MPa。出砂临界生产压差分布计算结果示意图如图1所示。横坐标水平段位置具体是指从水平井跟端水平段起始位置算起的井深坐标，单位为m。
S3：根据地层静压和标定井底流压，以及水平段流动压降计算模型，计算水平段各点的实际井底流压，以及水平段各点的实际生产压差(地层静压减去井底流压)；标定井底流压即水平段根端位置的井底流压；水平段各点的实际井底流压的计算公式如式(Ⅷ)所示，水平段各点的实际生产压差的计算公式如式(Ⅸ)所示：
Pwf(i)＝Pwf0+ΔPw(i) (Ⅷ)
ΔP(i)＝Pr-Pwf(i) (Ⅸ)
式(Ⅷ)、(Ⅸ)中，Pwf0—标定井底流压，MPa；ΔPw(i)—第i点的压力与根端压力之差，MPa；i是水平井水平段按照长度计算的任一分段，1≤i≤I；分段长度为1-3m，分段总数I等于水平井水平段总长度除以分段长度；根据井筒多相流模型计算得到；Pwf(i)—第i点实际井底流压，MPa；Pr—地层静压，MPa；ΔP(i)—第i点的实际生产压差，MPa。水平井水平段实际生产压差分布计算结果示意图如图2所示。
S4：定义出砂净压差为实际生产压差与出砂临界生产压差的差值，计算得到出砂净压差，进行无量纲化计算得到出砂强度指数；定义出砂净压差为实际生产压差与出砂临界生产压差的差值，计算得到出砂净压差，如式(Ⅹ)所示：
Δp(i)＝ΔP(i)-ΔPc(i),ifΔp(i)≤0thenΔp(i)＝0 (Ⅹ)
式(Ⅹ)中，ΔPc(i)—第i点的出砂临界生产压差，MPa；ΔP(i)—第i点的出砂净压差，MPa；
计算得到出砂强度指数，计算公式如式(Ⅺ)所示：
式(Ⅺ)中，△pmax—全井段出砂净压差的最大值，Δpmax＝max{p(i)}；△pmin—全井段出砂净压差的最小值，Δpmin＝min{Δp(i)},ifΔpmin>0.25thenΔpmin＝0.25；Js(i)—第i点的出砂强度指数，无量纲。预测得到的水平井水平段出砂强度指数分布示意图如图3所示。
S5：绘制水平井相对出砂剖面图，划分出砂亏空程度等级，并实现不同等级出砂程度的分段；水平井相对出砂剖面图是指用来表达沿水平井井轴方向上出砂严重程度和亏空程度相对大小的剖面图。根据计算得到的出砂强度指数绘制，是出砂强度指数沿井深的分布曲线。包括：
以10m为最小计量井段长度间隔划分水平井水平段，当求取的某水平段的出砂强度指数不大于0.1时，将该水平段划分为不出砂井段；当求取的某水平段的出砂强度指数为0.1-0.35时，将该水平段划分为轻微出砂井段；当求取的某水平段的出砂强度指数为0.35-0.70时，将该水平段划分为中等出砂井段；当求取的某水平段的出砂强度指数为0.70-1.0时，将该水平段划分为严重出砂井段。预测得到的水平井相对出砂剖面及出砂亏空程度等级划分结果示意图如图4所示。
S6：根据出砂程度等级决定分段分级防砂的分段分级方案。包括：
对于步骤S5中划为的不出砂井段，如果不出砂井段的长度(是指出砂等级被判别为不出砂的井段的长度)大于20m，则设计为零级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中防砂等级较低的防砂等级；
对于步骤S5中划为的轻微出砂井段，如果轻微出砂井段的长度大于20m，则设计为一级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中的防砂等级较低的防砂等级；
对于步骤S5中划为的中等出砂井段，如果中等出砂井段的长度大于20m，则设计为二级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中的防砂等级较低的防砂等级；
对于步骤S5中划为的严重出砂井段，如果严重出砂井段的长度大于20m，则设计为三级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中的防砂等级较低的防砂等级。
零级防砂井段、一级防砂井段、二级防砂井段、三级防砂井段是指根据出砂井段的出砂程度和需要防砂的强度所划分的防砂等级，等级越高，表示越需要重点防砂；等级越低，表示防砂的要求越低。其中，零级防砂是指不需要防砂，具体实施时零级防砂井段不需要管柱或仅使用打孔管(在光油管上打孔孔密50-120孔/m，孔径8-10mm)；一级防砂是指筛管挡砂介质厚度不低于5-8mm或砾石层厚度不低于20mm；二级防砂是指筛管挡砂介质厚度不低于10-15mm或砾石层厚度不低于30mm；三级防砂是指筛管挡砂介质厚度不低于20-25mm或砾石层厚度不低于50mm。
水平井分级防砂多段多级划分结果示意图如图5所示。
本发明通过相对出砂剖面预测指导分段分级防砂，解决水平井笼统防砂的盲目性强、防砂效果差的问题，降低了防砂作业的成本，保证水平井长期稳定生产，提高经济效益。
实施例2
采用实施例1所述的一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，该方法应用于大港油田某油藏水平井水平段，该水平井水平段长580m，含水85％，产液量约180t/d。利用实施例1的方法预测得到水平段的出砂临界压差、生产压差、出砂净压差分布示意图如图6所示。水平段的分级防砂多段多级划分结果示意图如图7所示。
根据评价结果，整个水平段的出砂和防砂分为8段，如图7所示，A段为轻微出砂井段，需要一级防砂；B、D、F、H段为中等出砂井段，需要二级防砂；C、E、G段为严重出砂井段，需要三级防砂。经测算，按照上述方案进行分段分级防砂，可节约防砂成本15-20万元，降低全井段表皮系数约0.23左右，对应提高产量15％以上。
实施例3
采用实施例1所述的一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，该方法应用于胜利油田某稠油油藏水平井水平段，该水平井水平段穿过三个层位，水平段长分别为25m、35m和400m。
利用实施例1方法预测得到水平段的出砂剖面及分段防砂评价结果如图8所示。根据评价结果，第一层段A为中等出砂井段，需要二级防砂；第二层位B段为轻微出砂段，需要一级防砂；第三层位中C、E、G为中等出砂井段，需要二级防砂、D、F段为严重出砂井段，需要三级防砂。经测算，按照上述方案进行分段分级防砂，可节约防砂成本20-30万元，降低全井段表皮系数约0.28左右，对应提高产量18％以上。
Claims (8)
1.一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，其特征在于，包括：
S1：计算储层岩石力学参数沿水平井水平段的分布，储层岩石力学参数包括动态泊松比、静态泊松比、动态弹性模量、单轴抗压强度、内聚力以及内摩擦角；
S2：根据储层岩石力学参数及地应力分布，预测水平井水平段出砂临界生产压差分布；
S3：根据地层静压和标定井底流压，以及水平段流动压降计算模型，计算水平段各点的实际井底流压，以及水平段各点的实际生产压差；标定井底流压即水平段根端位置的井底流压；
S4：定义出砂净压差为实际生产压差与出砂临界生产压差的差值，计算得到出砂净压差，进行无量纲化计算得到出砂强度指数；
S5：绘制水平井相对出砂剖面图，划分出砂亏空程度等级，并实现不同等级出砂程度的分段；
S6：根据出砂程度等级决定分段分级防砂的分段分级方案。
2.根据权利要求1所述的一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，其特征在于，所述步骤S1中，
动态泊松比的计算公式如式(Ⅰ)所示：
式(Ⅰ)中，μd—岩石的动态泊松比，无量纲；Δth—横波时差，μs/m；Δtv—纵波时差，μs/m；
静态泊松比的计算公式如式(Ⅱ)所示：
μ＝A·μd+B (II)
式(Ⅱ)中，μ—岩石的静态泊松比，无量纲；A、B—经验系数，通过实测数据拟合得到，无量纲；A＝0.38，B＝0.082；
动态弹性模量的计算公式如式(Ⅲ)所示：
式(Ⅲ)中，Ed—岩石的动态弹性模量，MPa；ρr—地层岩石密度，kg/m3；
单轴抗压强度的计算公式如式(Ⅳ)所示：
σc＝[C·(1-VCL)+D·VCL]·Ed (Ⅳ)
式(Ⅳ)中，σc—岩石单轴抗压强度，MPa；C、D—为经验拟合系数，C＝0.00459，D＝0.00816；VCL—泥质含量，小数；
内聚力的计算公式如式(Ⅴ)所示：
式(Ⅴ)中，C0—内聚力，c/MPa；
内摩擦角的计算公式如式(Ⅵ)所示：
φf＝36.545-0.4952C0 (Ⅵ)
式(Ⅵ)中，φf—岩石内摩擦角，rad。
3.根据权利要求1所述的一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，所述步骤S2，出砂临界生产压差的计算公式如式(Ⅶ)所示：
式(Ⅶ)中，△Pc—出砂临界生产压差，MPa；C0—岩石内聚力，MPa；μ—岩石的静态泊松比，无量纲；σze—外边界处的垂向应力，MPa；β—比奥特常数，无量纲；α—失效角，rad；Pr—地层静压，MPa。
4.根据权利要求1所述的一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，所述步骤S3，水平段各点的实际井底流压的计算公式如式(Ⅷ)所示，水平段各点的实际生产压差的计算公式如式(Ⅸ)所示：
Pwf(i)＝Pwf0+ΔPw(i) (Ⅷ)
ΔP(i)＝Pr-Pwf(i) (Ⅸ)
式(Ⅷ)、(Ⅸ)中，Pwf0—标定井底流压，MPa；ΔPw(i)—第i点的压力与根端压力之差，MPa；i是水平井水平段按照长度计算的任一分段，1≤i≤I；分段总数I等于水平井水平段总长度除以分段长度，分段长度为1-3m；根据井筒多相流模型计算得到；Pwf(i)—第i点实际井底流压，MPa；Pr—地层静压，MPa；ΔP(i)—第i点的实际生产压差，MPa。
5.根据权利要求1所述的一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，所述步骤S4中，
定义出砂净压差为实际生产压差与出砂临界生产压差的差值，计算得到出砂净压差Δp(i)，如式(Ⅹ)所示：
Δp(i)＝ΔP(i)-ΔPc(i),ifΔp(i)≤0thenΔp(i)＝0 (Ⅹ)
式(Ⅹ)中，ΔPc(i)—第i点的出砂临界生产压差，MPa；ΔP(i)—第i点的出砂净压差，MPa；
计算得到出砂强度指数，计算公式如式(Ⅺ)所示：
式(Ⅺ)中，△pmax—全井段出砂净压差的最大值，Δpmax＝max{p(i)}；△pmin—全井段出砂净压差的最小值，Δpmin＝min{Δp(i)},ifΔpmin>0.25thenΔpmin＝0.25；Js(i)—第i点的出砂强度指数，无量纲。
6.根据权利要求1所述的一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，所述步骤S5，包括：
以5-15m为最小计量井段长度间隔划分水平井水平段，当求取的某水平段的出砂强度指数不大于0.1时，将该水平段划分为不出砂井段；当求取的某水平段的出砂强度指数为0.1-0.35时，将该水平段划分为轻微出砂井段；当求取的某水平段的出砂强度指数为0.35-0.70时，将该水平段划分为中等出砂井段；当求取的某水平段的出砂强度指数为0.70-1.0时，将该水平段划分为严重出砂井段。
7.根据权利要求6所述的一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，以10m为最小计量井段长度间隔划分水平井水平段。
8.根据权利要求7所述的一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法，所述步骤S6，根据出砂程度等级决定分段分级防砂的分段分级方案，包括：
对于步骤S5中划为的不出砂井段，如果不出砂井段的长度大于20m，则设计为零级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中防砂等级较低的防砂等级；
对于步骤S5中划为的轻微出砂井段，如果轻微出砂井段的长度大于20m，则设计为一级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中的防砂等级较低的防砂等级；
对于步骤S5中划为的中等出砂井段，如果中等出砂井段的长度大于20m，则设计为二级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中的防砂等级较低的防砂等级；
对于步骤S5中划为的严重出砂井段，如果严重出砂井段的长度大于20m，则设计为三级防砂井段；如果长度不大于20m，则防砂等级设置为与该井段相邻两侧井段中的防砂等级较低的防砂等级。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810013090.6A CN108331574B (zh) | 2018-01-08 | 2018-01-08 | 一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810013090.6A CN108331574B (zh) | 2018-01-08 | 2018-01-08 | 一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108331574A true CN108331574A (zh) | 2018-07-27 |
| CN108331574B CN108331574B (zh) | 2019-03-05 |
Family
ID=62924744
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810013090.6A Active CN108331574B (zh) | 2018-01-08 | 2018-01-08 | 一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108331574B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109356566A (zh) * | 2018-12-18 | 2019-02-19 | 中海石油(中国)有限公司 | 一种针对深水挥发性油田中高含水阶段自喷生产井停喷时间预测的方法 |
| CN110145286A (zh) * | 2019-05-22 | 2019-08-20 | 西南石油大学 | 一种低渗透油藏或气藏完井工程的设计方法 |
| CN110644980A (zh) * | 2019-09-11 | 2020-01-03 | 中国石油天然气股份有限公司 | 一种超低渗透油藏储层综合分类评价方法 |
| CN111594101A (zh) * | 2020-03-19 | 2020-08-28 | 中国石油大学(华东) | 一种垂直井管外出砂亏空剖面预测及高饱和砾石充填施工方法 |
| CN112240195A (zh) * | 2019-07-16 | 2021-01-19 | 中国石油大学(华东) | 基于分布式光纤声音监测的油气井出砂监测模拟实验装置及工作方法 |
| CN113591421A (zh) * | 2021-10-08 | 2021-11-02 | 中国石油大学(华东) | 非均质储层水平井砂液产出剖面动态演化模拟与预测方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010116228A2 (en) * | 2009-04-10 | 2010-10-14 | Schlumberger Technology B.V. | In-situ evaluation of reservoir sanding and fines migration and related completion, lift and surface facilities design |
| CN101915090A (zh) * | 2010-07-29 | 2010-12-15 | 中国海洋石油总公司 | 一种油气井出砂量监测系统及监测方法 |
| CN102575510A (zh) * | 2009-09-17 | 2012-07-11 | 雪佛龙美国公司 | 控制地质力学油藏系统中的出砂的计算机实施的系统和方法 |
| CN104265273A (zh) * | 2014-09-22 | 2015-01-07 | 中国石油大学(华东) | 水平井分段完井入流测试装置及测试方法 |
| CN205135615U (zh) * | 2015-09-08 | 2016-04-06 | 中国石油大学(北京) | 一种油气合采裸眼完井不同含水率出砂临界压差实验装置 |
| CN106126866A (zh) * | 2016-08-09 | 2016-11-16 | 中国石油天然气股份有限公司 | 基于地质力学模型的油气井防砂方法及装置 |
| US20170123095A1 (en) * | 2015-10-30 | 2017-05-04 | Schlumberger Technology Corporation | Method for estimating reservoir fluid characteristics in laminated media |
-
2018
- 2018-01-08 CN CN201810013090.6A patent/CN108331574B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010116228A2 (en) * | 2009-04-10 | 2010-10-14 | Schlumberger Technology B.V. | In-situ evaluation of reservoir sanding and fines migration and related completion, lift and surface facilities design |
| CN102575510A (zh) * | 2009-09-17 | 2012-07-11 | 雪佛龙美国公司 | 控制地质力学油藏系统中的出砂的计算机实施的系统和方法 |
| CN101915090A (zh) * | 2010-07-29 | 2010-12-15 | 中国海洋石油总公司 | 一种油气井出砂量监测系统及监测方法 |
| CN104265273A (zh) * | 2014-09-22 | 2015-01-07 | 中国石油大学(华东) | 水平井分段完井入流测试装置及测试方法 |
| CN205135615U (zh) * | 2015-09-08 | 2016-04-06 | 中国石油大学(北京) | 一种油气合采裸眼完井不同含水率出砂临界压差实验装置 |
| US20170123095A1 (en) * | 2015-10-30 | 2017-05-04 | Schlumberger Technology Corporation | Method for estimating reservoir fluid characteristics in laminated media |
| CN106126866A (zh) * | 2016-08-09 | 2016-11-16 | 中国石油天然气股份有限公司 | 基于地质力学模型的油气井防砂方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 任闵燕: "水平井管外砾石充填防砂完井模拟研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109356566A (zh) * | 2018-12-18 | 2019-02-19 | 中海石油(中国)有限公司 | 一种针对深水挥发性油田中高含水阶段自喷生产井停喷时间预测的方法 |
| CN109356566B (zh) * | 2018-12-18 | 2022-02-08 | 中海石油(中国)有限公司 | 一种针对深水挥发性油田中高含水阶段自喷生产井停喷时间预测的方法 |
| CN110145286A (zh) * | 2019-05-22 | 2019-08-20 | 西南石油大学 | 一种低渗透油藏或气藏完井工程的设计方法 |
| CN110145286B (zh) * | 2019-05-22 | 2021-10-01 | 西南石油大学 | 一种低渗透油藏或气藏完井工程的设计方法 |
| CN112240195A (zh) * | 2019-07-16 | 2021-01-19 | 中国石油大学(华东) | 基于分布式光纤声音监测的油气井出砂监测模拟实验装置及工作方法 |
| CN112240195B (zh) * | 2019-07-16 | 2024-01-30 | 中国石油大学(华东) | 基于分布式光纤声音监测的油气井出砂监测模拟实验装置及工作方法 |
| CN110644980A (zh) * | 2019-09-11 | 2020-01-03 | 中国石油天然气股份有限公司 | 一种超低渗透油藏储层综合分类评价方法 |
| CN111594101A (zh) * | 2020-03-19 | 2020-08-28 | 中国石油大学(华东) | 一种垂直井管外出砂亏空剖面预测及高饱和砾石充填施工方法 |
| CN111594101B (zh) * | 2020-03-19 | 2021-03-12 | 中国石油大学(华东) | 一种垂直井管外出砂亏空剖面预测及高饱和砾石充填施工方法 |
| CN113591421A (zh) * | 2021-10-08 | 2021-11-02 | 中国石油大学(华东) | 非均质储层水平井砂液产出剖面动态演化模拟与预测方法 |
| CN113591421B (zh) * | 2021-10-08 | 2021-12-07 | 中国石油大学(华东) | 非均质储层水平井砂液产出剖面动态演化模拟与预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108331574B (zh) | 2019-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108331574B (zh) | 一种水平井水平段相对出砂亏空剖面预测及防砂分段分级方法 | |
| CN104504230B (zh) | 低渗气井采收率与极限动用半径估计方法 | |
| CN108561106B (zh) | 一种人工缝控储量提高采收率的油气开采方法 | |
| CN105243182B (zh) | 致密油压裂水平井的动态储量计算方法 | |
| CN105386751B (zh) | 一种基于油藏渗流模型的水平井测井产能预测方法 | |
| CN105093313B (zh) | 一种岩溶型油藏单井油气产能预测方法及装置 | |
| CN108301825B (zh) | 一种高温高压储层孔隙压力分析方法 | |
| CN106150477A (zh) | 一种确定单井控制储量的方法 | |
| Stein et al. | Estimating maximum sand-free production rates from friable sands for different well completion geometries | |
| CN108694254A (zh) | 一种变产变压生产气井产量经验递减曲线分析方法 | |
| CN106246158B (zh) | 超深低孔裂缝性砂岩气藏布井方法及装置 | |
| CN109025965A (zh) | 一种确定超低渗致密储层注水开发渗透率下限的方法 | |
| CN110984974B (zh) | 基于水侵速度确定有水气藏合理采气速度的方法 | |
| Pham | Estimation of sand production rate using geomechanical and hydromechanical models | |
| CN116050626B (zh) | 一种多级压裂水平井生产分级优化与效果预测方法 | |
| CN112949901A (zh) | 页岩油藏水平井准自然能量开发方法及系统 | |
| CN114169559A (zh) | 一种硫沉积储层伤害时空分布特征预测方法 | |
| ZHOU et al. | Application of multilateral wells with limited sand production to heavy oil reservoirs | |
| CN108166973A (zh) | 一种变压力场下非常规气井合理生产压差确定方法 | |
| Yong et al. | A method of production profile quick prediction based on typical curves: A case study of the Upper Shale large multi-layered sandstone reservoir, Rumaila Oilfield, Iraq | |
| Su et al. | A new method to calculate sweep efficiency of horizontal wells in heterogeneous reservoir | |
| Jing et al. | Productivity prediction of fractured horizontal wells with low permeability flow characteristics | |
| Wang et al. | Evaluation and Recognition of Fracturing Effect in Ultra-Low Permeability Volatile Reservoir | |
| CN111651848B (zh) | 应力敏感性油藏直井产液能力预测方法及装置 | |
| Jiang et al. | An Optimal Combination Method of Injection Interval Based on Seepage Resistance |
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 | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20180727 Assignee: Tianjin Kailei Oilfield Technology Co.,Ltd. Assignor: CHINA University OF PETROLEUM (EAST CHINA) Contract record no.: X2023980035116 Denomination of invention: A Method for Predicting the Relative Sand Production Deficiency Profile and Grading Sand Control Segments in Horizontal Wells Granted publication date: 20190305 License type: Common License Record date: 20230427 |