CN109858177A - 一种基于拟流线法的底水油藏水平井水驱数值模拟方法 - Google Patents
一种基于拟流线法的底水油藏水平井水驱数值模拟方法 Download PDFInfo
- Publication number
- CN109858177A CN109858177A CN201910134624.5A CN201910134624A CN109858177A CN 109858177 A CN109858177 A CN 109858177A CN 201910134624 A CN201910134624 A CN 201910134624A CN 109858177 A CN109858177 A CN 109858177A
- Authority
- CN
- China
- Prior art keywords
- quasi
- stream
- flow tube
- grid
- 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.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000004088 simulation Methods 0.000 title claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 238000006073 displacement reaction Methods 0.000 claims abstract description 7
- 230000003068 static effect Effects 0.000 claims abstract description 5
- 239000003921 oil Substances 0.000 claims description 15
- 230000035699 permeability Effects 0.000 claims description 8
- 239000010779 crude oil Substances 0.000 claims description 4
- 238000005325 percolation Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- ONUFESLQCSAYKA-UHFFFAOYSA-N iprodione Chemical compound O=C1N(C(=O)NC(C)C)CC(=O)N1C1=CC(Cl)=CC(Cl)=C1 ONUFESLQCSAYKA-UHFFFAOYSA-N 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 abstract description 4
- 238000011161 development Methods 0.000 abstract description 3
- 238000009472 formulation Methods 0.000 abstract 1
- 239000000203 mixture Substances 0.000 abstract 1
- 230000008901 benefit Effects 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000001764 infiltration Methods 0.000 description 2
- 230000008595 infiltration Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于拟流线法的底水油藏水平井水驱数值模拟方法,该方法包括:步骤1,收集油藏基本静态和动态信息,为数值模拟计算提供参数;步骤2,根据相渗实验数据,计算各含水饱和度下的Kro、Krw、fw、fw’、PV等参数表;步骤3,根据底水水脊形态,确定拟流线分布,沿流线进行等体积网格剖分;步骤4,模型初始化,根据油藏实际情况,设定网格初始饱和度、压力边界等;步骤5,根据水驱前缘驱替理论,计算饱和度场分布,渗流阻力等模拟参数;该发明针对底水油藏水驱油规律的认识难题,采用拟流线法建立底水驱理论模型,对底水油藏水平井开发动态进行仿真模拟,为底水油藏开发技术政策的制定提供理论和技术支持。
Description
技术领域
本发明涉及石油勘探与开发的油藏工程领域,特别是涉及到底水油藏水平井开采水驱油过程数值模拟计算方法。
背景技术
底水油藏水平井开采过程中原油主要受垂向底水的驱动向井筒汇聚,驱替范围主要在油层纵向区域,底水距油层纵向距离远远小于常规注采井距,生产动态上,底水油藏常常表现出含水上升快,递减大的特点。底水油藏水平井开发的关键就是确定水平井底水水脊波及半径,确定波及区域的饱和度分布,以指导水平井井距的设计,国内外确定底水水脊动态变化的方法主要是采用数值模拟方法,常规数值模拟方法具有严格的数学物理理论基础,但由于需要大量的数学计算,成本太高,同时成果交付周期长,为此,我们发明了一种基于拟流线法的底水油藏水平井水驱数值模拟方法,解决了以上问题。
发明内容
本发明的目的是提供一种基于拟流线法的底水油藏水平井水驱数值模拟方法,该方法基于拟流线法建立了底水油藏渗流数值模型,解决了常规数值模拟存在数值弥散的问题。
本发明的目的可通过如下技术措施来实现:
一种基于拟流线法的底水油藏水平井水驱数值模拟方法,包括如下步骤:
步骤1、收集油藏水平井基本静态和动态信息获取数值模拟计算参数;
步骤2、对数值模拟计算参数中的原始相渗数值进行归一化计算获得模拟特征参数;
步骤3、根据模拟特征参数按照拟流线方法对拟流区域进行等体积网格划分建立拟流区域网格模型;
步骤4、对拟流区域模型中每个单独拟流区域进行驱替计算;
步骤5、计算整个拟流区域模型二维流场含水率叠加生成井口含水率;
步骤6、判断井口含水率是否大于或等于截止条件;如果满足,则进入下一步;否则返回步骤4;
步骤7、输出整个拟流区域模型数值。
所述步骤4对拟流区域模型中每个拟流区域进行驱替计算过程:
4.1、根据拟流区域模型初始化,按照如下公式计算拟流区域模型初始状态下的渗流阻力,计算公式:
式中,Rl为流管l的渗流阻力,MPa﹒d﹒m-3;Ln为流管l的第n个网格的长度,m;An为流管l的第n个网格的截面积,m2;kron为流管l的第n个网格的油相相对渗透率,mD;krwn为流管l的第n个网格的水相相对渗透率,mD;μo为原油地层粘度,mPa﹒s;μw为水的地层粘度,mPa﹒s;
4.2、设定生产压差为控制条件,计算第一个时间步的流量,计算公式:
式中,为流管l在t1时刻的流量,m3/d;pe为流管入口端压力,MPa;pf为流管出口端压力,MPa;
4.3、根据水驱前缘追踪理论,代入当前时间步下的流量,计算水驱前缘沿流管推进的距离,计算公式:
式中,ξ为表示沿流管的网格标号;fw表示出口端含水率,小数;Sw表示出口端含水饱和度,小数;
4.4、根据等饱和度面移动方程,计算拟流区域起点与水驱前缘之间的饱和度分布,建立拟流区域的饱和度剖面;
4.5、根据拟流区域出口处网格的饱和度,根据相渗曲线分流量方程计算出口端含水率。
4.6、按照上面的计算步骤对每个拟流区域进行计算,得到各个拟流区域出口处的流量、含水率及流管的饱和度分布,按照如下公式并对拟流区域出口处指标进行求和,即可得到整管个二维流场的含水率;
式中,f为二维流场的总含水率,小数;fwl为流管l的出口端含水率,小数。
本发明与现有技术相比,具有如下的优点和有益效果:
(1)本发明基于拟流线法建立了底水油藏渗流数值模型,解决了常规数值模拟存在数值弥散的问题。
(2)本发明首次基于流管法模型,建立底水油藏水平井渗流数学模型,填补了现有数值模拟方法研究的空白。
(3)本发明所述的数值模拟方法简便,计算结果更准确更接近实际情况,为今后底水油藏数值模拟技术的发展奠定基础。
附图说明
图1,是本发明方法的流程图。
图2,是本发明方法涉及具体步骤的框图。
图3,是底水油藏水平井数值模拟网格剖分图。
图4,是本发明油藏数值模拟含油饱和度图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
如图1所示,本发明提供一种基于拟流线法的底水油藏水平井水驱数值模拟方法,该方法包括:
步骤1,收集油藏基本静态和动态信息,为数值模拟计算提供参数;在步骤1中,收集的油藏基本静态和动态信息包括原始相渗数据、地层压力、油气水高压物性参数、储层平面和垂向平均渗透率、油层平均孔隙度、油层平均含水饱和度、油井油气水产量数据。
步骤2,根据步骤1中收集到的原始相渗数据,进行归一化处理,在标准化相渗曲线上,将横坐标从0到1划分为n等分,利用插值法求取各含水饱和度下Kro、Krw,利用分流量公式求取fw、fw’、PV等参数。
步骤3,油藏水平井数值模型按照拟流线方法进行等体积网格划分建立拟流区域模型;根据底水水脊形态,确定拟流线分布,沿流线进行等体积网格剖分;如图3所示。在步骤3中,根据水脊剖面描述公式,沿油水界面等距离确定拟流线,在两流线之间的区域即为拟流管,沿流管方向进行等体积剖分。水脊剖面表达式为:
式中,a1、a2、a3为常量系数,可通过机理模型多因素回归得到。
其中,包括对拟流区域模型初始化,根据油藏实际情况,设定网格初始饱和度、压力边界等;根据步骤1中收集的地层压力、储层平面和垂向平均渗透率、油层平均孔隙度、油层平均含水饱和度对网格进行赋值。
步骤4,对拟流区域模型中每个拟流区域进行驱前计算;根据水驱前缘驱替理论,对拟流区域模块计算获得整个油藏水平井的二维流场的含水率,即饱和度场分布,渗流阻力等模拟参数。包含如下过程:如图2所示,
4.1、根据拟流区域模型初始化,按照如下公式计算拟流区域模型初始状态下的渗流阻力,计算公式:
式中,Rl为流管l的渗流阻力,MPa﹒d﹒m-3;Ln为流管l的第n个网格的长度,m;An为流管l的第n个网格的截面积,m2;kron为流管l的第n个网格的油相相对渗透率,mD;krwn为流管l的第n个网格的水相相对渗透率,mD;μo为原油地层粘度,mPa﹒s;μw为水的地层粘度,mPa﹒s;
4.2、设定生产压差为控制条件,计算第一个时间步的流量,计算公式:
式中,为流管l在t1时刻的流量,m3/d;pe为流管入口端压力,MPa;pf为流管出口端压力,MPa;
4.3、根据水驱前缘追踪理论,代入当前时间步下的流量,计算水驱前缘沿流管推进的距离,计算公式:
式中,ξ为表示沿流管的网格标号;fw表示出口端含水率,小数;Sw表示出口端含水饱和度,小数;
4.4、根据等饱和度面移动方程,计算拟流区域起点与水驱前缘之间的饱和度分布,建立拟流区域的饱和度剖面;
4.5、根据拟流区域出口处网格的饱和度,根据相渗曲线分流量方程计算出口端含水率。
4.6、按照上面的计算步骤对每个拟流区域进行计算,得到各个拟流区域出口处的流量、含水率及流管的饱和度分布,按照如下公式并对拟流区域出口处指标进行求和,即可得到整管个二维流场的含水率;
式中,f为二维流场的总含水率,小数;fwl为流管l的出口端含水率,小数。
步骤5,计算整个拟流区域模型二维流场的含水率,即拟流区域模型由沿轴向的二维流场微元构成,对二维流场的含水率沿水平井井筒方向叠加即可得到整个水平井井口的流量、含水率指标;
步骤6,判断井口含水率是否大于或等于98%;如果满足,则进入下一步;否则返回步骤1进行循环计算;即拟流区域模型的饱和度分布,计算当前时间步水驱前缘沿流管推进的距离、沿流管的饱和度分布、流管出口处网格饱和度、流管出口端含水率等指标,直至水平井井口含水率达到设定的截止条件,整个循环计算结束。
步骤7、输出整个拟流区域模型数值,即如图4所示。
Claims (2)
1.一种基于拟流线法的底水油藏水平井水驱数值模拟方法,其特征在于,包括如下步骤:
步骤1、收集油藏水平井基本静态和动态信息获取数值模拟计算参数;
步骤2、对数值模拟计算参数中的原始相渗数值进行归一化计算获得模拟特征参数;
步骤3、根据模拟特征参数按照拟流线方法对拟流区域进行等体积网格划分建立拟流区域网格模型;
步骤4、对拟流区域模型中每个单独拟流区域进行驱替计算;
步骤5、计算整个拟流区域模型二维流场含水率叠加生成井口含水率;
步骤6、判断井口含水率是否大于或等于截止条件;如果满足,则进入下一步;否则返回步骤4;
步骤7、输出整个拟流区域模型数值。
2.根据权利要求1所述一种基于拟流线法的底水油藏水平井水驱数值模拟方法,其特征在于,所述步骤4对拟流区域模型中每个拟流区域进行驱替计算过程:
4.1、根据拟流区域模型初始化,按照如下公式计算拟流区域模型初始状态下的渗流阻力,计算公式:
式中,Rl为流管l的渗流阻力,MPa﹒d﹒m-3;Ln为流管l的第n个网格的长度,m;An为流管l的第n个网格的截面积,m2;kron为流管l的第n个网格的油相相对渗透率,mD;krwn为流管l的第n个网格的水相相对渗透率,mD;μo为原油地层粘度,mPa﹒s;μw为水的地层粘度,mPa﹒s;
4.2、设定生产压差为控制条件,计算第一个时间步的流量,计算公式:
式中,qt1l为流管l在t1时刻的流量,m3/d;pe为流管入口端压力,MPa;pf为流管出口端压力,MPa;
4.3、根据水驱前缘追踪理论,代入当前时间步下的流量,计算水驱前缘沿流管推进的距离,计算公式:
式中,ξ为表示沿流管的网格标号;fw表示出口端含水率,小数;Sw表示出口端含水饱和度,小数;
4.4、根据等饱和度面移动方程,计算拟流区域起点与水驱前缘之间的饱和度分布,建立拟流区域的饱和度剖面;
4.5、根据拟流区域出口处网格的饱和度,根据相渗曲线分流量方程计算出口端含水率。
4.6、按照上面的计算步骤对每个拟流区域进行计算,得到各个拟流区域出口处的流量、含水率及流管的饱和度分布,按照如下公式并对拟流区域出口处指标进行求和,即可得到整管个二维流场的含水率;
式中,f为二维流场的总含水率,小数;fwl为流管l的出口端含水率,小数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910134624.5A CN109858177B (zh) | 2019-02-23 | 2019-02-23 | 一种基于拟流线法的底水油藏水平井水驱数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910134624.5A CN109858177B (zh) | 2019-02-23 | 2019-02-23 | 一种基于拟流线法的底水油藏水平井水驱数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109858177A true CN109858177A (zh) | 2019-06-07 |
CN109858177B CN109858177B (zh) | 2022-06-07 |
Family
ID=66898669
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910134624.5A Active CN109858177B (zh) | 2019-02-23 | 2019-02-23 | 一种基于拟流线法的底水油藏水平井水驱数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109858177B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110259444A (zh) * | 2019-07-15 | 2019-09-20 | 西南石油大学 | 基于流场诊断的水驱油藏渗流场可视化表征与评价方法 |
CN112302638A (zh) * | 2019-07-31 | 2021-02-02 | 中国石油天然气股份有限公司 | 一种油藏的流线场表征方法及装置 |
CN112685970A (zh) * | 2020-12-25 | 2021-04-20 | 中国石油大学(华东) | 一种水驱油藏流动单元渗流界面定量表征方法及系统 |
CN113356817A (zh) * | 2021-07-19 | 2021-09-07 | 中海油田服务股份有限公司 | 一种底水油藏水平井基于延长无水采油期的射孔优化设计方法 |
CN113554259A (zh) * | 2020-04-26 | 2021-10-26 | 中国石油化工股份有限公司 | 一种渗流阻力与驱替压力双场匹配优化评价方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103924966A (zh) * | 2013-01-10 | 2014-07-16 | 中国石油化工股份有限公司 | 基于储层物性时变模型的层系井网再建方法 |
CN104912551A (zh) * | 2015-05-08 | 2015-09-16 | 中国海洋石油总公司 | 一种油水相渗曲线和驱油效率的标定方法 |
CN105239976A (zh) * | 2015-11-04 | 2016-01-13 | 中国海洋石油总公司 | 考虑水驱油藏渗流参数动态变化的数值模拟方法 |
CN107038268A (zh) * | 2016-10-13 | 2017-08-11 | 中国石油化工股份有限公司 | 一种确定非均质储层五点井网水驱波及系数的方法 |
-
2019
- 2019-02-23 CN CN201910134624.5A patent/CN109858177B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103924966A (zh) * | 2013-01-10 | 2014-07-16 | 中国石油化工股份有限公司 | 基于储层物性时变模型的层系井网再建方法 |
CN104912551A (zh) * | 2015-05-08 | 2015-09-16 | 中国海洋石油总公司 | 一种油水相渗曲线和驱油效率的标定方法 |
CN105239976A (zh) * | 2015-11-04 | 2016-01-13 | 中国海洋石油总公司 | 考虑水驱油藏渗流参数动态变化的数值模拟方法 |
CN107038268A (zh) * | 2016-10-13 | 2017-08-11 | 中国石油化工股份有限公司 | 一种确定非均质储层五点井网水驱波及系数的方法 |
Non-Patent Citations (2)
Title |
---|
CHANGWEI LIU 等: "Experimental studies on production performance of oil reservoirs with bottom water", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 * |
吴斌: "多层油藏五点井网注水开发指标研究", 《内蒙古石油化工》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110259444A (zh) * | 2019-07-15 | 2019-09-20 | 西南石油大学 | 基于流场诊断的水驱油藏渗流场可视化表征与评价方法 |
CN110259444B (zh) * | 2019-07-15 | 2022-04-22 | 西南石油大学 | 基于流场诊断的水驱油藏渗流场可视化表征与评价方法 |
CN112302638A (zh) * | 2019-07-31 | 2021-02-02 | 中国石油天然气股份有限公司 | 一种油藏的流线场表征方法及装置 |
CN112302638B (zh) * | 2019-07-31 | 2023-07-25 | 中国石油天然气股份有限公司 | 一种油藏的流线场表征方法及装置 |
CN113554259A (zh) * | 2020-04-26 | 2021-10-26 | 中国石油化工股份有限公司 | 一种渗流阻力与驱替压力双场匹配优化评价方法 |
CN113554259B (zh) * | 2020-04-26 | 2024-04-16 | 中国石油化工股份有限公司 | 一种渗流阻力与驱替压力双场匹配优化评价方法 |
CN112685970A (zh) * | 2020-12-25 | 2021-04-20 | 中国石油大学(华东) | 一种水驱油藏流动单元渗流界面定量表征方法及系统 |
CN113356817A (zh) * | 2021-07-19 | 2021-09-07 | 中海油田服务股份有限公司 | 一种底水油藏水平井基于延长无水采油期的射孔优化设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109858177B (zh) | 2022-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109858177A (zh) | 一种基于拟流线法的底水油藏水平井水驱数值模拟方法 | |
CN103939066B (zh) | 一种一注多采井组定注水量确定油井产液量的方法 | |
CN104234677B (zh) | 一种注气垂直驱替提高凝析气藏凝析油采收率方法 | |
CN104989341B (zh) | 一种确定低渗透油藏有效驱替注采井距的方法 | |
CN105041307A (zh) | 一种碎屑岩油气储层优势渗流通道识别的工业化流程 | |
CN107437127B (zh) | 一种油井停喷地层压力预测方法 | |
CN103246820A (zh) | 一种油气藏数值模拟计算方法 | |
CN103899285A (zh) | 多层砂岩油藏近极限含水期轮替水驱方法 | |
CN109958413B (zh) | 一种特高含水期油藏动态流动单元划分方法 | |
CN113836695B (zh) | 一种基于无网格连接元的油藏数值模拟方法 | |
CN106909757A (zh) | 一种低渗透油藏超前注水地层合理压力水平确定方法 | |
CN109281642A (zh) | 一种改善聚合物驱后多层油藏开发效果的方法 | |
CN115270662A (zh) | 一种欠饱和储层煤层气井生产边界定量预测的计算方法 | |
CN116108572A (zh) | 页岩凝析气井体积压裂外区产能贡献分析方法 | |
CN113107475B (zh) | 一种气驱油藏的单井动态储量与可采储量确定方法及系统 | |
CN111608647B (zh) | 一种注水注聚井区地层压力预测方法 | |
CN109296363A (zh) | 特低渗透油藏co2驱初期产能预测方法 | |
CN114429085A (zh) | 一种用于分析缝洞型油藏流体势的方法及系统 | |
CN114757029B (zh) | 海上长水平井α-β波多级降排充填施工模拟方法及系统 | |
CN109915123A (zh) | 利用瞬时方程计算缝洞型油藏拟稳态水侵水侵量的方法 | |
CN113468476B (zh) | 一种聚合物驱增油量预测方法及装置 | |
CN115841083A (zh) | 确定注水井压驱配注量的方法 | |
CN117390980A (zh) | 二维纵向剖面层内差异化剩余油预测方法 | |
CN113756770B (zh) | 刻画非均相调驱累计降低含水率效果的方法 | |
Sazonov | Hydrodynamic simulation of reservoir fluids filtration at diverter technology conditions |
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 |