CN109684685B - 一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法 - Google Patents
一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法 Download PDFInfo
- Publication number
- CN109684685B CN109684685B CN201811500559.5A CN201811500559A CN109684685B CN 109684685 B CN109684685 B CN 109684685B CN 201811500559 A CN201811500559 A CN 201811500559A CN 109684685 B CN109684685 B CN 109684685B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- multilateral well
- multilateral
- hydrate
- pressure
- 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.)
- Active
Links
Classifications
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法,包括(1):建立多分支井水合物降压开采储层的地质模型;(2)将地质模型进行非结构网格离散化;(3)基于固体沉积物固体变形,建立水合物降压开采的多场流固耦合数学模型;(4)在建立的非结构网格***上,利用Galerkin有限元方法建立单元刚度方程,并对多场流固耦合数学模型进行离散求解,获得多分支井上的应力分布和压力分布;(5)将储层中的应力与水合物储层的摩尔‑库伦强度对比,分析储层稳定性;(6)、利用压力计算分支孔上的压力梯度,实现对多分支井开采产能的分析,为多分支井开采产能定量化及稳定性分析提供了有效的描述手段,具有较高的实际应用及参考价值。
Description
技术领域
本发明涉及海洋天然气水合物资源开采数值模拟技术领域,具体涉及一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法。
背景技术
天然气水合物储量巨大,是未来满足人类能源需求的高效清洁能源。目前天然气水合物的开采处于试验探索阶段。国际上的天然气水合物试开采表明:低效率、储层失稳和出砂是制约水合物开采的重要因素。
为提高水合物开采效率,许多学者提出了水平井、双水平井和水平井五点井网的开采方式。比如,授权公告号为【CN106761587B】的发明专利公开一种海洋粉砂质储层天然气水合物多分支孔有限防砂开采方法,该方法利用多分支井增加泄油面积,提供开采效率。水合物降压开采是通过降低井底压力在储层中形成压降漏斗,当储层压力减低到水合物相平衡压力时,水合物开始分解。理论上,井底降压幅度越大、压降漏斗影响范围越大,产气速率就越大。
由于海洋天然气水合物沉积物储层胶结差,压力减低会导致储层内应力改变和变形,可能引起储层失稳破坏;同时水合物又在沉积物孔隙介质间起胶结作用,水合物分解会降低储层的强度,又进一步加大储层失稳的风险。特别是对于多分支井方法,由于多分支井的存在降低了主井眼周围的储层强度,同时高的产气速率意味着更大的储层失稳风险。因此,多分支井开采方法必须平衡产气速率和储层稳定性两者之间的关系。
然而,目前多分支井开采缺乏产能定量化的描述手段,现有的数值模拟器无法适用于水合物多分支井开采的模拟。主要困难是多分支井井筒结构复杂,现有的模拟方法和数值模拟器无法建立准确的地质模型,从而无法分析多分支井开采产能和力学响应特征。
因此,为满足水合物多分支井开采的参数设计,定量获取多分支井开采的产能,分析其储层稳定性的特征,亟待提出一种针对多分支井降压开采产能和储层稳定性分析的技术发方案。
发明内容
本发明针对目前多分支井开采产能和储层稳定性方法缺乏的问题,提出一种多分支井水合物降压开采条件下产能和储层稳定性的分析方法。
本发明是采用以下的技术方案实现的:一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法,其特征在于,包括以下步骤:
步骤A、建立多分支井水合物降压开采储层的地质模型,所述地质模型包括多分支井筒内边界和储层外边界;
步骤B、将地质模型进行非结构网格离散化,得到地质模型的非结构网格***,具体包括:
(1)输入多分支井地质模型的多分支井内边界和储层外边界信息;
(2)在多分支井内边界和储层外边界上设置离散点;
(3)以边界离散点为约束,对整个地质模型进行网格离散,得到模型的非结构网格***;
步骤C、考虑固体沉积物固体变形因素,建立水合物降压开采的多场流固耦合数学模型;
所述多场流固耦合数学模型包括气相压力方程、水相饱和度方程、水合物分解方程、沉积物固体变形方程、初始条件、多分支井内边界条件和储层外边界条件;
步骤D、在步骤B建立的非结构网格***的基础上,利用Galerkin有限元方法建立单元刚度方程,对多场流固耦合数学模型进行离散求解,获得多分支井上的应力分布和压力分布,包括:
(1)在非结构网格上计算每一个单元的单元刚度方程;
(2)由结点位置关系,获得总体刚度方程;
(3)基于初始条件和多分支井内边界条件和储层外边界条件求解总体刚度方程,获得各时刻的压力分布、应力分布;
步骤E、将储层中的应力与水合物储层的摩尔-库伦强度对比,分析储层的稳定性;
步骤F、根据计算得到的储层压力分布,计算各储层时刻压力梯度,进而计算单个多分支井上各时刻的产水和产气速率,实现对水合物产能的分析。
进一步的,所述步骤A中在建立地质模型时具体通过以下方式实现:
(1)根据主井眼的尺寸,建立主井眼圆柱体模型;
(2)根据多分支井的长度和大小,建立单个多分支井圆柱体模型;
(3)根据多分支井的数量、分布以及与主井眼的位置关系,建立多分支井和主井眼的模型;
(4)根据储层的厚度和大小,建立储层圆柱体模型;
(5)采用布尔运算,用储层圆柱体减去主井眼和多分支井眼的圆柱体,获得多分支井开采储层的地质模型。
进一步的,所述步骤B中,在对地质模型进行离散化之前,先将井筒附近网格进行加密处理,即设置地质模型的井筒内边界上的取样点的间距小于储层外边界上的取样点的间距。
进一步的,所述步骤C中,所述的沉积物固体变形方程为:
其中,A3=Gm,Gm为水合物储层抗剪强度,MPa;νm为水合物储层泊松比,wx,wy,wz表示沉积物固体变形位移,peff为有效孔隙压力。
进一步的,所述步骤C中,所述的多分支井内边界条件为:
其中,pg为气体压力,Γi表示多分支井内边界,pwf为井底压力。
进一步的,所述步骤C中,所述的储层外边界条件:
其中,Γo表示储层外边界;no表示储层外边界的法线方向。
进一步的,所述步骤E对储层稳定性分析时,具体包括以下步骤:
(1)根据计算得到的储层应力,绘制多分支孔上最大有效主应力和最小有效主应力随时间的关系曲线;
(2)将储层沉积物的摩尔-库伦强度线与最大和最小有效主应力的曲线对比,如果最大有效主应力和最小有效主应力的曲线远离摩尔-库伦强度线,则说明储层是稳定的,否则储层不稳定。
进一步的,所述步骤F中,对于获得的单个多分支井上各时刻的产水和产气速率,将各多分支井的产水产气相加,进而得到总产气流量和产水流量,具体的:
(1)利用各时刻储层压力分布计算多分支井上的各时刻的压力梯度;
求得各网格结点压力后,利用总体刚度方程,提取其包含的边界梯度信息,得到:
其中,krg为气相水两相渗流的气相相对渗透率,k为沉积物多孔介质的绝对渗透率,ρg
为气相密度,pg为气体压力,μg为甲烷气体黏度,A为气体压力变化率矩阵,N为单元形函数,n为边界s的外法线,Q为刚度矩阵,{fg′}为去掉边界条件后的右端项;
求解上述方程(13)计算得到多分支井边界上的压力梯度:
(2)利用压力梯度计算单个多分支井上的各时刻的水流量和气流量;
其中Qg为单个多分支井的产气量,ri为多分支井的半径,L为多分支井的长度;
(3)将各时刻所有多分支井的水流量相加、气流量相加,获得多分支井水流量和气流量
与现有技术相比,本发明的优点和积极效果在于:
本发明所提出的方案基于非结构网格实现了复杂多分支井的地质建模和离散化,建立了考虑沉积物固体变形的多场流固耦合数学模型,并基于有限元方法对水合物多分支井开采的多场流固耦合数学模型进行求解,获得了多分支孔开采的产能,基于摩尔-库伦强度理论,实现了水合物多分支井开采的储层稳定性的分析,为多分支井开采产能定量化及稳定性分析提供了有效的描述手段,具有较高的实际应用及参考价值。
附图说明
图1为本发明实施例所述多分支井模型结构示意图;
图2为本发明实施例所述多分支井开采储层模型结构示意图;
图3为基于摩尔-库伦强度理论的多分支孔降压开采储层稳定性分析示意图;
其中:1、主井眼;2、多分支井;3、水合物储层;4、下伏层;5、上覆层;6、含水合物沉积物的摩尔-库伦强度线;7、某多分支井上的应力随时间变化曲线;8、另一多分支井上的应力随时间变化曲线。
具体实施方式
为了能够更加清楚地理解本发明的上述目的、特征和优点,下面结合附图及实施例对本发明做进一步说明。需要说明的是,在不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。
一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法,包括以下步骤:
步骤A、根据地质资料,建立多分支井水合物降压开采储层的地质模型,所述地质模型包括多分支井筒内边界和储层外边界,具体实现方式为:
(1)根据主井眼的尺寸,建立圆柱体表示的主井眼;
(2)根据多分支井的长度和大小,建立单个多分支井的圆柱体;
(3)根据多分支井的数量、分布、与主井眼的位置关系,建立多分支井和主井眼的模型,如图1所示,井筒半径为0.1米,多分支井数为3,分支井长度为2米,半径为0.01米;
(4)根据储层的厚度和大小,建立储层柱体模型;
(5)使用布尔运算,用储层圆柱体减去主井眼和多分支井眼的圆柱体,获得多分支井开采储层的地质模型,如图2所示。
步骤B、将地质模型进行非结构网格离散化,得到地质模型的非结构网格***,具体处理方式为:
(1)输入多分支井地质模型的多分支井内边界和储层外边界信息;
(2)在多分支井边界和储层外边界上设置离散点;
(3)以边界离散点为约束,对整个地质模型进行网格离散,得到模型的非结构网格***,非结构网格***有单元坐标和单元结点编号表示;
为了确保计算精度,离散化之前,先将多分支井附近网格进行加密处理,比如,在Netgen 开源的网格划分程序中实现的具体方式为:设置地质模型的井筒内边界上的取样点的间距小于储层外边界上的取样点的间距。
步骤C、考虑固体沉积物固体变形,建立水合物降压开采的多场流固耦合数学模型,所述多场流固耦合数学模型包括气相压力方程、水相饱和度方程、水合物分解方程、沉积物固体变形方程、初始条件、多分支井内边界条件和储层外边界条件;
具体的:
1)气相压力方程:
2)水相饱和度方程:
3)水合物分解方程:
4)沉积物固体变形方程:
5)初始条件:
pg(t=0)=pi,Sw(t=0)=Swi,T(t=0)=Ti (6)
6)多分支井内边界条件:
7)储层外边界条件:
其中,ρg和ρw分别为气相和水相密度,kg/m3;Sg和Sw分别为气相和水相饱和度;krg和 krw分别为气相水两相渗流的气相相对渗透率和水相相对渗透率;k为沉积物多孔介质的绝对渗透率,m2;μg和μw分别为甲烷气体和水的黏度,Pa·s;Z为气体偏差因子;Rg为气体常数,8.314J/(mol·K);T为温度,K;φ为含水合物沉积物孔隙度;ρh为水合物密度,kg/m3; Sh为水合物饱和度;mh为单位体积水合物分解速率,kg/(m3·s);mg为单位体积水合物分解的产气速率,kg/(m3·s);mw为单位体积水合物分解产水速率,kg/(m3·s);为甲烷的摩尔质量,0.016kg/mol;pe为水合物平衡压力,Pa;pg为气体压力,Pa;pc为毛管力,Pa;kreac为水合物分解速率常数,mol/(m2·Pa·s);Ars为单位体积储层水合物分解的表面积,m-1;pi为储层初始压力,Pa;Swi为储层初始含水饱和度;Ti为储层初始温度,K;pwf为井底压力, Pa;Γi表示井筒边界;Γo表示外边界;no表示外边界的法线方向;wx,wy,wz表示沉积物固体变形位移,peff为有效孔隙压力,Pa;α为Biot系数;A3=Gm, Gm为水合物储层抗剪强度,MPa;νm为水合物储层泊松比。
步骤D、在步骤B建立的非结构网格***的基础上,利用Galerkin有限元方法建立单元刚度方程,对多场流固耦合数学模型进行离散求解,获得多分支井上的应力分布和压力分布,其基本原理为:1)在非结构网格上计算每一个单元的单元刚度方程;2)由结点位置关系,形成总体刚度方程;3)代入初始条件和边界条件求解总体刚度方程,获得各时刻的压力分布、应力分布,具体包括:
步骤D1、气相压力的有限元求解:
(1)气相压力方程的积分降阶:
取非结构网格***中的某一个单元,将气相压力方程(1)两端同乘以压力变分并在单元上积分、降阶:
(2)代入压力的插值函数计算刚度方程:
压力在单元上表示为:
pg=N{pg}e,δpg=N{δpg}e (9)
式中,N为单元形函数,s为单元边界;为单元边界法线方向的压力梯度,则最终单元刚度方程为:
其中,
代入插值函数,可得气相压力的单元刚度方程为:
其中, Si为四面体表面的面积;
(3)总体刚度方程叠加:
根据单元结点编号叠加形成总体刚度方程:
(4)求解最终总体刚度方程:结合相应的边界条件求解线性方程组获得网格结点的气相压力值。
步骤D2、水相饱和度的有限元求解:
(1)将水相饱和度方程(2)两端同乘以饱和度变分并在单元上积分、降阶:
(2)代入饱和度的插值函数计算单元刚度方程:
饱和度在单元上的插值为:
Sw=N{Sw}e,δSw=N{δSw}e (13)
最终饱和度的单元刚度方程为:
其中,Be=∫NφρwNTdV,
代入插值函数,可得水相饱和度的单元刚度方程为:
(3)总体刚度方程叠加:
所有网格的单元刚度方程按照单元结点编号进行叠加,形成最终的总体刚度方程;
步骤D3、沉积物固体变形方程求解:公式(5)中的三个方程分别表示x,y,z三个方向上的固体变形方程,三个方程有相同的形式,下面仅以x方向的方程为例说明其求解方法;
(1)固体变形方程的积分降阶:
将沉积物固体变形方程(5)中x方向的方程两端同乘以位移的变分并在单元上积分,利用高斯公司降阶可得:
(2)代入插值函数计算单元刚度方程:
位移在单元上的插值为:
wx=N{wx}e,δwx=N{δwx}e (17)
带入方程(16)可得x方向沉积物固体变形的有限元单元刚度方程为:
[k11]{wx}+[k12]{wy}+[k13]{wz}+[k14]{p}={f1} (18)
其中,
{f1}=-∫(FxNj)ds
采取同样的方法可得y方向和z方向的沉积物固体变形位移有限元单元刚度方程;
(3)总体刚度方程叠加:
将所有网格的单元刚度方程按照单元结点编号进行叠加,形成最终的沉积物固体变形位移的总体刚度方程;
(4)求解最终总体刚度方程:
结合相应的边界条件求解总体刚度方程的线性方程组获得网格结点的固体变形位移;
(5)根据位移计算储层应力:
将位移乘以应力矩阵计算得到储层的应力分布:
{σ}=S{wx} (19)
其中σ为应力,S为应力矩阵。
步骤E、将储层中的应力与水合物储层的摩尔-库伦强度对比,分析储层的稳定性,偏离摩尔-库伦强度线表明储层稳定;应力与摩尔-库伦线相交表明储层不稳定,具体的处理步骤为:
(1)根据计算得到的储层应力,绘制多分支孔上最大有效主应力和最小有效主应力随时间的关系曲线,如图3所示,其中7和8为多分支井上的应力随时间变化曲线;
(2)将储层沉积物的摩尔-库伦强度线与最大和最小有效主应力的曲线对比,如果最大有效主应力和最小有效主应力的曲线远离摩尔-库伦强度线,则说明储层是稳定的,否则储层不稳定,如图3所示,曲线7与摩尔-库伦强度线相交,说明该曲线所在的储层点会出现失稳;曲线8偏离摩尔-库伦强度线,则该曲线所在的储层点不会出现失稳。
步骤F、根据计算得到的储层压力分布,计算各储层时刻压力梯度,进而计算单个多分支井上各时刻的产水和产气速率,将各多分支井的产水产气相加,得到总产气速度和产水速度随时间的关系,具体包括:
(1)利用各时刻储层压力分布计算多分支井上的各时刻的压力梯度;
求得各网格结点压力后,利用总体刚度方程(12),将其中{fg}包含的边界梯度信息提取出来,方程写为:
其中{fg′}为去掉边界条件后的右端项;
求解上述方程计算得到多分支井边界上的压力梯度:
(2)利用压力梯度计算单个多分支井上的各时刻的水流量和气流量;
其中Qg为单个多分支井的产气量ri为多分支井的半径,L为多分支井的长度;
(3)将各时刻所有多分支井的水流量相加、气流量相加,获得多分支井水流量和气流量。
综上,本方案所提出的方法可以描述多分支井的井身结构,获得储层压力、温度和沉积物应力的分布特征,基于摩尔-库伦强度理论实现多分支井开采储层稳定性的分析,获得多分支孔开采条件下的产能,为多分支井开采产能定量化及稳定性分析提供了有效的描述手段。
以上所述,仅是本发明的较佳实施例而已,并非是对本发明作其它形式的限制,任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例应用于其它领域,但是凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型,仍属于本发明技术方案的保护范围。
Claims (4)
1.一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法,其特征在于,包括以下步骤:
步骤A、建立多分支井水合物降压开采储层的地质模型,所述地质模型包括多分支井筒内边界和储层外边界,具体包括以下步骤:
(1)根据主井眼的尺寸,建立主井眼圆柱体模型;
(2)根据多分支井的长度和大小,建立单个多分支井圆柱体模型;
(3)根据多分支井的数量、分布以及与主井眼的位置关系,建立多分支井和主井眼的模型;
(4)根据储层的厚度和大小,建立储层圆柱体模型;
(5)采用布尔运算,用储层圆柱体减去主井眼和多分支井眼的圆柱体,获得多分支井开采储层的地质模型;
步骤B、将地质模型进行非结构网格离散化,得到地质模型的非结构网格***,具体包括:
(1)输入多分支井地质模型的多分支井内边界和储层外边界信息;
(2)在多分支井内边界和储层外边界上设置离散点;
(3)以边界离散点为约束,对整个地质模型进行网格离散,得到模型的非结构网格***;
步骤C、考虑固体沉积物固体变形因素,建立水合物降压开采的多场流固耦合数学模型;
所述多场流固耦合数学模型包括气相压力方程、水相饱和度方程、水合物分解方程、沉积物固体变形方程、初始条件、多分支井内边界条件和储层外边界条件,其中:
所述的沉积物固体变形方程为:
其中,A3=Gm,Gm为水合物储层抗剪强度,MPa;νm为水合物储层泊松比,wx,wy,wz表示沉积物固体变形位移,peff为有效孔隙压力;
所述的多分支井内边界条件为:
其中,pg为气体压力,Γi表示多分支井内边界,pwf为井底压力;
步骤D、在步骤B建立的非结构网格***的基础上,利用Galerkin有限元方法建立单元刚度方程,对多场流固耦合数学模型进行离散求解,获得多分支井上的应力分布和压力分布,包括:
(1)在非结构网格上计算每一个单元的单元刚度方程;
(2)由结点位置关系,获得总体刚度方程;
(3)基于初始条件和多分支井内边界条件和储层外边界条件求解总体刚度方程,获得各时刻的压力分布、应力分布;
步骤E、将储层中的应力与水合物储层的摩尔-库伦强度对比,分析储层的稳定性;
步骤F、根据计算得到的储层压力分布,计算各储层时刻压力梯度,进而计算单个多分支井上各时刻的产水和产气速率,实现对水合物产能的分析;
对于获得的单个多分支井上各时刻的产水和产气速率,将各多分支井的产水产气相加,进而得到总产气流量和产水流量,具体的:
(1)利用各时刻储层压力分布计算多分支井上的各时刻的压力梯度;
求得各网格结点压力后,利用总体刚度方程,提取其包含的边界梯度信息,得到:
其中,krg为气相水两相渗流的气相相对渗透率,k为沉积物多孔介质的绝对渗透率,ρg为气相密度,pg为气体压力,μg为甲烷气体黏度,A为气体压力变化率矩阵,N为单元形函数,n为边界s的外法线,Q为刚度矩阵,{fg′}为去掉边界条件后的右端项;
进而得到多分支井边界上的压力梯度:
(2)利用压力梯度计算单个多分支井上的各时刻的水流量和气流量;
其中Qg为单个多分支井的产气量,ri为多分支井的半径,L为多分支井的长度;
(3)将各时刻所有多分支井的水流量相加、气流量相加,获得多分支井水流量和气流量。
2.根据权利要求1所述的多分支井水合物降压开采条件下的产能及储层稳定性分析方法,其特征在于:所述步骤B中,在对地质模型进行离散化之前,先将井筒附近网格进行加密处理,即设置地质模型的井筒内边界上的取样点的间距小于储层外边界上的取样点的间距。
3.根据权利要求1所述的多分支井水合物降压开采条件下的产能及储层稳定性分析方法,其特征在于:所述步骤C中,所述的储层外边界条件:
其中,Γo表示储层外边界;no表示储层外边界的法线方向。
4.根据权利要求1所述的多分支井水合物降压开采条件下的产能及储层稳定性分析方法,其特征在于:所述步骤E对储层稳定性分析时,具体包括以下步骤:
(1)根据计算得到的储层应力,绘制多分支孔上最大有效主应力和最小有效主应力随时间的关系曲线;
(2)将储层沉积物的摩尔-库伦强度线与最大和最小有效主应力的曲线对比,如果最大有效主应力和最小有效主应力的曲线远离摩尔-库伦强度线,则说明储层是稳定的,否则储层不稳定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811500559.5A CN109684685B (zh) | 2018-12-10 | 2018-12-10 | 一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811500559.5A CN109684685B (zh) | 2018-12-10 | 2018-12-10 | 一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109684685A CN109684685A (zh) | 2019-04-26 |
CN109684685B true CN109684685B (zh) | 2019-12-17 |
Family
ID=66187406
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811500559.5A Active CN109684685B (zh) | 2018-12-10 | 2018-12-10 | 一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109684685B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110259443A (zh) * | 2019-07-12 | 2019-09-20 | 西南石油大学 | 一种基于3dec离散元的煤层井壁稳定性预测方法 |
CN111507024B (zh) * | 2020-06-05 | 2023-10-13 | 南京大学 | 基于gpu矩阵的离散元流固耦合数值模拟方法及*** |
CN113591338B (zh) * | 2021-05-28 | 2023-07-28 | 河海大学 | 一种荷载和地下水开采引发地面沉降三维变参数全耦合模拟计算方法 |
CN114548691B (zh) * | 2022-01-25 | 2024-05-31 | 华南理工大学 | 一种深海水合物降压开采时储层稳定性安全评估方法 |
CN117211763A (zh) * | 2023-09-28 | 2023-12-12 | 北京科技大学 | 开采天然气水合物的多分支水平井产量预测方法及*** |
CN117522173B (zh) * | 2024-01-04 | 2024-04-26 | 山东科技大学 | 基于深度神经网络的天然气水合物降压开采产能预测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105223315A (zh) * | 2015-10-12 | 2016-01-06 | 大连理工大学 | 一种天然气水合物储层-结构物相互作用的计算预测方法 |
CN106103887A (zh) * | 2014-03-17 | 2016-11-09 | 沙特***石油公司 | 在储层仿真中对交叉故障和复杂井眼进行建模 |
CN108241777A (zh) * | 2017-12-27 | 2018-07-03 | 青岛海洋地质研究所 | 基于非结构网格有限元法计算水合物沉积物中渗流速度场的方法 |
-
2018
- 2018-12-10 CN CN201811500559.5A patent/CN109684685B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106103887A (zh) * | 2014-03-17 | 2016-11-09 | 沙特***石油公司 | 在储层仿真中对交叉故障和复杂井眼进行建模 |
CN105223315A (zh) * | 2015-10-12 | 2016-01-06 | 大连理工大学 | 一种天然气水合物储层-结构物相互作用的计算预测方法 |
CN108241777A (zh) * | 2017-12-27 | 2018-07-03 | 青岛海洋地质研究所 | 基于非结构网格有限元法计算水合物沉积物中渗流速度场的方法 |
Non-Patent Citations (2)
Title |
---|
万义钊 等.南海神狐海域天然气水合物降压开采过程中储层的稳定性.《开发工程》.2018,第38卷(第4期), * |
程家望 等.天然气水合物降压开采储层稳定性模型分析.《新 能 源 进 展》.2016,第4卷(第1期), * |
Also Published As
Publication number | Publication date |
---|---|
CN109684685A (zh) | 2019-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109684685B (zh) | 一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法 | |
Wilson et al. | Optimization of shale gas field development using direct search techniques and reduced-physics models | |
US8515720B2 (en) | Determine field fractures using geomechanical forward modeling | |
US8682628B2 (en) | Multiphase flow in a wellbore and connected hydraulic fracture | |
CN104533370A (zh) | 压裂水平井油藏、裂缝、井筒全耦合模拟方法 | |
Durlofsky et al. | Uncertainty quantification for subsurface flow problems using coarse-scale models | |
CN105386751B (zh) | 一种基于油藏渗流模型的水平井测井产能预测方法 | |
CN108241777B (zh) | 基于非结构网格有限元法计算水合物沉积物中渗流速度场的方法 | |
CN114169263B (zh) | 基于压裂液返排数据的页岩油气藏裂缝参数反演方法 | |
Bourbiaux et al. | A rapid and efficient methodology to convert fractured reservoir images into a dual-porosity model | |
Li et al. | Estimation of relative permeability by assisted history matching using the ensemble Kalman filter method | |
EP2831804B1 (en) | System and method for automatic local grid refinement in reservoir simulation systems | |
CN110738001B (zh) | 一种非常规储层压裂增产改造区计算方法 | |
CN107066718B (zh) | 一种四维动态地应力模拟方法 | |
CN114427432B (zh) | 一种气藏剩余气开发潜力确定方法 | |
CN107130959B (zh) | 一种煤层气产量预测方法 | |
Mascarenhas et al. | Coarse scale simulation of horizontal wells in heterogeneous reservoirs | |
CN105089615A (zh) | 一种基于油藏模型的测井数据历史回归处理方法 | |
CN112394404A (zh) | 一种渐进式储层精细表征方法 | |
CN111706318B (zh) | 一种确定低渗储层剩余油分布状况的方法 | |
Yermekov et al. | Features of permeability anisotropy accounting in the hydrodynamic model | |
US20220010678A1 (en) | Optimization of discrete fracture network (dfn) using streamlines and machine learning | |
Nakashima et al. | Near-well upscaling for three-phase flows | |
Milliotte et al. | Well-data-based discrete fracture and matrix modelling and flow-based upscaling of multilayer carbonate reservoir horizons | |
CA3196445A1 (en) | A hybrid method for reservoir simulation |
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 |