CN113901378A - 压力梯度影响因素权重系数确定方法 - Google Patents
压力梯度影响因素权重系数确定方法 Download PDFInfo
- Publication number
- CN113901378A CN113901378A CN202010577036.1A CN202010577036A CN113901378A CN 113901378 A CN113901378 A CN 113901378A CN 202010577036 A CN202010577036 A CN 202010577036A CN 113901378 A CN113901378 A CN 113901378A
- Authority
- CN
- China
- Prior art keywords
- pressure gradient
- formula
- factor
- formation
- weight coefficient
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000004519 manufacturing process Methods 0.000 claims abstract description 38
- 238000002347 injection Methods 0.000 claims abstract description 12
- 239000007924 injection Substances 0.000 claims abstract description 12
- 238000004088 simulation Methods 0.000 claims abstract description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 42
- 230000035699 permeability Effects 0.000 claims description 27
- 239000010779 crude oil Substances 0.000 claims description 25
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 25
- 238000010606 normalization Methods 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 15
- 150000001875 compounds Chemical class 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000012216 screening Methods 0.000 abstract description 5
- 239000012530 fluid Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 1
- 230000000704 physical effect 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- 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
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供一种压力梯度影响因素权重系数确定方法,该压力梯度影响因素权重系数确定方法包括:步骤1,选取影响压力梯度的主要因素;步骤2,确定排列组合方案,并利用数值模拟计算不同方案的注采井间最小压力梯度;步骤3,对各方案参数进行无量纲化;步骤4,利用灰色关联度方法计算各影响因素的关联度;步骤5,将各影响因素的关联度转化为各因素对压力梯度影响的权重系数。该压力梯度影响因素权重系数确定方法针对不同影响因素参数变化水平,评判区域压力梯度分布状况,为后期井位筛选提供参考。
Description
技术领域
本发明涉及油田开发技术领域,特别是涉及到一种压力梯度影响因素权重系数确定方法。
背景技术
低渗透油藏物性差,渗透率级差大,平面动用极不均衡,而压力梯度作为现场难以获取的参数,在很大程度上决定着单井的有效动用范围,影响着单井动用界限,因此有必要确定影响压力梯度的各因素重要性。目前国内外关于注采井间压力梯度的计算多采用解析方法,在均质储层条件下得到注采井间压力梯度的变化规律,而未考虑多因素不同参数水平情况下的压力梯度变化规律,且没有对影响启动压力梯度的各因素的重要程度进行论证。
为此我们发明了一种新的压力梯度影响因素权重系数确定方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种可以明晰多因素影响下的压力梯度复杂分布状况,对影响压力梯度的各因素进行权重系数赋值,为后期新井位筛选提供参考的压力梯度影响因素权重系数确定方法。
本发明的目的可通过如下技术措施来实现:压力梯度影响因素权重系数确定方法,该压力梯度影响因素权重系数确定方法包括:步骤1,选取影响压力梯度的主要因素;步骤2,确定排列组合方案,并利用数值模拟计算不同方案的注采井间最小压力梯度;步骤3,对各方案参数进行无量纲化;步骤4,利用灰色关联度方法计算各影响因素的关联度;步骤5,将各影响因素的关联度转化为各因素对压力梯度影响的权重系数。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,选取影响压力梯度的主要因素:地层原油粘度、地层渗透率、地层厚度、含水饱和度和注采井距,根据目标区域取心井的岩心分析资料和地质研究工作,得到各参数取值范围,并根据各参数的数据集中度确定各参数的因素水平,以每个参数的数据集中度确定各参数基准值。
在步骤2中,依据各参数的因素水平,利用各参数的非基准值与其它参数的基准值确定排列组合方案,并利用数值模拟计算不同方案的注采井间最小压力梯度。
地层渗透率的无量纲化公式为:k标准=(k-k平均)/Sk
地层厚度的无量纲化公式为:h标准=(h-h平均)/Sh
注采井距的无量纲化公式为:L标准=(L-L平均)/SL
其中:
x标准=(xi-x平均)/S
式中,—地层原油粘度标准化后的变量值;μo—地层原油粘度,mPa·s;—地层原油粘度算数平均值,mPa·s;—地层原油粘度标准差;k标准—地层渗透率标准化后的变量值;k—地层渗透率,10-3μm2;k平均—地层渗透率算数平均值,10-3μm2;Sk—地层渗透率标准差;h标准—地层厚度标准化后的变量值;h—地层厚度,10-3μm2;h平均—地层厚度算数平均值,10-3μm2;Sh—地层厚度标准差;—含水饱和度标准化后的变量值;Sw—含水饱和度,%;—含水饱和度算数平均值,10-3μm2;—含水饱和度标准差;L标准—注采井距标准化后的变量值;L—注采井距,m;L平均—注采井距算数平均值,10-3μm2;SL—注采井距标准差;—最小压力梯度标准化后的变量值;Gp—最小压力梯度,MPa;—最小压力梯度算数平均值,MPa;—最小压力梯度标准差;x平均—x因素的算数平均值;xi—x因素值;S—x因素标准差;x标准—x因素标准化后的变量值。
在步骤4中,灰色关联系数的计算公式为:
x0(k)是各影响因素的值,xi(k)是压力梯度的值。
关联度的计算公式为:
根据灰色关联度的计算公式,计算可得到各影响因素的关联度;
式中,x0(k)—各影响因素的值,xi(k)—压力梯度的值;ξi(k)—各因素的灰色关联系数;ri—各因素的灰色关联度。
在步骤5中,权重系数的计算公式为:
ri—各因素的灰色关联度;αi—各因素对压力梯度影响的权重系数;
根据上述公式计算得到权重系数表。
本发明中的压力梯度影响因素权重系数确定方法,主要用于计算影响压力梯度的因素权重系数,可用于后期新井井位筛选,优化最佳井位。该方法选取影响压力梯度的主要因素:地层原油粘度、地层渗透率、地层厚度、含水饱和度和注采井距等,根据各因素实际区块的参数取值范围以及数据集中度,确定各因素的参数取值,在各参数取值不同排列组合情况下计算注采井距间最小压力梯度,并进行无量纲化,根据无量纲值计算各因素的关联度,得到最终的影响压力梯度的各因素的权重系数。压力梯度的分布对极限泄油面积起着至关重要的作用,因此需要一种压力梯度影响因素权重确定方法,针对不同影响因素参数变化水平,评判区域压力梯度分布状况,为后期井位筛选提供参考。
附图说明
图1为本发明的压力梯度影响因素权重系数确定方法的一具体实施例的流程图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的压力梯度影响因素权重系数确定方法的流程图。
步骤101,选取影响压力梯度的主要因素:地层原油粘度、地层渗透率、地层厚度、含水饱和度和注采井距等,根据目标区域取心井的岩心分析资料和地质研究工作,得到各参数取值范围。并根据各参数的数据集中度确定各参数的因素水平,以每个参数的数据集中度确定各参数基准值。
表1参数取值范围和基准值确定表
表中,—地层原油粘度分布范围,mPa·s;k1-k4—地层渗透率分布范围,10-3μm2;h1-h4—地层厚度分布范围,m;—含水饱和度分布范围,%;L1-L4—注采井距分布范围,m;k1、h2、L1—根据地层原油粘度、地层渗透率、地层厚度、含水饱和度和注采井距的数据集中度确定的各参数基准值。
步骤102,依据各参数的因素水平,利用各参数的非基准值与其它参数的基准值确定排列组合方案,并利用数值模拟计算不同方案的注采井间最小压力梯度。
表2各参数排列组合方案表
表中,—地层原油粘度排列值;k1、k2、k3、k4—地层渗透率排列值;h1、h2、h3、h4—地层厚度排列值;—含水饱和度排列值;L1、L2、L3、L4—注采井距排列值; —不同排列组合方案计算所得最小压力梯度。
步骤103,对各方案参数进行无量纲化。
无量纲化方法见表3。
表3无量纲化方法表
表中,—地层原油粘度标准化后的变量值;μo—地层原油粘度,mPa·s;—地层原油粘度算数平均值,mPa·s;—地层原油粘度标准差;k标准—地层渗透率标准化后的变量值;k—地层渗透率,10-3μm2;k平均—地层渗透率算数平均值,10-3μm2;Sk—地层渗透率标准差;h标准—地层厚度标准化后的变量值;h—地层厚度,10-3μm2;h平均—地层厚度算数平均值,10-3μm2;Sh—地层厚度标准差;—含水饱和度标准化后的变量值;Sw—含水饱和度,%;—含水饱和度算数平均值,10-3μm2;—含水饱和度标准差;L标准—注采井距标准化后的变量值;L—注采井距,m;L平均—注采井距算数平均值,10-3μm2;SL—注采井距标准差;—最小压力梯度标准化后的变量值;Gp—最小压力梯度,MPa;—最小压力梯度算数平均值,MPa;—最小压力梯度标准差。
其中:
x标准=(xi-x平均)/S
式中,x平均—x因素的算数平均值;xi—x因素值;S—x因素标准差;x标准—x因素标准化后的变量值。
步骤104,利用灰色关联度方法计算各影响因素的关联度。
灰色关联系数的计算公式为:
关联度的计算公式为:
根据灰色关联度的计算公式,计算可得到各影响因素的关联度。
式中,x0(k)—各影响因素的值,xi(k)—压力梯度的值;ξi(k)—各因素的灰色关联系数;ri—各因素的灰色关联度。
表4各影响因素关联度表
表中,r1—地层原油粘度的灰色关联度;r2—地层渗透率的灰色关联度;r3—地层厚度的灰色关联度;r4—含水饱和度的灰色关联度;r5—注采井距的灰色关联度。
步骤105,将各影响因素的关联度转化为各因素对压力梯度影响的权重系数。
权重系数的计算公式为:
ri—各因素的灰色关联度;αi—各因素对压力梯度影响的权重系数。
根据上述公式计算得到权重系数表。
表5各影响因素权重系数表
表中,α1—地层原油粘度对压力梯度影响的权重系数;α2—地层渗透率对压力梯度影响的权重系数;α3—地层厚度对压力梯度影响的权重系数;α4—含水饱和度对压力梯度影响的权重系数;α5—注采井距对压力梯度影响的权重系数。
在应用本发明的一具体实施例中,包括了以下步骤:
1.选取影响压力梯度的主要因素:原油粘度、地层渗透率、地层厚度、含水饱和度和注采井距等,根据目标区域取心井的岩心分析资料和地质研究工作,得到各参数取值范围(表6)。并根据各参数的数据集中度确定各参数的因素水平,以每个参数的数据集中度确定各参数基准值(表6)。
表6参数取值范围和基准值确定表
2)依据各参数的因素水平,利用各参数的非基准值与其它参数的基准值确定排列组合方案,并利用数值模拟计算不同方案的注采井间最小压力梯度。
表7各参数排列组合方案表
3)对各方案参数进行无量纲化。
无量纲化方法见表8。
表8无量纲化方法表
其中:
x标准=(xi-x平均)/S
式中,—地层原油粘度标准化后的变量值;μo—地层原油粘度,mPa·s;—地层原油粘度算数平均值,mPa·s;—地层原油粘度标准差;k标准—地层渗透率标准化后的变量值;k—地层渗透率,10-3μm2;k平均—地层渗透率算数平均值,10-3μm2;Sk—地层渗透率标准差;h标准—地层厚度标准化后的变量值;h—地层厚度,10-3μm2;h平均—地层厚度算数平均值,10-3μm2;Sh—地层厚度标准差;—含水饱和度标准化后的变量值;Sw—含水饱和度,%;—含水饱和度算数平均值,10-3μm2;—含水饱和度标准差;L标准—注采井距标准化后的变量值;L—注采井距,m;L平均—注采井距算数平均值,10-3μm2;SL—注采井距标准差;—最小压力梯度标准化后的变量值;Gp—最小压力梯度,MPa;—最小压力梯度算数平均值,MPa;—最小压力梯度标准差;x平均—x因素的算数平均值;xi—x因素值;S—x因素标准差;x标准—x因素标准化后的变量值。
无量纲化后的数据见表9。
表9无量纲化数据表
4)利用灰色关联度方法计算各影响因素的关联度。
灰色关联系数的计算公式为:
关联度的计算公式为:
式中,x0(k)—各影响因素的值,xi(k)—压力梯度的值;ξi(k)—各因素的灰色关联系数;ri—各因素的灰色关联度。
根据灰色关联度的计算公式,计算可得到各影响因素的关联度。
表10各影响因素关联度表
5)将各影响因素的关联度转化为各因素对压力梯度影响的权重系数。
权重系数的计算公式为:
式中,ri—各因素的灰色关联度;αi—各因素对压力梯度影响的权重系数。
根据上述公式计算得到权重系数表。
表11各影响因素权重系数表
地层原油粘度 | 地层渗透率 | 地层厚度 | 含水饱和度 | 注采井距 |
0.200 | 0.195 | 0.188 | 0.226 | 0.191 |
本发明综合考虑影响压力梯度的各种因素确定其权重系数,压力梯度作为现场难以获取的参数,在很大程度上决定着单井的有效动用范围,影响着单井动用界限,该方法根据实际区块的参数范围和数据集中度,确定不同方案得到各参数,进行无量纲化和关联度计算,最终得到影响压力梯度的各因素的权重系数。该方法对于多因素的复杂变化情况,可以根据权重系数的排序,为后期新井位筛选提供参考。
Claims (6)
1.压力梯度影响因素权重系数确定方法,其特征在于,该压力梯度影响因素权重系数确定方法包括:
步骤1,选取影响压力梯度的主要因素;
步骤2,确定排列组合方案,并利用数值模拟计算不同方案的注采井间最小压力梯度;
步骤3,对各方案参数进行无量纲化;
步骤4,利用灰色关联度方法计算各影响因素的关联度;
步骤5,将各影响因素的关联度转化为各因素对压力梯度影响的权重系数。
2.根据权利要求1所述的压力梯度影响因素权重系数确定方法,其特征在于,在步骤1中,选取影响压力梯度的主要因素:原油粘度、地层渗透率、地层厚度、含水饱和度和注采井距,根据目标区域取心井的岩心分析资料和地质研究工作,得到各参数取值范围,并根据各参数的数据集中度确定各参数的因素水平,以每个参数的数据集中度确定各参数基准值。
3.根据权利要求1所述的压力梯度影响因素权重系数确定方法,其特征在于,在步骤2中,依据各参数的因素水平,利用各参数的非基准值与其它参数的基准值确定排列组合方案,并利用数值模拟计算不同方案的注采井间最小压力梯度。
4.根据权利要求1所述的压力梯度影响因素权重系数确定方法,其特征在于,在步骤3中,
地层厚度的无量纲化公式为:h标准=(h-h平均)/Sh
注采井距的无量纲化公式为:L标准=(L-L平均)/SL
其中:
x标准=(xi-x平均)/S
式中,—地层原油粘度标准化后的变量值;μo—地层原油粘度,mPa·s;—地层原油粘度算数平均值,mPa·s;—地层原油粘度标准差;k标准—地层渗透率标准化后的变量值;k—地层渗透率,10-3μm2;k平均—地层渗透率算数平均值,10-3μm2;Sk—地层渗透率标准差;h标准—地层厚度标准化后的变量值;h—地层厚度,10-3μm2;h平均—地层厚度算数平均值,10-3μm2;Sh—地层厚度标准差;—含水饱和度标准化后的变量值;Sw—含水饱和度,%;—含水饱和度算数平均值,10-3μm2;—含水饱和度标准差;L标准—注采井距标准化后的变量值;L—注采井距,m;L平均—注采井距算数平均值,10-3μm2;SL—注采井距标准差;—最小压力梯度标准化后的变量值;Gp—最小压力梯度,MPa;—最小压力梯度算数平均值,MPa;—最小压力梯度标准差;x平均—x因素的算数平均值;xi—x因素值;S—x因素标准差;x标准—x因素标准化后的变量值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010577036.1A CN113901378A (zh) | 2020-06-22 | 2020-06-22 | 压力梯度影响因素权重系数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010577036.1A CN113901378A (zh) | 2020-06-22 | 2020-06-22 | 压力梯度影响因素权重系数确定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113901378A true CN113901378A (zh) | 2022-01-07 |
Family
ID=79186571
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010577036.1A Pending CN113901378A (zh) | 2020-06-22 | 2020-06-22 | 压力梯度影响因素权重系数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113901378A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116087053A (zh) * | 2022-12-15 | 2023-05-09 | 长江大学 | 一种研究稠油启动压力梯度对剩余油分布影响的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110766293A (zh) * | 2019-09-30 | 2020-02-07 | 煤炭科学技术研究院有限公司 | 基于灰色系统理论的覆岩裂隙带影响因素权重分析方法 |
-
2020
- 2020-06-22 CN CN202010577036.1A patent/CN113901378A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110766293A (zh) * | 2019-09-30 | 2020-02-07 | 煤炭科学技术研究院有限公司 | 基于灰色系统理论的覆岩裂隙带影响因素权重分析方法 |
Non-Patent Citations (3)
Title |
---|
倪淑敏;韩晓冰;赵肖;周晋冲;刘奎;: "启动压力梯度影响因素的灰色关联分析", 化工管理, no. 18, 21 June 2017 (2017-06-21), pages 49 * |
吴世东: "南七聚驱不同阶段合理压力梯度界限确定", 中国优秀硕士学位论文全文数据库工程科技Ⅰ辑, 15 February 2018 (2018-02-15), pages 019 - 363 * |
黄昕;刘思良;岳三琪;付玉;朱泽正;荐小纯;: "地下储气库地层压力影响因素灰色关联分析", 石油化工应用, no. 08, 25 August 2017 (2017-08-25), pages 50 - 53 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116087053A (zh) * | 2022-12-15 | 2023-05-09 | 长江大学 | 一种研究稠油启动压力梯度对剩余油分布影响的方法 |
CN116087053B (zh) * | 2022-12-15 | 2024-05-07 | 长江大学 | 一种研究稠油启动压力梯度对剩余油分布影响的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110644980B (zh) | 一种超低渗透油藏储层综合分类评价方法 | |
CN111425193B (zh) | 一种基于聚类分析测井岩石物理相划分的储层可压性评价方法 | |
CN108104806B (zh) | 剩余油分布规律定量分析方法 | |
CN107387051B (zh) | 低渗透非均质油藏多段压裂水平井重复压裂选井的方法 | |
CN106569288B (zh) | 基于储层裂缝有效性聚类分析的裂缝性储层质量评价方法 | |
CN106600443A (zh) | 动态含水饱和度的油井产量劈分方法 | |
CN111173507B (zh) | 高含水油田剩余油预测方法 | |
CN108150160B (zh) | 一种地层中欠压实作用超压的求取方法 | |
CN111456709B (zh) | 一种基于测井曲线的水平井多级压裂分段分簇方法 | |
CN108133087A (zh) | 一种气-水两相渗流应力敏感储层原始渗透率反演方法 | |
CN108444884A (zh) | 基于Harris模型的油相相对渗透率曲线表征方法 | |
CN110552691B (zh) | 一种油藏裂缝优势通道识别方法及其系统 | |
CN113901378A (zh) | 压力梯度影响因素权重系数确定方法 | |
CN113762614B (zh) | 一种基于静态油气田数据的页岩气井产能预测的方法 | |
CN111749688B (zh) | 一种优势渗流通道发育层位和方向的预测方法 | |
CN112257349B (zh) | 一种判断致密砂岩可动水气藏气井是否具有开发价值的方法 | |
CN109522579B (zh) | 水平井压裂施工破裂压力预测方法 | |
CN107992690A (zh) | 一种诱导应力干扰下多裂缝扩展形态均衡程度的评价方法 | |
CN110909466B (zh) | 提高差孔隙分选储层平均毛管压力曲线计算精度的方法 | |
CN106845685B (zh) | 一种降低水平井试油压裂作业成本的参数优化方法 | |
CN112796749B (zh) | 一种基于毛管数的油藏水驱数值模拟方法 | |
CN115705452A (zh) | 整装砂岩油藏开发中后期的新型采收率预测方法 | |
CN104747154A (zh) | 一种利用驱油效率比提高蒸汽驱剩余油研究精度的方法 | |
CN109869144A (zh) | 一种辨别油藏高压物性参数真伪的方法 | |
CN113094864B (zh) | 强非均质碳酸盐岩水平井分段设计方法 |
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 |