CN106246171B - 部分连通断层边界的无限大双重介质油藏数学建模方法 - Google Patents

部分连通断层边界的无限大双重介质油藏数学建模方法 Download PDF

Info

Publication number
CN106246171B
CN106246171B CN201610814107.9A CN201610814107A CN106246171B CN 106246171 B CN106246171 B CN 106246171B CN 201610814107 A CN201610814107 A CN 201610814107A CN 106246171 B CN106246171 B CN 106246171B
Authority
CN
China
Prior art keywords
area
pressure
well
fault boundary
boundary
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.)
Expired - Fee Related
Application number
CN201610814107.9A
Other languages
English (en)
Other versions
CN106246171A (zh
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 CN201610814107.9A priority Critical patent/CN106246171B/zh
Publication of CN106246171A publication Critical patent/CN106246171A/zh
Application granted granted Critical
Publication of CN106246171B publication Critical patent/CN106246171B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (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)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提供一种部分连通断层边界的无限大双重介质油藏数学建模方法，首先考虑部分连通断层边界的两区无限大双重介质油藏的不稳定压力特性；建立存在部分连通断层边界的两区无限大双重介质油藏的试井解释模型；然后对数学模型进行求解，并对无因次试井模型式进行Laplace变换，得到类似的渗流微分方程；最终得到实空间内的考虑井储和表皮效应影响的井底压力pwfD(xD,yD,tD)的数值解，从而绘制两区无限大双重介质油藏中存在部分连通断层边界的典型曲线；利用点源函数计算，得到不稳定试井解释模型和相关曲线，对部分连通断层边界的无限大双重介质油藏状况分析提供可靠依据。

Description

部分连通断层边界的无限大双重介质油藏数学建模方法
技术领域
本发明属于地质勘探领域，涉及一种部分连通断层边界的无限大双重介质油藏数学建模方法。
背景技术
目前国内外学者对断层边界的研究及相应的压力和压力导数曲线的分析都是针对封闭断层，而对非封闭断层的研究相对较少，对存在非封闭断层边界的试井解释模型研究也较少。对存在部分连通断层边界的双重介质油藏模型研究也很少。
下面是针对部分连通断层的试井模型方面的研究现状：
1.Bixel等人认为部分连通断层实际上是油层性质在平面上线性不连续的分布，最先对此类线性不连续边界进行研究。
2.1976年，Cinco等人基于点源函数的理论建立了存在无限导流裂缝或者断层边界的理论模型。
3.Stewart和Streltsova通过数值模拟研究了部分连通断层对干扰试井的影响。
4.Yaxely推导了无限大均质油藏中含有部分连通断层模型的解析解。
5.Ambastha进一步研究了复合油藏中存在部分连通断层边界的井底压力响应特征，他引入了无限小表皮的概念，将局部连通断层视为一个无限薄的表皮边界。
6.Abbaszadeh用点源函数方法建立了考虑断层内部流体流动的有限导流断层与储层耦合模型。
目前存在的针对部分连通断层边界的试井模型，是建立在均质地层基础上的；而在实际上，很多存在双重介质特性的地层中，也有部分连通断层边界的存在，现有方法无法正确反应部分连通断层边界的试井状况。
发明内容
本发明目的在于提供一种针对部分连通断层边界的两区无限大双重介质油藏不稳定试井模型的数学建模方法，利用点源函数计算，得到不稳定试井解释模型和相关曲线。
为得到上述目的，本发明采用如下方案：
部分连通断层边界的无限大双重介质油藏数学建模方法，其特征在于，包括以下步骤：
首先考虑部分连通断层边界的两区无限大双重介质油藏的不稳定压力特性；根据假设条件，基于渗流力学基本理论和点源函数计算方法，建立存在部分连通断层边界的两区无限大双重介质油藏的试井解释模型；
根据假设条件，激动井为定产量线源，且以激动井产量q生产，根据渗流力学的相关知识，建立两区渗流微分方程如下：
I区，x＞0：
II区，x＜0：
上式中，为单位体积下基岩向裂缝系统的窜流量，α为形状因子；pmi为基岩系统的压力；pfi为裂缝系统的压力；pf1为I区裂缝系统的压力；pm1为I区基岩系统的压力；kf2为II区裂缝系统的压力；μ2为II区流体的压力；km2为II区基岩系统的压力；pm2为II区基岩系统的压力；pf2为II区裂缝系统的压力；φf为裂缝孔隙度；Ctf为裂缝系统的压缩系数；φm为基质系统的孔隙度；
引入无因次变量，并将它们代入渗流微分方程及边界条件中，得到两区无限大双重介质油藏中存在部分连通断层的无因次试井解释模型；
然后对数学模型进行求解，并对无因次试井模型式进行Laplace变换，得到类似的渗流微分方程；
与渗流微分方程形式上相似，采用与之相同的求解方法，得到如下压力表达式：
其中，无因次的井底流压；xD为x方向的无因次坐标；yD为y方向的无因次坐标；α1为傅里叶变换的中间变量；α2为傅里叶变换的中间变量；M为I区II区的流度比；hD为无因次厚度比；为井到边界的距离；
根据Duhamel原理将井储和表皮叠加进去，再利用Stehfest数值反演对其进行拉斯逆变换，编制计算程序，最终得到实空间内的考虑井储和表皮效应影响的井底压力pwfD(xD,yD,tD)的数值解，从而绘制两区无限大双重介质油藏中存在部分连通断层边界的典型曲线；所述tD为无因次时间；
上述公式中符号说明如下：
q 激动井产量,m3/d
t 时间,hr
x x坐标,m
y y坐标,m
a 井的x坐标,m
b 井的y坐标,m
S 井表皮,无因次
δ δ函数
ηD 无因次扩散系数
z 拉普拉斯变换参数。
本发明对存在部分连通断层边界的无限大双重介质油藏进行数学建模，利用点源函数计算，得到不稳定试井解释模型和相关曲线，对部分连通断层边界的无限大双重介质油藏状况分析提供可靠依据。
附图说明
图1距离aD对典型曲线的影响示意图
图2界面表皮S对典型曲线的影响示意图
图3窜流系数λ1对典型曲线的影响示意图
图4窜流系数λ2对典型曲线的影响示意图
图5储容比ω1对典型曲线的影响示意图
图6储容比ω2对典型曲线的影响示意图
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
部分连通断层边界的无限大双重介质油藏数学建模方法，其特征在于，包括以下步骤：
首先考虑部分连通断层边界的两区无限大双重介质油藏的不稳定压力特性；根据假设条件，基于渗流力学基本理论和点源函数计算方法，建立存在部分连通断层边界的两区无限大双重介质油藏的试井解释模型；
根据假设条件，激动井为定产量线源，且以激动井产量q生产，根据渗流力学的相关知识，建立两区渗流微分方程如下：
I区，x＞0：
II区，x＜0：
上式中，为单位体积下基岩向裂缝系统的窜流量，α为形状因子；pmi为基岩系统的压力；pfi为裂缝系统的压力；pf1为I区裂缝系统的压力；pm1为I区基岩系统的压力；kf2为II区裂缝系统的压力；μ2为II区流体的压力；km2为II区基岩系统的压力；pm2为II区基岩系统的压力；pf2为II区裂缝系统的压力；φf为裂缝孔隙度；Ctf为裂缝系统的压缩系数；φm为基质系统的孔隙度；
引入无因次变量，并将它们代入渗流微分方程及边界条件中，得到两区无限大双重介质油藏中存在部分连通断层的无因次试井解释模型；
然后对数学模型进行求解，并对无因次试井模型式进行Laplace变换，得到类似的渗流微分方程；
与渗流微分方程形式上相似，采用与之相同的求解方法，得到如下压力表达式：
其中，无因次的井底流压；xD为x方向的无因次坐标；yD为y方向的无因次坐标；α1为傅里叶变换的中间变量；α2为傅里叶变换的中间变量；M为I区II区的流度比；hD为无因次厚度比；为井到边界的距离；
根据Duhamel原理将井储和表皮叠加进去，再利用Stehfest数值反演对其进行拉斯逆变换，编制计算程序，最终得到实空间内的考虑井储和表皮效应影响的井底压力pwfD(xD,yD,tD)的数值解，从而绘制两区无限大双重介质油藏中存在部分连通断层边界的典型曲线；所述tD为无因次时间；
上述公式中符号说明如下：
q 激动井产量,m3/d
t 时间,hr
x x坐标,m
y y坐标,m
a 井的x坐标,m
b 井的y坐标,m
S 井表皮,无因次
δ δ函数
ηD 无因次扩散系数
z 拉普拉斯变换参数。
图1所示为断层界面表皮S＝104时，井到断层的距离aD对典型曲线的影响。从图中可知，经过早期井储阶段后，当aD较大，如图aD＝1000，压力波还没有传到断层边界时，先出现I区裂缝系统的径向流阶段、基岩向裂缝系统的窜流阶段。当压力波到达断层边界，由于界面表皮S的存在，压力导数曲线上翘。当断层边界的影响结束后，出现II区系统的窜流阶段，由于在很短的时间内整个系统就达到了径向流阶段，因此II区双孔介质的反映不明显。
当aD较小，如图aD＝50时，I区系统的窜流发生在压力波传播断层边界之后。因此，经过井储阶段后，压力导数曲线先出现上翘，之后出现反映I区系统窜流的“凹子”，最后当压力波到达II区，出现II区的窜流及整个系统的径向流阶段。
当aD的值居中，如图5 0<aD<1000时，I区系统的窜流与断层界面表皮的影响同时进行，之后曲线的反映与前面讨论的情况类似。
界面表皮S的影响
图2所示为当I区系统的窜流发生在压力波传播到断层边界之前时，断层的界面表皮S对典型曲线的影响。从图中可知，经过井储阶段、I区裂缝系统的径向流阶段、基岩向裂缝系统的窜流阶段后，压力波到达断层边界，界面表皮的影响开始，压力导数曲线出现上升，并且S值越大，压力导数曲线上升越多，当S＝106时，表现出近似于封闭断层的特征，压力导数曲线出现值为1的水平线，并且界面表皮的影响掩盖了II区双孔介质的反映。
窜流系数的影响
图3和图4表示的是I、II区窜流系数λ1和λ2对典型曲线的影响。从图中可以看出，I区窜流系数λ1主要影响压力导数曲线上“凹子”出现时间的早晚，其它参数一定时，λ1越大，“凹子”出现时间越早，界面表皮的影响就发生在I区系统窜流之后，λ1越小，“凹子”出现时间越晚，相应的界面表皮的影响就发生在I区系统窜流之前。同样的，II区窜流系数λ2主要影响导数曲线上第二个“凹子”出现时间的早晚，λ2越大，“凹子”出现时间越早，但总的来说II区窜流影响不明显，这是因为压力波传播到II区后，在很短的时间内，就出现了整个系统的径向流阶段。
图5和图5-6表示的是I区系统的窜流发生在压力波传播到断层边界之前时，I、II区储容比ω1和ω2对典型曲线的影响。从图中可以看出，I区储容比ω1主要影响反映I区窜流的“凹子”形态，ω1越小，过渡段越长，“凹子”就越宽越深，之后出现部分连通断层界面表皮的影响，压力导数曲线上翘。II区储能比ω2主要影响第二个“凹子”形态，但总的来说ω2的影响不如ω1明显。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (1)

1.部分连通断层边界的无限大双重介质油藏数学建模方法，其特征在于，包括以下步骤：
首先考虑部分连通断层边界的两区无限大双重介质油藏的不稳定压力特性；根据假设条件，基于渗流力学基本理论和点源函数计算方法，建立存在部分连通断层边界的两区无限大双重介质油藏的试井解释模型；
根据假设条件，激动井为定产量线源，且以激动井产量q生产，根据渗流力学的相关知识，建立两区渗流微分方程如下：
I区，x＞0：
II区，x＜0：
上式中，为单位体积下基岩向裂缝系统的窜流量，α为形状因子；pmi为基岩系统的压力；pfi为裂缝系统的压力；pf1为I区裂缝系统的压力；pm1为I区基岩系统的压力；kf2为II区裂缝系统的压力；μ2为II区流体的压力；km2为II区基岩系统的压力；pm2为II区基岩系统的压力；pf2为II区裂缝系统的压力；φf为裂缝孔隙度；Ctf为裂缝系统的压缩系数；φm为基质系统的孔隙度；
引入无因次变量，并将它们代入渗流微分方程及边界条件中，得到两区无限大双重介质油藏中存在部分连通断层的无因次试井解释模型；
然后对数学模型进行求解，并对无因次试井模型式进行Laplace变换，得到类似的渗流微分方程；
与渗流微分方程形式上相似，采用与之相同的求解方法，得到如下压力表达式：
其中，无因次的井底流压；xD为x方向的无因次坐标；yD为y方向的无因次坐标；α1为傅里叶变换的中间变量；α2为傅里叶变换的中间变量；M为I区II区的流度比；hD为无因次厚度比；为井到边界的距离；
根据Duhamel原理将井储和表皮叠加进去，再利用Stehfest数值反演对其进行拉斯逆变换，编制计算程序，最终得到实空间内的考虑井储和表皮效应影响的井底压力pwfD(xD,yD,tD)的数值解，从而绘制两区无限大双重介质油藏中存在部分连通断层边界的典型曲线；所述tD为无因次时间；
上述公式中符号说明如下：
q激动井产量,m3/d
t时间,hr
x x坐标,m
y y坐标,m
a井的x坐标,m
b井的y坐标,m
S井表皮,无因次
δδ函数
ηD无因次扩散系数
z拉普拉斯变换参数。
CN201610814107.9A 2016-09-09 2016-09-09 部分连通断层边界的无限大双重介质油藏数学建模方法 Expired - Fee Related CN106246171B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610814107.9A CN106246171B (zh) 2016-09-09 2016-09-09 部分连通断层边界的无限大双重介质油藏数学建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610814107.9A CN106246171B (zh) 2016-09-09 2016-09-09 部分连通断层边界的无限大双重介质油藏数学建模方法

Publications (2)

Publication Number Publication Date
CN106246171A CN106246171A (zh) 2016-12-21
CN106246171B true CN106246171B (zh) 2019-09-20

Family

ID=57599541

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610814107.9A Expired - Fee Related CN106246171B (zh) 2016-09-09 2016-09-09 部分连通断层边界的无限大双重介质油藏数学建模方法

Country Status (1)

Country Link
CN (1) CN106246171B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108596392A (zh) * 2018-04-27 2018-09-28 中国石油大学(华东) 一种基于相似理论的水合物藏产能预测方法
CN110593860B (zh) * 2018-05-25 2021-08-13 中国科学技术大学 一种参数计算方法、装置及设备
CN109033519A (zh) * 2018-06-22 2018-12-18 中国石油天然气股份有限公司 一种异常高压碳酸盐岩挥发性油藏的试井解释方法及装置
CN109885894B (zh) * 2019-01-24 2023-05-05 江西理工大学 一种基于Trefftz多源点配点法的稀土矿层地下水渗流模型
CN110821486B (zh) * 2019-11-18 2022-04-01 西南石油大学 一种储层优势通道物性参数计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4797821A (en) * 1987-04-02 1989-01-10 Halliburton Company Method of analyzing naturally fractured reservoirs
CN101899972A (zh) * 2010-06-30 2010-12-01 中国石油大学(北京) 砂岩油藏注水开发可预测物理模型建立方法
CN103278436A (zh) * 2013-02-01 2013-09-04 西安石油大学 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法
CN105044283A (zh) * 2015-07-08 2015-11-11 中国海洋石油总公司 一种考虑聚合物粘弹性的多油层注聚井注入能力评价方法
CN105569641A (zh) * 2015-12-17 2016-05-11 西南石油大学 一种改善近井压降的水平井瞬态压力快速计算模型建立方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7200539B2 (en) * 2001-02-21 2007-04-03 Baker Hughes Incorporated Method of predicting the on-set of formation solid production in high-rate perforated and open hole gas wells

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4797821A (en) * 1987-04-02 1989-01-10 Halliburton Company Method of analyzing naturally fractured reservoirs
CN101899972A (zh) * 2010-06-30 2010-12-01 中国石油大学(北京) 砂岩油藏注水开发可预测物理模型建立方法
CN103278436A (zh) * 2013-02-01 2013-09-04 西安石油大学 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法
CN105044283A (zh) * 2015-07-08 2015-11-11 中国海洋石油总公司 一种考虑聚合物粘弹性的多油层注聚井注入能力评价方法
CN105569641A (zh) * 2015-12-17 2016-05-11 西南石油大学 一种改善近井压降的水平井瞬态压力快速计算模型建立方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《存在局部连通断层的复合油藏不稳定试井分析及典型曲线拟合方法》;张望明等;《石油勘探与开发》;19920228;第19卷(第2期);第69-76页 *
《非封闭断层的试井解释方法》;李传亮;《新疆石油地质》;19971231;第18卷(第4期);第370-376页 *

Also Published As

Publication number Publication date
CN106246171A (zh) 2016-12-21

Similar Documents

Publication Publication Date Title
CN106246171B (zh) 部分连通断层边界的无限大双重介质油藏数学建模方法
Haddad et al. XFEM-based CZM for the simulation of 3D multiple-cluster hydraulic fracturing in quasi-brittle shale formations
CN105840187B (zh) 致密性油藏水平井分段压裂产能计算方法
Wang et al. A review of stimulated reservoir volume characterization for multiple fractured horizontal well in unconventional reservoirs
Wang et al. Anisotropic permeability evolution of coal with effective stress variation and gas sorption: model development and analysis
Zeng et al. Numerical simulation of fracture network generation in naturally fractured reservoirs
Xu et al. Characterization of hydraulically-induced fracture network using treatment and microseismic data in a tight-gas formation: A geomechanical approach
Nie et al. Dual porosity and dual permeability modeling of horizontal well in naturally fractured reservoir
Hofmann et al. Numerical simulation of complex fracture network development by hydraulic fracturing in naturally fractured ultratight formations
CN108319756A (zh) 一种致密储层体积压裂缝网扩展模拟及表征方法
Yu et al. 3D analytical modeling of hydraulic fracturing stimulated reservoir volume
Zhao et al. Numerical investigation of shut-in time on stress evolution and tight oil production
Wang et al. Analytical solutions for a quad-linear flow model derived for multistage fractured horizontal wells in tight oil reservoirs
Haddad et al. Integration of XFEM and CZM to model 3D multiple-stage hydraulic fracturing in quasi-brittle shale formations: solution-dependent propagation direction
CN111927420A (zh) 一种任意形状气藏有限导流非对称裂缝井压力模拟方法
Luo et al. A new fracture-unit model and its application to a z-fold fracture
Brown Analytical trilinear pressure transient model for multiply fractured horizontal wells in tight shale reservoirs
CN110094196A (zh) 一种碳酸盐岩裸眼水平井分段酸压效果评价方法
Liu et al. A simplified and efficient method for water flooding production index calculations in low permeable fractured reservoir
Wang Asymptotic solutions to the diffusivity equation: validation and field applications
Li et al. Extended finite element method for analysis of multi-scale flow in fractured shale gas reservoirs
CN110390154A (zh) 一种提高复杂断块油气田油藏数值模拟效率的方法
He et al. Numerical investigation of interaction between hydraulic and natural fractures under various geologic and engineering parameters using particle flow code method
Lu et al. Numerical simulation on the initiation and propagation of synchronous perforating fractures in horizontal well clusters
CN108920824A (zh) 一种窄河道储层的产量递减分析方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Feng Guoqing

Inventor after: Guo Yu

Inventor before: Feng Guoqing

Inventor before: Kong Bing

GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190920

Termination date: 20200909